[HNOI2007][bzoj1187] 神奇游乐园 [插头dp]

题面：

传送门

给定一个四联通棋盘图，每个格子有权值，求一条总权值最大的回路

思路：

插头dp基础教程

棋盘？

回路？

n,m<=10？

当然是插头dp啦~\(≧▽≦)/~

然后发现这道题并不是哈密顿回路了......需要考虑一下终止条件的不同，对比URAL1519

诚然，回路问题中依然可以用括号表示法来压缩状态，但是此时形成回路、统计答案的就不仅仅在最后一个非障碍格子才可以，而是任何时候，只要状态中只剩下一对左右插头满足条件（看不懂这里的可以去上面那道题看看）

同时请注意，这道题统计的是最值、不是和值

同时请注意，这题可能会爆longlong......我丧病地开了__int128，避免了高精度的地狱

Code：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define ll long long
 #define hash deep_dark_fantasy
 #define inf 1e9
 using namespace std;
 inline int read(){
     int re=,flag=;char ch=getchar();
     while(ch>''||ch<''){
         if(ch=='-') flag=-;
         ch=getchar();
     }
     while(ch>=''&&ch<='') re=(re<<)+(re<<)+ch-'',ch=getchar();
     return re*flag;
 }
 int n,m,x[][],cur,pre,ex,ey;
 int st[][];ll ans[][],re;
 int tot[],bit[],state[],st_tot,hash=;
 struct edge{
     int to,next;
 }a[];
 void insert(int sta,ll val){
     int p=sta%hash,i;
     for(i=state[p];i;i=a[i].next){
         if(st[cur][a[i].to]==sta){
             ans[cur][a[i].to]=max(ans[cur][a[i].to],val);return;
         }
     }
     tot[cur]++;
     a[++st_tot].to=tot[cur];
     a[st_tot].next=state[p];
     state[p]=st_tot;st[cur][tot[cur]]=sta;ans[cur][tot[cur]]=val;
 }
 void dp(){
     int i,j,k,l,now,down,right;ll val;re=-inf;
     cur=;tot[cur]=;ans[cur][]=;st[cur][]=;
     for(i=;i<=n;i++){
         for(j=;j<=tot[cur];j++) st[cur][j]<<=;
         for(j=;j<=m;j++){
             pre=cur;cur^=;tot[cur]=;st_tot=;memset(state,,sizeof(state));
             for(k=;k<=tot[pre];k++){
                 now=st[pre][k];val=ans[pre][k];
                 right=(now>>bit[j-])%;down=(now>>bit[j])%;
                 if(!down&&!right){
                     insert(now,val);
                     if(j!=m)
                         insert(now+(<<bit[j-])+((<<bit[j])<<),val+x[i][j]);
                 }
                 if(down&&!right){
                     insert(now-down*(<<bit[j])+down*(<<bit[j-]),val+x[i][j]);
                     if(j!=m)insert(now,val+x[i][j]);
                 }
                 if(right&&!down){
                     insert(now,val+x[i][j]);
                     if(j!=m)
                         insert(now+right*(<<bit[j])-right*(<<bit[j-]),val+x[i][j]);
                 }
                 if(right==&&down==){
                     int cnt=;
                     for(l=j+;l<=m;l++){
                         if((now>>bit[l])%==) cnt++;
                         if((now>>bit[l])%==) cnt--;
                         if(!cnt){
                             insert(now-(<<bit[l])-(<<bit[j])-(<<bit[j-]),val+x[i][j]);
                             break;
                         }
                     }
                 }
                 if(right==&&down==){
                     int cnt=;
                     for(l=j-;l>=;l--){
                         if((now>>bit[l])%==) cnt--;
                         if((now>>bit[l])%==) cnt++;
                         if(!cnt){
                             insert(now+(<<bit[l])-((<<bit[j])<<)-((<<bit[j-])<<),val+x[i][j]);
                             break;
                         }
                     }
                 }
                 if(right==&&down==){
                     insert(now-((<<bit[j-])<<)-(<<bit[j]),val+x[i][j]);
                 }
                 if(right==&&down==){
                     if((now==(<<bit[j-])+((<<bit[j])<<))&&(val+x[i][j]>re)){
                         re=val+x[i][j];
                     }
                 }
             }
         }
     }
 }
 int main(){
     int i,j;
     n=read();m=read();
     for(i=;i<=;i++) bit[i]=(i<<);
     for(i=;i<=n;i++) for(j=;j<=m;j++) x[i][j]=read();
     dp();
     printf("%lld",re);
 }