HDU1575 Tr A（矩阵快速幂）

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

做题要认真认真再认真，我是小菜鸡

用了k,还定义k去做别的参数，傻子，有了模板不会做，傻子

模板链接：https://blog.csdn.net/qq_40293424/article/details/102160987

注意模运算时的运算律

 

import java.util.Arrays;
import java.util.Scanner;
class y_{
	int a[][]=new int[12][12];//定义矩阵，存储初始矩阵Tr
	public y_(){
		for(int i=0;i<12;i++)
			Arrays.fill(a[i],0);//初始化
	}
}
public class Main {
	static int T,n,k;
	static y_ I=new y_();

	static y_ quick(y_ s,y_ t){		//矩阵相乘+模运算函数
		y_ res=new y_();	
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int k=0;k<n;k++)
					res.a[i][j]=(res.a[i][j]+s.a[i][k]*t.a[k][j]%9973)%9973;//这里注意不要写成+=，否则最后应该还要模一次，一般答案大了，就是少运算了一次，反正多模运算，答案又错不了。							
		return res;
	}
	
	static y_ quickm(y_ mac,int n){		//快速幂函数
		y_ ans=new y_();
		ans.a=I.a;
		while(n>0){
			if((n&1)==1) ans=quick(ans,mac);								
			mac=quick(mac,mac);			
			n>>=1;
		}		
		return ans;
	}
	
	public static void main(String[] args) {
		// TODO Auto-generated method stub
		Scanner sc=new Scanner(System.in);
		T=sc.nextInt();
		while(T>0){
			n=sc.nextInt();
			k=sc.nextInt();
			int kk=0;
			while(kk<n){//单位矩阵初始化
				I.a[kk][kk]=1;
				kk++;
			}			
			y_ mac =new y_();
			y_ res =new y_();
			
			for(int i=0;i<n;i++)
				for(int j=0;j<n;j++)
					mac.a[i][j]=sc.nextInt();							
			res=quickm(mac,k);
			kk=0;
			int sum=0;
			while(kk<n){
				sum=(sum+res.a[kk][kk])%9973;//求和
				kk++;
			}	
			System.out.println(sum);
			T--;
		}
	}

}