如果没有题里的“不会立刻沿着刚刚走来的路走回”限制,那么直接矩乘计算k步的方案数
但加了这个限制,就不能以点来矩乘了,考虑边数<=60,如果以边建邻接矩阵呢??
先拆边,再把每一个边和以其终点为起点的边相连,注意不能是拆前的同一条边。
然后矩阵乘喽。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define mod 45989
using namespace std;
int a[122][122],b[122],c[122][122],d[122];
int ans,n,m,t,S,T;
bool bo;
int e=1,head[22];
struct edge{
int u,v,next;
}ed[122];
void add(int u,int v){
ed[e].u=u; ed[e].v=v;
ed[e].next=head[u];
head[u]=e++;
}
void work(int A[122][122],int B[122][122],int C[122][122]){
int D[122][122];
for(int i=1;i<=2*m;i++)
for(int j=1;j<=2*m;j++){
D[i][j]=0;
for(int k=1;k<=2*m;k++)
D[i][j]=(D[i][j]+(A[i][k]*B[k][j])%mod)%mod;
}
for(int i=1;i<=2*m;i++)
for(int j=1;j<=2*m;j++)
C[i][j]=D[i][j];
}
int main()
{
int u,v;
scanf("%d%d%d%d%d",&n,&m,&t,&S,&T); S++; T++;
for(int i=1;i<=m;i++){
scanf("%d%d",&u,&v);
u++; v++;
add(u,v); add(v,u);
}
for(int i=1;i<=2*m;i++)
for(int j=1;j<=2*m;j++){
if(ed[i].v==ed[j].u)
a[i][j]=1;
}
for(int i=1;i<=2*m;i+=2) a[i][i+1]=a[i+1][i]=0;
for(int i=1;i<=2*m;i++) c[i][i]=1;
for(int i=head[S];i;i=ed[i].next) b[i]=1;
t--;
while(t){
if(t&1) work(a,c,c);
work(a,a,a);
t>>=1;
}
for(int i=1;i<=2*m;i++)
for(int j=1;j<=2*m;j++)
d[i]=(d[i]+(b[j]*c[j][i])%mod)%mod;
for(int i=1;i<=2*m;i++)
if(ed[i].v==T) ans=(ans+d[i])%mod;
printf("%d\n",ans);
return 0;
}