poj1185状压dp

传送门](http://poj.org/problem?id=1185)
汉语题，就不解释了。
m<=10,且炮的距离是2，那么仅仅考虑一行，0000000000—-1111111111，（1表示放，0表示不放）这些状态中最多有60个能放下炮的状态。称为可行状态S；
放置4个情况1种，
3个情况(4+3+2+1)+(3+2+1)+(2+1)+1=20,
2个情况为7+6+5+4+3+2+1=28,
放置1个为10种，
放置0个为1种，其和为60种。
1，我们预先处理出可以放的这60个状态。
2，我们输入时处理出每一行的状态，称为实际状态W，PHHPH用01101表示，因为只有P能放炮，让P为0的话，S&W只要不是0就不可以放，这里可以自己简单思考一下，然后我们用dp[i][j][k]一个三维数组进行状态转移，dp[i][j][k]的值表示第i行的状态为第j个状态，第i-1行的状态为第k个状态时 第i行可以放的炮的个数，那么dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][L]+num[i]);（num[i]表示第i个状态可以放的炮的个数，L表示i-2行的第L个状态。这些都好处理的）。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
const int MAXN=105;
const int MAXM=65;
int satisfy[MAXM],tot,Mapstatus[MAXN];
int n,m;
char str[MAXN][12];
int dp[MAXN][65][65],num_one[MAXM];
bool satis(int num)
{
    if(num&(num<<1))return false;
    if(num&(num<<2))return false;
    return true;
}
int one_num(int num)
{
    int ans=0;
    while(num)
    {
        if((num&1))
            ans++;
        num=num/2;
    }
    return ans;
}
void init()
{
    tot=0;
    for(int i=0;i<((1<<m));i++)
    {
        if(satis(i))
        {
            satisfy[tot]=i;
            num_one[tot++]=one_num(i);
        }
    }
}
void solve()
{
    memset(dp,-1,sizeof(dp));
    for(int i=0;i<tot;i++)
    {
        if(!(Mapstatus[0]&satisfy[i]))
             dp[0][i][0]=num_one[i];
    }

    for(int i=1;i<n;i++)
    {
        for(int j=0;j<tot;j++)///第i行状态
        {
            if(Mapstatus[i]&satisfy[j])///第j个状态根本无法放在第i行上
                continue;
            for(int k=0;k<tot;k++)///第i-1行状态
            {
                if(satisfy[k]&satisfy[j])///检测第i-1行的状态与第i行的状态是否冲突
                continue;
                for(int l=0;l<tot;l++)///第i-2行状态
                {
                    if(satisfy[l]&satisfy[j])///检测第i-2行的状态与第i行的状态是否冲突
                    continue;
                    ///上行的状态不冲突，且与环境状态不冲突，就继续向下执行
                    if(dp[i-1][k][l]==-1)///该情况是要绝对有结果的，不然不能作为下面的转移式
                        continue;
                    dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+num_one[j]);
                }
            }
        }
    }
}
int main()
{
    while(scanf("%d%d",&n,&m)!=-1)
    {
        init();
        memset(Mapstatus,0,sizeof(Mapstatus));
        for(int i=0;i<n;i++)
        {
            scanf("%s",str[i]);
            for(int j=m-1;j>=0;j--)
            {
                if(str[i][j]=='P')///P为0，H为1,
                {
                    Mapstatus[i]*=2;
                }
                else
                {
                    Mapstatus[i]=Mapstatus[i]*2+1;
                }
            }
        }
        solve();
        int ans=0;
        for(int i=0;i<tot;i++)
        {
            for(int j=0;j<tot;j++)
                ans=max(ans,dp[n-1][i][j]);
        }
        cout<<ans<<endl;
    }
    return 0;
}

/**

5 4
PHPP
PPHH
PPPP
PHPP
PHHP
*/