HDU 2157

http://acm.hdu.edu.cn/showproblem.php?pid=2157

求A到B经过K个点的方案数

http://www.matrix67.com/blog/archives/276  这里面的经典问题8

存图的邻接矩阵自乘k次，得到的新矩阵A行B列就是答案

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std ;

#define Matr 105 //矩阵大小,注意能小就小  

struct mat//矩阵结构体，a表示内容，r行c列 矩阵从1开始
{
    int a[Matr][Matr];
    int r,c;
    mat()
    {
        r=c=;
        memset(a,,sizeof(a));
    }
};
void print(mat m)
{
    int i,j;
    //printf("%d\n",m.size);
    for(i=;i<m.r;i++)
    {
        for(j=;j<m.c;j++)printf("%d ",m.a[i][j]);
        printf("\n");
    }
}  

mat mul(mat m1,mat m2,int mod)
{
    mat ans=mat();ans.r=m1.r;ans.c=m2.c;
    for(int i=;i<=m1.r;i++)
        for(int j=;j<=m2.r;j++)
            if(m1.a[i][j])
                for(int k=;k<=m2.c;k++)
                    ans.a[i][k]=(ans.a[i][k]+m1.a[i][j]*m2.a[j][k])%mod;
    return ans;
}
mat quickmul(mat m,int n,int mod)
{
    mat ans=mat();
    int i;
    for(i=;i<=m.r;i++)ans.a[i][i]=;
    ans.r=m.r;ans.c=m.c;
    while(n)
    {
        if(n&)ans=mul(m,ans,mod);
        m=mul(m,m,mod);
        n>>=;
    }
    return ans;
}
/*
ans^=n ->
mat ans=mat();
ans.r=R;ans.c=C;
初始化ans矩阵
ans=quickmul(ans,n,mod);
*/  

int main()
{
    int n,m ;
    while(~scanf("%d%d",&n,&m))
    {
        if(!n && !m)break ;
        mat M=mat() ;
        M.r=n ;M.c=n ;
        while(m--)
        {
            int a,b ;
            scanf("%d%d",&a,&b) ;
            M.a[a+][b+]= ;
        }
        int T ;
        scanf("%d",&T) ;
        while(T--)
        {
            int a,b,k ;
            scanf("%d%d%d",&a,&b,&k) ;
            mat res=quickmul(M,k,) ;
            printf("%d\n",res.a[a+][b+]) ;
        }
    }
    return  ;
}