POJ 1185 炮兵阵地

一道非常典型的状态压缩题，对于熟练者来说，绝对算是个水题，但是对于状态压缩入门者，这算是一个好题，能让我们对状态压缩有深刻理解。

解题：

1.先找出一行的符合状态，进行01转换。

2.因为是和前2行有关的所以建立一个三位数组dp[i][j][k] 表示第i行 第j个状态 前一个状态是为k的 最大的炮兵数。

          状态转移方程dp[i][j][k]=max(dp[i][j][k],(num[j]+dp[i-1][k][t]));

3.预先处理第一第二行的状态，方便处理。

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int m,n;
int dp[105][100][100],cur[105],state[105],num[105];
int top;
int max(int a,int b)
{
	return a>b?a:b;
}
int num1_number(int a)
{
    int count=0;
    while(a){
        if(a%2!=0) count++;
        a=a/2;
    }
    return count;
}
void Init()
{
	int total=1<<n;
	for(int i=0;i<total;i++){
		if(i&(i<<1)) continue;
		if(i&(i<<2)) continue;
		state[top]=i;
		num[top++]=num1_number(i);
	}
}
bool fit(int a,int b)
{
	if(a&b) return false;
	return true;
}
void DP()
{
	int i,j,k,t;
	for(i=0;i<top;i++){
		if(fit(state[i],cur[1]))
			dp[1][i][0]=num[i];
	}
	for(i=0;i<top;i++){
		if(fit(state[i],cur[2])){
			for(j=0;j<top;j++){
				if(!fit(cur[1],state[j])) continue;
				if(fit(state[j],state[i]))
					dp[2][i][j]=max(dp[2][i][j],dp[1][j][0]+num[i]);
			}
		}
	}
	for(i=3;i<=m;i++){
		for(j=0;j<top;j++){
			if(!fit(cur[i],state[j])) continue;
			for(k=0;k<top;k++){
				if(!fit(cur[i-1],state[k])) continue;
				if(fit(state[k],state[j]))
					for(t=0;t<top;t++){
						if(!fit(cur[i-2],state[t])) continue;
						if(fit(state[t],state[j])&&fit(state[t],state[k]))
							dp[i][j][k]=max(dp[i][j][k],(num[j]+dp[i-1][k][t]));
					}
			}
		}
	}
	int ans=0;
	for(i=1;i<=m;i++)
		for(j=0;j<top;j++)
			for(k=0;k<top;k++)
				ans=max(ans,dp[i][j][k]);
	printf("%d\n",ans);
}
int main()
{
	int i,j;
	char s;
	scanf("%d%d",&m,&n);
	getchar();
	memset(num,0,sizeof(num));
	for(i=1;i<=m;i++){
		for(j=1;j<=n;j++){
			s=getchar();
			if(s=='H') cur[i]+=1<<(n-j);
			else num[i]++;
		}
		getchar();
	}
	top=0;
	Init();
	DP();
}