hdu2157：How many ways??

n<=20个点m<=100条边有向图不带权，t个询问问Ai到Bi的经过k<=20条边方案数多少。

f[i][j]--i到j的方案数，，初始化成初邻接矩阵，这样做一次就得到2条路最短路，做两次就是4条……这东西不就是矩乘吗？

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
//#include<iostream>
using namespace std;

int n,m,t;
#define maxn 23
typedef int mat[maxn][maxn];
mat a,ans;
const int mod=1000;
void copy(mat a,mat b)
{
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            a[i][j]=b[i][j];
}
void mul(mat a,mat b,mat ans)
{
    mat t;memset(t,0,sizeof(t));
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            for (int k=1;k<=n;k++)
                t[i][j]=(t[i][j]+a[i][k]*b[k][j]%mod)%mod;
    copy(ans,t);
}
void init(mat t)
{
    memset(t,0,sizeof(mat));
    for (int i=1;i<=n;i++) t[i][i]=1;
}
void pow(mat a,int b,mat ans)
{
    mat t,tmp;init(t);copy(tmp,a);
    while (b)
    {
        if (b&1) mul(t,tmp,t);
        mul(tmp,tmp,tmp);
        b>>=1;
    }
    copy(ans,t);
}
int main()
{
    while (scanf("%d%d",&n,&m),n||m)
    {
        int x,y,v;memset(a,0,sizeof(a));
        for (int i=1;i<=m;i++)
        {
            scanf("%d%d",&x,&y);
            x++;y++;
            a[x][y]=1;
        }
        scanf("%d",&t);
        for (int i=1;i<=t;i++)
        {
            scanf("%d%d%d",&x,&y,&v);
            x++;y++;
            pow(a,v,ans);
            printf("%d
",ans[x][y]);
        }
    }
    return 0;
}