Description
HH有个一成不变的习惯,喜欢饭后百步走。所谓百步走,就是散步,就是在一定的时间 内,走过一定的距离。 但
是同时HH又是个喜欢变化的人,所以他不会立刻沿着刚刚走来的路走回。 又因为HH是个喜欢变化的人,所以他每
天走过的路径都不完全一样,他想知道他究竟有多 少种散步的方法。 现在给你学校的地图(假设每条路的长度都
是一样的都是1),问长度为t,从给定地 点A走到给定地点B共有多少条符合条件的路径
数据保证Ai != Bi,但不保证任意两个路口之间至多只有一条路相连接。
路口编号从0到N -1。
同一行内所有数据均由一个空格隔开,行首行尾没有多余空格。没有多余空行。
答案模45989。
N ≤ 20,M ≤ 60,t ≤ 2^30,0 ≤ A,B
来自 https://www.lydsy.com/JudgeOnline/problem.php?id=1875
Solution
仔细审题发现并不能来回走,因此直接上矩阵乘法是会wa的
考虑dp,设f[I,j]表示走了i步到达了第j条边的终点的方案数,那么f[I+1,j]=Σf[I,k] (满足e[k].y=e[j].x)
第一维可以去掉,那么转移矩阵就很显然了。注意边是双向的
Code
#include <stdio.h>
#include <string.h>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))
const int MOD=45989;
const int N=125;
int rc[N][N],f[N];
int n,m,t;
int ls[N],edCnt=1;
struct edge {int x,y,next;} e[N];
struct mat {
int a[N][N];
mat operator *(mat b) {
mat c; fill(c.a,0);
rep(i,1,edCnt) rep(j,1,edCnt) {
rep(k,1,edCnt) c.a[i][j]=(c.a[i][j]+a[i][k]*b.a[k][j]%MOD)%MOD;
}
return c;
}
} A,ret;
void add_edge(int x,int y) {
e[++edCnt]=(edge) {x,y,ls[x]}; ls[x]=edCnt;
e[++edCnt]=(edge) {y,x,ls[y]}; ls[y]=edCnt;
}
mat ksm(mat x,int dep) {
ret=x; dep-=1;
while (dep) {
if (dep&1) ret=ret*x;
x=x*x; dep/=2;
}
return ret;
}
int main(void) {
int st,ed;
scanf("%d%d%d%d%d",&n,&m,&t,&st,&ed);
st++; ed++;
rep(i,1,m) {
int x,y; scanf("%d%d",&x,&y);
x++; y++;
add_edge(x,y);
}
rep(i,2,edCnt) rep(j,2,edCnt) {
if (e[i].x==e[j].y&&i!=(j^1)) A.a[j][i]=1;
}
A=ksm(A,t-1);
for (int i=ls[st];i;i=e[i].next) f[i]=1;
int ans=0;
rep(i,1,edCnt) if (e[i].y==ed) {
rep(j,1,edCnt) ans=(ans+f[j]*A.a[j][i]%MOD)%MOD;
}
printf("%d\n", ans);
return 0;
}