斐波那契数列 矩阵乘法解法

博客围绕斐波那契数列求和展开,通过矩阵运算进行分析。定义了Fn和Fn+1矩阵,计算出运算矩阵A,得出Fn+1=F1*A^(n - 1),即只需乘上A矩阵的(n - 1)次方。还给出输入输出样例,并附上C++代码。

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

分析


Fn=[fn(第n项)fn+1(第n+1项)Sn(前n项和)]\begin{array}{ll} &F_{n}= &\begin{bmatrix} f_{n}(第n项) & f_{n+1}(第n+1项) & S_{n}(前n项和)\\ \end{bmatrix} \end{array}Fn=[fn(n)fn+1(n+1)Sn(n)]

Fn+1=[fn+1fn+2Sn+1]\begin{array}{ll} &F_{n+1}= &\begin{bmatrix} f_{n+1} & f_{n+2} & S_{n+1}\\ \end{bmatrix} \end{array}Fn+1=[fn+1fn+2Sn+1]
由这两个矩阵可以计算出所要进行运算的矩阵A:

A=[010111001]\begin{array}{ll} &A= &\begin{bmatrix} 0 & 1& 0\\ 1 & 1& 1\\ 0 & 0& 1\\ \end{bmatrix} \end{array}A=010110011

由此可计算出矩阵的和:Fn+1=F1∗An−1 由此可计算出矩阵的和:F_{n+1}=F_{1}*A^{n-1} Fn+1=F1An1
所以只用乘上A矩阵的(n-1)次方

输入样例

25000 100000

输出样例

84375

(计算第25000项的斐波那契数列的总和对100000取模后的结果)

C++ 代码
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
int m,n;
LL a[3][3]={    //预定义A矩阵
  {0,1,0},
  {1,1,1},
  {0,0,1}
};
void mul1(LL c[3],LL a[3],LL b[3][3])   //计算F矩阵乘A矩阵
{
    LL temp[3]={0};
    for(int i=0;i<3;i++)
        for(int j=0;j<3;j++)
            temp[i]=(temp[i]+a[j]*b[j][i])%m;
            
    memcpy(c,temp,sizeof temp);
}
void mul2(LL c[][3],LL a[][3],LL b[][3])    //计算A矩阵乘A矩阵
{
    LL temp[3][3]={0};
    for(int i=0;i<3;i++)
        for(int j=0;j<3;j++)
            for(int k=0;k<3;k++)
                temp[i][j]=(temp[i][j]+a[i][k]*b[k][j])%m;  

    memcpy(c,temp,sizeof temp);
}
int main()
{
    LL f[3]={1,1,1};
    cin>>n>>m;
    n--;    //计算第n项的和只需要乘n-1次A矩阵,所以n--
    while(n)
    {
        if(n&1) mul1(f,f,a); 
        mul2(a,a,a);
        n>>=1;
    }
    cout<<f[2];
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

Jay_fearless

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

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

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

打赏作者

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

抵扣说明:

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

余额充值