POJ - 1185 - 炮兵阵地（dp & 状压）

#POJ - 1185 - 炮兵阵地
司令部的将军们打算在NM的网格地图上部署他们的炮兵部队。一个NM的地图由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的好题，我刚开始实在不懂，就各种看题解，这里先分享一个题解，写的很棒。

这里记一下自己的理解，所谓状压，就是把一个状态进行压缩成为一个数字能够以数字的形式进行存储，这里是一个地图，有平原，有山地两种地形，平原地方可以放炮兵，那么就定为平原为0，山地为1，那么每一行都可以用一个数字来存储。又因为每一个炮兵都有攻击范围，所以我们用一个数组来记录一下一行内能够使得炮兵两两不能攻击的排布方式。然后就是每一行枚举并且dp了，这里要用一个三维数组，即dp[i][j][k]，其中 i 是指当前处于哪一行，j是当前这一行所处的状态，也就是我们刚刚预处理过的使得炮兵两两不能攻击的排布方式，k是 i - 1 行，即上一行所处的状态（可行的分布方式）。我们还需要记录每一种状态下的炮兵数量。这样我们需要枚举当先行，上一行，上上行的所有状态，然后找出所有可行状态下最多的炮兵数。
AC代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
//行列数和状态数
const int maxr = 110, maxc = 15, maxm = 70;
int row, col, nums; //行列，两个炮兵互不攻击的条件下符合条件的状态个数
int base[maxr]; //将第i行的地图压缩成一个状态
int state[maxm];    //两个炮兵互不攻击的条件下，符合条件的状态
int soldier[maxm];  //对应state[i]状态下能放多少炮兵
int dp[maxr][maxm][maxm];   //dp[i][j][k] 表示第i行状态为state[j]，第i-1行状态为state[k]时的最优解
char g[maxr][maxc];

int main()
{
    memset(base, 0, sizeof(base));
    memset(state, 0, sizeof(state));
    memset(soldier, 0, sizeof(soldier));
    memset(dp, 0, sizeof(dp));
    nums = 0;

    scanf("%d%d", &row, &col);
    for(int i = 0; i < row; i++)
    {
        scanf("%s", &g[i]);
        for(int j = 0; j < col; j++)
            if(g[i][j] == 'H')  base[i] += 1 << j;  //计算原地图的状态
    }

    for(int i = 0; i < (1 << col); i++) //仅是两个炮兵不互相攻击的条件下计算所有状态
    {
        if(i & (i << 1) || i & (i << 2))    continue ;  //i这个状态出现了士兵两两攻击
        int k = i;
        while(k) //这个循环计算状态i的二进制形式里面有多少个1，也就是放了多少个士兵
        {
            soldier[nums] += k & 1; //等价于k%2，相当于判断k的二进制形式里面有多少个1
            k = k >> 1;
        }
        state[nums++] = i; //保存这个合法的状态
    }

    for(int i = 0; i < nums; i++)   //先初始化dp[0][i][0]，即初始化第1行的情况
    {
        if(state[i] & base[0]) continue;
        //在state[i]的基础上，还要满足士兵不能放在山上，这个判断就是处理这个问题的
        dp[0][i][0] = soldier[i];
    }

    for(int r = 1; r < row; r++)
        for(int i = 0; i < nums; i++)   //枚举第r行的状态
        {
            if(state[i] & base[r])  continue ;
            for(int j = 0; j < nums; j++) //枚举第r-1行的状态
            {
                if(state[j] & base[r-1]) continue;
                if(state[i] & state[j]) continue;   //第r行的士兵和第r-1行的士兵相互攻击
                for(int k = 0; k < nums; k++) //枚举第r-2行的状态
                {
                    if(state[k] & base[r-2]) continue;
                    if(state[j] & state[k]) continue;//第r-1行的士兵和第r-2行的士兵相互攻击
                    if(state[i] & state[k]) continue;//第r行的士兵和第r-1
                    dp[r][i][j] = max(dp[r][i][j], dp[r-1][j][k] + soldier[i]);
                }
            }
        }
    int ans = 0;
    for(int i = 0; i < nums; i++)
        for(int j = 0; j < nums; j++) //枚举dp[row-1][i][j]
           ans = max(ans, dp[row - 1][i][j]);   //找最后一行所能达到的最大值
    printf("%d\n", ans);
    return 0;
}