[caioj 1486及hdu 2157，利用矩阵乘法解决的经典题目八]多少条路呢？？

给定一个有向图，有n个点，m条边。求A点到B点恰好经过k条边的方案数（可走重复边）。

这道题表面上是dp，dfs什么的，但实际上是一道矩阵乘法加快速幂的裸题。为什么呢，其实是因为矩阵乘法它那相乘的特殊。首先把有向图转换成邻接矩阵X（可达矩阵），即如果x点能到达y点，那么a[x][y]=1，如果x,y两点之间不能到达，那么a[x][y]=0。假设A矩阵为X的平方，那么A[i][j]=ΣX[i][k]*X[k][j]，这其实就是乘法原理与加法原理的结合，那么这个A矩阵的A[i][j]就表示着i到j经过了2条边的方案数，懂了吗。那么经过k条边，只需用快速幂将X平方k次就可以了，这题就解决了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
struct node
{
    int a[25][25];
    node()
    {
        memset(a,0,sizeof(a));
    }
};
int n;
node chengfa(node a,node b)
{
    node c;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            for(int k=1;k<=n;k++)
            {
                c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j])%1000;
            }
        }
    }
    return c;
}
int main()
{
    int m,t;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n==0 && m==0)break;
        node pre,ans;
        for(int i=1;i<=m;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            pre.a[x+1][y+1]=1;
        }
        for(int i=1;i<=n;i++)ans.a[i][i]=1;
        scanf("%d",&t);
        while(t--)
        {
            int x,y,k;
            scanf("%d%d%d",&x,&y,&k);
            x++;y++;
            node s=ans,ss=pre;
            while(k>0)
            {
                if(k%2==1)s=chengfa(ss,s);
                ss=chengfa(ss,ss);
                k/=2;
            }
            printf("%d\n",s.a[x][y]%1000);
        }
    }
    return 0;
}