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;
}