POJ-1185 炮兵阵地（状压dp)

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

炮兵阵地

Time Limit: 2000MS	Memory Limit: 65536K
Total Submissions: 37531	Accepted: 14352

Description

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N *M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示：

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。

Input

第一行包含两个由空格分割开的正整数，分别表示N和M；
接下来的N行，每一行含有连续的M个字符(‘P’或者’H’)，中间没有空格。按顺序表示地图中每一行的数据。N <= 100；M <= 10。

Output

仅一行，包含一个整数K，表示最多能摆放的炮兵部队的数量。
Sample Input

5 4
PHPP
PPHH
PPPP
PHPP
PHHP
Sample Output

6
思路：一道状压dp，和我刚写的那道题很类似方格取数（1）但是我感觉比那个要难一点点，昨天下午弄了一下午也没弄好，一个小细节卡死我了，我用的是(tot[j]&tot[k]&tot[u])=0,但是之前的循环中我已经保证了（tot[j]&tot[k])=0,对于任何的tot[u],0&tot[u]一定为0。（一下午也没看出来，都快看傻逼了都）
因为这道题牵扯到前两行的状态，所有我把一二行单独判断，对于每个m位的二进制数，我们要判断((i&(i>>1))&&(i&(i>>2)))=0,因为不能与前两位的炮发生冲突，还要判断，这些位置能不能放，因为有的是"H",不能放炮车，详情看代码。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
using namespace std;
char mp[105][15];
int dp[105][200][200];//dp[i][j][k]表示第i行状态为j，第i-1行状态为k,符合要求的数量。
int tot[1000];
int n,m;
int pan(int i,int x)
{
    int t=1;
    int sum=0;
    while(x)
    {
        if(x&1)
        {
            if(mp[i][m-t]=='H')
            {
                sum=1;
                break;
            }
        }
        x/=2;
        t++;
    }
    return sum;
}
int jis(int i,int x)
{

    int sum=0;
    while(x)
    {
        if(x&1)
        {
            sum++;
        }
        x/=2;

    }
    return sum;
}

int main()
{
    scanf("%d%d",&n,&m);
    int cut=0;
    for(int i=0; i<(1<<m); i++)//判断符合要求的数字
    {
        if(((i&(i>>2))==0)&&((i&(i>>1))==0))
        {
            tot[++cut]=i;
        }
    }
    int ma=-1;
    for(int i=0; i<n; i++)
    {
        scanf("%s",mp[i]);
    }
    for(int j=1; j<=cut; j++) //单独判断第一行
    {
        if(!pan(0,tot[j]))
        {
            int val=jis(0,tot[j]);
           dp[1][j][1]=val;

        }
    }
    for(int i=2; i<=n; i++)
    {
        for(int j=1; j<=cut; j++)
        {
            if(!pan(i-1,tot[j]))
            {
                int val=jis(i-1,tot[j]);
                for(int k=1; k<=cut; k++)
                {
                    if((tot[j]&tot[k])==0&&!pan(i-2,tot[k]))
                    {
                        if(i<3)
                        {

                                    dp[i][j][k]=dp[i-1][k][1]+val;

                        }
                        else
                        {
                            for(int u=1; u<=cut; u++)
                            {
                                if((tot[j]&tot[u])==0&&(!pan(i-3,tot[u])))
                                {
                                    dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][u]+val);
                                     //printf("%d#%d**%d**%d**%d\n",i,tot[j],val,tot[k],tot[u]);
                                }
                            }
                        }

                    }

                }

            }
        }
    }
    for(int j=1; j<=cut; j++)
    {
        for(int k=1; k<=cut; k++)
        {
            ma=max(ma,dp[n][j][k]);
        }
    }
    printf("%d\n",ma);
    return 0;
}