POJ 1185--炮兵阵地（状压dp）

来源：

题意：给出一个n*m的矩阵，矩阵的每个方格标有P/H p表示可以安置大炮，H表示不能安置大炮，当大炮安置于(i,j)点时，其左右两个单位以及上下两个单位都在攻击范围，求在两支大炮不会相互攻击的前提下，最多能够安置大炮的数量。

分析：

当前行大炮的的安置要受其前两行的影响，所以状态转移方程有：

dp[i][j][k] = max(dp[i][j][k],dp[i - 1][k][l] + sum[j]) dp[i][j][k]表示第i行的状态为j第i - 1的状态为k sum[j]表示该行取j状态时可能增加的数量；

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

int n,m,t;
int map[105],s[105],num_1[105];//s[]存储了所有的合法状态；num_1[]是每个合法状态中炮兵的个数
int dp[105][105][105];

int dfsdp(int n,int now,int last)//记忆化搜索，也就是有记录的深搜
{
	int i,k,max=0;
	int judge=s[now]|s[last];//judge的每位记录了下行的该位置是否能放炮兵，因为本行放过和上行放过的地方显然不能放

	if(dp[n][now][last]!=-1)//如果被搜过，直接返回，记忆化搜索的思想
		return dp[n][now][last];

	if(n==1)
		return dp[n][now][last]=num_1[now];

	for(i=0;i<t;i++)//对所有可能的方法往下继续搜
	{
		if(!(s[i]&judge||map[n-2]&s[i]))
		{
			k=dfsdp(n-1,last,i);
			if(max<k)
			{
				max=k;
			}
		}
	}

	dp[n][now][last]=max+num_1[now];
	return dp[n][now][last];
}

int main()
{
	int i,j;
	char temp[11];

	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++)
	{
		scanf("%s",temp);
		map[i]=0;

		for(j=0;j<m;j++)//将地图转化成2进制数(H=1,P=0)，并用十进制数存储
			if(temp[j]=='H')
				map[i]+=1<<j;
	}

	memset(dp,-1,sizeof(dp));

    t=0;
	for(i=0;i<(1<<m);i++)//计算合法状态以及合法状态中的1的个数
	{
	    //cout<<((i<<1)&i)<<" "<<((i<<2)&i)<<" "<<((i>>1)&i)<<" "<<((i>>2)&i)<<" i "<<i<<endl;
	    if(!((i<<1)&i)&&!((i<<2)&i)&&!((i>>1)&i)&&!((i>>2)&i))
		{
		    //cout<<i<<" iii"<<endl;
			s[t]=i;
			int tmp=i;
//			num_1[t]=tmp&1;
//			while(tmp=tmp>>1)
//			{
//			    num_1[t]+=tmp&1;
//			}
//			cout<<num_1[t]<<" num "<<endl;
			for(j=i;j>0;j-=j&-j)
				num_1[t]++;
			t++;
		}
	}
    printf("%d\n",dfsdp(n+2,0,0));//开搜，从n+2开始，但是那额外2行什么都不放
	return 0;
}