炮兵阵地 POJ - 1185 状压dp

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

这是提交十三天后写的 

打开代码一声卧槽 我当初都写了啥？的懵逼感。花了数分钟理清了关系

首先初步枚举横向的合法状态存在s数组里面

bitc是用来计算一个状态位的二进制1的个数

judge是用来判断 这块地 能不能合法的放下sta状态

dp【i】【j】【k】方程的意思是  第i行 的状态为j 第i-1行状态为j时最优的个数

dp[i][t][j]=max(dp[i][t][j],dp[i-1][k][t]+num[j]);

由i-1层向i层转移 但携带了很多状态 相对于Corn Fields要6很多

其他的在代码里说明

 

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#define de cout<<"debug"<<endl;
using namespace std;
int land[110];
int s[500],num[500];
int dp[110][100][100];
char map[110][100];
int m,n;
int top;

void init()
{
	top=0;
	for(int i=0;i<(1<<n);i++)
	{
		if((i&(i<<1))==0&&(i&(i<<2))==0) s[++top]=i;
	}
}

int bitc(int x)
{
	int cnt=0;
	while(x)
	{
		cnt++;
		x&=(x-1);
	}
	return cnt;
}

int judge(int sta,int x)
{
	if(land[x]&sta) return 0;
	return 1;
}

int main()
{
	while(scanf("%d%d",&m,&n)!=EOF)
	{
		init();
		//cout<<top<<endl;
		memset(dp,-1,sizeof(dp));
		for(int i=1;i<=m;i++) scanf("%s",map[i]+1);		
		for(int i=1;i<=m;i++)
		{
			land[i]=0;
			for(int j=1;j<=n;j++)
			{
				if(map[i][j]=='H') land[i]+=(1<<(j-1));//不能放的存为1 这样在judge里面 和状态sta按位与一下 出现1就可以表示不合法辣
			}
		}
		
		for(int i=1;i<=top;i++)
		{
			num[i]=bitc(s[i]);//每个状态放的炮数
		}		
		for(int i=1;i<=top;i++)
		{
			if(judge(s[i],1)) dp[1][1][i]=num[i];
		}
		//for(int i=1;i<=top;i++) cout<<dp[1][1][i]<<endl;
		for(int i=2;i<=m;i++)
		{
			for(int j=1;j<=top;j++)
			{
				if(!judge(s[j],i)) continue;//找到合法的i层状态
				for(int k=1;k<=top;k++)
				{
					if(s[j]&s[k]) continue;//找到合法的i-2层状态
					for(int t=1;t<=top;t++)
					{
						if(s[j]&s[t]) continue;//合法的i-1层状态
						if(dp[i-1][k][t]==-1) continue;
						dp[i][t][j]=max(dp[i][t][j],dp[i-1][k][t]+num[j]);
					}
				}								
			}
		}
		int ans=0;		
		for(int i=1;i<=m;i++)
		for(int j=1;j<=top;j++)
		for(int k=1;k<=top;k++) 
		ans=max(ans,dp[i][j][k]);
		printf("%d\n",ans);
	}
}