Luogu P2704 [NOI2001] 炮兵阵地

题目描述

司令部的将军们打算在 $N\times M$ 的网格地图上部署他们的炮兵部队。

一个 $N\times M$ 的地图由 $N$ 行 $M$ 列组成，地图的每一格可能是山地（用 $\text{H}$ 表示），也可能是平原（用 $\text{P}$ 表示），如下图。

在每一格平原地形上最多可以布置一支炮兵部队（山地上不能够部署炮兵部队）；一支炮兵部队在地图上的攻击范围如图中黑色区域所示：

如果在地图中的灰色所标识的平原上部署一支炮兵部队，则图中的黑色的网格表示它能够攻击到的区域：沿横向左右各两格，沿纵向上下各两格。

图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。

现在，将军们规划如何部署炮兵部队，在防止误伤的前提下（保证任何两支炮兵部队之间不能互相攻击，即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内），在整个地图区域内最多能够摆放多少我军的炮兵部队。

输入格式

第一行包含两个由空格分割开的正整数，分别表示 $N$ 和 $M$。

接下来的 $N$ 行，每一行含有连续的 $M$ 个字符，按顺序表示地图中每一行的数据。

输出格式

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

样例 #1

样例输入 #1

5 4
PHPP
PPHH
PPPP
PHPP
PHHP

样例输出 #1

6

提示

对于 $100\%$ 的数据，$N\le 100$，$M\le 10$，保证字符仅包含 P 与 H。

状压dp - AC

定义$dp[i][j][k]$，前i行，第i行状态j，第i-1行k，最大炮兵数量。后面不讲。

时间复杂度看似是$O(n{2}^{3m})$，过不了。事实上，先测一下极端数据PPPPPPPPPP会发现，它只有60种可能的合法状态。也就是说，每行的状态最多60种，而非$2^m$。
时间复杂度$O(60^{3}n)$。

int n, m, dp[MAXN][MAXSTA][MAXSTA], ans;
bool g[MAXN][MAXM];
vector<int> s[MAXN];

void dfs(int h, int k, int sta) {
	if (k >= m) {
		s[h].push_back(sta);
		return;
	}
	dfs(h, k + 1, sta);
	if (g[h][k] && (k < 1 || !(sta & (1 << k - 1))) && (k < 2 || !(sta & (1 << k - 2)))) dfs(h, k + 1, sta | (1 << k));
}

int main() {
	setIO("");

	cin >> n >> m;
	for (int i = 1; i <= n; ++i) {
		string tmp; cin >> tmp;
		for (int j = 0; j < m; ++j) {
			if (tmp[j] == 'P') g[i][j] = 1;
		}
	}
	
	s[0].push_back(0);
	for (int i = 1; i <= n; ++i) {
		dfs(i, 0, 0);
//		cout << s[i].size() << "\n";
	}
	
	for (int i = 0; i < s[1].size(); ++i) {
		dp[1][i][0] = __builtin_popcount(s[1][i]);
		if (i == n) ckmax(ans, dp[1][i][0]);//不要忘记n=1的情况
	}
	for (int i = 2; i <= n; ++i) {
		for (int j = 0; j < s[i].size(); ++j) {
			for (int k = 0; k < s[i - 1].size(); ++k) {
				if (!(s[i][j] & s[i - 1][k])) {
					for (int l = 0; l < s[i - 2].size(); ++l) {
						if (!(s[i][j] & s[i - 2][l])) {
							ckmax(dp[i][j][k], dp[i - 1][k][l] + __builtin_popcount(s[i][j]));
						}
					}
					if (i == n) ckmax(ans, dp[i][j][k]);
				}
			}
		}
	}
	cout << ans << "\n";
	
	return 0;
}