解题报告：HDU_6185 Covering (轮廓线DP+高斯消元+矩阵快速幂)

题目链接

题意：

给一个4*n的表格，你有两种矩阵(1*2)，(2*1)，询问放满的方案数。

n<=1e9

思路：

显然公式应该是一个线性递推方程，知道后可以用矩阵快速幂在O( log(n) * m^3 )求得答案（m为方程的项数）

为了求这个方程，我们可以用轮廓线DP求的方程的前k项

然后假设一个k>m，用高斯消元求k*k的矩阵的秩，从而求得m

再用高斯消元求得方程即可

求的方程为

代码：

#include<bits/stdc++.h>

using namespace std;

const int mod = 1e9+7;



void init(){
   int dp[100][16]={0};
   dp[1][0]=dp[1][3]=dp[1][6]=dp[1][12]=dp[1][15]=1;
   for(int x=2;x<=16;x++){
      for(int i=0;i<16;i++){
         dp[x][i] += dp[x-1][i^15];
         for(int j=0;j<3;j++){
            if((i&(1<<j))&&(i&(1<<(j+1)))){
               int t = i^(1<<j)^(1<<(j+1));
               dp[x][i] += dp[x-1][t^15];
               //if(x==2)printf("%d is from %d\n",i,i^(1<<j)^(1<<(j+1)));
            }
         }//printf("dp[%d][%d]--->%d\n",x,i,dp[x][i]);
      }dp[x][15]+=dp[x-1][15];
      printf("dp[%d]--->%d\n",x,dp[x][15]);
   }
}


namespace fast_Matrix{

int n=4;
const int MAX_E = 4;

inline void debug(long long a[MAX_E][MAX_E]){
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            printf("%I64d%c",a[i][j],j!=n-1?' ':'\n');
        }
    }printf("\n");
}

long long ans[5]={1,1,5,11,36};

long long fast_mat(long long y)
{

long long tmp[MAX_E][MAX_E]={0};
long long mut[MAX_E][MAX_E]={
1,5,1,mod-1,
1,0,0,0,
0,1,0,0,
0,0,1,0,
};
long long res[MAX_E][MAX_E]={0};

for(int i=0;i<n;i++){
    res[i][i] = 1;
}


    while(y){
        if(y & 1){
            for(int i = 0; i < n; i++)
                for(int j = 0; j < n; j++)
                    tmp[i][j] = res[i][j];
            memset(res,0,sizeof(res));
            for(int i = 0; i < n; i++)
                for(int j = 0; j < n; j++)if(tmp[i][j])
                    for(int k = 0; k < n; k++)
                        res[i][k] = (res[i][k] + tmp[i][j] * mut[j][k])%mod ;
        }

        for(int i = 0; i < n; i++)
            for(int j = 0; j < n; j++)
                tmp[i][j] = mut[i][j];
        memset(mut,0,sizeof(mut));
        for(int i = 0; i < n; i++)
            for(int j = 0; j < n; j++)if(tmp[i][j])
                for(int k = 0; k < n; k++)
                    mut[i][k] = (mut[i][k] + tmp[i][j] * tmp[j][k])%mod ;
        y >>= 1;
    }
    long long val = 0;
    for(int i=0;i<n;i++)val = (val+res[0][i]*ans[n-i-1])%mod;
    return val;
}

}
using namespace fast_Matrix;

int main()
{
   init();
   long long t;
   while(scanf("%lld",&t)==1){
      if(t>=3)printf("%lld\n",fast_mat(t-3));
      else printf("%lld\n",ans[t]);
   }return 0;
}