luogu 1939 【模板】矩阵加速（数列）

题目大意：

a[1]=a[2]=a[3]=1

a[x]=a[x-3]+a[x-1] (x>3)

求a数列的第n项%1000000007

思路：

使用矩阵快速幂进行加速

在草稿纸上填了填数

然后就A了

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 #include<cstdlib>
 #include<cstring>
 #include<algorithm>
 #include<map>
 #include<vector>
 #include<queue>
 #define inf 2147483611
 #define ll long long
 #define MAXN 10010000
 #define MOD 1000000007
 using namespace std;
 inline ll read()
 {
     ll x=,f=;char ch=getchar();
     while(!isdigit(ch)) {if(ch=='-') f=-;ch=getchar();}
     while(isdigit(ch)) {x=x*+ch-'';ch=getchar();}
     return x*f;
 }
 ll T,n;
 struct mat {ll num[][];};
 mat mul (mat x,mat y)
 {
     mat res;
     memset(res.num,,sizeof(res.num));
     for(ll i=;i<;i++)
         for(ll j=;j<;j++)
             for(ll k=;k<;k++) (res.num[i][j]+=x.num[i][k]*y.num[k][j])%=MOD;
     return res;
 }
 ll q_pow(ll n)
 {
     if(n<=) return ;
     n-=;
     mat t,res;
     t.num[][]=t.num[][]=t.num[][]=t.num[][]=,t.num[][]=t.num[][]=t.num[][]=t.num[][]=t.num[][]=;
     res.num[][]=res.num[][]=res.num[][]=res.num[][]=res.num[][]=res.num[][]=res.num[][]=res.num[][]=,res.num[][]=;
     while(n)
     {
         //cout<<n<<" "<<res.num[0][0]<<" "<<res.num[0][1]<<" "<<res.num[0][2]<<endl;
         if(n&) res=mul(res,t);
         t=mul(t,t);
         n>>=;
     }
     return res.num[][];
 }
 int main()
 {
     T=read();
     while(T--) {n=read();printf("%lld\n",q_pow(n));}
 }