Problem Description
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的内容。
每组数据的第一行有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<stdio.h> #include<string.h> #include<iostream> using namespace std; int c[20][20] ,ans[20][20],ans1[20][20],n,cc[20][20],b[20][20]; void fastpow(int B)// 矩阵快速幂 { int i ,j,l,sum ; while(B) { if(B%2==1) { for( l=1 ; l<=n ; l++ ) { for( i = 1; i<=n; i++) { sum=0; for( j=1; j<=n; j++) { sum=sum+(ans[l][j]*cc[j][i])%9973; } ans1[l][i]=sum%9973; } } for(i=1; i<=n; i++) for(j=1; j<=n; j++) { ans[i][j]=ans1[i][j]; } } for( l=1 ; l<=n ; l++ ) { for( i = 1; i<=n; i++) { int sum=0; for( j=1; j<=n; j++) { sum=sum+(cc[l][j]*cc[j][i])%9973; } b[l][i]=sum%9973; } } for(i=1; i<=n; i++) for(j=1; j<=n; j++) cc[i][j]=b[i][j]; B=B>>1; } } int main() { int a,sum; int t,k,i,j; scanf("%d",&t); while(t--) { scanf("%d%d",&n,&k); sum=0; for( i = 1 ; i<=n ; i++) for(j=1 ; j<=n ; j++) { scanf("%d",&c[i][j]); cc[i][j]=c[i][j]; // cout<<cc[i][j]; } for(i=1; i<=20; i++) for(j=1;j<=20;j++) //初始asn为单位矩阵 { if(i==j) ans[i][j]=1; else ans[i][j]=0; } fastpow(k); for( i = 1; i<=n; i++) { for( j=1; j<=n; j++) { if(i==j) sum+=ans[i][j]%9973; } } cout<<sum%9973<<endl; } return 0; }