题目描述
Description
N个节点的有向图, 求从start到finish刚好经过时间time的总方案数 mod 502630.
输入描述
Input Description
第一行包含一个整数n, 所有点是从0到N-1编号.
接下来n行,每行包含n个字符. 第i行第j个字符表示i到j需要的时间. 字符只可能是’1’到’5’, 或者是’.’表示i不能到达j. 保证主对角线都是’.’.
接下来一行3个整数start, finish, time.
输出描述
Output Description
输出总方案数.
样例输入
Sample Input
3
.12
2.1
12.
0 2 5
样例输出
Sample Output
8
数据范围及提示
Data Size & Hint
对于20%的数据, 输入的字符不是’1’就是’.’;
对于100%的数据, 1 <= n <= 10; 1 <= start,finish <= n; 1 <= time <= 10^9.
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
/* 矩阵乘法快速幂 本题路径长度较小,所以将路径长度为n的边拓展成n条路径长度为1的边来做 */ #include<cstdio> #include<iostream> #include<cstring> #define M 60 #define mod 502630 using namespace std; struct node { int map[M][M]; node() { memset(map,0,sizeof(map)); } };node yong; int n,s,t,step; node jv(node a,node b) { node c; for(int k=1;k<=n*5;k++) for(int i=1;i<=n*5;i++) for(int j=1;j<=n*5;j++) c.map[i][j]+=((long long)a.map[i][k]*b.map[k][j]%mod)%mod, c.map[i][j]%=mod; return c; } void work() { node ans; for(int i=1;i<=5*n;i++) ans.map[i][i]=1; while(step) { if(step&1)ans=jv(ans,yong); yong=jv(yong,yong); step/=2; } printf("%d",ans.map[s+1][t+1]); } int main() { scanf("%d",&n); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { char c; cin>>c; if(c!='.') { int x=c-'0'; for(int k=1;k<x;k++) yong.map[(k-1)*n+i][k*n+i]=1; yong.map[(x-1)*n+i][j]=1; } } scanf("%d%d%d",&s,&t,&step); work(); return 0; }