POJ 1185 炮兵布阵 （经典状压DP）

题意：一个n*m矩阵(n<=100， m<=10)，每个方格可以放大炮用0表示，不可以放大炮用1表示（原题用字母），让放最多的大炮，大炮与大炮间不会

互相攻击。每个炮的攻击范围是上下左右两个距离。

思路：（点击打开链接）可以发现，对于每一行放大炮的状态，只与它上面一行和上上一行的状态有关，每一行用状态压缩的表示方法，0表示不放大

炮，1表示放大炮，同样的，先要满足硬件条件，即有的地方不能放大炮，然后就是每一行中不能有两个1的距离小于2（保证横着不互相攻击），这些

要预先处理一下。然后就是状态表示和转移的问题了，因为是和前两行的状态有关，所以要开个三维的数组来表示状态，当前行的状态可由前两行的状

态转移而来。即如果当前行的状态符合前两行的约束条件（不和前两行的大炮互相攻击），则当前行的最大值就是上一个状态的值加上当前状态中1的

个数（当前行放大炮的个数） 

dp[i][curStatus][preStatus]:第i行状态为curStatus，i-1行状态为preStatus状态时的最多炮兵数

状态转移方程：dp[i][j][k] =max(dp[i][j][k],dp[i-1][k][t]+num[j]); num[j]为j状态中1的个数

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn = 105;
int dp[maxn][maxn][maxn], maze[maxn], status[maxn], num[maxn];
//status:同一行两个炮兵不相互攻击的状态
//maze:地形的限制
//dp[i][curStatus][preStatus]:第i行状态为curStatus，i-1行状态为preStatus状态时的最多炮兵数
//num[i]:一行状态为status[i]的炮兵数（1的数量）
char ss[maxn];
int n, m;

int main(void)
{
    while(cin >> n >> m)
    {
        memset(dp, 0, sizeof(dp));
        memset(maze, 0, sizeof(maze));
        for(int i = 1; i <= n; i++)
        {
            scanf(" %s", ss);
            for(int j = 0; j < m; j++)
                if(ss[j] == 'H')
                    maze[i] += (1<<j);
        }
        int len = 0;
        for(int i = 0; i < (1<<m); i++)
            if(!(i&(i<<1)) && !(i&(i<<2)))
            {
                int t = i, cnt = 0;
                while(t)
                {
                    if(t%2) cnt++;
                    t /= 2;
                }
                num[len] = cnt;
                status[len++] = i;
            }
        for(int i = 0; i < len; i++)
            if(!(maze[1]&status[i]))
                dp[1][i][0] = num[i];   //第一行赋初值
        for(int i = 2; i <= n; i++)
        {
            for(int j = 0; j < len; j++)    //枚举第i行状态
            {
                if(status[j]&maze[i]) continue;
                for(int k = 0; k < len; k++)    //枚举第i-1行状态
                {
                    if(status[j]&status[k]) continue;//i-1与i行上下是否互攻
                    if(status[k]&maze[i-1]) continue;//i-1地形是否符合
                    for(int t = 0; t < len; t++)    //枚举第i-2行状态
                    {
                        if(status[t]&status[j]) continue;//与第i行不能冲突
                        if(status[t]&status[k]) continue;//与第i-1行不能冲突
                        if(status[t]&maze[i-2]) continue;//与地形不能冲突
                        dp[i][j][k] = max(dp[i][j][k], dp[i-1][k][t]+num[j]);
                    }
                }
            }
        }
        int ans = 0;
        for(int i = 0; i < len; i++)
            for(int j = 0; j < len; j++)
                ans = max(ans, dp[n][i][j]);
        printf("%d\n", ans);
    }
    return 0;
}