司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成，地图的每一格可能是山地（用"H" 表示），也可能是平原（用"P"表示），如下图。在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示： 

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。
现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。 

第一行包含两个由空格分割开的正整数，分别表示N和M； 
接下来的N行，每一行含有连续的M个字符('P'或者'H')，中间没有空格。按顺序表示地图中每一行的数据。N <= 100；M <= 10。

仅一行，包含一个整数K，表示最多能摆放的炮兵部队的数量。

5 4
PHPP
PPHH
PPPP
PHPP
PHHP

6

#include <iostream>
#include <bitset>
#include <algorithm>
using namespace std;int N, M;
const int MAXN = 110;
int map[MAXN];
int state[MAXN][MAXN], ones[MAXN][MAXN], nums[MAXN];
int dp[MAXN][MAXN][MAXN];
const int INF = 0X3F3F3F3F;bool conflict(const int &s1, const int &s2) {return (s1 & s2);
}
/** 已知当前状态为 s, 能否在 pos 处摆放一门大炮* 由递推关系知, s 中摆放的大炮位置都在 pos 的前面*   s == ..1001, 表示第 1, 4 个位置上有炮台*/
bool underAttack(const int &pos, const int s) { if(pos == 0) {//首列, 肯定能放return false;}else if(pos == 1) {if(s&1)//若首列已放则第二列不能再放return true;}else if(pos >= 2) {if(s & (1<<(pos-1)))//倒数第一列有没有放return true;if(s & (1<<(pos-2)))//倒数第二列有没有放return true;}return false;
}/**    枚举第 line 行摆放大炮的可行方案并统计每种方案下的大炮数目*   cs, line 是全局变量, 就不需要通过参数传入函数了 *      cs 记录已经有多少种摆放方案*/
int cs, line;
void enumCannon(const int &st, const int s) {if(st == M) {state[line][cs++] = s;return ;}if(!(map[line] & (1<<st)) && !underAttack(st, s)) { // 是平原且不在其他炮台的工具范围内int tmp = s;tmp = (s|1<<st); enumCannon(st+1, s);enumCannon(st+1, tmp);}else{enumCannon(st+1, s);}
}/**    dp[1] 单独初始化*    dp[r][i][j], i, j 都是排列的 index, 而不是状态本身, 这样就能保证 i, j < 60*/
void pre_process() {// 初始化 state 数组for(line = 1; line <= N; line ++) {cs = 0;enumCannon(0, 0);sort(state[line], state[line]+cs);nums[line] = cs;for(int i = 0; i < cs; i ++) {bitset<10> bs(state[line][i]);ones[line][i] = bs.count(); // ones 记录有多少个炮台}}// 对第 0 行的炮台进行初始化memset(state[0], 0, sizeof(state[0]));memset(ones[0], 0, sizeof(ones[0]));nums[0] = 60;//    dp[1]初始化, 没有任何限制, 随便摆for(int i = 0; i < nums[1]; i ++) {for(int j = 0; j < nums[0]; j++) {dp[1][i][j] = ones[1][i]; // 能放多少放多少, 不需要考虑前几行对当前行的影响}                       }
}
/** i, j, k 分别是指 dp[i][j][k], i 表示行, j 比较第 i 行有可能可行的炮台摆放方案索引, k 是第 i-1 行的摆放方案*  dp[i-1][k][s] 是已知*/
int mainFunc() {for(int i = 2; i <= N; i ++) {for(int j = 0; j < nums[i]; j ++) {for(int k = 0; k < nums[i-1]; k ++) {if(conflict(state[i][j], state[i-1][k])) // i, i-1 行方案冲突continue;for(int s = 0; s < nums[i-2]; s ++) {    // i, i-2 行方案是否冲突if(conflict(state[i][j], state[i-2][s]))continue;dp[i][j][k] = max(dp[i][j][k], dp[i-1][k][s]+ones[i][j]);}}}}int ans = 0;if(N == 1) {for(int i = 0; i < nums[N]; i ++) {ans = max(ans, dp[N][i][0]);}return ans;}for(int i = 0; i < nums[N]; i ++) {for(int j = 0; j < nums[N-1]; j++) {ans = max(ans, dp[N][i][j]);}}return ans;
}
int main() {freopen("E:\\Copy\\ACM\\poj\\1185\\in.txt", "r", stdin);while(cin >> N >> M) {for(int i = 1; i <= N; i ++) {char c;getchar();for(int j = 0; j < M; j ++) { // 从左往右压scanf("%c", &c);if(c == 'H')map[i] = (map[i]|(1<<j)); // 技巧}}pre_process();cout << mainFunc() << endl;}return 0;
}