poj1185 炮兵布阵(状压dp)

思路来源

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

心得

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

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

正确性显然

 

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

自己水平不够，没法合并

 

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

代码

#include <iostream>
#include <algorithm> 
#include <cstring>
#include <cstdio>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <stack>
#include <queue>
#include <functional>
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<int> 
#define si set<int>
#define pii pair<int,int> 
#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<<m)-1)
   	{
   		if(adjacant(i))
   		{
   			int k=i;
   			while(k)e[cnt].num+=(k&1),k>>=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<iostream>
#include<cstdio>
#include<cmath> 
#include<cstring>
#include<vector>
#include<algorithm>
#include<map>
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<<n);++i)
	{
		if(!(i&(i<<1)))
		{
			if(!(i&(i<<2)))
			{
				state[cnt]=i;
				sum[cnt]=__builtin_popcount(i);
				cnt++;
            }
	    }
	}
} 
int main()
{
	while(~scanf("%d%d",&m,&n))
	{
	init();
	for(int i=1;i<=m;++i)
	{
	 scanf("%s",a);
	 int len=strlen(a);
	 mp[i]=0;
	 for(int j=0;j<len;++j)
	 mp[i]=(mp[i]<<1)+(a[j]=='H');
        }
	for(int j=0;j<cnt;++j)//这行 
	{
		if(state[j]&mp[1])continue;
	    dp[1][j][0]=sum[j];
	}
	for(int j=0;j<cnt;++j)//这行 
	{
			if(state[j]&mp[2])continue;
			for(int k=0;k<cnt;++k)//上行 
			{
				if(state[k]&mp[1])continue;
				if(state[j]&state[k])continue;
				dp[2][j][k]=max(dp[2][j][k],dp[1][k][0]+sum[j]);
				//printf("%d:%lld:%lld\n",i,state[j],dp[i][state[j]]);
			}
	}
	for(int i=3;i<=m;++i)
	{
		for(int j=0;j<cnt;++j)//这行 
	    {
			if(state[j]&mp[i])continue;
			for(int k=0;k<cnt;++k)//上行 
			{
				if(state[k]&mp[i-1])continue;
				if(state[j]&state[k])continue;
				for(int l=0;l<cnt;++l)
				{
					if(state[l]&mp[i-2])continue;
					if(state[j]&state[l])continue;
					if(state[k]&state[l])continue;
					dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+sum[j]);
				}
			}
        }
	}
	for(int i=0;i<cnt;++i)//注意到m<3可能成立 
	{
		for(int j=0;j<cnt;++j)
		res=max(res,dp[m][i][j]);
	}
	printf("%lld\n",res);
    }
	return 0;
}