春天到了, HDU校园里开满了花, 姹紫嫣红, 非常美丽. 葱头是个爱花的人, 看着校花校草竞相开放, 漫步校园, 心情也变得舒畅. 为了多看看这迷人的校园, 葱头决定, 每次上课都走不同的路线去教室, 但是由于时间问题, 每次只能经过k个地方, 比方说, 这次葱头决定经过2个地方, 那他可以先去问鼎广场看看喷泉, 再去教室, 也可以先到体育场跑几圈, 再到教室. 他非常想知道, 从A 点恰好经过k个点到达B点的方案数, 当然这个数有可能非常大, 所以你只要输出它模上1000的余数就可以了. 你能帮帮他么?? 你可决定了葱头一天能看多少校花哦
Input
输入数据有多组, 每组的第一行是2个整数 n, m(0 < n <= 20, m <= 100) 表示校园内共有n个点, 为了方便起见, 点从0到n-1编号,接着有m行, 每行有两个整数 s, t (0<=s,t<n) 表示从s点能到t点, 注意图是有向的.接着的一行是两个整数T,表示有T组询问(1<=T<=100),
接下来的T行, 每行有三个整数 A, B, k, 表示问你从A 点到 B点恰好经过k个点的方案数 (k < 20), 可以走重复边。如果不存在这样的走法, 则输出0
当n, m都为0的时候输入结束
接下来的T行, 每行有三个整数 A, B, k, 表示问你从A 点到 B点恰好经过k个点的方案数 (k < 20), 可以走重复边。如果不存在这样的走法, 则输出0
当n, m都为0的时候输入结束
把给定的图转为邻接矩阵,即A(i,j)=1当且仅当存在一条边i->j。令C=A*A,那么C(i,j)=ΣA(i,k)*A(k,j),实际上就等于从点i到点j恰好经过2条边的路径数(枚举k为中转点)。类似地,C*A的第i行第j列就表示从i到j经过3条边的路径数。同理,如果要求经过k步的路径数,我们只需要二分求出A^k即可。
#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <vector>
#include <set>
#include <list>
#include <queue>
#include <map>
#include <stack>
using namespace std;
#define L(i) i<<1
#define R(i) i<<1|1
#define INF 0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-9
#define maxn 100010
#define MOD 1000000007
int N,n,m;
struct Matrix
{
long long a[22][22];
Matrix()
{
memset(a,0,sizeof(a));
}
};
Matrix operator*(Matrix A,Matrix B)
{
Matrix ans;
for(int i = 0; i < N; i++)
for(int j = 0; j < N; j++)
{
ans.a[i][j] = 0;
for(int k = 0; k < N; k++)
ans.a[i][j] += (A.a[i][k] * B.a[k][j]) % 1000;
ans.a[i][j] %= 1000;
}
return ans;
}
Matrix Pow(Matrix A,int n)
{
Matrix ans;
for(int i = 0; i < N; i++)
ans.a[i][i] = 1;
while(n)
{
if(n % 2)
ans = ans * A;
n /= 2;
A = A * A;
}
return ans;
}
int main()
{
int t;
//scanf("%d",&t);
while(scanf("%d%d",&n,&m) && (n+m))
{
N = n;
Matrix A;
int x,y;
for(int i = 0; i < m; i++)
{
scanf("%d%d",&x,&y);
A.a[x][y] = 1;
}
scanf("%d",&t);
while(t--)
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
Matrix B = Pow(A,c);
printf("%lld\n",B.a[a][b]);
}
}
}