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