一发非常猎奇的AC求大佬解答

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=3,mod=1e9+7;
ll n;
struct matrix{
    int a[N][N];
    matrix(){
        for(int i=1;i<N;i++)
            for(int j=1;j<N;j++)
                a[i][j]=0;
    }
    void init(){
        for(int i=1;i<N;i++)
            for(int j=1;j<N;j++)
                a[i][j]=1;
    }
};
matrix operator*(const matrix& a,const matrix& b){
    matrix c;
    for(int i=1;i<N;i++)
        for(int j=1;j<N;j++)
            for(int k=1;k<N;k++)
                c.a[i][j]=(1ll*a.a[i][k]*b.a[k][j]%mod+c.a[i][j])%mod;
    return c;
}
matrix qpow(matrix d,ll z){
    matrix ret;
    ret.init();
    while(z){
        if(z&1)ret=ret*d;
        z>>=1;
        d=d*d;
    }
    return ret;
}
matrix a;
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin>>n;
    if(n==1){
        cout<<1<<endl;
        return 0;
    }
    a.a[1][1]=1;a.a[1][2]=1;
    a.a[2][1]=1;a.a[2][2]=0;
    a=qpow(a,n-2);
    cout<<a.a[1][1]<<endl;
    return 0;
}