hdu-4539-郑厂长系列故事――排兵布阵 状态压缩dp

http://blog.youkuaiyun.com/catalyst1314/article/details/19015401

http://blog.youkuaiyun.com/z309241990/article/details/8742321

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

int e[220],s[220],c[220];
int dp[110][220][220];
int count(int x)
{
	int sum=0;
	while(x)
	{
		x=x&(x-1);
		sum++;
	}
	return sum;
}
int main()
{
	int n,m;
	int i,j,k,l;
	while(cin>>n>>m)
	{
		memset(e,0,sizeof(e));
		for(i=0;i<n;i++)
			for(j=0;j<m;j++)
			{
				scanf("%d",&l);
				if(!l)
					e[i]|=(1<<j);
			}
		int tot=(1<<m);
		int cnt=0;
		for(i=0;i<tot;i++)
		{
			if(i&(i<<2)) continue;
			s[cnt]=i;
			c[cnt++]=count(i);
		}
		memset(dp,-1,sizeof(dp));
		for(i=0;i<cnt;i++){
			if(s[i]&e[0]) continue;
			dp[0][i][0]=c[i];
		}
		for(i=1;i<n;i++)
			for(j=0;j<cnt;j++)
			{
				if(s[j]&e[i]) continue;
				for(k=0;k<cnt;k++)
				{
					if(s[j]&(s[k]<<1)||s[j]&(s[k]>>1)) continue;
					for(l=0;l<cnt;l++)
					{
						if(s[j]&s[l]) continue;
						if(s[k]&(s[l]<<1)||s[k]&(s[l]>>1)) continue;
						if(dp[i-1][k][l]==-1) continue;
						dp[i][j][k]=max(dp[i][j][k],dp[i-1][k][l]+c[j]);
					}
				}
			}
		int ans=0;
		for(i=0;i<cnt;i++)
			for(j=0;j<cnt;j++)
				ans=max(ans,dp[n-1][i][j]);
		printf("%d\n",ans);
	}
}