标签：10 mat int 矩阵 base ans 快速 模板 matrix

矩阵快速幂与正常的整数快速幂完全一致，只不过重定义了一下乘号罢了。

但注意此代码下的ans应初始化为矩阵乘法的幺元——即单位矩阵。

联动#F.Tr A

#include<bits/stdc++.h>
using namespace std;
int n, moyn = 9973, t, s;
struct matrix{
	long long mat[10][10];
	matrix() {memset(mat,0,sizeof(mat));}
	matrix operator *(const matrix& T) const {
		matrix res;
		int tmp;
		for(int i = 0;i < 10;++i) {
			for(int k = 0;k < 10;++k) {
				tmp = mat[i][k];
				for(int j = 0;j < 10;++j) {
					res.mat[i][j] += tmp*T.mat[k][j];
					res.mat[i][j] %= moyn;
				}
			}
		}
		return res;
	}
} ans, base;
void init() {
	memset(ans.mat,0,sizeof(ans.mat));
	for(int i = 0;i < 10;++i)
		ans.mat[i][i] = 1;
	memset(base.mat,0,sizeof(base.mat));
	return;
}
void quickpow(int x) {
	while(x) {
		if(x&1) ans = ans*base;
		base = base*base;
		x >>= 1;
	}
	return;
}
int main() {
	scanf("%d",&t);
	for(int i = 1;i <= t;++i) {
		scanf("%d %d",&n,&s);
		init();
		for(int j = 0;j < n;++j)
			for(int k = 0;k < n;++k) 
				scanf("%lld",&base.mat[j][k]);
		quickpow(s);
		int res = 0;
		for(int i = 0;i < 10;++i) {
			res += ans.mat[i][i];
			res %= moyn;
		}
		printf("%d\n",res);
	}
	return 0;
}

标签：10,mat,int,矩阵,base,ans,快速,模板,matrix
来源： https://www.cnblogs.com/bikuhiku/p/Matrix_quick-pow.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。