矩阵问题&&http://acm.hdu.edu.cn/showproblem.php?pid=1575

一道矩阵水题，，，，主要是二分法，，，终于1A了，，，

#include<iostream> #include<cstdio> #include<string.h> #define N 11 #define MOD 9973 using namespace std; typedef struct { int num[N][N]; }Node; Node a,b; int n,k; Node mul(Node x,Node y) { Node c; memset(c.num,0,sizeof(c.num)); for(int i=0;i<n;++i) for(int j=0;j<n;++j) { for(int t=0;t<n;++t) c.num[i][j]=(c.num[i][j]+x.num[i][t]*y.num[t][j])%MOD; } return c; } Node ceil() { Node p=a,q=b; while(k) { if(1&k) p=mul(p,q); q=mul(q,q); k=k>>1; } return p; } int main() { int Case; scanf("%d",&Case); while(Case--) { scanf("%d%d",&n,&k); for(int i=0;i<n;++i) for(int j=0;j<n;++j) scanf("%d",&b.num[i][j]),a.num[i][j]=(i==j); Node c=ceil(); int sum=0; for(int i=0;i<n;++i) sum+=c.num[i][i],sum%=MOD; printf("%d\n",sum); }return 0; }