bzoj 1187: [HNOI2007]神奇游乐园

人生第一次写插头DP，从此患上了插头DP恐惧症

昨天晚上写的时候以为是轮廓线动规，然后发现没法确定联通块的数目，于是乎设计出了一种诡异的分层图DP，然后发现不对。。。。。。跪了……

遂翻题解，发现时插头DP。

膜了下论文，表示好坑爹，怎么这么多分类讨论啊。

然后经过了艰苦卓绝的调试过程，终于A掉了

话说这代码写这么烂我自己都看不懂了

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int inf=1e9;
int f[2][18000],bin[20];
int state[10000],tot,n,m;
int get(int s,int pos){
	if(s&bin[pos*2])return 2;
	else if(s&bin[pos*2-1])return 1;
	return 0; 
}
#define refresh(x,y) (x)=max((x),(y))
#define cl(x) bin[2*(x)-1]
#define cr(x) bin[2*(x)]
#define c(x) ((x)<<2)
void dfs(int pos,int s){
	if(pos>m+1)state[++tot]=s;
	else{
		dfs(pos+1,s);
		dfs(pos+1,s|cl(pos));
		dfs(pos+1,s|cr(pos));
	}
}
void print(int s){
	for(int i=1;i<=m+1;i++)
	if(get(s,i)==1)putchar('(');
	else if(get(s,i)==2)putchar(')');
	else putchar('#');
	putchar('\n');
}
int main(){
	//freopen("a.in","r",stdin);
	//freopen("a.out","w",stdout); 
	scanf("%d%d",&n,&m);
	bin[1]=1;
	for(int i=2;i<=20;i++)bin[i]=bin[i-1]<<1;
	dfs(1,0);
	int v;
	memset(f[0],-0x3f,sizeof(f[0]));
	int now=0,last=1;
	f[now][0]=0;
	int ans=-inf;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++){
		now^=1;last^=1;scanf("%d",&v);
		memset(f[now],-0x3f,sizeof(f[now]));
		for(int k=1;k<=tot;k++){
			int s=state[k];
			if(f[last][s]<-inf)continue;
			int x=get(s,j),y=get(s,j+1);
			//if(i==2&&j==2)printf("%d %d %d %d\n",i,j,f[last][s],s),print(s);
			if(!x&&!y){
				if(j!=m)refresh(f[now][s],f[last][s]);
				else refresh(f[now][c(s)],f[last][s]);
				if(i!=n&&j!=m)refresh(f[now][s|cl(j)|cr(j+1)],f[last][s]+v);
			}else if(x==1&&y==1){
				for(int t=j+2;t<=m+1;t++)
				if(get(s,t)==2){
					int si=s^cr(t)^cl(j)^cl(j+1);
					if(j!=m)refresh(f[now][si|cl(t)],f[last][s]+v);
					else refresh(f[now][c(si|cl(t))],f[last][s]+v);
					break;
				}
			}else if(x==2&&y==2){
				for(int t=j-1;t>=1;t--)
				if(get(s,t)==1){
					int si=s^cl(t)^cr(j)^cr(j+1);
					if(j!=m)refresh(f[now][si|cr(t)],f[last][s]+v);
					else refresh(f[now][c(si|cr(t))],f[last][s]+v);
					break;
				}
			}else if(x==2&&y==1){
				if(j!=m)refresh(f[now][s^cr(j)^cl(j+1)],f[last][s]+v);
			}else if(x==1&&y==2){
				bool flag=true;
				for(int t=1;t<=m+1;t++)
				if(t!=j&&t!=j+1&&get(s,t))flag=false; 
				if(flag)ans=max(ans,f[last][s]+v);
			}else if(x==1&&!y){
				if(i!=n){
					if(j!=m)refresh(f[now][s],f[last][s]+v);
					else refresh(f[now][c(s)],f[last][s]+v);
				}
				if(j!=m)refresh(f[now][s^cl(j)|cl(j+1)],f[last][s]+v);
			}else if(y==1&&!x){
				if(j!=m)refresh(f[now][s],f[last][s]+v);
				if(i!=n){
					if(j!=m)refresh(f[now][s^cl(j+1)|cl(j)],f[last][s]+v);
					else refresh(f[now][c(s^cl(j+1)|cl(j))],f[last][s]+v);
				}
			}else if(x==2&&!y){
				if(i!=n){
					if(j!=m)refresh(f[now][s],f[last][s]+v);
					else refresh(f[now][c(s)],f[last][s]+v);
				}
				if(j!=m)refresh(f[now][s^cr(j)|cr(j+1)],f[last][s]+v);
			}else if(y==2&&!x){
				if(j!=m)refresh(f[now][s],f[last][s]+v);
				if(i!=n){
					if(j!=m)refresh(f[now][s^cr(j+1)|cr(j)],f[last][s]+v);
					else refresh(f[now][c(s^cr(j+1)|cr(j))],f[last][s]+v);
				}
			}
		}
	}
	printf("%d\n",ans);
	return 0;
}