P3390 【模板】矩阵快速幂

题意：斐波那契数列，但是n的范围到1e18；

思路： 矩阵快速幂裸题；

#include<cstdio>
#include<algorithm>
#include<math.h>
#include<string.h>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
struct node
{
    ll a[10][10];
}ans,A,B;
node mat(node x,node y)
{
    node c;
    for(int i=0;i<=1;i++)
    for(int j=0;j<=1;j++)
        c.a[i][j]=0;
    for(int i=0;i<=1;i++)
        for(int j=0;j<=1;j++)
        for(int k=0;k<=1;k++)
            c.a[i][j]=(c.a[i][j]+x.a[i][k]*y.a[k][j]%mod)%mod;
    return c;
}
void init()
{
    A.a[0][0]=1;A.a[0][1]=1;
    A.a[1][0]=1;A.a[1][1]=0;
}
void quick_mod(ll n)
{
    for(int i=0;i<=1;i++)
    for(int j=0;j<=1;j++)
        if(i==j) B.a[i][j]=1;
        else B.a[i][j]=0;
    init();
    while(n){
        if(n&1) B=mat(B,A);
        A=mat(A,A);
        n>>=1;
    }
}
int main()
{
    ll n;
    scanf("%lld",&n);
    if(n==1||n==2){
        printf("1\n");
        return 0;
    }
    quick_mod(n-2);
    ans.a[0][0]=ans.a[0][1]=1;
    ans=mat(ans,B);
    printf("%lld\n",ans.a[0][0]);
    return 0;
}