hdu 2157 How many ways??

这是一道经典的矩阵快速幂的题（可参考acm题集方向中的矩阵快速幂第8条），只需套一下模板，对应求一下数据就可以了，矩阵快速幂有点像disjkstra或者FLOYD 懂模板会用解题很方便（不太像万恶的背包或者说dp hh）

AC：

#include<iostream>  
#include<cstdio>  
#include<cstdlib>  
#include<cstring>  
#include<string>    
#include<algorithm>      
using namespace std;    
const int MAX=20+10;  
const int mod=1000;  
int f[MAX][MAX],sum[MAX][MAX],array[MAX][MAX];  
int n,m;  
  
void MatrixMult(int a[MAX][MAX],int b[MAX][MAX]){//矩阵相乘 
    int c[MAX][MAX]={0};  
    for(int i=0;i<n;++i){  
        for(int j=0;j<n;++j){  
            for(int k=0;k<n;++k){  
                c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod;  
            }  
        }  
    }  
    for(int i=0;i<n;++i){  
        for(int j=0;j<n;++j)a[i][j]=c[i][j];  
    }  
}  
  
int Matrix(int a,int b,int k){//矩阵快速幂   
    for(int i=0;i<n;++i){//初始化矩阵sum,使sum*b=b(b是矩阵)   
        for(int j=0;j<n;++j){  
            if(i == j)sum[i][j]=1;  
            else sum[i][j]=0;  
        }  
    }  
    for(int i=0;i<n;++i){//将最初矩阵f复制给array   
        for(int j=0;j<n;++j)array[i][j]=f[i][j];  
    }  
    while(k){  
        if(k&1)MatrixMult(sum,array);  
        MatrixMult(array,array);  
        k>>=1;  
    }  
    return sum[a][b];  
}
int main()
{
	int s,t,T,a,b,k;
	while(scanf("%d %d",&n,&m)==2&&(n||m))
	{
		memset(f,0,sizeof(f));
		for(int i=0;i<m;i++)
		{
			scanf("%d %d",&s,&t);
			f[s][t]=1;
		}
		scanf("%d",&T);
		while(T--)
		{
			scanf("%d %d %d",&a,&b,&k);
			cout<<Matrix(a,b,k)%1000<<endl;
		}
	}
	return 0;
}  

其中非main函数部分是模板可以收集一下嘿嘿