[poj1185][NOI2001]炮兵阵地 状压dp

                                                        炮兵阵地

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions:0		Accepted:0

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

预处理每一行可以有的状态

dp[i][j][k]:前i行，i行状态为j，i-1行状态为k的最大解（注意状态是编号不是二进制押位，因为这道题就算全是平原也没有多少种所以要开编号不然目测会mle（然而我开的很大））

然后这不是水题吗？（开始看成求方案数想问怎么没有模数写完样例输出400多的zz）

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <bitset>
#include <queue>
#include <ctime>
#include <map>
#include <set>
using namespace std;
typedef long long ll;
int q[105][1025], a[105], n, m, ans, dp[105][505][505];
char ch[105][15];
void init(){
	for( int j = 1; j <= n; j++ )
		for( int i = 0; i < (1<<m); i++ )
			if( !(i&(i<<1)) && !(i&(i<<2)) )
				if( !(i&a[j]) )
					q[j][++q[j][0]] = i;
}
int count( int x ){ int cnt = 0; while(x){ cnt++; x -= x&-x; } return cnt; }
int main(){
	scanf( "%d%d", &n, &m );
	for( int i = 1; i <= n; i++ ) scanf( "%s", ch[i] );
	for( int i = 1; i <= n; i++ ) for( int j = 0; j < m; j++ ) if( ch[i][j] == 'H' ) a[i] += (1<<j);
	init();
	if( n == 1 ){
		for( int i = 1; i <= q[1][0]; i++ ) ans = max( ans, count(q[1][i]) );
		cout<<ans<<endl;
		return 0;
	}
	for( int i = 1; i <= q[2][0]; i++ )
		for( int j = 1; j <= q[1][0]; j++ )
			if( (q[1][j]&q[2][i]) == 0 )
				dp[2][i][j] = count(q[1][j]) + count(q[2][i]);
	for( int i = 3; i <= n; i++ )
		for( int j = 1; j <= q[i][0]; j++ )
			for( int k = 1; k <= q[i - 1][0]; k++ )
				if( (q[i][j]&q[i - 1][k]) == 0 )
					for( int l = 1; l <= q[i - 2][0]; l++ )
						if( (q[i][j]&q[i - 2][l]) == 0 && (q[i - 1][k]&q[i - 2][l]) == 0 )
							dp[i][j][k] = max( dp[i][j][k], dp[i - 1][k][l] + count(q[i][j]) );
	for( int i = 1; i <= q[n][0]; i++ )
		for( int j = 1; j <= q[n - 1][0]; j++ )
			ans = max( ans, dp[n][i][j] );
	cout<<ans<<endl;
	return 0;
}
/*
5 4
PHPP
PPHH
PPPP
PHPP
PHHP
*/