bzoj 1187: [HNOI2007]神奇游乐园 插头dp

1187: [HNOI2007]神奇游乐园

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 668  Solved: 337
[Submit][Status][Discuss]

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&
lt;=m<=6。接下来的n行，每行有m个整数，第i行第j列表示游乐场的第i行第j列的小格子中的娱乐项目的满意度，这个满意度的范围是
[-1000，1000]。同一行的两个整数之间用空格隔开。

Output

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

Sample Input

4 4
100 300 -400 400
-100 1000 1000 1000
-100 -100 -100 -100
-100 -100 -100 1000

Sample Output

4000

HINT

大家测下这个数据 5 5 1 1 -100 3 3 1 1 -100 3 3 1 1 -100 3 3 1 1 -100 3 3 1 1 -100 3 3 结果是30?

Source

　　普通的单回路插头dp，注意结束时判断无插头是在两个地方（行中，行末）。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define MAXN 120
#define MAXM 7
#define MAXV 1<<MAXM<<MAXM
#define MAXL 400
#define INF 0x3f3f3f3f
#define getp(v,p) (v>>(p)>>(p)&3)
#define clrp(v,p) (v&(~(3<<(p)<<(p))))
#define update1(v,d) dp[i+1][0][plst[v]]=max(dp[i+1][0][plst[v]],dp[i][j][k]+mp[i][j]*d);
#define update2(v,d) dp[i][j+1][plst[v]]=max(dp[i][j+1][plst[v]],dp[i][j][k]+mp[i][j]*d);
int mp[MAXN][MAXN];
int dp[MAXN][MAXM][MAXV];
int lst[MAXL];
int lstp[MAXL][MAXM];
int plst[MAXV];
int totl=;
int stack[MAXN],tops=-;
int n,m;
void dfs(int now,int sum,int status)
{
        if (now==m+)
        {
                if (!sum)
                {
                        plst[status]=totl;
                        lst[totl]=status;
                        for (int i=;i<=m;i++)
                        {
                /*                switch(status%4)
                                {
                                        case 0:printf(".");break;
                                        case 1:printf("(");break;
                                        case 2:printf(")");break;
                                }*/
                                if (status%==)
                                {
                                        stack[++tops]=i;
                                }else if (status%==)
                                {
                                        lstp[totl][i]=stack[tops];
                                        lstp[totl][stack[tops]]=i;
                                        tops--;
                                }
                                status/=;
                        }
                //        printf("%d\n",totl);
                        totl++;
                }
                return ;
        }
        dfs(now+,sum,status);
        dfs(now+,sum+,status+(<<now<<now));
        if (sum)
                dfs(now+,sum-,status+(<<now<<now));
}

int main()
{
        freopen("input.txt","r",stdin);
        scanf("%d%d",&n,&m);
        for (int i=;i<n;i++)
        {
                for (int j=;j<m;j++)
                {
                        scanf("%d",&mp[i][j]);
                }
        }
        dfs(,,);
        for (int i=;i<n;i++)
                for (int j=;j<m;j++)
                        for (int k=;k<totl;k++)
                                dp[i][j][k]=-INF;
        register int p,q,v;
        int ans=-INF;
        dp[][][]=;
        for (register int i=;i<n;i++)
        {
                for (register int j=;j<m;j++)
                {
                        for (register int k=;k<totl;k++)
                        {
                                if (dp[i][j][k]==-INF)continue;
                                p=getp(lst[k],j);
                                q=getp(lst[k],j+);
                                v=clrp(lst[k],j);
                                v=clrp(v,j+);
                                if (j==m-)
                                {
                                        if (p== && q==)
                                        {
                                                update1((v)<<,);
                                        }else if ((p== && q==) || (p== && q==))
                                        {
                                                update1((v+(<<j<<j))<<,);
                                        }else if (p== && q== && !v)
                                        {
                                                ans=max(ans,dp[i][j][k]+mp[i][j]);
                                        }else if (p== && q==)
                                        {
                                                update1((v+(<<lstp[k][j]<<lstp[k][j]))<<,);
                                        }
                                }else
                                {
                                        if (p== && q==)
                                        {
                                                update2(v,);
                                                update2(v+(<<j<<j)+(<<(j+)<<(j+)),);
                                        }else if ((p== && q==) || (p== && q==))
                                        {
                                                update2(v+(<<j<<j),);
                                                update2(v+(<<(j+)<<(j+)),);
                                        }else if ((p== && q==) || (p== && q==))
                                        {
                                                update2(v+(<<j<<j),);
                                                update2(v+(<<(j+)<<(j+)),);
                                        }else if (p== && q==)
                                        {
                                                update2(v-(<<lstp[k][j+]<<lstp[k][j+]),);
                                        }else if (p== && q== && !v)
                                        {
                                                ans=max(ans,dp[i][j][k]+mp[i][j]);
                                        }else if (p== && q==)
                                        {
                                                update2(v,);
                                        }else if (p== && q==)
                                        {
                                                update2(v+(<<lstp[k][j]<<lstp[k][j]),);
                                        }
                                }
                        }
                }
        }
        printf("%d\n",ans);
}