bzoj1151: [CTSC2007]动物园zoo 状压dp

Description

新建的圆形动物园是亚太地区的骄傲。圆形动物园坐落于太平洋的一个小岛上,包含一大圈围栏,每个围栏里有一
种动物。如下图所示:

你是动物园的公共主管。你要做的是,让每个来动物园的人都尽可能高兴。今天有一群小朋友来动物园参观,你希
望能让他们在动物园度过一段美好的时光。但这并不是一件容易的事——有的动物有一些小朋友喜欢,有的动物有
一些小朋友害怕。如,Alex 喜欢可爱的猴子和考拉,而害怕拥牙齿锋利的狮子。而Polly 会因狮子有美丽的鬃毛
而喜欢它,但害怕有臭味的考拉。你可以选择将一些动物从围栏中移走以使得小朋友不会害怕。但你不能移走所有
的动物,否则小朋友们就没有动物可看了。每个小朋友站在大围栏圈的外面,可以看到连续的 5 个围栏。你得到
了所有小朋友喜欢和害怕的动物信息。当下面两处情况之一发生时,小朋友就会高兴:
至少有一个他害怕的动物被移走
至少有一个他喜欢的动物没被移走
例如,考虑下图中的小朋友和动物:

假如你将围栏 4 和 12 的动物移走。Alex 和 Ka-Shu 将很高兴,因为至少有一个他们害怕的动物被移走了。这也
会使 Chaitanya 高兴,因为他喜欢的围栏 6 和8 中的动物都保留了。但是,Polly 和 Hwan 将不高兴,因为他们
看不到任何他们喜欢的动物,而他们害怕的动物都还在。这种安排方式使得三个小朋友高兴。现在,换一种方法,
如果你将围栏 4 和 6 中的动物移走,Alex 和 Polly 将很高兴,因为他们害怕的动物被移走了。Chaitanya 也会
高兴,虽然他喜欢的动物 6被移走了,他仍可以看到围栏 8 里面他喜欢的动物。同样的 Hwan 也会因可以看到自
己喜欢的动物 12 而高兴。唯一不高兴的只有 Ka-Shu。如果你只移走围栏 13 中的动物,Ka-Shu 将高兴,因为有
一个他害怕的动物被移走了,Alex, Polly, Chaitanya 和 Hwan 也会高兴,因为他们都可以看到至少一个他们喜
欢的动物。所以有 5 个小朋友会高兴。这种方法使得了最多的小朋友高兴。

Input

输入的第一行包含两个整数N, C,用空格分隔。
N是围栏数(10≤N≤10 000),C是小朋友的个数(1≤C≤50 000)。
围栏按照顺时针的方向编号为1,2,3,…,N。
接下来的C行,每行描述一个小朋友的信息,
以下面的形式给出: E F L X1 X2 … XF Y1 Y2 … YL
其中: E表示这个小朋友可以看到的第一个围栏的编号(1≤E≤N),
换句话说,该小朋友可以看到的围栏为E, E+1, E+2, E+3, E+4。
注意,如果编号超过N将继续从1开始算。
如:当N=14, E=13时,这个小朋友可以看到的围栏为13,14,1, 2和3。
F表示该小朋友害怕的动物数。
L表示该小朋友喜欢的动物数。
围栏X1, X2, …, XF 中包含该小朋友害怕的动物。
围栏Y1, Y2, …, YL 中包含该小朋友喜欢的动物。
X1, X2, …, XF, Y1, Y2, …, YL是两两不同的整数,
而且所表示的围栏都是该小朋友可以看到的。
小朋友已经按照他们可以看到的第一个围栏的编号从小到大的顺序排好了
(这样最小的E对应的小朋友排在第一个,最大的E对应的小朋友排在最后一个)。
注意可能有多于一个小朋友对应的E是相同的。

Output

仅输出一个数,表示最多可以让多少个小朋友高兴

题解

看数据范围,只有小朋友看到的栏杆范围可以状压。
然后我们用\(num[i]][j]\)表示观察范围为\(i-(i+4)\)的小朋友在状态为\(j\)的情况下的高兴人数。
接着用\(f[i][j]\)表示观察范围为\(i-(i+4)\),状态为\(j\)最多的小朋友高兴人数,然后由第\(i-1\)个围栏去不去掉转移
\[
f[i][j]=max(f[i-1][(j\&15)<<1],f[i-1][(j\&15)<<1|1])+num[i][j]
\]
又因为它是一个环,我们枚举初始的状态,只有结尾处的状态与初始状态相同的才合法,来更新答案。
然后记得每次枚举答案时给别的状态致负值,防止影响答案。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1e4+10;
const int maxm=(1<<5)+10;
const int INF=0x3f3f3f3f;
int n,c,ans;
int num[maxn][maxm];
int f[maxn][maxm];
int id[maxn];
int main(){
    ios::sync_with_stdio(false);
    cin>>n>>c;
    for(int i=1;i<=c;i++){
        int e,a,f,l,x=0,y=0,ht,lv;
        cin>>e>>f>>l;a=e;
        for(int i=1;i<=5;i++){
            if(e>n)e=1;
            id[e++]=i;
        }
        for(int i=1;i<=f;i++){
            cin>>ht;ht=id[ht];
            x|=1<<(ht-1);
        }
        for(int i=1;i<=l;i++){
            cin>>lv;lv=id[lv];
            y|=1<<(lv-1);
        }
        for(int i=0;i<(1<<5);i++)
            if((i&y)||((i^((1<<5)-1))&x))
                num[a][i]++;
    }
    for(int i=0;i<(1<<5);i++){
        memset(f[0],-INF,sizeof(f[0]));
        f[0][i]=0;
        for(int j=1;j<=n;j++) for(int k=0;k<(1<<5);k++)
            f[j][k]=max(f[j-1][(k&15)<<1],f[j-1][(k&15)<<1|1])+num[j][k];
        ans=max(f[n][i],ans);
    }
    printf("%d\n",ans);
}

原文地址:https://www.cnblogs.com/Nan-Cheng/p/9748777.html

时间: 2024-11-11 19:08:09

bzoj1151: [CTSC2007]动物园zoo 状压dp的相关文章

ZOJ3305Get Sauce 状压DP,

状压DP的题目留个纪念,首先题意一开始读错了,搞了好久,然后弄好了,觉得DFS可以,最后超时,修改了很久还是超时,没办法看了一下n的范围,然后觉得状压可以,但是没有直接推出来,就记忆化搜索了一下,可是一直错,莫名奇妙,然后没办法看了一下题解,发现了下面这个比较好的方法,然后按照这个方程去推,然后敲,也是WA了好多把,写的太搓了,没人家的清楚明了,唉~也算是给自己留个纪念,状压一直做的都不太好~唉~还好理解了, 参考了  http://blog.csdn.net/nash142857/articl

poj 2411 Mondriaan&#39;s Dream(状压DP)

Mondriaan's Dream Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 12232   Accepted: 7142 Description Squares and rectangles fascinated the famous Dutch painter Piet Mondriaan. One night, after producing the drawings in his 'toilet series

(状压dp)uva 10817 Headmaster&#39;s Headache

题目地址 1 #include <bits/stdc++.h> 2 typedef long long ll; 3 using namespace std; 4 const int MAX=1e5+5; 5 const int INF=1e9; 6 int s,m,n; 7 int cost[125]; 8 //char sta[MAX]; 9 string sta; 10 int able[125]; 11 int dp[125][1<<8][1<<8]; 12 in

HDU5816 Hearthstone(状压DP)

题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5816 Description Hearthstone is an online collectible card game from Blizzard Entertainment. Strategies and luck are the most important factors in this game. When you suffer a desperate situation an

HDU 4336 容斥原理 || 状压DP

状压DP :F(S)=Sum*F(S)+p(x1)*F(S^(1<<x1))+p(x2)*F(S^(1<<x2))...+1; F(S)表示取状态为S的牌的期望次数,Sum表示什么都不取得概率,p(x1)表示的是取x1的概率,最后要加一因为有又多拿了一次.整理一下就可以了. 1 #include <cstdio> 2 const int Maxn=23; 3 double F[1<<Maxn],p[Maxn]; 4 int n; 5 int main() 6

Travel(HDU 4284状压dp)

题意:给n个城市m条路的网图,pp在城市1有一定的钱,想游览这n个城市(包括1),到达一个城市要一定的花费,可以在城市工作赚钱,但前提有工作证(得到有一定的花费),没工作证不能在该城市工作,但可以走,一个城市只能工作一次,问pp是否能游览n个城市回到城市1. 分析:这个题想到杀怪(Survival(ZOJ 2297状压dp) 那个题,也是钱如果小于0就挂了,最后求剩余的最大钱数,先求出最短路和 Hie with the Pie(POJ 3311状压dp) 送披萨那个题相似. #include <

BZOJ 1087: [SCOI2005]互不侵犯King( 状压dp )

简单的状压dp... dp( x , h , s ) 表示当前第 x 行 , 用了 h 个 king , 当前行的状态为 s . 考虑转移 : dp( x , h , s ) = ∑ dp( x - 1 , h - cnt_1( s ) , s' ) ( s and s' 两行不冲突 , cnt_1( s ) 表示 s 状态用了多少个 king ) 我有各种预处理所以 code 的方程和这有点不一样 ------------------------------------------------

BZOJ 1072 排列 状压DP

题意:链接 方法:状压DP? 题解:这题其实没啥好写的,不算很难,推一推就能搞出来. 首先看到这个问题,对于被d整除这个条件,很容易就想到是取余数为0,所以想到可能状态中刚开始含有取余数. 先说我的第一个想法,f[i][j]表示选取i个数且此时的mod为j,这样的思想是第一下蹦出来的,当时想的就是在线来搞最终的答案.不过转瞬即发现,这TM不就是暴力吗魂淡!并没有什么卵用,于是开始想这个状态可不可以做什么优化. 显然第二维的j并不需要太大的优化,暂且先将其搁置一边,来考虑第一维的i怎么优化. 把滚

COdeforces#417D Cunning Gena(状压DP)

A boy named Gena really wants to get to the "Russian Code Cup" finals, or at least get a t-shirt. But the offered problems are too complex, so he made an arrangement with his n friends that they will solve the problems for him. The participants