给出一个n个点有向图,求从s点到f点恰好经过时间t的路径总数。不能
,可以重复的走各点。
【输入数据】
第一行包含一个整数n, 所有点是从0到n-1编号.
接下来n行,每行包含n个字符.
第i行第j个字符表示i到j需要的时间,字符只可能是’1’到’5’,
或者是’.’表示i不能到达j, 保证主对角线都是’.’。
接下来一行3个整数s, f, t。
【输出数据】
输出总方案数mod 502630的值。
【样例输入】 route.in
3
.12
2.1
12.
0 2 5
【样例输出】 route.out
8
【数据范围及提示】
对于20%的数据, 输入的字符不是’1’就是’.’;
对于100%的数据, 1 <= n <= 10; 1 <= s,f <= n; 1 <= t <= 10^9
思路:
这道题关键在于拆点。由于题中的点很少,时间很长,所以可以用矩阵乘法来做。
由于每个边有权值,所以我们可以把每个点拆成5个点。例如,如果有一条从i到j权值为5的边,可以将之
拆成i->i+n->i+2*n->i+3*n->i+4*n->j。
拆了之后,就可以用矩阵乘法求出该矩阵k次自乘之后的矩阵,而矩阵中i->j的数的大小即为方案数
(本题由于从0开始标号,所以(i+1)->(j+1)的值为方案数。)
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define mod 502630
int n,s,t,k;
struct node
{
int map[100][100];
node()
{
memset(map,0,sizeof(map));
}
}sai;
void read()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("\n");
for(int j=1;j<=n;j++)
{
char x;
scanf("%c",&x);
if(x!='.')
{
int yy=x-'0';
for(int l=1;l<yy;l++)
{
sai.map[(l-1)*n+i][l*n+i]=1;
}
sai.map[(yy-1)*n+i][j]=1;
}
}
}
scanf("%d%d%d",&s,&t,&k);
}
node mut(node x,node y)
{
node z;
for(int i=1;i<=5*n;i++)
{
for(int j=1;j<=5*n;j++)
{
for(int l=1;l<=5*n;l++)
{
z.map[i][j]=(z.map[i][j]+(long long)x.map[i][l]*y.map[l][j]%mod)%mod;
}
z.map[i][j]%=mod;
}
}
return z;
}
void work()
{
node e;
for(int i=1;i<=5*n;i++)
{
e.map[i][i]=1;
}
while(k)
{
if(k&1)e=mut(e,sai);
k=k>>1;
sai=mut(sai,sai);
}
s++;
t++;
printf("%d",e.map[s][t]);
}
int main()
{
read();
work();
return 0;
}