标签：sizeof int ble memset 快速 set 注意 eof for
题目简洁明了（一点都不好伐）
照例，化简题目
给一张图，每一个时间点有一些点不能走，（有周期性），求从起点第k秒恰好在终点的方案数，可重复，不可停留。
额dp实锤
于是就被打脸了....
有一种东西叫做邻接矩阵，还有一种东西叫做矩阵乘法，bk201大仙曾经讲过，能用邻接矩阵的k次方求这个东西。
那，难度下降了很多了。
但是，对于那周期是3,4,6的鳄鱼怎么办呢？
答案就是：
观察3,4,6，lcm（3,4,6）=12，所以，当k%12==0，鳄鱼们都回到了原状态。
所以，我们就暴力地对于每一种情况建一个邻接矩阵，把不能走的点刨掉，然后把他们依次乘起来（因为不满足交换律）
然后把它给k次方就行了。。。
值得注意的地方：
1、因为点是从0开始，所以根据我的习惯，都加上了1；
2、不满足交换律不满足交换律不满足交换律
3、k不一定整除12，所以在剩余的那一部分要手动乘一下初始矩阵
代码：
#include<iostream> #include<cstdio> #include<cstring> using namespace std; const int maxn=55; const int mod=10000; int m,n,st,en,K,fn; int map[maxn][maxn]; int unable[maxn][maxn]; struct node { int a[maxn][maxn]; node() { memset(a,0,sizeof(a)); } }q,w[55]; inline node operator * (const node& x,const node& y) { node z; for(int i=1;i<=n;i++) { for(int k=1;k<=n;k++) { for(int j=1;j<=n;j++) { z.a[i][j]=(z.a[i][j]+x.a[i][k]*y.a[k][j])%mod; } } } return z; } node ksm(node x,int y) { node z; for(int i=1;i<=n;i++) z.a[i][i]=1; while(y) { if(y&1) z=z*x; x=x*x; y>>=1; } return z; } int main() { scanf("%d%d%d%d%d",&n,&m,&st,&en,&K); st++;en++; for(int i=1;i<=m;i++) { int x,y; scanf("%d%d",&x,&y); x++;y++; map[x][y]=map[y][x]=1; } scanf("%d",&fn); while(fn--) { int x; scanf("%d",&x); for(int j=0;j<x;j++) { int y; scanf("%d",&y); y++; for(int k=j;k<12;k+=x) { unable[k][y]=true; } } } for(int i=1;i<=n;i++) { q.a[i][i]=1; } for(int i=0;i<12;i++) { for(int j=1;j<=n;j++) { for(int k=1;k<=n;k++) { w[i].a[j][k]=(map[j][k]&&!unable[i][k]); } } } for(int i=1;i<12;i++) { q=q*w[i]; } q=q*w[0]; int k=K/12,t=K-k*12; node ans; for(int i=1;i<=n;i++) { ans.a[i][i]=1; } ans=ans*ksm(q,k); for(int i=1;i<=t;i++) ans=ans*w[i]; printf("%d ",ans.a[st][en]); return 0; }
（完）
P2579 [ZJOI2005]沼泽鳄鱼（邻接矩阵，快速幂）
标签：sizeof int ble memset 快速 set 注意 eof for
原文地址：https://www.cnblogs.com/ajmddzp/p/11374728.html
