[ZJOI 2005] 沼泽鳄鱼

题目传送-Luogu2579

题目传送-BZOJ1898

题意：

有一个无向图。
有许多条鱼在若干个点之间循环出现，循环节长度小于等于\(4\)
要求从\(A\)花费\(K\)个单位时间到达\(B\),中途不能和鱼相碰,问有多少方案。
\(n\le50,K \le 2000000000\)

题解：

这个循环节长度太小了，它们的\(lcm\)只有\(12\)
意味着只有12张不同的图
我们以12为一单位跑矩阵快速幂就行了
最后几幅图暴力

过程：

乘法的时候一不小心乘反了

代码：

const int N=60,U=12;
const int P=10000;
int n,m,s,S,T,K;
struct MAT {
	int mat[N][N];
	inline void clear() {
		mem(mat,0);
	}
	inline void print() {
		for(int i=1;i<=n;i++) {
			for(int j=1;j<=n;j++)
				printf("%d ",mat[i][j]);
			puts("");
		}puts("");
	}
	inline MAT operator * (const MAT &a)const {
		MAT c; c.clear();
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				for(int k=1;k<=n;k++) {
					c.mat[i][j]+=mat[i][k]*a.mat[k][j]%P;
					if(c.mat[i][j]>=P) c.mat[i][j]-=P;
				}
		return c;
	}
}unit,u,ma[12],st,trans,ans;
inline MAT Pow(MAT x,int y) {
	MAT ret=unit;
	for(;y;y>>=1,x=x*x)
		if(y&1) ret=ret*x;
	return ret;
}
signed main() {
	// freopen("2.in","r",stdin);
	// freopen("my.out","w",stdout);
	read(n); read(m); read(S); read(T); read(K); ++S; ++T;
	for(int i=1;i<=n;i++) unit.mat[i][i]=1;
	for(int i=1;i<=m;i++) {
		int x,y; read(x); read(y); ++x; ++y;
		u.mat[x][y]=u.mat[y][x]=1;
	}
	for(int i=0;i<12;i++) ma[i]=u;
	read(s);
	for(int i=1;i<=s;i++) {
		int loop; read(loop);
		for(int j=0;j<loop;j++) {
			int x; read(x); ++x;
			for(int k=j;k<12;k+=loop) {
				for(int p=1;p<=n;p++)
					ma[k].mat[p][x]=0;
			}
		}
	}
	// printf("%d\n",(ma[1]*ma[2]*ma[3]).mat[S][T]);
	// for(int i=0;i<4;i++) ma[i].print();
	for(int i=2;i<12;i++) ma[i]=ma[i-1]*ma[i]; ma[0]=ma[11]*ma[0];
	int t=K/12,mod=K%12;
	if(mod==0) st=unit;
	else st=ma[mod];
	trans=Pow(ma[0],t);
	// ma[1].print();
	// st.print(); trans.print();
	ans=trans*st;
	printf("%d\n",ans.mat[S][T]);
	return 0;
}
/*
20 19 1 15 20
3 14
17 16
7 10
0 2
17 19
12 15
5 7
19 8
19 12
8 0
14 10
14 13
17 6
1 3
3 15
11 17
9 15
9 0
16 15
15
3 12 15 18
4 9 10 11 12
3 16 13 6
2 19 16
3 15 10 5
4 2 0 15 11
2 0 8
3 9 5 0
3 0 4 8
2 8 4
3 7 16 6
2 3 7
4 7 19 8 2
4 3 5 12 18
2 9 13
*/

用时：1h