【HNOI2007】bzoj1187 神奇游乐园

Description

经历了一段艰辛的旅程后，主人公小P乘坐飞艇返回。在返回的途中，小P发现在漫无边际的沙漠中，有一块狭长的绿地特别显眼。往下仔细一看，才发现这是一个游乐场，专为旅途中疲惫的人设计。娱乐场可以看成是一块大小为n×m的区域，且这个n×m的区域被分成n×m个小格子，每个小格子中就有一个娱乐项目。然而，小P并不喜欢其中的所有娱乐项目，于是，他给每个项目一个满意度。满意度为正时表示小P喜欢这个项目，值越大表示越喜欢。为负时表示他不喜欢，这个负数的绝对值越大表示他越不喜欢。为0时表示他对这个项目没有喜恶。小P决定将飞艇停在某个小格中，然后每步他可以移动到相邻的上下左右四个格子的某个格子中。小P希望找一条路径，从飞艇所在格出发，最后又回到这个格子。小P有一个习惯，从不喜欢浪费时间。因此，他希望经过每个格子都是有意义的：他到一个地方后，就一定要感受以下那里的惊险和刺激，不管自己是不是喜欢那里的娱乐项目。而且，除了飞艇所在格，其他的格子他不愿意经过两次。小P希望自己至少要经过四个格子。在满足这些条件的情况下，小P希望自己玩过的娱乐项目的满意度之和最高。你能帮他找到这个最高的满意度之和吗？
Input

输入文件中的第一行为两个正整数n和m，表示游乐场的大小为n×m。因为这个娱乐场很狭窄，所以n和m满足：2<=n<=100，2<=m<=6。接下来的n行，每行有m个整数，第i行第j列表示游乐场的第i行第j列的小格子中的娱乐项目的满意度，这个满意度的范围是[-1000，1000]。同一行的两个整数之间用空格隔开。
Output

输出文件中仅一行为一个整数，表示最高的满意度之和。

括号序列维护插头dp。注意以下几点：
如果左边和上边没有插头可以不选。
更新答案的时候必须保证没有其他插头，而且更新的值不能放到dp数组里，防止选到了回路以外的元素。
换行需要特殊处理。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int oo=0x3f3f3f3f;
int dp[2][20010],que[2][20010],last[20010],a[110][10],
n,m,x,i,j,t[2],cur,ans=-oo,flag=1;
void upd(int p)
{
    if (j==m-1)
    {
        if ((p>>2*m)&3) return;
        p<<=2;
    }
    if (last[p]<cur)
    {
        last[p]=cur;
        dp[cur&1][p]=-oo;
        que[cur&1][++t[cur&1]]=p;
    }
    dp[cur&1][p]=max(dp[cur&1][p],dp[cur&1^1][x]+flag*a[i][j]);
}
int main()
{
    int p,q,x1,u,tem;
    scanf("%d%d",&n,&m);
    for (i=0;i<n;i++)
        for (j=0;j<m;j++)
            scanf("%d",&a[i][j]);
    t[0]=1;
    for (i=0;i<n;i++)
        for (j=0;j<m;j++)
        {
            cur++;
            t[cur&1]=0;
            for (int k=1;k<=t[cur&1^1];k++)
            {
                x=que[cur&1^1][k];
                p=(x>>2*j)&3;
                q=(x>>2*(j+1))&3;
                x1=x^(p<<2*j)^(q<<2*(j+1));
                if (p+q==0)
                {
                    upd(x|(1<<2*j)|(2<<2*(j+1)));
                    flag=0;
                    upd(x1);
                    flag=1;
                }
                else
                {
                    if (p*q==0)
                    {
                        upd(x1|(p+q<<2*(j+1)));
                        upd(x1|(p+q<<2*j));
                    }
                    else
                    {
                        if (p==1)
                        {
                            if (q==1)
                            {
                                for (int k=j+2,cnt=0;;k++)
                                {
                                    tem=(x>>2*k)&3;
                                    if (tem==1) cnt++;
                                    if (tem==2)
                                    {
                                        if (cnt) cnt--;
                                        else
                                        {
                                            u=k;
                                            break;
                                        }
                                    }
                                }
                                upd(x1^(3<<u*2));
                            }
                            else if (!x1) ans=max(ans,dp[cur&1^1][x]+a[i][j]);
                        }
                        else
                        {
                            if (q==1) upd(x1);
                            else
                            {
                                for (int k=j-1,cnt=0;;k--)
                                {
                                    tem=(x>>2*k)&3;
                                    if (tem==1)
                                    {
                                        if (cnt) cnt--;
                                        else
                                        {
                                            u=k;
                                            break;
                                        }
                                    }
                                    if (tem==2) cnt++;
                                }
                                upd(x1^(3<<2*u));
                            }
                        }
                    }
                }
            }
        }
    printf("%d\n",ans);
}