【HNOI2007/Luogu3190】神奇游乐园 插头DP

原题走这里

状态和模板一样，转移也几乎和模板一样。
但是由于本题并不要求覆盖所有点，因此要允许“##”到“##”的转移
此外由于本题是单回路，所以“()”到“##”的转移也只能发生一次，为了方便，我们可以干脆把形如“####()###”的状态，也就是只有一对相邻括号的状态作为最终解。
取所有这些最终解的最大值即可。

代码如下：

#include <bits/stdc++.h>
#define LL long long
#define chkmax(x,y) x=max(x,y)
#define clean(x) for(int k=0;k<s;k++)x[k]=-1e16;
using namespace std;
int n,m,s,H[200000],a[102][7];
map<int,int> IH;
LL d[2][200000],ret=-1e16;
void getHash(int k,int x,int y)
{
    if(y<0||y>m-k+1)return;
    if(k==m+1)IH[x]=s,H[s++]=x;
    getHash(k+1,x<<2,y);
    getHash(k+1,x<<2|2,y-1);
    getHash(k+1,x<<2|3,y+1);
}
LL DP()
{
    clean(d[1]);
    d[1][0]=0;
    for(int i=1,cur=1;i<=n;i++,cur^=1)
    {
        for(int j=0;j<m;j++,cur^=1)
        {
            clean(d[cur^1]);
//          d[cur^1][0]=0;
            for(int k=0;k<s;k++)
            {
//              cout<<i<<" "<<j<<' '<<k<<' '<<d[cur][k]<<endl;
                LL temp=d[cur][k];
                if(temp>1e16)continue;
                int t=H[k],x=t>>(j<<1)&3,y=t>>((j+1)<<1)&3;
                if(x+y==3||x+y==2)
                {
                    chkmax(d[cur^1][k],temp+a[i][j]);
                    chkmax(d[cur^1][IH[t^((x|y)<<(j<<1))^((x|y)<<((j+1)<<1))]],temp+a[i][j]);
                }   
                if(x==3&&y==2) 
                {
                    chkmax(d[cur^1][IH[t^(x<<(j<<1))^(y<<((j+1)<<1))]],temp+a[i][j]);
                }
                if(x==2&&y==3&&(H[k]^(x<<(j<<1))^(y<<((j+1)<<1)))==0)//终结状态 
                {
                    chkmax(ret,temp+a[i][j]); 
                }
                if((x|y)==0)
                {
                    if(k==0&&temp)ret=min(ret,temp);
                    chkmax(d[cur^1][IH[t^(2<<(j<<1))^(3<<((j+1)<<1))]],temp+a[i][j]);
                    chkmax(d[cur^1][k],temp);
                }
                if(x==2&&y==2)
                {
                    for(int ii=j+2,jj=0;ii<=m;ii++)
                    {
                        if((t>>(ii<<1)&3)==2)jj++;
                        if((t>>(ii<<1)&3)==3)jj--;
                        if(jj<0)
                        {
                            chkmax(d[cur^1][IH[t^(2<<(j<<1))^(2<<((j+1)<<1))^(1<<(ii<<1))]],temp+a[i][j]);
                            break;
                        }   
                    } 
                }
                if(x==3&&y==3)
                {
                    for(int ii=j-1,jj=0;ii>=0;ii--)
                    {
                        if((t>>(ii<<1)&3)==3)jj++;
                        if((t>>(ii<<1)&3)==2)jj--;
                        if(jj<0)
                        {
                            chkmax(d[cur^1][IH[t^(3<<(j<<1))^(3<<((j+1)<<1))^(1<<(ii<<1))]],temp+a[i][j]);
                            break;
                        }   
                    } 
                }
            }
        }
        clean(d[cur^1]);
        for(int j=0;H[j]<(1<<(m<<1));j++)
        {
            d[cur^1][IH[H[j]<<2]]=d[cur][j];
        }
    }
    return ret;
}
int main() {
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<m;j++)
        {
            cin>>a[i][j];
        }
    }
    getHash(0,0,0);
    cout<<DP()<<endl;
    return 0;
}