POJ1185 状态压缩DP 又是一道有意义的水题

之前做过POJ3285，那道题下一行的状态只和上一行有关，所以每一行只需要记录一行的状态。

这一个炮兵阵地，炮可以打两行，所以不仅仅与上一行有关，还与上上一行有关，所以每一行的状态都由上两行状态推出，因为要记录两个状态，所以用了一个三维的dp，

dp[i][j][k]，表示第i行状态为state[k]，第i-1行状态为state[j]时前I行最多有几个炮兵。

写的时候一开始在初始化的时候没有注意硬件条件，挂了一发。

写完这个程序，100+行，编译了一下，还是有一个for循环中int忘记定义，除了这个之外就没有错误了。这种感觉真不错，不过还是有逻辑错。

第一个逻辑错是混淆了POJ3285和这个题的所求的东西，在特判行数为1的时候统计成了此时的状态数量。

第二个逻辑错是初始化的时候没加上对硬件条件的判断。

真是的，要谨慎啊。

不能疏忽。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;

int state[62];
int num[62];
int dp[102][62][62];//dp[i][j][k] 表示前i行中 i行为k种状态 i-1行为j种状态时所能放置的最大炮兵数
int tot=0;
int row[200];
int cal(int i)
{
    int count=0;
    while(i>0)
    {
        count++;
        i-=(i&(-i));
    }
    return count;
}
void init(int n)
{
    int k=1<<n;
    for(int i=0;i<k;i++)
    {
        if((i&(i<<1))==0&&(i&(i<<2))==0)
        {
            num[tot]=cal(i);
            state[tot++]=i;
        }
    }
}

char tmp[100];
int main()
{
    memset(dp,0,sizeof(dp));
    int m,n;
    scanf("%d%d",&m,&n);
    init(n);
    for(int i=0;i<m;i++)
    {
        scanf("%s",tmp);
        for(int j=n-1;j>=0;j--)
        {
            if(tmp[n-1-j]=='P')
                row[i]+=1<<j;
        }
    }
    if(m==1)
    {
        int count=-1;
        for(int i=0;i<tot;i++)
        {
            if((row[0]&state[i])==state[i])
                count=max(num[i],count);
        }
        printf("%d\n",count);
        return 0;
    }
    for(int i=0;i<tot;i++)
    {
        for(int j=0;j<tot;j++)
        {
            if((row[0]&state[i])==state[i]&&(row[1]&state[j])==state[j]&&(state[i]&state[j])==0)
            {
                dp[1][i][j]=num[i]+num[j];
            }
        }
    }
    if(m==2)
    {
        int count=-1;
        for(int i=0;i<tot;i++)
        {
            for(int j=0;j<tot;j++)
            {
                count=max(dp[1][i][j],count);
            }
        }
        printf("%d\n",count);
        return 0;
    }
    for(int i=2;i<m;i++)
    {
        for(int j=0;j<tot;j++)
        {
            if((state[j]&row[i])==state[j])
            {
                for(int k=0;k<tot;k++)
                {
                    for(int l=0;l<tot;l++)
                    {
                        if(dp[i-1][k][l]&&((state[j]&state[k])==0)&&((state[j]&state[l])==0))
                        {
                            dp[i][l][j]=max(dp[i][l][j],dp[i-1][k][l]+num[j]);
                        }
                    }
                }
            }
        }
    }
    int max_=-1;
    for(int i=0;i<tot;i++)
    {
        for(int j=0;j<tot;j++)
        {
            if(dp[m-1][i][j]>max_)
                max_=dp[m-1][i][j];
        }
    }
    printf("%d\n",max_);
    return 0;
}