POJ1185 炮兵阵地 状压dp

题目描述：

司令部的将军们打算在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，与HDU4539非常类似。题目大意是在一块n*m的地上，P表示可以放置大炮，H表示障碍。一个大炮能攻击其前后左右距离为2的位置（包括自己的位置），求在不会互相攻击的前提下最多放置的大炮数量。
具体看代码 代码注释有解释。

代码如下：

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <string>
#include <math.h>
#include <stdlib.h>
#define max(a,b) (a)>(b)?(a):(b)
using namespace std;
typedef long long ll;

int n,m,k;
int sum[105];//sum[i]表示在第i状态下的大炮放置数量
int dp[105][60][60];//dp[i][j][k]表示的是在i行j状态，i-1行k状态的值
//60这个数字的意思是在一行中，所可能出现的状态个数不多于60个
int mp[105];//mp[i]表示第i行的输入状态
int state[60];//state[i]表示第i个状态的表示（二进制转化十进制存储）

bool ok(int x)//去掉横向干扰的
{
    if (x & (x<<1)) return false;//排除有连续两个1
    if (x & (x<<2)) return false;//排除101,111这些间隔一个有同有1的
    return true;
}

int get1(int x)//1的个数
{
    int ret=0;
    while(x)
    {
        if (x&1) ret++;
        x>>=1;
    }
    return ret;
}

void find()//无障碍情况下找出所有合法状态
{
    memset(state,0,sizeof(state));
    for(int i=0; i<(1<<m); i++)
    {
        if (ok(i))
        {
            state[k]=i;
            sum[k]=get1(i);
            k++;
        }
    }
}

void DP()
{
    for(int i=0; i<k; i++)//i=0时 dp值和sum值相同
    {
        if (!(state[i] & mp[0])) dp[0][i][0]=sum[i];
    }
    for(int i=1; i<n; i++) //1到n-1行
    {
        for(int j=0; j<k; j++) //i行状态
        {
            if (state[j] & mp[i]) continue;//非法状态
            for(int p=0; p<k; p++) //i-1行状态
            {
                if (state[j] & state[p]) continue;//i行状态与i-1行状态是否合理
                for(int q=0; q<k; q++)//i-2行状态
                {
                    if (state[q] & state[j]) continue;//i行状态与i-2行状态是否合理
                    if (dp[i-1][p][q]==-1) continue;//排除非法状态
                    dp[i][j][p]=max(dp[i][j][p],dp[i-1][p][q]+sum[j]);
                    //i行的j状态，i-1行p状态dp值由i-1行p状态，i-2行q状态dp值更新
                }
            }
        }
    }
}

int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        memset(mp,0,sizeof(mp));
        memset(dp,-1,sizeof(dp));
        char tmp;
        for(int i=0; i<n; i++)
            for(int j=0; j<m; j++)
            {
                cin>>tmp;
                if (tmp=='H') mp[i]=mp[i]*2+1;//H置为1,P为0
                else mp[i]*=2;
            }
        k=0;
        find();
        DP();
        int ans=0;
        for(int i=0; i<k; i++)
            for(int j=0; j<k; j++)
                ans=max(ans,dp[n-1][i][j]);
        printf("%d\n",ans);
    }
    return 0;
}
//再着重解释一下dp[105][60][60]的60
//因为其攻击范围的原因，一行最多10格，最多放置4个。
//放置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种。
//所以为了减少空间消耗，60种即可。