矩阵乘法模板[P1939]

本文深入探讨了矩阵快速幂算法的实现细节,通过具体的代码示例展示了如何利用矩阵乘法和快速幂运算解决复杂的问题。文章首先定义了一个矩阵结构体,并实现了矩阵乘法和矩阵快速幂函数。随后,通过主函数中的示例输入,演示了如何应用这些函数来求解特定数学问题。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

传送门

#include<bits/stdc++.h>
#define Mod 1000000007
#define LL long long
using namespace std;
struct Matrix{
	LL x[4][4];
	Matrix(){memset(x,0,sizeof(x));}
}; 
int T,n;
Matrix mul(Matrix a,Matrix b){
	Matrix c;
	for(int i=1;i<=3;i++)
		for(int j=1;j<=3;j++)
			for(int k=1;k<=3;k++)
				c.x[i][j] = (c.x[i][j] + a.x[i][k] * b.x[k][j]) % Mod;
	return c;
}
Matrix power(Matrix a,int k){
	Matrix ans;
	for(int i=1;i<=3;i++) ans.x[i][i]=1;
	for(;k;k>>=1){
		if(k&1) ans = mul(ans,a);
		a = mul(a,a);
	}return ans; 
}
int main(){
	scanf("%d",&T); while(T--){
		scanf("%d",&n); if(n<=3){printf("1\n"); continue;}
		Matrix a,b;
		a.x[1][1]=a.x[1][2]=a.x[1][3]=1;
		b.x[1][1]=b.x[1][2]=b.x[2][3]=b.x[3][1]=1;
		b = power(b,n-3); a = mul(a,b);
		printf("%lld\n",a.x[1][1]);
	}
} 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值