poj1185 炮兵布阵(状压dp)

poj1185 炮兵布阵(状压dp)_第1张图片

思路来源

https://blog.csdn.net/qq_34374664/article/details/54900865

心得

判断当前枚举的放1的位置a是否合法,既可以和原地图对应的十进制数b作a|b==b的判断

也可以对地图取反,做a&(!b)==0的判断

正确性显然

 

对于第一行和第二行就特殊讨论一下吧

自己水平不够,没法合并

 

mdzz这种题真是wa得查不出错心态爆炸

poj1185 炮兵布阵(状压dp)_第2张图片

代码

#include 
#include  
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
const int INF=0x3f3f3f3f;
const int maxn=1e2+10; 
const int mod=1e8;
const int MOD=998244353;
const double eps=1e-7;
typedef long long ll;
#define vi vector 
#define si set
#define pii pair 
#define pi acos(-1.0)
#define pb push_back
#define lowbit(x) (x&(-x))
#define sci(x) scanf("%d",&(x))
#define scll(x) scanf("%lld",&(x))
#define sclf(x) scanf("%lf",&(x))
#define pri(x) printf("%d",(x))
#define rep(i,j,k) for(int i=j;i<=k;++i)
#define per(i,j,k) for(int i=j;i>=k;--i)
#define mem(a,b) memset(a,b,sizeof(a)) 
using namespace std;
int n,m;
char s[15];
int mp[105],cnt;//平原置0 山地H置1 A|B==B等价于A&(!B)==0 
int dp[110][105][105];
struct node
{
	int state;//行合法状态对应十进制数 
	int num;//炮个数 
	node():num(0){
	}
}e[105];
bool disjoint(int a,int b)
{
	return (a&b)==0;
}
bool adjacant(int now)
{
	return disjoint(now,now<<1)&&disjoint(now,now<<2);//显然包括(now<<1)&(now<<2)==0的情况 
}
int main()
{ 
   while(~scanf("%d%d",&n,&m))
   {
   	rep(i,0,n-1)
   	{
   		scanf("%s",s);
   		rep(j,0,m-1)
   		{
   		 mp[i]=(mp[i]<<1)+(s[j]=='H');
   		}
   	}
   	rep(i,0,(1<>=1;
   			e[cnt++].state=i;
   		}
   	}
   	rep(j,0,cnt-1)//特判第一行 
   	{
   		int pre=e[j].state;
   		if(!disjoint(pre,mp[0]))continue;//和山地相邻 
   		dp[0][j][0]=e[j].num;
   	}
   	rep(j,0,cnt-1)//特判第二行 
	{
		int pre=e[j].state;
   		if(!disjoint(pre,mp[1]))continue;
   		rep(k,0,cnt-1)
		{
			int now=e[k].state;
			if(!disjoint(now,mp[0]))continue;
			if(!disjoint(pre,now))continue;
			dp[1][j][k]=max(dp[1][j][k],dp[0][k][0]+e[j].num);
		} 
	}
   	rep(i,2,n-1)
   	{
   		rep(j,0,cnt-1)
   		{
   			int now=e[j].state;
   			if(!disjoint(now,mp[i]))continue; 
   			rep(k,0,cnt-1)
   			{
   				int pre=e[k].state;
   				if(!disjoint(pre,mp[i-1]))continue;
   				if(!disjoint(pre,now))continue;
   				rep(l,0,cnt-1)
   				{
   					int pp=e[l].state;
   					if(!disjoint(pp,mp[i-2]))continue;
   					if(!disjoint(pp,pre))continue;
   					if(!disjoint(pp,now))continue;
   					dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+e[j].num);
   				}
   			}
   		}
   	}
   	int ans=0;
   	rep(i,0,cnt-1)
   	{
   		rep(j,0,cnt-1)
   		{
   			ans=max(ans,dp[n-1][i][j]);
   		}
   	}
   	printf("%d\n",ans);
   }
   return 0;
}

后记(20190212)

好叭我又做了一遍……

又温习了一下这个两个状态的转移(毕竟cf场前两天刚出过一个牌的三维dp)……

又温习了一下这个把山地状态记录下来,和山地状态相与为0的判断……

毕竟又是T了两发WA了一发……真菜啊……

ans得最后判,m可能小于3,不然每个循环里都得加……

还是不会把三个合成一个的操作(应该是i=2开始输入然后把i=1和i=0全赋成平原???)

不想写了……太菜了……

#include
#include
#include 
#include
#include
#include
#include
using namespace std;
typedef long long ll;
int m,n;
char a[12];
ll state[1<<11],sum[1<<11],cnt;
ll res,mp[111],dp[111][111][111];//这行 上行状态标号 上上行状态标号 三维滚动 
void init()
{
	cnt=res=0;
	memset(dp,0,sizeof(dp));
	for(int i=0;i<(1<

 

你可能感兴趣的:(状压dp)