BZOJ1187 [HNOI2007]神奇游乐园(插头dp)
麻麻我会写插头dp了!
推荐陈丹琦论文:https://wenku.baidu.com/view/3e90d32b453610661ed9f4bd.html
破题调一年
#include <cstdio> #include <cstring> const int N=1<<14; int n,m,tt,a1=0xcfcfcfcf,st[7],a[105][7],b[N],g[N][7],f[105][7][N]; void up(int &x,int y) {if(x<y) x=y;} void dfs(int x,int s,int k) { if(x==m+2) { if(!k) { for(int i=0,tp=0;i<m;i++) { int t=(s>>(i*2))&3; if(t==2) st[++tp]=i; else if(t==1) g[s][i]=st[tp],g[s][st[tp]]=i,tp--; } b[++tt]=s; } return; } dfs(x+1,s<<2,k),dfs(x+1,s<<2|1,k+1); if(k) dfs(x+1,s<<2|2,k-1); } int main() { scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&a[i][j]); memset(f,0xcf,sizeof f),f[1][1][0]=0,f[1][1][(1<<(m*2))|(2*(1<<(m-1)*2))]=a[1][1],dfs(1,0,0); for(int i=1;i<=n;i++) { for(int j=1;j<m;j++) { for(int k=1;k<=tt;k++) if(f[i][j][b[k]]!=0xcfcfcfcf) { int t1=(b[k]>>((m-j)*2))&3,t2=(b[k]>>((m-j-1)*2))&3,t3=1<<((m-j)*2),t4=1<<((m-j-1)*2); int s=b[k]&((1<<((m+1)*2))-1-t3*3-t4*3),t=f[i][j][b[k]]+a[i][j+1]; if(!t1&&!t2) up(f[i][j+1][s],t-a[i][j+1]),up(f[i][j+1][s|t3|(t4*2)],t); else if(!t1) up(f[i][j+1][s|(t3*t2)],t),up(f[i][j+1][s|(t4*t2)],t); else if(!t2) up(f[i][j+1][s|(t3*t1)],t),up(f[i][j+1][s|(t4*t1)],t); else if(t1==1&&t2==1) up(f[i][j+1][s-(1<<(g[b[k]][m-j-1]*2))],t); else if(t1==1&&t2==2) { if(b[k]==t3+t4*2&&a1<t) a1=t; } else if(t1==2&&t2==1) up(f[i][j+1][s],t); else up(f[i][j+1][s+(1<<(g[b[k]][m-j]*2))],t); } } for(int k=1;k<=tt;k++) if(f[i][m][b[k]]!=0xcfcfcfcf) { int t1=(b[k]>>(m*2))&3,t3=1<<(m*2),t4=1<<((m-1)*2),s=(b[k]&((1<<(m*2))-1))>>2,t=f[i][m][b[k]]+a[i+1][1]; if(!t1) up(f[i+1][1][s],t-a[i+1][1]),up(f[i+1][1][s|t3|(t4*2)],t); else up(f[i+1][1][s|(t3*t1)],t),up(f[i+1][1][s|(t4*t1)],t); } } printf("%d",a1); return 0; }