P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)
Posted ajmddzp
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)相关的知识,希望对你有一定的参考价值。
题目简洁明了(一点都不好伐)
照例,化简题目
给一张图,每一个时间点有一些点不能走,(有周期性),求从起点第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]沼泽鳄鱼(邻接矩阵,快速幂)的主要内容,如果未能解决你的问题,请参考以下文章
bzoj 1898 [Zjoi2005]Swamp 沼泽鳄鱼——矩阵快速幂
bzoj1898 [Zjoi2005]Swamp 沼泽鳄鱼——矩阵快速幂
BZOJ 1898: [Zjoi2005]Swamp 沼泽鳄鱼 [矩阵乘法]