【NOI2001】炮兵阵地

【NOI2001】炮兵阵地

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

Solution

=状压dp的经典题目之一。
这道题和其他一些普通的题目不同的是每一个格子对周围的影响范围有两格。
我们先预处理出每一行有可能的情况，并且用十进制来表示集合的状态，然后
对于每一行的每一种状态，枚举前两行的每一种状态，如果状态可以兼容，那
么这一行的累计个数就是上一行累计的炮兵个数加上当前状态的炮兵个数。
最后得出的答案就是最后一行的各种情况中最大的一种。

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
inline int read(){
	char c;int rec=0;
	while((c=getchar())<'0'||c>'9');
	while(c>='0'&&c<='9')rec=rec*10+c-'0',c=getchar();
	return rec;
}
int f[105][105][105];
int a[105][105];
int cnt[105][105];
int st[105],num[105];
int n,m,ans;
char ch[15];
void ST(){
	for(int i=0;i<(1<<m);i++){
		if(i&(i<<1)||i&(i<<2)||i&(i>>1)||i&(i>>2))continue;
		st[++st[0]]=i;
	}
	for(int i=1;i<=st[0];i++)
		for(int j=0;j<m;j++)
		    if(st[i]&(1<<j))num[i]++;
	return ;
}
inline void Pre(int v,int temp){
	for(int i=1;i<=st[0];i++){
		if(temp&st[i])continue;
		a[v][++a[v][0]]=st[i];
		cnt[v][a[v][0]]=num[i];
	}return ;
}
int main(){
	n=read();m=read();ST();
	for(int i=1;i<=n;i++){
		scanf("%s",ch+1);
		int temp=0;
		for(int j=1;j<=m;j++)
			temp=(temp<<1)+(ch[j]=='H');
		Pre(i,temp);
	}
	a[0][0]=1;
	for(int i=1;i<=a[1][0];i++)f[1][i][1]=cnt[1][i],ans=max(ans,f[1][i][1]);
	for(int i=2;i<=n;i++)
		for(int j=1;j<=a[i][0];j++)
			for(int k=1;k<=a[i-1][0];k++)
				if(!(a[i][j]&a[i-1][k]))
					for(int t=1;t<=a[i-2][0];t++)
						if(!(a[i][j]&a[i-2][t])&&!(a[i-1][k]&a[i-2][t])) { 
                            f[i][j][k]=max(f[i][j][k],f[i-1][k][t]+cnt[i][j]);
                            ans=max(ans,f[i][j][k]);
                    }
	cout<<ans;
	return 0;
}