poj 1185 炮兵阵地 【状压dp】

题目链接：http://poj.org/problem?id=1185

---

题目大意：一个方格组成的矩阵，每个方格可以放大炮用P表示，不可以放大炮用H表示\，让放最多的大炮，大炮与大炮间不会互相攻击。

/*
dp[i][j][k] 表示第i行状态为k，第i-1状态为j时的最大炮兵个数。 
递归方程：dp[i][l][j]=max(dp[i][l][j],dp[i-1][k][l]+num[j]); num[j]为j状态中1的个数 
dp[1][1][i]=num[i] 状态i能够满足第一行的硬件条件（注意：这里的i指的是第i个状态，
不是一个二进制数，开一个数组保存二进制状态） 
*/
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<math.h>
#include<map>
#include<set>
#include<queue>
#include<vector>
#include<stack>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const double pi=acos(-1.0);
const ll mod=1e8;
const int N=70+10;

int cur[105],e,m,n;
int num[105];//满足每个状态的二进制1的个数
int dp[105][N][N];//一维是第i行，二维前一行的状态，三维当前行的状态
int state[N];//满足“同一行互打不到”情况的状态

void init()//保存满足情况的状态
{
    e=0;
    for(int i=0;i<(1<<m);i++)
    {
        if((i&(i<<1))==0&&(i&(i<<2))==0)
            state[e++]=i;
    }
}

bool judge(int x,int k)//判断该状态是否满足当前行的情况
{
    if(x&cur[k]) return false;
    return true;
}

int cal(int x)//就算二进制中1的个数
{
    int s=0;
    while(x)
    {
        s++;
        x&=(x-1);
    }
    return s;
}

int main()
{
    while(~scanf("%d %d",&n,&m))
    {
        if(n==0&&m==0) break;
        init();
        char str[N];
        int ans=0;
        for(int i=1;i<=n;i++)
        {
            cur[i]=0;
            scanf("%s",str);
            for(int j=0;j<m;j++)
            {
                if(str[j]=='H')
                    cur[i]+=(1<<j);
            }
        }
        memset(dp,-1,sizeof(dp));
        for(int i=0;i<e;i++)
        {
            num[i]=cal(state[i]);
            if(judge(state[i],1))
                dp[1][0][i]=num[i];
        }
        for(int i=2;i<=n;i++)
        {
            for(int j=0;j<e;j++)//遍历当前行的状态
            {
                if(!judge(state[j],i)) continue;
                for(int k=0;k<e;k++)//遍历上两行的状态
                {
                    if(state[j]&state[k]) continue;
                    for(int l=0;l<e;l++)//遍历前一行的状态
                    {
                        if(state[j]&state[l]) continue;
                        if(dp[i-1][k][l]==-1) continue;
                        dp[i][l][j]=max(dp[i][l][j],dp[i-1][k][l]+num[j]);
                    }   //找最大值
                }
            }
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<e;j++)
            {
                for(int k=0;k<e;k++)
                    ans=max(ans,dp[i][j][k]);
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}