您的位置:首页 > 其它

BZOJ1187 [HNOI2007]神奇游乐园(插头dp)

2017-06-15 19:23 417 查看
麻麻我会写插头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
,g
[7],f[105][7]
;
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;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: