【洛谷P1962 斐波那契数列】矩阵快速幂+数学推导
来提供两个正确的做法:
- 斐波那契数列双倍项的做法(附加证明)
- 矩阵快速幂
一、双倍项做法
在偶然之中,在百度中翻到了有关于斐波那契数列的词条(传送门),那么我们可以发现一个这个规律$ \frac{F_{2n}}{F_{n}}=F_{n-1}+F_{n+1} $,那么我就想到了是不是可以用这个公式实现类似于快速幂之类的东西:power(n,m)=power(n*n,m/2) m mod 2=0 power(n,m)=power(n*n,m/2)*n m mod 2=1
快速幂这个东西,是分成偶数情况和奇数情况,所以我们只是知道偶数想的计算公式,所以我们接下来要推导一下奇数项的递归式
\[F_{2n}=F_{n}\times(F_{n-1}+F_{n+1})
\]
\[F_{2n+2}=F_{n+1}\times(F_{n}+F_{n+2})
\]
那么我们就是要从\(F_{2n}\)和\(F_{2n+2}\)推导求出\(F_{2n+1}\)
\[F_{2n+1}=F_{2n+2}-F_{2n}
\]
\[F_{2n+1}=F_{n+1}\times(F_{n}+F_{n+2})-F_{n}*(F_{n-1}+F_{n+1})
\]
\[F_{2n+1}=F_{n+1}\times F_{n}+F_{n+1}\times F_{n+2} - F_{n}\times F_{n-1}-F_{n}\times F_{n+1}
\]
\[F_{2n+1}=F_{n+1}\times F_{n+2}-F_{n}\times F_{n-1}
\]
\[F_{2n+1}=F_{n+1}\times(F_{n+1}+F_{n})-F_{n}\times(F_{n+1}-F_{n})
\]
\[F_{2n+1}={F_{n+1}}^2+{F_{n}}^2
\]
以上就是我们对于这个公式的推导
那么我们就得到了
**F[2n] = F[n+1]² - F[n-1]² = (2F[n-1] + F[n]) · F[n] **
F[2n+1] = F[n+1]² + F[n]²
那么,我们在写一个map,那么就可以不用全部都递归到底了,优化一下。
用map映射一下大数,映射到我们的答案上。
#include <bits/stdc++.h>
using namespace std;
const int Mod=1e9+7;//mod数
long long n;
map<long long,long long> ma;//搞映射
inline long long work(long long x){
if(x==1||x==0)return 1;//边界
if(ma.count(x))return ma[x];//count如果是返回1那么就是这个答案已经在map中映射过了,0就是没有
long long res=0,t=x/2;
if(x&1) res=work(t)*(work(t-1)+work(t+1))%Mod;//公式2
else res=work(t)*work(t)%Mod+work(t-1)*work(t-1)%Mod;//公式1
return ma[x]=res;
}
int main() {//主程序
cin>>n;
long long res=work(n-1)%Mod;
cout<<res<<endl;
return 0;
}
注:这个程序的复杂度是也差不多是log(n),也是非常优的解法
二、矩阵乘法解法
这个解法应该是这一道题的正解。
我是一个蒟蒻,还是只是初懂矩阵乘法的小白。
我就贴一下自己的代码,详细的题解还是看一下别的大佬的题解。
#include <bits/stdc++.h>
using namespace std;
#define mod 1000000007 //Mod数
struct Matrix{//这个是矩阵的结构体
long long ma[2][2];
};
Matrix mul(Matrix A,Matrix B)//矩阵乘法
{
Matrix C;//答案矩阵
C.ma[0][0]=C.ma[0][1]=C.ma[1][0]=C.ma[1][1]=0;//初始化
for(int i=0;i<2;i++)
{
for(int j=0;j<2;j++)
{
for(int k=0;k<2;k++)
{
C.ma[i][j]=(C.ma[i][j]+A.ma[i][k]*B.ma[k][j])%mod;
}
}
}
return C;
}
Matrix pow_mod(Matrix A,long long n)//卡苏米+矩阵乘法优化
{
Matrix B;
B.ma[0][0]=B.ma[1][1]=1;
B.ma[0][1]=B.ma[1][0]=0;
while(n) {
if(n&1) B=mul(B,A);
A=mul(A,A);
n>>=1;
}
return B;
}
int main()
{
long long n;
Matrix A;
A.ma[0][0]=1;A.ma[0][1]=1;
A.ma[1][0]=1;A.ma[1][1]=0;//初始的数组
Matrix ans=pow_mod(A,n);
printf("%lld\n",ans.ma[0][1]);//输出答案
return 0;
}
注:关于矩阵乘法加速,这个矩阵有多种写法,这个只是其中的一种,不需要纠结于矩阵是否唯一。
黎明的朝阳,会为苦难中最坚强的信念升起
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 开发中对象命名的一点思考
· .NET Core内存结构体系(Windows环境)底层原理浅谈
· C# 深度学习:对抗生成网络(GAN)训练头像生成模型
· .NET 适配 HarmonyOS 进展
· .NET 进程 stackoverflow异常后,还可以接收 TCP 连接请求吗?
· 本地部署 DeepSeek:小白也能轻松搞定!
· 基于DeepSeek R1 满血版大模型的个人知识库,回答都源自对你专属文件的深度学习。
· 在缓慢中沉淀,在挑战中重生!2024个人总结!
· 大人,时代变了! 赶快把自有业务的本地AI“模型”训练起来!
· Tinyfox 简易教程-1:Hello World!