铺砖问题

最猥琐的一道题。
分两块做：
①前50%
很显然是状压dp，我们设f[i][j]表示前i行，第i行状态为j的方案数。
状态中，1表示下凸，0表示其他。
然后我们预处理出所有搭配起来合法的状态，就可以直接转移。
什么叫合法呢？就是这两个状态i和j，i&j=0因为不可能两个都下凸，这样就重叠了。
而且i|j以后，连续的0必须是偶数个，不然就不会合法，这很显然。

#include<cstdio>
#include<iostream>
#include<cstring>
#define mo 1000000007
using namespace std;
long long n,m;
long long f[105][2050];
int b[2050][1005];
bool bz[2050]; 
char s[110]; 
int main()
{
    scanf("%s%lld",s+1,&m);
    int i,j,x,k;
    memset(bz,true,sizeof(bz));
    for (i=0;i<=(1<<m)-1;++i)
    {
        x=i;
        j=0;
        k=1;
        while (x||k<=m+1)
        {
            if ((x&1)==1||k>m)
            {
                if(j%2) 
                {
                    bz[i]=false;
                    break;
                }
                j=0;    
            }
            if ((x&1)==0)
                ++j;
            ++k;    
            x>>=1;
        }           
    } 
    for (i=0;i<=(1<<m)-1;++i)
        for (j=0;j<=(1<<m)-1;++j)
            if((i&j)==0&&bz[i|j]) b[i][++b[i][0]]=j;
    if(strlen(s+1)<=3)
    {   
        for (i=1;i<=strlen(s+1);++i)
            n=n*10+(s[i]-48); 
        f[0][0]=1;
        for (i=1;i<=n;++i)
            for (j=0;j<=(1<<m)-1;++j)
                for (k=1;k<=b[j][0];++k)
                    (f[i][j]+=f[i-1][b[j][k]])%=mo;
        printf("%lld\n",f[n][0]%mo);
    }
}

②后50%
我们发现每次我们预处理的东西p[i][j]=1表示从i这个状态转移到j这个状态合法。
那么每次转移一行都是这个所以直接就可以用一个矩阵乘法，也就是你预处理那东西，直接快速幂上n次，答案就是p[0][0]了。

#include<cstdio>
#include<iostream>
#include<cstring>
#define mo 1000000007
using namespace std;
long long n,m;
long long f[105][2050];
int b[2050][1005];
bool bz[2050]; 
char s[210]; 
long long matr[35][35],ans[35][35],matrix[35][35],ans2[35][35];
int n1[205];
int main()
{
    scanf("%s%lld",s+1,&m);
    int i,j,x,k;
    memset(bz,true,sizeof(bz));
    for (i=0;i<=(1<<m)-1;++i)
    {
        x=i;
        j=0;
        k=1;
        while (x||k<=m+1)
        {
            if ((x&1)==1||k>m)
            {
                if(j%2) 
                {
                    bz[i]=false;
                    break;
                }
                j=0;    
            }
            if ((x&1)==0)
                ++j;
            ++k;    
            x>>=1;
        }           
    } 
    int l=0;
    if(m<=5) l=1;
    for (i=0;i<=(1<<m)-1;++i)
        for (j=0;j<=(1<<m)-1;++j)
            if((i&j)==0&&bz[i|j]) 
            {
                b[i][++b[i][0]]=j;
                if(l) matr[i][j]=1;
            }
    if(strlen(s+1)<=3)
    {   
        for (i=1;i<=strlen(s+1);++i)
            n=n*10+(s[i]-48); 
        f[0][0]=1;
        for (i=1;i<=n;++i)
            for (j=0;j<=(1<<m)-1;++j)
                for (k=1;k<=b[j][0];++k)
                    (f[i][j]+=f[i-1][b[j][k]])%=mo;
        printf("%lld\n",f[n][0]%mo);
    }
    else
    {
        for (i=strlen(s+1);i>=1;--i)
            n1[i]=s[strlen(s+1)-i+1]-48;
        n1[0]=strlen(s+1);
        ans[0][0]=1;
        while(n1[0]>1||(n1[0]==1&&n1[1]>0))
        {
            if(n1[1]&1) 
            {
                memset(ans2,0,sizeof(ans2));
                for (k=0;k<=(1<<m)-1;++k)
                    for (i=0;i<=(1<<m)-1;++i)
                        for (j=0;j<=(1<<m)-1;++j)
                            ans2[i][j]=(ans2[i][j]+ans[i][k]*matr[k][j])%mo;
                memcpy(ans,ans2,sizeof(ans2));
            }
            memset(matrix,0,sizeof(matrix));
            for (k=0;k<=(1<<m)-1;++k)
                for (i=0;i<=(1<<m)-1;++i)
                    for (j=0;j<=(1<<m)-1;++j)
                        (matrix[i][j]+=matr[i][k]*matr[k][j])%=mo;
            memcpy(matr,matrix,sizeof(matrix));
            if (n1[1]%2==1) n1[1]--;
            for (i=1;i<n1[0];i++)
            {
                if (n1[i+1]%2==1) n1[i+1]--,n1[i]+=10;
                n1[i]/=2;
            }
            n1[n1[0]]/=2;
            if (n1[n1[0]]==0) n1[0]--;
        }   
        printf("%lld\n",ans[0][0]);
    }
}