HDU1575

Keep_Trying_Go

于 2021-03-01 12:30:55 发布

阅读量150

点赞数 1

分类专栏：矩阵连乘文章标签：算法

本文链接：https://blog.csdn.net/Keep_Trying_Go/article/details/114260107

版权

矩阵连乘专栏收录该内容

8 篇文章 0 订阅

订阅专栏

矩阵快速幂的使用（模版）：

#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const ll mod=9973;
const ll maxx=1e2+7;
ll ax,bx;
struct node{
	ll p[11][11];
}ans,res;
node mul(node a,node b){
	node temp;
	for(int i=0;i<ax;i++){
		for(int j=0;j<ax;j++){
			temp.p[i][j]=0;
		}
	}
	for(int i=0;i<ax;i++){
		for(int j=0;j<ax;j++){
			for(int k=0;k<ax;k++){
				temp.p[i][j]+=a.p[i][k]*b.p[k][j]%mod;
				temp.p[i][j]%=mod;
			}
		}
	}
	return temp;
}
void quick_mod(ll b){
	for(int i=0;i<ax;i++){
		for(int j=0;j<ax;j++){
			if(i==j){
				ans.p[i][j]=1;
			}else{
				ans.p[i][j]=0;
			}
		}
	}
	while(b){
		if(b&1){
			ans=mul(ans,res);
		}
		res=mul(res,res);
		b>>=1;
	}
}
int main(){
	int t;
	cin>>t;
	node tx,ty;
	while(t--){
		cin>>ax>>bx;
		for(int i=0;i<ax;i++){
			for(int j=0;j<ax;j++){
				cin>>res.p[i][j];
			}
		}
		quick_mod(bx);
		ll sum=0;
		for(int i=0;i<ax;i++){
			sum+=ans.p[i][i];
		} 
		cout<<sum%mod<<endl;
	}
	return 0;
}