C - 炮兵阵地（压状dp+三维数组）

题目

司令部的将军们打算在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

解释

使用maze数组记录每行地图的二进制地形形态， 使用status和num数组记录m长度下每一种符合1相邻两个0的地形信息,stauts记录符合条件二进制地形形态， num数组记录相对应的status中1的个数（即该地形下的炮台放置位置）。
使用三维数dp[i] [j] [k],第一维为第几行， 第二维代表第i行状态， 第三维代表第i-1行状态。
转移方程为 dp[i] [j] [k] = max(dp[i ][j] [k], dp[i-1] [k ][t] + num[j]);
其中j状态要符合第i行地形， 即maze[i]&status[j] == 0, 说明j状态下的1 都对应在maze地形上的0位置。
重点是dp[i-1] [k ][t]， k代表i的上一行状态， t代表i的上一行的上一行的状态，我们的目的就是找出所有符合条件的k,t， 然后更新 dp[i] [j] [k] 。条件为(!(status[j]&status[k]) && !(maze[i-1] & status[k])) 和 f(!(status[t]&status[k])&&!(status[t]&status[j]) && !(status[t]& maze[i-2]))。 我们要先更新dp[1] [i][0], 更新图中第一行， 由于dp[0][][]图中其实是不存在的， 但是初始化都为0， 对第二行按照假装存在运行，是不会产生影响的。

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#define mod 100000000
#define mem(x) memset(x, 0, sizeof(x))
using namespace std;
int n, m;
int dp[105][70][70];
int maze[105];
int status[70];
int num[70];
int main(){
    int i, j;
    char s[15];
    while(~scanf("%d %d", &n, &m)){
        mem(dp);
        mem(maze);
        mem(status); mem(num);
        for(i = 1; i <= n; i++){
                scanf("%s", s);
            for(j = 0; j < m; j++)
                if(s[j] == 'H')
                    maze[i] += (1<<j);
        }

        int len = 0, k, sumn;
        for(i = 0; i < (1<<m); i++){
            sumn= 0;
            if(!(i & (i<<1)) && !(i &(i<<2))){
                k = i;
                while(k){
                    if(k&1)
                        sumn++;
                    k >>= 1;
                }
                num[len] = sumn;
                status[len++] = i;
            }
        }

        int t;
        for(i = 0; i < len; i++)
            if(!(maze[1]&status[i]))
                dp[1][i][0] = num[i];

        for(i = 2; i <= n; i++){
            for(j = 0; j < len; j++){
                if(!(maze[i] & status[j])){
                    for(k = 0; k < len; k++)
                        if(!(status[j]&status[k]) && !(maze[i-1] & status[k]))
                           for(t = 0; t < len; t++)
                                if(!(status[t]&status[k])&&!(status[t]&status[j]) && !(status[t]& maze[i-2]))
                                    dp[i][j][k] = max(dp[i][j][k], dp[i-1][k][t] + num[j]);
                }
            }
        }
        int maxx = 0;
        for(i = 0; i < len; i++)
            for(j = 0; j < len; j++)
                maxx = max(maxx, dp[n][i][j]);
        printf("%d\n", maxx);
    }
    return 0;
}