1875: [SDOI2009]HH去散步
Time Limit: 20 Sec Memory Limit: 64 MBSubmit: 823 Solved: 373
[ Submit][ Status][ Discuss]
Description
HH有个一成不变的习惯,喜欢饭后百步走。所谓百步走,就是散步,就是在一定的时间 内,走过一定的距离。 但是同时HH又是个喜欢变化的人,所以他不会立刻沿着刚刚走来的路走回。 又因为HH是个喜欢变化的人,所以他每天走过的路径都不完全一样,他想知道他究竟有多 少种散步的方法。 现在给你学校的地图(假设每条路的长度都是一样的都是1),问长度为t,从给定地 点A走到给定地点B共有多少条符合条件的路径
Input
第一行:五个整数N,M,t,A,B。其中N表示学校里的路口的个数,M表示学校里的 路的条数,t表示HH想要散步的距离,A表示散步的出发点,而B则表示散步的终点。 接下来M行,每行一组Ai,Bi,表示从路口Ai到路口Bi有一条路。数据保证Ai = Bi,但 不保证任意两个路口之间至多只有一条路相连接。 路口编号从0到N − 1。 同一行内所有数据均由一个空格隔开,行首行尾没有多余空格。没有多余空行。 答案模45989。
Output
一行,表示答案。
Sample Input
4 5 3 0 0
0 1
0 2
0 3
2 1
3 2
0 1
0 2
0 3
2 1
3 2
Sample Output
4
HINT
对于30%的数据,N ≤ 4,M ≤ 10,t ≤ 10。
对于100%的数据,N ≤ 20,M ≤ 60,t ≤ 230,0 ≤ A,B<n,0 ≤="" ai,bi="" <n。<="" p="">
Source
矩阵乘法应用。
《十个利用矩阵乘法解决的经典题目》from matrix67
矩阵乘法可以用来统计路径:
把原图建成邻接矩阵,另C=A*A,C(i,j)=sigma(A(i,k)*A(k,j)),C(i,j)就是从i到j且经过两条边的方案数。那么经过k条边就是A^(k-1)。
但是这道题有个限制:
不会沿着刚刚走来的路立刻走回。
上面的方法是会计算到这种情况的,因此需要改进一下:
把原来点的邻接矩阵变成边的邻接矩阵,如果A(i,j)=1,表示第i条边和第j条边相连,那么就不会立刻走刚刚走的边了(因为A(i,i)一定是0)
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#define mod 45989
using namespace std;
struct matrix
{
int f[125][125];
}a,b;
int n,m,s,t;
struct edge
{
int x,y;
}e[500];
matrix Mult(matrix a,matrix b)
{
matrix ans;
for (int i=0;i<2*m;i++)
for (int j=0;j<2*m;j++)
{
ans.f[i][j]=0;
for (int k=0;k<2*m;k++)
ans.f[i][j]=(ans.f[i][j]+a.f[i][k]*b.f[k][j]%mod)%mod;
}
return ans;
}
matrix Pow(int n)
{
matrix ans;
int f=0;
while (n)
{
if (n&1)
{
if (f) ans=Mult(ans,b);
else ans=b,f=1;
}
n>>=1;
b=Mult(b,b);
}
return ans;
}
int main()
{
int l;
scanf("%d%d%d%d%d",&n,&m,&l,&s,&t);
for (int i=0;i<2*m;i+=2)
{
scanf("%d%d",&e[i].x,&e[i].y);
e[i+1].x=e[i].y,e[i+1].y=e[i].x;
}
for (int i=0;i<2*m;i++)
{
for (int j=0;j<2*m;j++)
if (j!=i&&j!=(i^1)&&e[j].x==e[i].y)
b.f[i][j]=1;
if (e[i].x==s)
a.f[0][i]=1;
}
matrix ans=Mult(a,Pow(l-1));
int x=0;
for (int i=0;i<2*m;i++)
if (e[i].y==t)
x=(x+ans.f[0][i])%mod;
cout<<x<<endl;
return 0;
}