Poj 1185&Poj 2411 状压dp（from Poj2663）

鱼头最近查刷题有点紧，事情又多，红了。不得已敲了几发水题，罪恶感爆棚。
不过敲水题过程中发现了这样一个题：
POJ2663：问3*n的方格用1*2的domino填充有多少种方法。
瞎几把YY了之后得出递推公式如下
a[i]=3*a[i-1]+2*b[i-1];
b[i]=a[i-1]+b[i-1];
其中a[i]是n=2*i的答案啦（显然n为奇数的结果为0）
b[i]是这种图形的答案
。。。111
。。。111
。。。011
然后就xjb搞完了，代码：

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
ll a[30];
ll b[30];
int main()
{
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    a[1]=3;b[1]=1;
    for(int i=2;i<=20;i++)
    {
        a[i]=3*a[i-1]+2*b[i-1];
        b[i]=a[i-1]+b[i-1];
    }
    int n;
    while(scanf("%d",&n),n!=-1)
    {
        if(!n){printf("1\n");continue;}
        else if(n%2){printf("0\n");continue;}
        else printf("%lld\n",a[n/2]);
    }
    return 0;
}

然后考虑m*n的情况，找到POJ2411，记起来宝哥报告讲的状压dp什么的，然而当时只是理论ac了，遂开敲。
题意：就把上面的3换成m，还是问domino填充种数。
做法就是把每行保存成一个2进制数，0表示空，1表示放，注意：竖着放的第一行先空着，视为0。
举个例子：

1–11—
1–11–1
–1—–1
–1——

对应
01100111
11111110
11011111
11111111

状态转移就是dp[i][j]=sum(dp[i-1][j1]),其中j1为与j不冲突的状态。
因为懒，没高兴想通过j构造j1,直接枚举所有j1判断是否和j冲突了，所以1300ms过的（菜如我）;
代码：

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
int m,n;
ll dp[12][4000];
bool judge(int p,int q)
{
    int pos=0;
    for(int i=0;i<n;)
    {
        if(p&(1<<i))
        {
            if(q&(1<<i))
                if(i==n-1||!(p&(1<<(i+1)))||!(q&(1<<(i+1))))return false;
                else i+=2;
            else i++;
        }
        else
        {
            if(q&(1<<i))i++;
            else return false;
        }
    }
    return true;
}
ll dfs(int f,int s)
{
    if(dp[f][s]!=-1)return dp[f][s];
    ll ans=0;
    for(int i=0;i<(1<<n);i++)
        if(judge(s,i))ans+=dfs(f-1,i);
    return dp[f][s]=ans;
}
bool jj1(int num)
{
    for(int i=0;i<n;)
    {
        if(num&(1<<i))
            if(i==n-1||!(num&(1<<(i+1))))return false;
            else i+=2;
        else i++;
    }
    return true;
}
int main()
{
    while(~scanf("%d %d",&m,&n))
    {
        if(!m&&!n)break;
        memset(dp,-1,sizeof(dp));
        for(int i=0;i<(1<<n);i++)
            if(jj1(i))dp[1][i]=1;
            else dp[1][i]=0;
        ll ans=dfs(m,(1<<n)-1);
        printf("%lld\n",ans);
        //printf("%lld\n",dp[2][1]);
    }
    return 0;
}

然后是POJ1185，就是那个炮兵阵地的题，题是中文的，略题意。
与上题不同的是要保存两行的状态，dp[i][now][past]=max(dp[i-1][past][pastspast]),其中pastspast和now和past和地图都不冲突，bingo！
因为每行每个之间至少距离2，所以根本没有1024种，只有六十几吧，先预处理一下就ok了。
还要讨论只有一行的状况。
代码：

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
int m,n;
ll dp[12][4000];
bool judge(int p,int q)
{
    int pos=0;
    for(int i=0;i<n;)
    {
        if(p&(1<<i))
        {
            if(q&(1<<i))
                if(i==n-1||!(p&(1<<(i+1)))||!(q&(1<<(i+1))))return false;
                else i+=2;
            else i++;
        }
        else
        {
            if(q&(1<<i))i++;
            else return false;
        }
    }
    return true;
}
ll dfs(int f,int s)
{
    if(dp[f][s]!=-1)return dp[f][s];
    ll ans=0;
    for(int i=0;i<(1<<n);i++)
        if(judge(s,i))ans+=dfs(f-1,i);
    return dp[f][s]=ans;
}
bool jj1(int num)
{
    for(int i=0;i<n;)
    {
        if(num&(1<<i))
            if(i==n-1||!(num&(1<<(i+1))))return false;
            else i+=2;
        else i++;
    }
    return true;
}
int main()
{
    while(~scanf("%d %d",&m,&n))
    {
        if(!m&&!n)break;
        memset(dp,-1,sizeof(dp));
        for(int i=0;i<(1<<n);i++)
            if(jj1(i))dp[1][i]=1;
            else dp[1][i]=0;
        ll ans=dfs(m,(1<<n)-1);
        printf("%lld\n",ans);
        //printf("%lld\n",dp[2][1]);
    }
    return 0;
}