做p2p网站案例/朝阳区搜索优化seosem

A为一个方阵，则Tr A表示A的迹（就是主对角线上各项的和），现要求Tr(A^k)%9973。 

Input

数据的第一行是一个T，表示有T组数据。 
每组数据的第一行有n(2 <= n <= 10)和k(2 <= k < 10^9)两个数据。接下来有n行，每行有n个数据，每个数据的范围是[0,9]，表示方阵A的内容。 

Output

对应每组数据，输出Tr(A^k)%9973。

Sample Input

2
2 2
1 0
0 1
3 99999999
1 2 3
4 5 6
7 8 9

Sample Output

2
2686

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<vector>
#include<cmath>const int maxn =1e5+5;
const int mod =9973;
typedef long long ll;
using namespace std;
struct mat
{int a[15][15];
};
int n;
mat Mul(mat a,mat b)
{mat ans;memset(ans.a,0,sizeof(ans.a));for(int t=0;t<n;t++){for(int j=0;j<n;j++){for(int k=0;k<n;k++){ans.a[t][j]=(ans.a[t][j]+a.a[t][k]*b.a[k][j])%mod;}}}return ans;
}
mat ans;
ll quickPow(ll x)
{mat res;memset(res.a,0,sizeof(res));for(int t=0;t<n;t++){res.a[t][t]=1;}while(x){if(x&1){res=Mul(ans,res);}ans=Mul(ans,ans);x>>=1;}ll ss=0;for(int t=0;t<n;t++){ss=(ss+res.a[t][t])%mod;}return ss;
}int main()
{int T;cin>>T;int k;while(T--){cin>>n>>k;for(int t=0;t<n;t++){for(int j=0;j<n;j++){scanf("%d",&ans.a[t][j]);}}ll s=quickPow(k);cout<<s<<endl;}return 0;
}