又一道矩阵乘法加速的题目。看到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);
}