标签：int YBTOJ ll ans 矩阵 波拉 res trans mod

文章目录

• 题目描述
• 题目描述
• 代码

题目描述

题目描述

关键在于如何转化为本题的题目。。。
设

y=(1-根号5）/2$$

再令：

A(n)=xⁿ + + + yⁿ

通过尝试可以发现，A其实就是一个1,3为前两项的斐波拉契数列
则

xⁿ=A(n)-yⁿ

A的值可以用矩阵快速幂来求
而y是在(-1,0)的区间
所以分一下n的奇偶性，就可以求出答案了

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int mod=1e9+7;
ll a,b,c,k;
ll ans[4][4],res[4][4],trans[4][4];
void mul1(){//res*res
	memset(trans,0,sizeof(trans));
	for(int i=1;i<=2;i++){
		for(int j=1;j<=2;j++){
			for(int p=1;p<=2;p++){
				trans[i][j]+=res[i][p]*res[p][j];
				trans[i][j]%=mod;
			}
		}
	}
	for(int i=1;i<=2;i++){
		for(int j=1;j<=2;j++){
			res[i][j]=trans[i][j];
		}
	}
}
void mul2(){//res*ans
	memset(trans,0,sizeof(trans));
	for(int i=1;i<=2;i++){
		for(int j=1;j<=2;j++){
			for(int p=1;p<=2;p++){
				trans[i][j]+=ans[i][p]*res[p][j];
				trans[i][j]%=mod;
			}
		}
	}
	for(int i=1;i<=2;i++){
		for(int j=1;j<=2;j++){
			ans[i][j]=trans[i][j];
		}
	}
}
void ksm(ll w){
	while(w){
		if(w&1) mul2();
		mul1();
		w>>=1;
	}
	return;
}
void print(){
	for(int i=1;i<=2;i++){
		for(int j=1;j<=2;j++){
			printf("%d ",ans[i][j]);
		}
		printf("\n");
	}
}
int main()
{
	ans[1][1]=ans[2][2]=1;
	res[1][1]=res[1][2]=res[2][1]=1;
	scanf("%lld",&a);
	if(a<=2) {printf("1");return 0;}
	ksm(a-2);
//	print();
	printf("%lld",(ans[1][1]+ans[2][1])%mod);
    return 0;
}

标签：int,YBTOJ,ll,ans,矩阵,波拉,res,trans,mod
来源： https://blog.csdn.net/BUG_Creater_jie/article/details/119301502

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