bzoj1875: [SDOI2009]HH去散步

又一道矩阵乘法加速的题目。看到N ≤ 20，M ≤ 60，t ≤ 2^30的数据范围，很容易想到用矩阵乘法。本题不能立刻走回头路，所以可以以边为状态，把无向边拆成2条有向边a[i]表示目前在第i条边终点，b[i][j]表示第i条边是否接着j，先把起点为A的a[1]=1，然后a*b^t。统计终点为B的a[i]值。别忘%。。。

#include<iostream>
#include<cstdio>
#define S 200
using namespace std;
int n,m,t,A,B,a[S][S],b[S][S],u[S],map[S][S],Ans;
inline int read()
{  
    int x=0; char ch=getchar();  
    while (ch<'0' || ch>'9') ch=getchar();  
    while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }  
    return x;  
}  
void mul(int a[S][S],int b[S][S],int ans[S][S])
{
	int tmp[S][S];
	for(int i=0;i<2*m;i++)
	for(int j=0;j<2*m;j++)
	{
	    tmp[i][j]=0;
	    for(int l=0;l<2*m;l++)
			tmp[i][j]=(tmp[i][j]+a[i][l]*b[l][j]%45989)%45989;
	}
	for(int i=0;i<2*m;i++)
	for(int j=0;j<2*m;j++)
	    ans[i][j]=tmp[i][j];
}
void ksm(int a[S][S],int b[S][S],int t)
{		
	for (;t>0;mul(b,b,b),t>>=1) if (t&1) mul(a,b,a);
}
int main()
{ 
	n=read();m=read();t=read();A=read();B=read();
	if (t==0)
	{
		printf("%d\n",A==B); return 0;
	}
	for (int i=0;i<m;i++)
		u[i<<1|1]=read(),u[i<<1]=read();
	for (int i=0;i<2*m;i++)
	for (int j=0;j<2*m;j++)
	if (u[i]==u[j^1]&&i>>1!=j>>1)
		b[i][j]=1;

	for (int i=0;i<2*m;i++)
	if (u[i]==A) a[0][i^1]=1;
	ksm(a,b,t-1);
	for (int i=0;i<2*m;i++)
	if (u[i]==B)
		Ans=(Ans+a[0][i])%45989;
	printf("%d\n",Ans);
}