Mondriaan's Dream POJ

Mondriaan's Dream POJ - 2411

可以用状压dp，但是要打一下表。暴力枚举行、这一行的状态、上一行的状态，判断如果上一行的状态能转移到这一行的状态就转移。

状态定义：ans[i][S]表示i行前已经全部填满，i行已经填上的列为集合S。如果有竖着的，全部当做用这一行的去补满上一行缺的。

（貌似还是插头dp的入门题）

#include<cstdio>
#include<cstring>
typedef long long LL;
LL f[20][20];
LL h,w;
/*LL ans[12][2050];
LL h,w;
bool judge(LL a,LL b)
{
    LL i,num=0,t1,t2;
    for(i=1;i<=w;i++)
    {
        t1=a&1;
        t2=b&1;
        a>>=1;
        b>>=1;
        if(t1==1&&t2==1)
            num++;
        else if(t1==0&&t2==0)
            return false;
        else
        {
            if(num&1)
                return false;
            num=0;
        }
    }
    if(num&1)
        return false;
    else
        return true;
}
LL get(LL h,LL w)
{
    LL i,j,k;
    memset(ans,0,sizeof(ans));
    ans[0][(1<<w)-1]=1;
    for(i=1;i<=h;i++)
        for(j=0;j<(1<<w);j++)
            for(k=0;k<(1<<w);k++)
                if(judge(k,j))
                {
                    ans[i][j]+=ans[i-1][k];
                }
    return ans[h][(1<<w)-1];
}
int main()
{
    for(h=1;h<=11;h++)
        for(w=1;w<=11;w++)
            printf("f[%lld][%lld]=%lld;
",h,w,get(h,w));
    return 0;
}*/
int main()
{
    f[1][1]=0;
    f[1][2]=1;
    f[1][3]=0;
    f[1][4]=1;
    f[1][5]=0;
    f[1][6]=1;
    f[1][7]=0;
    f[1][8]=1;
    f[1][9]=0;
    f[1][10]=1;
    f[1][11]=0;
    f[2][1]=1;
    f[2][2]=2;
    f[2][3]=3;
    f[2][4]=5;
    f[2][5]=8;
    f[2][6]=13;
    f[2][7]=21;
    f[2][8]=34;
    f[2][9]=55;
    f[2][10]=89;
    f[2][11]=144;
    f[3][1]=0;
    f[3][2]=3;
    f[3][3]=0;
    f[3][4]=11;
    f[3][5]=0;
    f[3][6]=41;
    f[3][7]=0;
    f[3][8]=153;
    f[3][9]=0;
    f[3][10]=571;
    f[3][11]=0;
    f[4][1]=1;
    f[4][2]=5;
    f[4][3]=11;
    f[4][4]=36;
    f[4][5]=95;
    f[4][6]=281;
    f[4][7]=781;
    f[4][8]=2245;
    f[4][9]=6336;
    f[4][10]=18061;
    f[4][11]=51205;
    f[5][1]=0;
    f[5][2]=8;
    f[5][3]=0;
    f[5][4]=95;
    f[5][5]=0;
    f[5][6]=1183;
    f[5][7]=0;
    f[5][8]=14824;
    f[5][9]=0;
    f[5][10]=185921;
    f[5][11]=0;
    f[6][1]=1;
    f[6][2]=13;
    f[6][3]=41;
    f[6][4]=281;
    f[6][5]=1183;
    f[6][6]=6728;
    f[6][7]=31529;
    f[6][8]=167089;
    f[6][9]=817991;
    f[6][10]=4213133;
    f[6][11]=21001799;
    f[7][1]=0;
    f[7][2]=21;
    f[7][3]=0;
    f[7][4]=781;
    f[7][5]=0;
    f[7][6]=31529;
    f[7][7]=0;
    f[7][8]=1292697;
    f[7][9]=0;
    f[7][10]=53175517;
    f[7][11]=0;
    f[8][1]=1;
    f[8][2]=34;
    f[8][3]=153;
    f[8][4]=2245;
    f[8][5]=14824;
    f[8][6]=167089;
    f[8][7]=1292697;
    f[8][8]=12988816;
    f[8][9]=108435745;
    f[8][10]=1031151241;
    f[8][11]=8940739824;
    f[9][1]=0;
    f[9][2]=55;
    f[9][3]=0;
    f[9][4]=6336;
    f[9][5]=0;
    f[9][6]=817991;
    f[9][7]=0;
    f[9][8]=108435745;
    f[9][9]=0;
    f[9][10]=14479521761;
    f[9][11]=0;
    f[10][1]=1;
    f[10][2]=89;
    f[10][3]=571;
    f[10][4]=18061;
    f[10][5]=185921;
    f[10][6]=4213133;
    f[10][7]=53175517;
    f[10][8]=1031151241;
    f[10][9]=14479521761;
    f[10][10]=258584046368;
    f[10][11]=3852472573499;
    f[11][1]=0;
    f[11][2]=144;
    f[11][3]=0;
    f[11][4]=51205;
    f[11][5]=0;
    f[11][6]=21001799;
    f[11][7]=0;
    f[11][8]=8940739824;
    f[11][9]=0;
    f[11][10]=3852472573499;
    f[11][11]=0;
    scanf("%lld%lld",&h,&w);
    while(h!=0&&w!=0)
    {
        printf("%lld
",f[h][w]);
        scanf("%lld%lld",&h,&w);
    }
}