炮兵阵地 POJ1185

炮兵阵地

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions: 10789		Accepted: 3878

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

分析：

第i行的放置情况仅与i-1和i-2行有关，与i-2之前行的放置状态无关，且每次分析第i行的放置过程一样，满足动态规划算法的条件。对于N*M的阵地每个位置有两种状态（部署、不部署），利用状态压缩的DP算法。借助位运算来判断冲突与否的问题。

要注意 &和<<运算符的优先级比较低，一般都加上括号。

源代码如下：

#include <iostream>
using namespace std;

#define MAX 60			//M=10时 所有可能方案数为60
int dp[103][MAX][MAX];     //dp[i][j][k]表示第i行状态为s[j],第i-1行状态为s[k]
int s[MAX];				//对于M列 得出的所有可能状态
int num[MAX];			//对应s[] 每种状态中1的个数

int map[100];			//阵地地图
int realMax;			//对于给定的M，得到的实际方案最大值

int getMaxScheme(int m) //由M获得一行的最大放置方案
{
	int top= 1<<m;
	int i,j=0;
	for(i=0;i<top;i++)
	{
		if((i & (i<<1)) || (i & (i<<2)))continue;
		s[j]=i;
		j++;
	}
	return j;
}

void setOneNumber()
{
	int i,count,t;
	for(i=0;i<realMax;i++)
	{
		count=0;
		t=s[i];
		while(t)
		{
			count++;
			t&=(t-1);
		}
		num[i]=count;
	}
}

inline int max(int a,int b)
{
	if(a>b)return a;
	else return b;
}

int main()
{
	freopen("in.txt","r",stdin);
	int N,M;
	cin>>N>>M;
	memset(map,0,sizeof(map));
	memset(dp,-1,sizeof(dp));

	int i,j,k,t;
	char charRow[11];
	for(i=0;i<N;i++)
	{
		cin>>charRow;
		for(j=0;j<M;j++)
		{
			if(charRow[j]=='H')map[i]=map[i]+ (1<<j);	//P P H P 对应 0 1 0 0 即4 反向倒置的，不过不影响
		}
//		cout<<map[i]<<endl;
	}

	realMax=getMaxScheme(M);
	setOneNumber();

	for(i=0;i<realMax;i++)
	{
		if((s[i]&map[0])==0)dp[0][i][0]=num[i];
	}

	for(i=1;i<N;i++)
	{
		for(j=0;j<realMax;j++)		//枚举第i行
		{
			if(map[i]&s[j])continue;
			for(k=0;k<realMax;k++)	//枚举第i-1行
			{
				if(s[k]&s[j])continue;
				for(t=0;t<realMax;t++)	//枚举第i-2行
				{
					if(s[t]&s[j])continue;
					if(dp[i-1][k][t]==-1)continue;
					dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][t]+num[j]);
				}
			}

		}
	}

	int res=0;

	for(i=0;i<realMax;i++)
	{
		for(j=0;j<realMax;j++)
		{
			res=max(res,dp[N-1][i][j]);
		}
	}

	cout<<res<<endl;

	return 0;
}