AC日记——斐波那契数列 洛谷 P1962

斐波那契数列

思路：

　　矩阵快速幂；

来，上代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

#define mod 1000000007

struct MatrixType {
    long long ai[][];

    void mem()
    {
        for(int i=;i<;i++)
        {
            for(int j=;j<;j++) ai[i][j]=;
        }
    }
};

long long n;

MatrixType mul(MatrixType a,MatrixType b)
{
    MatrixType res;
    res.mem();
    for(int i=;i<=;i++)
    {
        for(int j=;j<=;j++)
        {
            for(int v=;v<=;v++) res.ai[i][j]=(res.ai[i][j]+(a.ai[i][v]*b.ai[v][j])%mod)%mod;
        }
    }
    return res;
}

long long poww(long long mi)
{
    MatrixType pos,mii;
    pos.mem(),mii.mem();
    pos.ai[][]=,pos.ai[][]=;
    pos.ai[][]=,pos.ai[][]=;
    mii.ai[][]=,mii.ai[][]=;
    mii.ai[][]=,mii.ai[][]=;
    while(mi>)
    {
        if(mi&) pos=mul(pos,mii);
        mi=mi>>,mii=mul(mii,mii);
    }
    return pos.ai[][];
}

int main()
{
    scanf("%lld",&n);
    cout<<poww(--n)%mod;
    return ;
}