HDU 6470 /// 矩阵快速幂

题目大意：

f[1]=1 f[2]=2

f[n]=f[n-1]+2*f[n-2]+n^3

在某博客截的图 现在忘记原博位置了 抱歉

根据递推式1和递推式3构造出两个矩阵

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define LL long long
#define inc(i,j,k) for(int i=j;i<=k;i++)
#define dec(i,j,k) for(int i=j;i>=k;i--)
#define gcd(i,j) __gcd(i,j)
#define mem(i,j) memset(i,j,sizeof(i))
const int N=1e5+;
const int M=;
const int mod=;

LL n,m;
struct MAT {
    LL a[M][M];
    MAT(){ mem(a,); }
    MAT operator*(MAT p) {
        MAT res;
        for(int i=;i<M;i++)
        for(int j=;j<M;j++)
        for(int k=;k<M;k++)
            res.a[i][j]=(res.a[i][j]+a[i][k]*p.a[k][j])%mod;
        return res;
    }
};
MAT mod_pow(MAT A,LL x) {
    MAT res;
    res.a[][]=;
    while(x) {
        if(x&) res=res*A;
        A=A*A; x>>=;
    }
    return res;
}
void init(MAT& A,MAT& B) {
    A.a[][]=,A.a[][]=,A.a[][]=,
    A.a[][]=,
                            A.a[][]=,A.a[][]=,A.a[][]=,A.a[][]=,
                                        A.a[][]=,A.a[][]=,A.a[][]=,
                                                    A.a[][]=,A.a[][]=,
                                                                A.a[][]=;
    B.a[][]=,B.a[][]=,B.a[][]=,B.a[][]=,B.a[][]=,B.a[][]=;
}

int main()
{
    int _; scanf("%d",&_);
    while(_--) {
        LL n; scanf("%lld",&n);
        MAT A; MAT ans;
        init(A,ans);
        ans=mod_pow(A,n-2LL)*ans;
        printf("%lld\n",ans.a[][]);
    }

    return ;
}