[HNOI2007]神奇游乐园

题面在这里

题意

四通格中每格都有权值\(a_{ij}\),求单回路上路径权值和的最大值

数据范围

\[2\le n\le 100,2\le m\le 6,-10^3\le a_{ij}\le10^3\]

sol

插头DP,最后合并单回路的时候就不要记状态了,直接更新答案就好

代码

#include<bits/stdc++.h>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<iomanip>
#include<cstring>
#include<complex>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define mp make_pair
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define RG register
#define il inline
using namespace std;
typedef unsigned long long ull;
typedef vector<int>VI;
typedef long long ll;
typedef double dd;
const dd eps=1e-10;
const int mod=1e8;
const int N=1000010;
const int inf=2147483647;
il ll read(){
    RG ll data=0,w=1;RG char ch=getchar();
    while(ch!=‘-‘&&(ch<‘0‘||ch>‘9‘))ch=getchar();
    if(ch==‘-‘)w=-1,ch=getchar();
    while(ch<=‘9‘&&ch>=‘0‘)data=data*10+ch-48,ch=getchar();
    return data*w;
}

il void file(){
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
}

int n,m,mz[105][8],ans=-inf;
map<int,int>M[2];
map<int,int>::iterator tmp;
il void print(int x){putchar(48+(x&3));if(x>>2)print(x>>2);}
il void init(){
    n=read();m=read();
    memset(mz,129,sizeof(mz));
    for(RG int i=1;i<=n;i++)
        for(RG int j=1;j<=m;j++)
            mz[i][j]=read();
}

int p1,p2,p3,q1,q2,q3,k,d,pos,cal;
il void dpblank(int i,int j,int cur){
    p1=1<<(2*(j-1));q1=1<<(2*j);
    p2=2<<(2*(j-1));q2=2<<(2*j);
    p3=3<<(2*(j-1));q3=3<<(2*j);
    for(tmp=M[cur^1].begin();tmp!=M[cur^1].end();tmp++){
        k=tmp->first;d=tmp->second+mz[i][j];cal=0;
        if(!(k&p3)&&!(k&q3)&&i+1<=n&&j+1<=m){//新建,##
            k=k|p1|q2;
            if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
            else M[cur][k]=max(M[cur][k],d);
        }
        else if((k&p2)&&(k&q1)){//合并,)(
            k=k^p2^q1;
            if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
            else M[cur][k]=max(M[cur][k],d);
        }
        else if((k&p1)&&(k&q2)&&!(k^p1^q2)){//合并,()
            ans=max(ans,d);//直接统计答案
        }
        else if((k&p2)&&(k&q2)){//合并,))
            for(pos=2*(j-2);cal||((k>>pos)&3)!=1;pos-=2)
                if(((k>>pos)&3)==2)cal++;
                else if(((k>>pos)&3)==1)cal--;
            k=k^p2^q2^(3<<pos);if(j==m)k<<=2;
            if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
            else M[cur][k]=max(M[cur][k],d);
        }
        else if((k&p1)&&(k&q1)){//合并,((
            for(pos=2*(j+1);cal||((k>>pos)&3)!=2;pos+=2)
                if(((k>>pos)&3)==1)cal++;
                else if(((k>>pos)&3)==2)cal--;
            k=k^p1^q1^(3<<pos);
            if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
            else M[cur][k]=max(M[cur][k],d);
        }
        else if((k&p3)&&!(k&q3)){//保持, )# / (#
            if(i!=n){//下
                if(j==m)k<<=2;
                if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
                else M[cur][k]=max(M[cur][k],d);
            }
            if(j!=m){//右
                k=k^(k&p3)|((k&p3)<<2);
                if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
                else M[cur][k]=max(M[cur][k],d);
            }
        }
        else if(!(k&p3)&&(k&q3)){//保持, #) / #(
            if(j!=m){//右
                if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
                else M[cur][k]=max(M[cur][k],d);
            }
            if(i!=n){//下
                k=k^(k&q3)|((k&q3)>>2);if(j==m)k<<=2;
                if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
                else M[cur][k]=max(M[cur][k],d);
            }
        }
    }
}

il void dpblock(int i,int j,int cur){
    p3=3<<(2*(j-1));q3=3<<(2*j);
    for(tmp=M[cur^1].begin();tmp!=M[cur^1].end();tmp++){
        k=tmp->first;d=tmp->second;
        if(!(k&p3)&&!(k&q3)){
            if(j==m)k<<=2;
            if(M[cur].find(k)==M[cur].end())M[cur][k]=d;
            else M[cur][k]=max(M[cur][k],d);
        }
    }
}

il void solve(){
    M[0][0]=0;RG int cur=0;
    for(RG int i=1;i<=n;i++)
        for(RG int j=1;j<=m;j++)
            M[cur^=1].clear(),dpblank(i,j,cur),dpblock(i,j,cur);
    printf("%d\n",ans);
}

int main()
{
    init();solve();
    return 0;
}

原文地址:https://www.cnblogs.com/cjfdf/p/8692147.html

时间: 2024-10-13 07:21:32

[HNOI2007]神奇游乐园的相关文章

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

1187: [HNOI2007]神奇游乐园 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 668  Solved: 337[Submit][Status][Discuss] Description 经 历了一段艰辛的旅程后,主人公小P乘坐飞艇返回.在返回的途中,小P发现在漫无边际的沙漠中,有一块狭长的绿地特别显眼.往下仔细一看,才发现这是一个游 乐场,专为旅途中疲惫的人设计.娱乐场可以看成是一块大小为n×m的区域,且这个n×m的区域被分成n×m个

【BZOJ 1187】 [HNOI2007]神奇游乐园

1187: [HNOI2007]神奇游乐园 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 706 Solved: 358 [Submit][Status][Discuss] Description 经历了一段艰辛的旅程后,主人公小P乘坐飞艇返回.在返回的途中,小P发现在漫无边际的沙漠中,有一块狭长的绿地特别显眼.往下仔细一看,才发现这是一个游乐场,专为旅途中疲惫的人设计.娱乐场可以看成是一块大小为n×m的区域,且这个n×m的区域被分成n×m个小格

【BZOJ1187】[HNOI2007]神奇游乐园 插头DP

[BZOJ1187][HNOI2007]神奇游乐园 Description 经历了一段艰辛的旅程后,主人公小P乘坐飞艇返回.在返回的途中,小P发现在漫无边际的沙漠中,有一块狭长的绿地特别显眼.往下仔细一看,才发现这是一个游乐场,专为旅途中疲惫的人设计.娱乐场可以看成是一块大小为n×m的区域,且这个n×m的区域被分成n×m个小格子,每个小格子中就有一个娱乐项目.然而,小P并不喜欢其中的所有娱乐项目,于是,他给每个项目一个满意度.满意度为正时表示小P喜欢这个项目,值越大表示越喜欢.为负时表示他不喜欢

bzoj:1187: [HNOI2007]神奇游乐园

Description 经历了一段艰辛的旅程后,主人公小P乘坐飞艇返回.在返回的途中,小P发现在漫无边际的沙漠中,有一块狭长的绿地特别显眼.往下仔细一看,才发现这是一个游乐场,专为旅途中疲惫的人设计.娱乐场可以看成是一块大小为n×m的区域,且这个n×m的区域被分成n×m个小格子,每个小格子中就有一个娱乐项目.然而,小P并不喜欢其中的所有娱乐项目,于是,他给每个项目一个满意度.满意度为正时表示小P喜欢这个项目,值越大表示越喜欢.为负时表示他不喜欢,这个负数的绝对值越大表示他越不喜欢.为0时表示他对

【插头DP】BZOJ1187- [HNOI2007]神奇游乐园

[题目大意] 在n*m的网格中选一条回路,使权值和最大. [思路] 和之前裸的插头DP差不多,只不过现在回路不需要经过所有的格子.所以有以下几个注意点(具体看注释): (1)left和up插头相等的时候,直接更新答案; (2)left和up插头不存在的时候,还要考虑当前格子不取的情况. orz写了半天,感觉被掏空.jpg 1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long ll; 4 const int MA

【插头DP】 HYSBZ 1187 神奇游乐园

通道:http://www.lydsy.com/JudgeOnline/problem.php?id=1187 题意:单回路,权值,可选可不选,权值和最大. 代码: 1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 5 using namespace std; 6 7 const int MAX_N = 107; 8 const int MAX_M = 7; 9 const int HAS

【Bzoj1187】神奇游乐园(插头DP)

Description 题意给定一个矩阵,每个格子有权值,在[-1000.1000]内,求一条回路使得回路经过权值和最大,每个格子最多经过一次 2≤n≤100,2≤m≤6 Code #include <cstdio> #include <algorithm> #define N 10 #define Inf 0x7fffffff using namespace std; int A[N],now,op[N],k[N],L[N],Ans=-Inf; int n,m,g[110][N]

[入门向选讲] 插头DP:从零概念到入门 (例题:HDU1693 COGS1283 BZOJ2310 BZOJ2331)

转载请注明原文地址:http://www.cnblogs.com/LadyLex/p/7326874.html 最近搞了一下插头DP的基础知识……这真的是一种很锻炼人的题型…… 每一道题的状态都不一样,并且有不少的分类讨论,让插头DP十分锻炼思维的全面性和严谨性. 下面我们一起来学习插头DP的内容吧! 插头DP主要用来处理一系列基于连通性状态压缩的动态规划问题,处理的具体问题有很多种,并且一般数据规模较小. 由于棋盘有很特殊的结构,使得它可以与“连通性”有很强的联系,因此插头DP最常见的应用要数

计数类问题专题

主要是前两天被uoj的毛爷爷的题虐的不轻,心里很不爽啊,必须努力了,, 计数类问题分为:1.组合数学及数论计数 2.dp:状态压缩dp,插头轮廓线dp,树形dp,数位dp,普通dp 3.容斥原理 4.polya原理 5.图论计数 6.生成函数 7.其它(生成树计数等等) 本文主要研究前3个内容 考虑基本计数原理:加法原理,减法原理,乘法原理,除法原理 计数的基本原则:结果不重不漏 加法原理比较自然,中间过程有时减法原理 考虑到无向,有向图的各种量值(生成树之类)计数,状态压缩dp解决 论文:ht