hdu1575矩阵快速幂模板

什么时候用：        1.有递推公式         2.n比较大

（注意：全用ll，以防万一数太大）

#include<bits/stdc++.h>
using namespace std;
const int mod=9973;
typedef long long ll;
typedef struct jz{
	ll m[13][13];
}jz;//结构体定义矩阵 作函数参数等等 好写点？
int n;
//对矩阵里面的数运算，所以要取余 
jz mul(jz a,jz b){
	jz ans;
	for(int i=1;i<=n;i++)
    	for(int j=1;j<=n;j++){
    		ans.m[i][j]=0;
			for(int k=1;k<=n;k++)
				ans.m[i][j]=(ans.m[i][j]+a.m[i][k]*b.m[k][j]%mod)%mod;
				//加起来之后整体有个mod否则相加后可能过大。相乘类似 		
		}
	return ans; 
} 
//对矩阵的运算，矩阵没法取余，所以不用  
jz ksm(jz a,ll b){
	jz ans;
	for(int i=1;i<=n;i++)
    	for(int j=1;j<=n;j++)
    		ans.m[i][j]=1*(i==j);
	while(b){
		if(b%2) ans=mul(ans,a);
		a=mul(a,a);
		b/=2;
	}
	return ans;
}
int main(){
	ll T,k,ans;
	jz A,B;
	cin>>T;
	while(T--){
		ans=0;//ans初始化为0别忘了  
		cin>>n>>k;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				scanf("%lld",&A.m[i][j]);
		B=ksm(A,k);
		for(int i=1;i<=n;i++)
			ans=(ans+B.m[i][i])%mod;
		printf("%lld\n",ans);	
	}
}