H - Log Concave Sequences —— 矩阵快速幂

This way

题意：

你现在有无穷个0,1，2，现在让你构造一个长度为n的串，使得满足以下条件
对于第i位（2<=i=a[i-1]*a[i+1]
问你能构造多少种这个串

题解：

很明显它首先是dp，于是我写出这样一个式子：
dp[i][j]表示第i位取j时的情况数
但是这个式子是要从dp[i-1][k]dp[i-2][l]这里转移过来的，同时n的范围是1e18，所以这个很明显是矩阵快速幂，但是这种方法不行。
于是想到dp[i][j][k]表示在第i位取j，第i+1位取k的时候的状态数。这样转移就变成了一维的。
那么转移矩阵就是枚举ijk三个位置，j和i表示后面两位，k和j表示前面两位。在i*k<=j*j成立的时候，$b.m[i\ast 3+j][j\ast 3+k]++;$
因为每个值都有3种，所以i3+j表示dp[i+1]的状态，j*3+k表示dp[i]的状态。
好久没做矩阵快速幂了啊，对于
转移矩阵*初始值
的时候，初始值的矩阵是m*1的，也就是：

	for(int i=0;i<9;i++)
        a.m[i][0]=1;
    n-=2;
    while(n){
        if(n&1)
            a=b*a;
        b=b*b;
        n>>=1;
    }

但是对于
初始值*转移矩阵
的时候，初始矩阵是m*m的，也就是

	for(int i=0;i<9;i++)
        a.m[i][i]=1;
    n-=2;
    while(n){
        if(n&1)
            a=a*b;
        b=b*b;
        n>>=1;
    }

呕，我要吐了

#include 
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
int maxn=9;
struct Matrix
{
    ll m[9][9];
    Matrix(){
        memset(m,0,sizeof(m));
    }
    Matrix  operator +(const Matrix &b)const
    {
        Matrix c;
        for(int i=0; i>=1;
    }
    ll ans=0;
    for(int i=0;i<9;i++)
            ans=(ans+a.m[i][0])%mod;
    printf("%lld\n",ans);
    return 0;
}
/**/