解题思路：我的第一题矩阵快速幂，一定要敲熟模板！！！
AC代码如下：
#include #include #include #include #include #include #include #include #include #include #include using namespace std; #define MAX 11 typedef struct { int a[MAX][MAX]; }Det; int ksm(int a,int b,int c) { int ans=1; a%=c; while(b) { if(b&1) ans=(ans*a)%c; b>>=1; a=(a*a)%c; } return ans; } Det mul(Det a,Det b,int n,int MOD) { Det c; int i,j,k; memset(c.a,0,sizeof(c.a)); for(i=1;i<=n;i++) for(j=1;j<=n;j++) for(k=1;k<=n;k++) { c.a[i][j]+=a.a[i][k]*b.a[k][j]; c.a[i][j]%=MOD; } return c; } Det mulksm(Det a,int n,int m,int MOD) { Det b; int i,j,k; memset(b.a,0,sizeof(b.a)); for(i=1;i<=n;i++) b.a[i][i]=1; while(m) { if(m&1) b=mul(b,a,n,MOD); m>>=1; a=mul(a,a,n,MOD); } return b; } int main() { int n,i,j,k,t,sum; Det a; cin>>t; while(t--) { cin>>n>>k; sum=0; for(i=1;i<=n;i++) for(j=1;j<=n;j++) cin>>a.a[i][j]; a=mulksm(a,n,k,9973); for(i=1;i<=n;i++) sum+=a.a[i][i]; cout<
