斐波那契数列与矩阵

#include <iostream>
#include <math.h>
#include <algorithm>
#include <stack>
#include <string.h>
#include <vector>
#include <cstring>
using namespace std;
#define ll long long
ll read() {
    ll f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9') {
        if(ss=='-')
            f=-1;
        ss=getchar();
    }
    while(ss>='0'&&ss<='9') {
        x=x*10+ss-'0';
        ss=getchar();
    }
    return x*f;//符号乘输的大小
}
const int mod=1000000007;
struct Matrix {//矩阵类
    ll a[3][3];
    Matrix() {
        memset(a,0,sizeof(a));//当使用Matrix m1的时候m1.a[3][3]数组里面的值将清空为0；
    }
    Matrix operator*(const Matrix &b)const {
    //operator*：操作符重载，即乘重新定义的意思，普通的*对正常规定的类型有效，无法对自己定义的结构体进行操作。
    //(const Matrix &b) 即m1.*(m2),类似于把m2当作参数传给m1；
        Matrix res;//两个矩阵相乘的结果；即结果矩阵。
        for(int i=1; i<=2; ++i)
            for(int j=1; j<=2; ++j)
                for(int k=1; k<=2; ++k)
                    res.a[i][j]=(res.a[i][j]+a[i][k]*b.a[k][j])%mod;//行和列相乘相加取mod
        return res;
    }
} ans,base;
void init() {
    base.a[1][1]=base.a[1][2]=base.a[2][1]=1;
    ans.a[1][1]=ans.a[1][2]=1;
}
void qpow(int b) {//矩阵的快速幂
    while(b) {
        if(b&1)ans=ans*base;
        base=base*base;//因为重载了*，所以可以直接乘；
        b>>=1;
    }
}
int main() {
    int n=read();
    if(n<=2)return puts("1");
    init();
    qpow(n-2);
    printf("%lld\n",ans.a[1][1]%mod);
    return 0;
}