hdu 2855 Fibonacci Check-up 矩阵+斐波那契通项公式

最新推荐文章于 2022-06-24 12:37:59 发布

knownothing

最新推荐文章于 2022-06-24 12:37:59 发布

阅读量1.6k

点赞数

分类专栏：矩阵数论* 文章标签：矩阵斐波那契数列

本文链接：https://blog.csdn.net/a601025382s/article/details/10195939

版权

数论* 同时被 2 个专栏收录

113 篇文章

订阅专栏

矩阵

38 篇文章

订阅专栏

本文介绍了一种利用矩阵快速幂算法高效计算斐波那契数列的方法，并通过具体的C++代码实现展示了该算法的过程。适用于解决斐波那契数列中大数值计算的问题。

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

#include <cstdio>
#include <cstring>
#include <cmath>
#include <vector>
#include <algorithm>
#include <iostream>
using namespace std;
#define LL __int64
int mod;
struct matrix{
    int f[2][2];
};
matrix mul(matrix a,matrix b)
{
    matrix c;
    int i,j,k;
    memset(c.f,0,sizeof(c.f));
    for(k=0;k<2;k++)
    {
        for(i=0;i<2;i++)
        {
            if(!a.f[i][k])continue;
            for(j=0;j<2;j++)
            {
                if(!b.f[k][j])continue;
                c.f[i][j]=(c.f[i][j]+a.f[i][k]*b.f[k][j])%mod;
            }
        }
    }
    return c;
}
matrix pow_mod(matrix a,int b)
{
    matrix s;
    memset(s.f,0,sizeof(s.f));
    for(int i=0;i<2;i++)
        s.f[i][i]=1;
    while(b)
    {
        if(b&1)
            s=mul(s,a);
        a=mul(a,a);
        b=b>>1;
    }
    return s;
}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        int n;
        cin>>n>>mod;
        if(n==0){cout<<0<<endl;continue;}
        n=n*2;
        matrix e;
        e.f[0][0]=e.f[0][1]=e.f[1][0]=1;e.f[1][1]=0;
        e=pow_mod(e,n);
        cout<<e.f[0][1]<<endl;
    }
    return 0;
}
/*
    矩阵:
    斐波那契数列
    f[n]=f[n-1]+f[n-2];
    |f[n+1] f[n]|=|f[n] f[n-1]|*|1 1|
                                |1 0|
     方法：
     通项公式：f[n]=(1/sqrt(5))*(((1+sqrt(5))/2)^n-((1-sqrt(5))/2)^n)
     二项式：(1+a)^n=∑(C(n,k)*a^k)(0<=k<=n)
     
     ∑C(n,k)*f[k]
     =(1/sqrt(5))*∑C(n,k)*(((1+sqrt(5))/2)^k-((1-sqrt(5))/2)^k)
     =(1/sqrt(5))*(∑C(n,k)*((1+sqrt(5))/2)^k-∑C(n,k)*((1-sqrt(5))/2)^k)
     =(1/sqrt(5))*((3+sqrt(5))/2)^n-((3-sqrt(5))/2)^n)
     =(1/sqrt(5))*((1+sqrt(5))/2)^(2*n)-((1-sqrt(5))/2)^(2*n))
     =f[2*n]
    
     还可化开下。
     又g[n]=∑C(n,k)*f[k];
     g[n]=[2*n]=f[2*n-1]+f[2*n-2]=2*f[2*n-2]+f[2*n-3]=3*f[2*n-2]-f[2*n-2]+f[2*n-3]
                =3*f[2*n-2]-f[2*n-4]=3*g[n-1]-g[n-2];
    矩阵：
        |g[n+1] g[n]|=|g[n] g[n-1]|*| 3 1|(g[0]=0,g[1]=3)
                                    |-1 0|
    结果:0 1 3 8 21 55...如果你思维够活跃，也可以直接看出来g[n]=3*g[n-1]-g[n-2]的吧。。
*/