Description
HH有个一成不变的习惯,喜欢饭后百步走。所谓百步走,就是散步,就是在一定的时间
内,走过一定的距离。
但是同时HH又是个喜欢变化的人,所以他不会立刻沿着刚刚走来的路走回。
又因为HH是个喜欢变化的人,所以他每天走过的路径都不完全一样,他想知道他究竟有多
少种散步的方法。
现在给你学校的地图(假设每条路的长度都是一样的都是1),问长度为t,从给定地
点A走到给定地点B共有多少条符合条件的路径。
Solution
经典的矩阵乘法路径数问题
但是加了一个限制,就是不能走刚走过的路。
原来的邻接矩阵就不能用了
先将边拆成有向的
考虑将点与点间转移转化成边与边转移
将原先邻接矩阵
a[i][j]
转化成
b[i][j]
表示从第
i
条边的起点走到第
然后原来从一条无向边之间拆出来的两条有向边之间不连即可。
Code
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define N 22
#define M 122
#define mo 45989
using namespace std;
int n,m,t,st,ed,al[2*M][2],fr[N][2*M];
struct node
{
int a[M][M];
};
node ti(node a,node b)
{
node c;
fo(i,1,2*m)
fo(j,1,2*m)
{
c.a[i][j]=0;
fo(k,1,2*m)
c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j])%mo;
}
return c;
}
node s;
node ksm(node k,int n)
{
if(n==1) return k;
node s=ksm(k,n/2);
return (n%2)?ti(s,ti(s,k)):ti(s,s);
}
int main()
{
cin>>n>>m>>t>>st>>ed;
fo(i,1,m)
{
int x,y;
scanf("%d%d",&x,&y);
al[2*i-1][0]=al[2*i][1]=x;
al[2*i-1][1]=al[2*i][0]=y;
fr[x][++fr[x][0]]=2*i-1;
fr[y][++fr[y][0]]=2*i;
}
fo(i,1,2*m)
{
int x=al[i][0],y=al[i][1];
fo(j,1,fr[y][0])
{
if(abs(fr[y][j]/2-i/2)==1&&abs(fr[y][j]-i)==1) continue;
s.a[i][fr[y][j]]=1;
}
}
s=ksm(s,t-1);
int ans=0;
fo(i,1,2*m)
{
fo(j,1,2*m)
{
if(al[i][0]==st&&al[j][1]==ed) ans=(ans+s.a[i][j])%mo;
}
}
cout<<ans;
}