POJ 1185炮兵阵地（状压DP+优化）

最新推荐文章于 2021-07-15 11:32:37 发布

FainTwinkle

最新推荐文章于 2021-07-15 11:32:37 发布

阅读量277

点赞数

分类专栏：状压DP

本文链接：https://blog.csdn.net/sinat_30126425/article/details/52068147

版权

状压DP 专栏收录该内容

4 篇文章 0 订阅

订阅专栏

这题是个十分经典的状态压缩DP了

如果直接状态压缩，每一行都是 2^10的复杂度，爆炸。

所以可以根据题目的一些性质来进行优化，这里他说每一个炮弹周围两个格子不能放炮弹了。所以就可以把一些没有用的状态处理掉。

只留下的有用的状态，这里可以用一个初始化来进行。stu这个数组来记录所有可能出现的状态，这个预处理一下就好了，num表示每个状态中1的个数。

然后就暴力dp就好辣！

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

这个式子的意思是第i行的状态为 stu[j] 上一行的状态为stu[ k ] ，所对应的放的炮弹的个数。

AC代码：

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
char mm[110][20];
int pp[110];
int dp[101][70][70];
int stu[70],num[70];
int tot;
void init(int m){
	tot=0;
	int flag=0,nu=0;
	memset(num,0,sizeof(num));
	for(int s=0;s<(1<<m);s++){
		flag=0,nu=0;
		for(int i=0;i<=9;i++){
			int dd=(1<<i)&s;
			if((1<<i)&s){
				nu++;
				if(1<<(i-1)&s  && (i-1)>=0) {flag=1;break;}
				if(1<<(i-2)&s  && (i-2)>=0) {flag=1;break;}
				if(1<<(i+1)&s  && (i+1)<=9) {flag=1;break;}
				if(1<<(i+2)&s  && (i+2)<=9) {flag=1;break;}
			}
		}
		if(flag==1)	continue;
		num[tot]=nu;
		stu[tot++]=s;
	}
}

int main()
{
	
	int n,i,j;
	int m;
	scanf("%d%d",&n,&m);
	init(m);
	for(int i=1;i<=n;i++)
		scanf("%s",mm[i]+1);
	memset(pp,0,sizeof(pp));
	memset(dp,0,sizeof(dp));
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(mm[i][j]=='H')
				pp[i]=pp[i]^(1<<(m-j));
		}
	}
	for(i=0;i<tot;i++){
		int w=stu[i]&pp[1];
		if(w==0) {dp[1][i][0]=num[i];}
	}
	if(n>=2)
	for(int i=0;i<tot;i++){
		if(stu[i]&pp[2])	continue;
		for(int j=0;j<tot;j++){
			if(stu[j]&pp[1])	continue;
			if(stu[j]&stu[i])	continue;
			dp[2][i][j]=dp[1][j][0]+num[i];
		}
	}
	for(int i=3;i<=n;i++){
		for(int j=0;j<tot;j++){
			if(pp[i]& stu[j])	continue;
			for(int k=0;k<tot;k++){
				if(pp[i-1]&stu[k])	continue;
				if(stu[k]&stu[j])	continue;
				for(int l=0;l<tot;l++){
					if(pp[i-2]&stu[l])	continue;
					if(stu[l]&stu[k])	continue;
					if(stu[l]&stu[j])	continue;
					dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+num[j]);
				}
			}
		}
	}
	int ans=0;
	for(int i=0;i<tot;i++)
		for(int j=0;j<tot;j++)
			{ans=max(dp[n][i][j],ans);}
	cout<<ans<<endl;
}