Poj1185状态压缩

因为行数很少，所以我们可以按行来定义状态。那么每一行的状态就只与它的上两行有关。由于只有放于不放，所以可以用二进制数表示状态。

那么对于列呢，我们可以先预处理出一个数组表示对于列数，有多少种组合,保存在stk[]中，（不考虑平地与高原，只考虑两两不互轰），用num[]表示对应的组合里有多少个1，即部署了多少个炮兵。用cur[]记录下图的原始状态。

       状态转移为:dp[r][i][j]=max(dp[r][i][j],dp[r-1][j][k]+num[i]); (对于第r行，它的r-1行为状态i，它的r-2行为状态j时对应的布置炮兵的最大数量）

       状态判断是否可行，直接两个状态&就行。

#include<iostream>
using namespace std;
int stk[70],num[70],cur[105],dp[105][70][70];
int n,m,tot;
bool judge(int a,int b)
{
	return a&b;
}
int max(int a,int b)
{
	return a>b?a:b;
}
void deal()
{
	for(int i=0;i<(1<<m);i++)
	{
		if(judge(i,i<<1) || judge(i,i<<2)) continue;
		int k = i;
		while(k)
		{
			num[tot]+=k&1;
			k=k>>1;
		}
		stk[tot++]=i;
	}
}
int main()
{
	char c;
	int ans;
	cin>>n>>m;
	tot=0;
	for(int i=0;i<n;i++)
	 for(int j=0;j<m;j++)
	 {
	 	cin>>c;
	 	if(c=='H') cur[i]+=(1<<j);
	 }
	 
	 deal();
	 
	 //初始化 边界 
	 for(int i=0;i<tot;i++)
	 {
	 	if(judge(stk[i],cur[0])) continue;
	 	dp[0][i][0]=num[i];	
	 }
	 for(int i=0;i<tot;i++)
	 {
	 	if(judge(stk[i],cur[1])) continue;
	 	for(int j=0;j<tot;j++)
	 	{
	 		if(judge(stk[j],cur[0])) continue;
	 		if(judge(stk[i],stk[j])) continue;
	 		dp[1][i][j]=max(dp[1][i][j],dp[0][j][0]+num[i]);
	 	}
	 }
	 
	 for(int r=2;r<n;r++)
	 {
	  	for(int i=0;i<tot;i++)
	  	{
	  		if(judge(stk[i],cur[r])) continue;
	   	 	for(int j=0;j<tot;j++)
	   	 	{
	   	 		if(judge(stk[j],cur[r-1])) continue;
	   	 		if(judge(stk[i],stk[j])) continue;
		  	 	for(int k=0;k<tot;k++)
		  		{
		  		 	if(judge(stk[k],cur[r-2])) continue;
		  		 	if(judge(stk[k],stk[i])) continue;
		  		 	if(judge(stk[k],stk[j])) continue;
		  		 	dp[r][i][j]=max(dp[r][i][j],dp[r-1][j][k]+num[i]);	 	
		  		}
		  	}
		}
	}
	ans=-100;
	for(int i=0;i<tot;i++)
	 for(int j=0;j<tot;j++)
	 	ans=max(ans,dp[n-1][i][j]);
	cout<<ans<<endl;
	return 0;
}