Codeforces 781D Axel and Marston in Bitland 矩阵 bitset
Posted zhouzhendong
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Codeforces 781D Axel and Marston in Bitland 矩阵 bitset相关的知识,希望对你有一定的参考价值。
原文链接https://www.cnblogs.com/zhouzhendong/p/CF781D.html
题目传送门 - CF781D
题意
有一个 n 个点的图,有 m 条有向边,边有两种类型:0 和 1 。
有一个序列,其构造方案为:初始情况为 0 ;对于当前串,将当前串的 0 变成 1 ,1 变成 0 ,接到当前串的后面,得到一个新串;不断进行这个操作。
最终得到一个形如 0110100110010110…… 的无限串。
问从节点 1 出发,依次经过上述串对应的 0/1 边,最多能走多少步。
如果答案大于 $10^{18}$ ,输出 -1 。
$nleq 500,mleq 2n^2$
题解
首先考虑处理出两个矩阵,分别处理 0/1 两种类型的边的转移。
于是我们可以直接矩阵倍增一下得到一个 $O(n^3 log 10^{18})$ 的做法。
再注意到我们只关系这些矩阵元素是否为 0 ,那么我们只需要用 bitset 优化一下矩阵乘法就好了。
时间复杂度 $O(n^3 log 10^{18} / 32)$ 。
代码
#include <bits/stdc++.h> using namespace std; typedef long long LL; LL read(){ LL x=0,f=1; char ch=getchar(); while (!isdigit(ch)&&ch!=‘-‘) ch=getchar(); if (ch==‘-‘) f=-1,ch=getchar(); while (isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar(); return x*f; } const int N=505; const LL INF=1000000000000000000LL; int n,m; struct Mat{ bitset <N> v[N]; Mat(){} Mat(int x){ for (int i=1;i<=n;i++) v[i].reset(); for (int i=1;i<=n;i++) v[i][i]=x; } friend Mat operator * (Mat a,Mat b){ Mat ans(0); for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) if (a.v[i][j]) ans.v[i]|=b.v[j]; return ans; } }M[2][35]; int main(){ n=read(),m=read(); M[0][0]=M[1][0]=Mat(0); for (int i=1;i<=m;i++){ int a=read(),b=read(),t=read(); M[t][0].v[a][b]=1; } int k=0; while (M[0][k].v[1].count()&&(1LL<<(k<<1))<=INF){ k++; M[0][k]=M[0][k-1]*M[1][k-1]*M[1][k-1]*M[0][k-1]; M[1][k]=M[1][k-1]*M[0][k-1]*M[0][k-1]*M[1][k-1]; } LL s=0,p=0; Mat now(0); now.v[1][1]=1; while (k>0&&s<=INF){ k--; if ((now*M[p][k]).v[1].count()){ now=now*M[p][k]; s+=1LL<<(k<<1); p^=1; if ((now*M[p][k]).v[1].count()){ now=now*M[p][k]; s+=1LL<<(k<<1); if ((now*M[p][k]).v[1].count()){ now=now*M[p][k]; s+=1LL<<(k<<1); p^=1; } } } } if (s>INF) puts("-1"); else cout << s << endl; return 0; }
以上是关于Codeforces 781D Axel and Marston in Bitland 矩阵 bitset的主要内容,如果未能解决你的问题,请参考以下文章