2017江西icpc邀请赛快速幂板子

最新推荐文章于 2023-09-21 22:56:57 发布

原创最新推荐文章于 2023-09-21 22:56:57 发布 · 163 阅读

0 ·

CC 4.0 BY-SA版权

模板同时被 2 个专栏收录

10 篇文章

订阅专栏

快速幂数学题

2 篇文章

订阅专栏

该博客围绕2017江西ICPC邀请赛，提供了快速幂的代码模板。快速幂是信息技术领域常用算法，可高效计算幂运算，在竞赛等场景有重要应用。

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

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MD=1000000007;
int p[5]= {0,1,5,11,36};
struct Matrix {
	ll a[4][4];
	Matrix() {
		memset(a,0,sizeof(a));
	}
	Matrix operator*(const Matrix  y) {
		Matrix  ans;
		for(int i=0; i<4; i++)
			for(int j=0; j<4; j++)
				for(int k=0; k<4; k++)
					ans.a[i][j]=(ans.a[i][j]+a[i][k]*y.a[k][j]+MD)%MD;
		return ans;
	}
	void operator=(const Matrix b) {
		for(int i=0; i<4; i++)
			for(int j=0; j<4; j++)
				a[i][j]=b.a[i][j];
	}
};
ll mat_pow(ll n) {
	Matrix c,res;
	c.a[0][0]=c.a[0][1]=c.a[1][2]=c.a[2][3]=c.a[2][0]=1;
	c.a[3][0]=-1,c.a[1][0]=5;
	res.a[0][0]=36,res.a[0][1]=11,res.a[0][2]=5,res.a[0][3]=1;
	while(n) {
		if(n&1) res=res*c;
		c=c*c;
		n=n>>1;
	}
	//for(int i=0;i<4;i++){
	//	cout<<res.a[0][i]<<endl;
	//}
	return res.a[0][0];
}
int main() {
	ll n;
	while(~scanf("%lld",&n)) {
		if(n<=4) {
			printf("%d\n",p[n]);
		} else printf("%lld\n",mat_pow(n-4));
	}
	return 0;
}