HDU 2157 How many ways??（矩阵快速幂）

中文题就不说题意了，其中 i 到 j 经过k个点的所有方案就是两个邻接矩阵相乘k-1次之后的结果。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<queue>
#include<cstring>
using namespace std;
const int mo = 1000; 
int n,m;
struct node{
	int m[31][31];
	node(){
		memset(m,0,sizeof(m));
	}
	
};

node cmp(node a, node b){
		node te;
	for(int i = 0; i < n; i++)
		for(int j = 0; j < n; j++){
		te.m[i][j]=0;
			for(int k = 0; k< n; k++){
				te.m[i][j] += a.m[i][k] * b.m[k][j]; 
				te.m[i][j] %= mo; 
			}
		} 
	return te;
}

node powermod(node base,int k,int a,int b){
		node ans;
		for(int i = 0; i <n; i++ )
			ans.m[i][i] = 1;
	    
	    while(k){
	    	if(k&1) ans = cmp(ans, base);
	    	k >>= 1;
	    	base=cmp(base, base);
		}
		cout<<ans.m[a][b]<<endl;
		return ans;
	}

 
int main(){

	int a,b,T,k;
	
	while(scanf("%d%d",&n,&m)&&(n || m)){
				node te;//这个定义在while函数外，导致wa 
			while(m--){
				scanf("%d%d",&a,&b);
				te.m[a][b] = 1;
			}
			scanf("%d",&T);
			while(T--){
				scanf("%d%d%d",&a,&b,&k);
				 powermod(te,k,a,b);		
			}
	}
		
	return 0;
}