codevs 1482 路线统计

【问题描述】
给出一个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;
}