POJ 1185 炮兵阵地

题目链接：http://poj.org/problem?id=1185

炮兵阵地

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions: 26021		Accepted: 10042

Description

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

Source

Noi 01

思路：算是一道入门的状压DP题目吧。我是看完我们老师的ppt后写的，如下所示。

多斟酌几次就会明白了。

附上AC代码：

#include <cstdio>
#include <cstring>
//#pragma comment(linker, "/STACK:102400000, 102400000")
using namespace std;
const int maxn = 105;
const int maxm = 60;
int mat[maxn], sta[maxm], one[maxm];
int dp[maxn][maxm][maxm], lmt[10];
char str[15];
int n, m;

int cnt(int x){
	int ans = 0;
	while (x){
		ans += x&1;
		x >>= 1;
	}
	return ans;
}

bool check(int x){
	if (x & (x<<1))
		return false;
	if (x & (x<<2))
		return false;
	return true;
}

void init(){
	int j=1, k=0;
	for (int i=0; i<(1<<10); ++i)
		if (check(i)){
			sta[k] = i;
			one[k] = cnt(i);
			if (i > ((1<<j)-1)){
				lmt[j-1] = k;
				++j;
			}
			++k;
		}
	lmt[9] = maxm;
}

int main(){
	#ifdef LOCAL
	freopen("input.txt", "r", stdin);
	freopen("output.txt", "w", stdout);
	#endif
	init();
	while (~scanf("%d%d", &n, &m)){
		for (int i=1; i<=n; ++i){
			scanf("%s", str);
			mat[i] = 0;
			for (int j=0; j<m; ++j)
				mat[i] = mat[i]*2+(str[j]=='H');
		}
		m = lmt[m-1];
		memset(dp, 0, sizeof(dp));
		for (int j=0; j<m; ++j)
			if (!(mat[1] & sta[j]))
				for (int k=0; k<m; ++k)
					dp[1][j][k] = one[j];
		for (int i=2; i<=n; ++i)
			for (int j=0; j<m; ++j)
				for (int k=0; k<m; ++k)
					for (int s=0; s<m; ++s){
						int t = sta[k]|sta[s]|mat[i];
						if (t & sta[j])
							continue;
						t = dp[i-1][k][s]+one[j];
						if (t > dp[i][j][k])
							dp[i][j][k] = t;
					}
		int ans = 0;
		for (int j=0; j<m; ++j)
			for (int k=0; k<m; ++k)
				if (dp[n][j][k] > ans)
					ans = dp[n][j][k];
		printf("%d\n", ans);
	}
	return 0;
}