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=;
 void copy(mat a,mat b)
 {
     for (int i=;i<=n;i++)
         for (int j=;j<=n;j++)
             a[i][j]=b[i][j];
 }
 void mul(mat a,mat b,mat ans)
 {
     mat t;memset(t,,sizeof(t));
     for (int i=;i<=n;i++)
         for (int j=;j<=n;j++)
             for (int k=;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,,sizeof(mat));
     for (int i=;i<=n;i++) t[i][i]=;
 }
 void pow(mat a,int b,mat ans)
 {
     mat t,tmp;init(t);copy(tmp,a);
     while (b)
     {
         if (b&) mul(t,tmp,t);
         mul(tmp,tmp,tmp);
         b>>=;
     }
     copy(ans,t);
 }
 int main()
 {
     while (scanf("%d%d",&n,&m),n||m)
     {
         int x,y,v;memset(a,,sizeof(a));
         for (int i=;i<=m;i++)
         {
             scanf("%d%d",&x,&y);
             x++;y++;
             a[x][y]=;
         }
         scanf("%d",&t);
         for (int i=;i<=t;i++)
         {
             scanf("%d%d%d",&x,&y,&v);
             x++;y++;
             pow(a,v,ans);
             printf("%d\n",ans[x][y]);
         }
     }
     return ;
 }