1643：Fibonacci 前 n 项和

 

 s[n]=s[n-1]*1+f[n]+f[n-1]*0

f[n+1]=s[n-1]*0+f[n]+f[n-1]

f[n]=s[n-1]*0+f[n-1]*0+f[n]

按照这三个方程构造矩阵即可

#include <iostream>
#include <cstdio>
#include <cstring>
#define N 3+5

using namespace std;
typedef long long ll;
    ll n,m;
    int i,j,k;
    //int a[N];
    ll mod;

    struct node
    {
        ll g[N][N];
        node(){ memset(g,0,sizeof g); }
        node operator*(node a){
            node ans;
            for(int i=1;i<=3;i++){
                for(int j=1;j<=3;j++){
                    for(int k=1;k<=3;k++){
                        ans.g[i][j]=(ans.g[i][j]+g[i][k]*a.g[k][j]%mod)%mod;
                    }
                }
            }
            return ans;
        }
    }ans;

node pow_mod(node a,int x)
{
    node ans;
    for(int i=1;i<=3;i++) ans.g[i][i]=1;
    while(x){
        if(x&1) ans=ans*a;
        a=a*a;
        x>>=1;
    }
    return ans;
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    while(cin>>n>>mod){
        node a; //转移矩阵
        a.g[1][1]=a.g[1][2]=a.g[2][2]=a.g[2][3]=a.g[3][2]=1;
        a=pow_mod(a,n-1);
        cout<<(a.g[1][1]+a.g[1][2]+a.g[1][3])%mod<<endl;
    }
    return 0;
}