JZOJ 2288. 【ZJOI2005】沼泽鳄鱼【矩阵乘法】

---

题目：

传送门

---

题意：

有$n$个柱子，$m$个桥连通这些柱子
还有些食人鱼来限制我们的行动，食人鱼的行动是有周期性的，而当食人鱼处于某根柱子时，我们就无法前去
现在给出起点和终点，问有多少种合法的方案

---

分析：

我们需要一个前置知识：对于一张图，$g[i][j{]}^n$表示在图中从$i$开始走$n$步到达$j$的方案
那么这题就有点显然了吧
唯一需要留意的就是鱼的周期了，我们看到，它们的周期都是$2,3,4$，也就是$12$步肯定是一个大周期
我们可以把每$12$步看成一个整体，再用快速幂进行求解
对于那些剩余的步数，我们可以单独用矩阵乘法顺次相乘解决

---

代码：

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define LL long long 
#define LZX 10000
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
int n=read(),m=read(),b=read()+1,e=read()+1,t=read();
struct QWQ{
	int f[55][55];
}ans,iv,g[15],p;
int nd[5];
QWQ operator * (QWQ a,QWQ b)
{
	QWQ tem;
	memset(tem.f,0,sizeof(tem.f));
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++)
	    for(int k=1;k<=n;k++)
	      tem.f[i][j]=(tem.f[i][j]+a.f[i][k]*b.f[k][j])%LZX;
	return tem;
}
QWQ ksm(QWQ x,int y)
{
	QWQ tem=iv;
	while(y)
	{
		if(y&1) tem=tem*x;
		x=x*x;y>>=1;
	}
	return tem;
}
int main()
{
	for(int i=1;i<=m;i++) 
	{
		int x=read()+1,y=read()+1;
		for(int j=1;j<=12;j++)
		  g[j].f[x][y]=g[j].f[y][x]=1;
	}
	int fish=read();
	while(fish--)
	{
		int z=read();
		for(int i=1;i<=z;i++) nd[i]=read()+1;
		for(int i=1;i<=n;i++)
		  for(int j=0;j<=12;j++)
		    g[j].f[i][nd[j%z+1]]=0;
	}
	for(int i=1;i<=n;i++) iv.f[i][i]=p.f[i][i]=1;
	for(int i=1;i<=12;i++) p=p*g[i];
	ans=ksm(p,t/12);
	for(int i=1;i<=t%12;i++) ans=ans*g[i];
	cout<<ans.f[b][e];
	return 0;
}