bzoj2310 ParkII（插头dp）

类似bzoj1187
这回要求一条简单路的最大权值和。
简单路与简单回路的区别就在于，我们可以有独立插头。
于是多一维状态表示独立插头即可。
转移时多讨论一些情况就好了qaq
此时p=1,q=2的插头不管其实也是可以的，因为一条简单回路一定也能视作一条简单路，在其他情况被讨论到。
注意到我们这样无法考虑到只有一个点的简单路，于是乎单独处理掉。
复杂度$O(nm4^mm)$

#include <bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define N 100003
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,a[110][10],ans=-inf,pp,c[]={0,-1,1,0};
struct Hash_table{
    int next,key,val;
};
inline void Max(int &x,int y){x=max(x,y);}
struct Icefox{
    int h[N],num;
    Hash_table data[N];
    inline void init(){num=0;memset(h,0,sizeof(h));}
    inline void ins(int key,int val){
        int x=key%N;
        for(int i=h[x];i;i=data[i].next)
            if(key==data[i].key){Max(data[i].val,val);return;}
        data[++num].key=key;data[num].val=val;data[num].next=h[x];h[x]=num;
    }
}dp[2];
inline int get1(int S,int x){x<<=1;return S>>x&3;}
inline int set1(int S,int x,int val){x<<=1;return S&~(3<<x)|val<<x;}
inline int getr(int S,int r){
    int cnt=-1;while(cnt) cnt+=c[get1(S,++r)];return r;
}
inline int getl(int S,int l){
    int cnt=1;while(cnt) cnt+=c[get1(S,--l)];return l;
}
inline void update(int x,int y,int S,int val){
    int p=get1(S,y),q=get1(S,y+1),ss;
    if(!p&&!q){
        ss=set1(S,y,0);ss=set1(ss,y+1,0);
        dp[pp^1].ins(ss,val-a[x][y]);//不放插头 
        if(x<n-1){//放下插头 
            ss=set1(S,y,3);ss=set1(ss,y+1,0);
            dp[pp^1].ins(ss,val);
        }if(y<m-1){//放右插头 
            ss=set1(S,y,0);ss=set1(ss,y+1,3);
            dp[pp^1].ins(ss,val);
        }if(x<n-1&&y<m-1){//右+下 
            ss=set1(S,y,1);ss=set1(ss,y+1,2);
            dp[pp^1].ins(ss,val);
        }return;
    }if(!p||!q){
        if(x<n-1){//对应左/上+下 
            ss=set1(S,y,p+q);ss=set1(ss,y+1,0);
            dp[pp^1].ins(ss,val);
        }if(y<m-1){//对应左/上+右 
            ss=set1(S,y,0);ss=set1(ss,y+1,p+q);
            dp[pp^1].ins(ss,val);
        }ss=set1(S,y,0);ss=set1(ss,y+1,0);//不放插头 
        if(p+q==1) ss=set1(ss,getr(S,y+1),3);
        if(p+q==2) ss=set1(ss,getl(S,y),3);
        if(p+q==3){if(!ss) Max(ans,val);return;}
        dp[pp^1].ins(ss,val);return;
    }ss=set1(S,y,0);ss=set1(ss,y+1,0);
    if(p==3&&q==3){if(!ss) Max(ans,val);return;}
    if(p==3&&q==1) ss=set1(ss,getr(S,y+1),3);
    if(p==3&&q==2) ss=set1(ss,getl(S,y),3);
    if(p==1&&q==3) ss=set1(ss,getr(S,y+1),3);
    if(p==2&&q==3) ss=set1(ss,getl(S,y),3);
    if(p==1&&q==1) ss=set1(ss,getr(S,y+1),1);
    if(p==2&&q==2) ss=set1(ss,getl(S,y),2);
    if(p==1&&q==2){if(!ss) Max(ans,val);return;}
    dp[pp^1].ins(ss,val);
}

int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=0;i<n;++i)
        for(int j=0;j<m;++j) a[i][j]=read(),Max(ans,a[i][j]);
    pp=0;dp[pp].init();dp[pp].ins(0,0);
    for(int i=0;i<n;++i){
        for(int j=1;j<=dp[pp].num;++j) dp[pp].data[j].key<<=2;
        for(int j=0;j<m;++j){
            dp[pp^1].init();
            for(int k=1;k<=dp[pp].num;++k)
                update(i,j,dp[pp].data[k].key,dp[pp].data[k].val+a[i][j]);
            pp^=1;
        }
    }printf("%d\n",ans);
    return 0;
}