来自FallDream的博客，未经允许，请勿转载，谢谢，
经历了一段艰辛的旅程后，主人公小P乘坐飞艇返回。在返回的途中，小P发现在漫无边际的沙漠中，有一块狭长的绿地特别显眼。往下仔细一看，才发现这是一个游乐场，专为旅途中疲惫的人设计。娱乐场可以看成是一块大小为n×m的区域，且这个n×m的区域被分成n×m个小格子，每个小格子中就有一个娱乐项目。然而，小P并不喜欢其中的所有娱乐项目，于是，他给每个项目一个满意度。满意度为正时表示小P喜欢这个项目，值越大表示越喜欢。为负时表示他不喜欢，这个负数的绝对值越大表示他越不喜欢。为0时表示他对这个项目没有喜恶。小P决定将飞艇停在某个小格中，然后每步他可以移动到相邻的上下左右四个格子的某个格子中。小P希望找一条路径，从飞艇所在格出发，最后又回到这个格子。小P有一个习惯，从不喜欢浪费时间。因此，他希望经过每个格子都是有意义的：他到一个地方后，就一定要感受以下那里的惊险和刺激，不管自己是不是喜欢那里的娱乐项目。而且，除了飞艇所在格，其他的格子他不愿意经过两次。小P希望自己至少要经过四个格子。在满足这些条件的情况下，小P希望自己玩过的娱乐项目的满意度之和最高。你能帮他找到这个最高的满意度之和吗？
n<=100 m<=6
很显然是插头dp
因为不要求全部走完，所以和回路有点区别。
如果两个插头是()，且其他都没有插头，可以直接更新答案。
如果两个插头都没有，那么可以不选择这一格进行转移。
可以预处理状态，最多127种，复杂度127nm
#include<iostream>
#include<cstdio>
#include<cstring>
#define INF 1000000000
using namespace std;
inline int read()
{
int x = , f = ; char ch = getchar();
while(ch < '' || ch > ''){ if(ch == '-') f = -; ch = getchar();}
while(ch >= '' && ch <= ''){x = x * + ch - '';ch = getchar();}
return x * f;
} int s[],q[],n,m,top=,ans=-INF,tot=,a[][],c[][],f[][<<];
inline void R(int&x,int y){y>x?x=y:;}
int main()
{
n=read();m=read();
for(int i=;i<=n;++i)
for(int j=;j<=m;++j)
a[i][j]=read();
for(int i=;i<<<(m<<)+;++i)
{
s[++tot]=i;top=;
for(int j=;j<=m;++j)
{
int x=i>>(j<<);
if((x&)==) {top=-;break;}
if((x&)==) q[++top]=j;
if((x&)==)
{
if(!top) {top=-;break;}
c[tot][j]=q[top];c[tot][q[top]]=j;
--top;
}
}
if(top) --tot;
}
for(int i=;i<;++i)
for(int j=;j<<<;++j)
f[i][j]=-INF;
f[m][]=;
for(int i=;i<=n;++i)
{
for(int j=;j<=tot;++j)
{
if(s[j]&) f[][s[j]]=-INF;
else f[][s[j]]=f[m][s[j]>>];
}
for(int j=;j<=m;++j)
{
int x=(j-)<<;
for(int k=;k<=tot;++k) f[j][s[k]]=-INF;
for(int k=;k<=tot;++k)
{
int p=(s[k]>>x)&,q=(s[k]>>(x+))&;
if(!p&&!q)
{
R(f[j][s[k]],f[j-][s[k]]);
R(f[j][s[k]^(<<x)],f[j-][s[k]]+a[i][j]);
}
else if(p&&q)
{
if(p==&&q==)
R(f[j][s[k]^(<<x)^(<<(c[k][j]<<))],f[j-][s[k]]+a[i][j]);
if(p==&&q==)
if(s[k]==(<<x)) R(ans,f[j-][s[k]]+a[i][j]);
if(p==&&q==)
R(f[j][s[k]^(<<x)],f[j-][s[k]]+a[i][j]);
if(p==&&q==)
R(f[j][s[k]^(<<x)^(<<(c[k][j-]<<))],f[j-][s[k]]+a[i][j]);
}
else
{
R(f[j][s[k]],f[j-][s[k]]+a[i][j]);
R(f[j][s[k]^(p<<x)^(q<<x+)|(p<<x+)|(q<<x)],f[j-][s[k]]+a[i][j]);
}
}
}
}
cout<<ans;
return ;
}