【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);}