炮兵阵地 （状压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

题意：如图是炮的攻击范围，在H处不能放置炮，问在每一个炮都不冲突的情况下，最多能放多少个炮。（如果两个炮在互相的攻击范围则冲突）

题解：状压dp，预处理每一个行能放置炮的所有状态，并记录每一个状态对应的炮的数量。

dp方程：dp[row][i][j] = min(dp[row-1][j][k] + Num[row][i] , dp[row][i][j]);

dp[row][i][j] 表示在第row行状态为i， 上一行状态为j 的最大数量， 需要row，i，j，k四重循环。

注意：前两行要单独处理。

最后附上AC代码

//#include<bits/stdc++.h>
//#include <unordered_map>
//#include<unordered_set>
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>
#include<climits>
#include<queue>
#include<cmath>
#include<stack>
#include<map>
#include<string>
using namespace std;
#define LL long long
#define ULL unsigned long long
#define MT(a,b) memset(a,b,sizeof(a))
const int INF  =  0x3f3f3f3f;
const int O    =  1e5;
const int mod  =  1e9+7;
const int maxn =  1e2+5;
const double PI  =  3.141592653589;
const double E   =  2.718281828459;

int n, m ;
int dp[maxn][maxn][maxn];
vector<int>NextSta[maxn];
vector<int>NextNum[maxn];
int getH[maxn];

int getNum(int sta)
{
    int ans = 0 ;
    while(sta)
    {
        if(sta&1) ans ++;
        sta >>= 1;
    }
    return ans ;
}

void pret(string s, int row)
{
    for(int sta=0;sta<=(1<<m)-1;sta++)
    {
        if(sta & getH[row]) continue; //判断是否与H冲突
        if(sta & sta<<1)    continue;
        if(sta & sta<<2)    continue; //判断是否与自己的状态冲突
        NextSta[row].push_back(sta);  //记录每一行能布置炮的所有状态
        NextNum[row].push_back(getNum(sta)); //记录每种能布置炮的状态所对应的炮的数量
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    MT(getH, 0);
    for(int row = 0;row<n;row++)
    {
        char s[maxn]; scanf("%s",s);
        NextNum[row].clear();
        NextSta[row].clear();
        for(int i=0;i<m;i++) if(s[i]=='H') getH[row] |= 1<<i;//记录每一行H的状态
        pret(s,row); //预处理
    }
    
    MT(dp, -1);
    for(int i=0;i<NextSta[0].size();i++) dp[0][i][0] = NextNum[0][i]; //特殊处理第一排
    for(int i=0;i<NextSta[1].size();i++) for(int j=0;j<NextSta[0].size();j++) //特殊处理第二排
    {
        if(NextSta[1][i] & NextSta[0][j]) continue;
        //if(dp[0][j][0] == -1) continue;
        dp[1][i][j] = max(dp[1][i][j], dp[0][j][0] + NextNum[1][i]);
    }
    
    for(int row=2; row<n; row++)
    for(int i=0; i<NextSta[row].size(); i++)
    for(int j=0; j<NextSta[row-1].size(); j++)
    {
        if(NextSta[row-1][j] & NextSta[row][i]) continue;//判段是否与上一行冲突
        for(int k=0; k<NextSta[row-2].size(); k++)
        {
            if(NextSta[row-2][k] & NextSta[row][i]) continue;//判断是否与上面第二行冲突
            //if(dp[row-1][j][k] == -1) continue;
            dp[row][i][j] = max(dp[row][i][j], dp[row-1][j][k] + NextNum[row][i]);//重要dp方程
        }
    }
    
    int ans = 0;
    if(n<=1) for(int i=0;i<NextNum[0].size();i++) ans = max( ans , dp[0][i][0]);
    else for(int i=0;i<NextSta[n-1].size();i++)  for(int j=0;j<NextSta[n-2].size();j++)
        ans = max(ans, dp[n-1][i][j]);
    
    printf("%d\n",ans);
    return  0;
}