HDU 6185 Covering （递推+矩阵快速幂）*

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6185

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))

#define pii pair<ll,ll>
#define mk(x,y) make_pair(x,y)

const int  maxn =4;
const int mod=1e9+7;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意：给定两种砖块格式，和4*n的地板，
问用这两种砖块铺满地板的方案数是多少（答案模1e9+7).

首先我们不难推出可以用五种情况来到达当前位置。
第一种和第二种都是直接完全铺满（无分割）,
其余三种都是可以延扩或者直接封满的，有个递推的展开式子在里面，
f(n)=f(n-1)+f(n-2)+2(f(n-2)+f(n-3)+f(n-4)+....)+f(n-2)+f(n-4)+f(n-6)+....f(1)/f(0).
用上述式子可以得到f(n)关于f(n-1),f(n-2),f(n-3),f(n-4)的递推式，用矩阵快速幂搞一搞就行了。
懒得画图了，贴下有图的博客：
https://blog.csdn.net/elbadaernu/article/details/77825979
时间复杂度：O（logn）.
*/

const ll mi[maxn][maxn]={{0,1,0,0},{0,0,1,0},{0,0,0,1},{-1,1,5,1}};
struct mat{
    ll a[maxn][maxn];
    mat(){
        mst(a,0);rep(i,0,maxn) a[i][i]=1;
    }
    mat operator*(const mat& y) const{
        mat ret;
        rep(i,0,maxn) rep(j,0,maxn){
            ret.a[i][j]=0;
            rep(k,0,maxn) ret.a[i][j]=(ret.a[i][j]+a[i][k]*y.a[k][j]%mod+mod)%mod;
        }
        return ret;
    }
};
mat qucik_pow(ll y){
    mat ret,x;mst(x.a,0);
    rep(i,0,maxn) rep(j,0,maxn) x.a[i][j]=mi[i][j];
    for(;y;y>>=1,x=x*x) if(y&1) ret=x*ret;return ret;
}
int main()
{
    ll n,res,ans[4]={1,5,11,36};
    while(scanf("%lld",&n)!=EOF){
        if(n<=4) printf("%lld\n",ans[n-1]);
        else{
            mat ret=qucik_pow(n-4);  res=0;
            rep(i,0,maxn) res=(res+ret.a[maxn-1][i]*ans[i]%mod+mod)%mod;
            printf("%lld\n",res);
        }
    }
    return 0;
}