人生第一次写插头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; }