hdu 2157 How many ways?? （ 经典矩阵 ）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2157

题解链接：http://blog.csdn.net/rowanhaoa/article/details/21134509

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<stdlib.h>

using namespace std;

const int mod=1e3;

typedef struct matrix
{
    int  tr[25][25];
    matrix (){ memset(tr,0,sizeof(tr)); }
}S;
int n,m,map[25][25];

S mul( S A, S B)
{
    S C;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            for(int k=1;k<=n;k++)
            C.tr[i][j]=(C.tr[i][j]+A.tr[i][k]*B.tr[k][j])%mod;
    return C;
}

S mulpow( S A, int k )
{
    S B;
    if( k == 1 || k== 0 ) return A;
    B = mulpow( A, k/2 );
    if( k%2 == 0 )
        return mul( B, B );

    return mul( A, mul( B, B ));
}
//
//S mulpow( S A, int k )
//{
//    S B;
//    for(int i=1;i<=n;i++)
//        B.tr[i][i]=1;
//    while( k )
//    {
//        if( k&1 ) B = mul( B, A );
//        A = mul( A, A );
//        k >>=1;
//    }
//    return B;
//}

int main( )
{
    int a,b,c;
    while(scanf("%d%d",&n,&m),n||m)
    {
    S t,t1;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&a,&b);
            t.tr[a+1][b+1]=1;
        }
        scanf("%d",&m);
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d%d",&a,&b,&c);
            if( a==b && c==0 )
            {
                printf("1
");
                continue;
            }
            if( a!=b && c==0 )
            {
                printf("0
");
                continue;
            }
            t1 = mulpow( t, c );
            printf("%d
",(t1.tr[a+1][b+1])%mod);
        }
    }
    return 0;
}