斐波那契数列相关整理

定义：
斐波那契数列（Fibonacci Sequence）：斐波那契数列（Fibonacci Sequence），又称黄金分割数列、因数学家列昂纳多·斐波那契（Leonardoda Fibonacci）以兔子繁殖为例子而引入，故又称为“兔子数列”，指的是这样一个数列：1、1、2、3、5、8、13、21、34、……在数学上，斐波那契数列以如下被以递推的方法定义：F(1)=1，F(2)=1, F(n)=F(n-1)+F(n-2)（n>=3，n∈N*）在现代物理、准晶体结构、化学等领域，斐波纳契数列都有直接的应用，为此，美国数学会从1963年起出版了以《斐波纳契数列季刊》为名的一份数学杂志，用于专门刊载这方面的研究成果。

实现方法：
fib[1]=1 , fib[2]=1
fib[n]=fib[n-1]+fib[n-2] n>=3

  1）递归法（效率低，复杂度 O(2^n) )
  2）递推法（效率可以，复杂度 O(n) )
  3）矩阵快速幂法 （效率很高，复杂度 log(n) ) 

递归

ll fibonacci(int n)
{
    if(n==1||n==2){
        return 1;
    }
    else{
        return fibonacci(n-1)+fibonacci(n-2);
    }
}

递推

    fib[1]=1;
    fib[2]=1;
    for(int i=3;i<=n;i++){
        fib[i]=fib[i-2]+fib[i-1];
    }

快速幂

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll MOD=1e7+7;
//通常情况下会取模

struct mat
{
    ll a[4][4];
};
mat mat_mul(mat x,mat y)
{
    mat res;
    memset(res.a,0,sizeof(res.a));
    for(int i=0;i<2;i++){
        for(int j=0;j<2;j++){
            for(int k=0;k<2;k++){
                res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j]%MOD)%MOD;
            }
        }
    }
    return res;
}
ll qpow(mat ans,ll p)
{
    mat res;
    res.a[0][0]=res.a[1][1]=1;
    res.a[0][1]=res.a[1][0]=0;
    while(p){
        if(p&1){
            res=mat_mul(res,ans);
        }
        ans=mat_mul(ans,ans);
        p>>=1;
    }
    return (res.a[0][0]+res.a[0][1])%MOD;
}

int main()
{
    mat ans;
    ans.a[0][0]=ans.a[0][1]=ans.a[1][0]=1;
    ans.a[1][1]=0;
    int n;
    cin>>n;//第n项
    if(n<=2){
        printf("1\n");
    }
    else{
        printf("%lld\n",qpow(ans,n-2));
    }
    return 0;
}

一些相关性质拓展

1. 前n项Fibonacci数列和
   1：sum(n),sum[n]=1sum[n-1]+1fib[n]+0*fib[n-1],由该等式我们转化出对应矩阵一般式：
   

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll MOD;
//通常情况下会取模

struct mat
{
    ll a[4][4];
};
mat mat_mul(mat x,mat y)
{
    mat res;
    memset(res.a,0,sizeof(res.a));
    for(int i=0;i<3;i++){
        for(int j=0;j<3;j++){
            for(int k=0;k<3;k++){
                res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j]%MOD)%MOD;
            }
        }
    }
    return res;
}
ll qpow(mat ans,ll p)
{
    mat res;
    for(int i=0;i<3;i++){
        for(int j=0;j<3;j++){
            res.a[i][j]=(i==j?1:0);
        }
    }
    while(p){
        if(p&1){
            res=mat_mul(res,ans);
        }
        ans=mat_mul(ans,ans);
        p>>=1;
    }
    return (res.a[0][0]+res.a[0][1]+res.a[0][2])%MOD;
}

int main()
{
    mat ans;
    ans.a[0][2]=ans.a[1][0]=ans.a[2][0]=ans.a[2][2]=0;
    ans.a[0][0]=ans.a[0][1]=ans.a[1][1]=ans.a[1][2]=ans.a[2][1]=1;
    int n,m;
    scanf("%d%lld",&n,&MOD);
    if(n>=2){
        printf("%lld\n",qpow(ans,n-1));
    }
    else{
        printf("1\n");
    }
    return 0;
}

  2：sum[n]=fib[n+2]-1,证明如下：
  sum[n]=fib[1]+fib[2]+fib[3]+.....+fib[n],等号两边同时加上fib[2];
  fib[2]+sum[n]=fib[2]+fib[1]+fib[2]+fib[3]+....+fib[n];
  等式右边进行化简 得 fib[n+2]
  移项得 sum[n]=fib[n+2]-fib[2]=fib[n+2]-1;证毕
  代码同上类似

2）前n项feibonacci的平方和 ,看图更易理解，代码同上类似

3） fib[1]=A,fib[2]=B,fib[i]=fib[i-1]+fib[i-2] (i>2)
gcd（fib[N],fib[N+1]）=gcd（A，B）