bzoj 1187: [HNOI2007]神奇遊樂園【插頭dp】
阿新 • • 發佈:2018-12-14
要判邊界!!要判邊界!!要判邊界!!if(j!=m)!!!
我真是zz橫著轉移要判斷到底能不能向右邊出邊……
然後剩下的和1814差不多,九十因為不要求經過所有格子,所以左右括號隨時可以合併,但是注意合併的時候輪廓線上不能有別的括號,然後還有是(0,0)可以轉移到(0,0),相當一不經過這個格子
#include<iostream> #include<cstdio> #include<cstring> using namespace std; const int N=300005,mod=299989; int n,m,nw,la,h[N],b[20],c[2]; long long a[105][105],ans=-1e18; struct qwe { int ne,to[2]; long long va[2]; }f[N]; int read() { int r=0,f=1; char p=getchar(); while(p>'9'||p<'0') { if(p=='-') f=-1; p=getchar(); } while(p>='0'&&p<='9') { r=r*10+p-48; p=getchar(); } return r*f; } void add(int x,long long v) {//cerr<<v<<endl; int u=x%mod+1; for(int i=h[u];i;i=f[i].ne) if(f[i].to[nw]==x) { f[i].va[nw]=max(f[i].va[nw],v); return; } c[nw]++; f[c[nw]].ne=h[u]; f[c[nw]].to[nw]=x; f[c[nw]].va[nw]=v; h[u]=c[nw]; } int main() { n=read(),m=read(); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) a[i][j]=read(); b[0]=1; for(int i=1;i<=10;i++) b[i]=b[i-1]<<2; c[0]=1,f[1].va[0]=0,f[1].to[0]=0; for(int i=1;i<=n;i++) { for(int j=1;j<=c[nw];j++) f[j].to[nw]<<=2; for(int j=1;j<=m;j++) {//cerr<<i<<" "<<j<<endl; la=nw,nw^=1; memset(h,0,sizeof(h)); c[nw]=0; for(int k=1;k<=c[la];k++) { int x=f[k].to[la],b1=(x>>(j*2-2))%4,b2=(x>>(j*2))%4; long long v=f[k].va[la]; if(!b1&&!b2) { add(x,v); if(j!=m) add(x+b[j-1]+2*b[j],v+a[i][j]); } else if(!b1&&b2) { add(x+b2*b[j-1]-b2*b[j],v+a[i][j]); if(j!=m) add(x,v+a[i][j]); } else if(b1&&!b2) { add(x,v+a[i][j]); if(j!=m) add(x-b1*b[j-1]+b1*b[j],v+a[i][j]); } else if(b1==1&&b2==1) { for(int t=1,l=j+1;l<=m;l++) { if((x>>(l*2))%4==1) t++; if((x>>(l*2))%4==2) t--; if(!t) { add(x-b[j-1]-b[j]-b[l],v+a[i][j]); break; } } } else if(b1==2&&b2==2) { for(int t=1,l=j-2;l>=0;l--) { if((x>>(l*2))%4==1) t--; if((x>>(l*2))%4==2) t++; if(!t) { add(x-2*b[j-1]-2*b[j]+b[l],v+a[i][j]); break; } } } else if(b1==2&&b2==1) add(x-2*b[j-1]-b[j],v+a[i][j]); else if(x==b[j-1]+b[j]*2) ans=max(ans,v+a[i][j]); } } } printf("%lld\n",ans); return 0; }