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

题目链接: Covering.

大致题意：

用1 * 2规格的地毯去铺4 * n规格的地面，已知n，问有多少种不同的方案使得地面恰好被铺满且地毯不重叠。答案对1000000007取模。
（1≤n≤10¹⁸)

思路：

考虑全部铺满时，前2列的放法。有如下5种情况：

故可导出递推式

---

AC代码：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#define mod 1000000007
using namespace std;

typedef long long ll;

struct M
{
    ll m[4][4];
};

M Multiply(M a,M b)
{
    M tmp;
    memset(tmp.m,0,sizeof(tmp.m));
    for(int i=0; i<4; i++)
        for(int j=0; j<4; j++)
            for(int k=0; k<4; k++)
            {
                tmp.m[i][j] += a.m[i][k] * b.m[k][j];
                tmp.m[i][j] %= mod;
            }
    return tmp;
}

M Fastpower(M a,ll n)
{
    M tmp;
    memset(tmp.m,0,sizeof(tmp.m));
    for(int i=0; i<4; i++)
        tmp.m[i][i]=1;
    while(n)
    {
        if(n & 1)
        {
            tmp = Multiply(tmp,a);
        }
        a = Multiply(a,a);
        n >>= 1;
    }
    return tmp;
}

int main()
{
    ll n;
    //freopen("D:\\1.txt","r",stdin);
    while(~scanf("%lld",&n))
    {
        if(n==1)
            printf("1\n");
        else if(n==2)
            printf("5\n");
        else if(n==3)
            printf("11\n");
        else if(n==4)
            printf("36\n");
        else
        {
            M st;
            memset(st.m,0,sizeof(st.m));
            st.m[0][0]=1;
            st.m[1][0]=5;
            st.m[2][0]=1;
            st.m[3][0]=-1;
            st.m[0][1]=1;
            st.m[1][2]=1;
            st.m[2][3]=1;

            M a;
            memset(a.m,0,sizeof(a.m));
            a.m[0][0]=36;
            a.m[0][1]=11;
            a.m[0][2]=5;
            a.m[0][3]=1;

            st = Fastpower(st,n-4);
            a = Multiply(a,st);
            printf("%lld\n",(a.m[0][0]+mod)%mod);

        }
    }
    return 0;
}