POJ

传送门

Solution:

1.矩阵分块 题解在这里

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define R register
#define go(i,a,b) for(R int i=a;i<=b;i++)
#define yes(i,a,b) for(R int i=a;i>=b;i--)
#define ll long long
#define db double
using namespace std;
ll n,m;int mod;
struct node{ll mt[3][3];};
node calc(node x,node y,int a,int b,int c)
{
    node z;memset(z.mt,0,sizeof(z.mt));
    go(i,1,a)
        go(j,1,b)
        go(k,1,c)
        z.mt[i][j]=(z.mt[i][j]+x.mt[i][k]*y.mt[k][j])%mod;
    return z;
}
node a,b,c;
void sol()
{
    while(m)
    {
        if(m&1) b=calc(b,a,2,2,2);
        m>>=1;a=calc(a,a,2,2,2);
    }
}
int main()
{
    scanf("%lld",&n);m=n-2;mod=1000000007;
    if(n<=2){printf("1");return 0;}
    a.mt[1][1]=1,a.mt[1][2]=1,a.mt[2][1]=1,a.mt[2][2]=0;
    b.mt[1][1]=1,b.mt[1][2]=0,b.mt[2][1]=0,b.mt[2][2]=1;
    c.mt[1][1]=1;c.mt[2][1]=1;
    sol();
    c=calc(b,c,2,1,2);
    printf("%lld",c.mt[1][1]);
    return 0;
}

但是这个代码有点问题 就是k=1和k=2的时候会错，是因为递推式子的原因，详见题解

但还是AC了 因为没有k=1,2这么水的数据 这里说一下我就不改了

2.分治 题解

没啦qwq