[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=0,flag=1;char ch=getchar();
    while(ch>'9'||ch<'0'){
        if(ch=='-') flag=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9') re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
    return re*flag;
}
int n,m,x[150][150],cur,pre,ex,ey;
int st[2][300010];ll ans[2][300010],re;
int tot[2],bit[20],state[300010],st_tot,hash=300000;
struct edge{
    int to,next;
}a[300010];
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=0;tot[cur]=1;ans[cur][1]=0;st[cur][1]=0;
    for(i=1;i<=n;i++){
        for(j=1;j<=tot[cur];j++) st[cur][j]<<=2;
        for(j=1;j<=m;j++){
            pre=cur;cur^=1;tot[cur]=0;st_tot=0;memset(state,0,sizeof(state));
            for(k=1;k<=tot[pre];k++){
                now=st[pre][k];val=ans[pre][k];
                right=(now>>bit[j-1])%4;down=(now>>bit[j])%4;
                if(!down&&!right){
                    insert(now,val);
                    if(j!=m) 
                        insert(now+(1<<bit[j-1])+((1<<bit[j])<<1),val+x[i][j]);
                }
                if(down&&!right){
                    insert(now-down*(1<<bit[j])+down*(1<<bit[j-1]),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*(1<<bit[j])-right*(1<<bit[j-1]),val+x[i][j]);
                }
                if(right==1&&down==1){
                    int cnt=1;
                    for(l=j+1;l<=m;l++){
                        if((now>>bit[l])%4==1) cnt++;
                        if((now>>bit[l])%4==2) cnt--;
                        if(!cnt){
                            insert(now-(1<<bit[l])-(1<<bit[j])-(1<<bit[j-1]),val+x[i][j]);
                            break;
                        }
                    }
                }
                if(right==2&&down==2){
                    int cnt=1;
                    for(l=j-2;l>=0;l--){
                        if((now>>bit[l])%4==1) cnt--;
                        if((now>>bit[l])%4==2) cnt++;
                        if(!cnt){
                            insert(now+(1<<bit[l])-((1<<bit[j])<<1)-((1<<bit[j-1])<<1),val+x[i][j]);
                            break;
                        }
                    }
                }
                if(right==2&&down==1){
                    insert(now-((1<<bit[j-1])<<1)-(1<<bit[j]),val+x[i][j]);
                }
                if(right==1&&down==2){
                    if((now==(1<<bit[j-1])+((1<<bit[j])<<1))&&(val+x[i][j]>re)){
                        re=val+x[i][j];
                    }
                }
            }
        }
    }
}
int main(){
    int i,j;
    n=read();m=read();
    for(i=1;i<=10;i++) bit[i]=(i<<1);
    for(i=1;i<=n;i++) for(j=1;j<=m;j++) x[i][j]=read();
    dp();
    printf("%lld",re);
}