洛谷 P4660 & bzoj 1168 [ Baltic OI 2008 ] 手套 —— 分析+单调栈
题目:https://www.luogu.org/record/show?rid=12702916
https://www.lydsy.com/JudgeOnline/problem.php?id=1168
一眼不可做...即使数据范围很小...
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const xn=;
int n;
ll a[xn],b[xn],inf=1e17;
bool cmp(int x,int y){return b[x]>b[y];}
int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++)scanf("%lld",&a[i]);
for(int i=;i<=n;i++)scanf("%lld",&b[i]);
ll ansa=,ansb=; int tot=;
for(int i=;i<=n;i++)
{
if(a[i]&&!b[i])ansa+=a[i];
else if(!a[i]&&b[i])ansb+=b[i];
else if(a[i]&&b[i])a[++tot]=a[i],b[tot]=b[i];
}
int mx=(<<tot); ll pa=inf,pb=inf;
for(int s=;s<mx;s++)
{
ll tmpa=,tmpb=,mxa=inf,mxb=inf; int cnta=;
for(int i=;i<=tot;i++)
{
if(s&(<<(i-)))tmpa+=a[i],mxa=min(mxa,a[i]-),cnta++;
else tmpb+=b[i],mxb=min(mxb,b[i]-);
}
if(cnta==tot)tmpa-=mxa; if(cnta==)tmpb-=mxb;
if(!cnta)tmpa++; else tmpb++;
// if(cnta==tot)tmpa-=mxa,tmpb++;
// else if(cnta==0)tmpb-=mxb,tmpa++;
// else tmpa-=mxa,tmpb++;
if(tmpa+tmpb<pa+pb)pa=tmpa,pb=tmpb;
}
printf("%lld\n%lld\n",pa+ansa,pb+ansb);
return ;
}
还写过10分的错误思路
首先,设每种颜色左右手套个数为 l[i],r[i],如果有 l[i]=0 或 r[i]=0,显然最劣情况要选,所以提前处理好;
设一个状态 ans(L),表示左手套取 L 个时的答案;
显然,这个答案是唯一的而且最劣的,现在要找到这个答案;
取 L 个左手套,根据取法的不同,可以取出许多种颜色集合,把每一个可能的集合记作 T;
而这些颜色集合又组成了一个集合,记作 A(L),表示选 L 个左手套会产生的颜色集合的集合;
考虑每个集合 T,在右边都对应一种最劣的选法,就是把不在 T 中的颜色都选了,最后再+1;
设这个为 R(T),即 R(T) = (∑r[i]) + 1,其中 i 不在集合 T 中;
所以,ans(L) = max{ R(T) },其中 T 是 A(L) 中的一个集合,这样做很妙地没有限制什么而找到了最劣答案,保证了正确性;
现在,考虑 A(L) 中有哪些 T,发现:
A(L) = { T | |T| <= L <= S(T) },其中 |T| 是 T 这个集合包含的颜色数量,S(T) 是 T 这个集合包含的那些颜色的手套个数的和;
由于我们找的是 max{ R(T) },而 |T| 越小,R(T) 越大,所以最终 ans(L) 的取值会是合法范围内某个 |T| 最小的,所以可以暂时忽略 |T| <= L 这个条件!!
所以,现在就变成若 S(T) >= L,则 R(T) 可以更新 ans(L);
然后,可以发现,会成为最终答案的 L,应该是某些颜色的手套全选,最后再+1,否则...应该不是最优的(感性理解...);
所以处理出这种 L,找 S(T) >= L 中 R(T) 最大的,可以用单调栈实现;
用结构体存一种集合的 S(T) 和 R(T),S(T) 同时也是上面所说的 L,按 S(T) 从小到大排序,不断弹栈令栈内 R(T) 是单调递减的,就能找到一个结构体后面(S(T) >= L)的 R(T) 最大值,对于一个 R(T) 也能去更新前面同等条件最小的 L;
还有些疑惑之处...就是一些 L 相同的结构体,排序时按 r<y.r 或 r>y.r 都可以过...但是感觉应该是 r<y.r 才对!后面就可以把前面都弹掉,则栈里后面元素的 S(T) 一定 >= L+1;
但如果栈里最后只留下一个值...难道因为这样构造所以不会出现这种情况?
总之勉强A了,思路还是很妙的。
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int const xn=;
int n,a[xn],b[xn],ansl,ansr,sta[<<],top;
struct N{
int l,r;
bool operator < (const N &y) const
{return l==y.l?r<y.r:l<y.l;}//
}f[<<];
int rd()
{
int ret=,f=; char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
while(ch>=''&&ch<='')ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
return f?ret:-ret;
}
void update(int x,int y)
{
if(ansl+ansr>x+y||(ansl+ansr==x+y&&ansl>x))
ansl=x,ansr=y;
}
int main()
{
n=rd(); int pl=,pr=,tot=;
for(int i=;i<=n;i++)a[i]=rd();
for(int i=;i<=n;i++)b[i]=rd();
for(int i=;i<=n;i++)
{
if(a[i]&&!b[i])pl+=a[i];
else if(!a[i]&&b[i])pr+=b[i];
else if(a[i]&&b[i])a[++tot]=a[i],b[tot]=b[i];
}
int mx=(<<tot);
for(int i=;i<mx;i++)
for(int j=;j<=tot;j++)
{
if(i&(<<(j-)))f[i].l+=a[j];
else f[i].r+=b[j];
}
sort(f,f+mx);
for(int i=;i<mx;i++)
{
while(top&&f[sta[top]].r<f[i].r)top--;
sta[++top]=i;
}
ansl=1e9; ansr=1e9;
for(int i=;i<top;i++)update(f[sta[i]].l+,f[sta[i+]].r+);
printf("%d\n%d\n",ansl+pl,ansr+pr);
return ;
}
洛谷 P4660 & bzoj 1168 [ Baltic OI 2008 ] 手套 —— 分析+单调栈的更多相关文章
- 洛谷10月月赛Round.1| P3400 仓鼠窝[单调栈]
题目描述 萌萌哒的Created equal是一只小仓鼠,小仓鼠自然有仓鼠窝啦. 仓鼠窝是一个由n*m个格子组成的行数为n.列数为m的矩阵.小仓鼠现在想要知道,这个矩阵中有多少个子矩阵!(实际上就是有 ...
- 洛谷P2178 品酒大会【后缀数组】【单调栈】
题目描述 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战 两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个奖项,吸引了众多品酒师参加. 在大会的晚餐上,调酒师 Rainb ...
- 【洛谷5294】[HNOI2019] 序列(主席树维护单调栈+二分)
点此看题面 大致题意: 给你一个长度为\(n\)的序列\(A\),每次询问修改一个元素(只对当前询问有效),然后让你找到一个不下降序列\(B\),使得这两个序列相应位置之差的平方和最小,并输出这个最小 ...
- 洛谷 P4198 楼房重建 线段树维护单调栈
P4198 楼房重建 题目链接 https://www.luogu.org/problemnew/show/P4198 题目描述 小A的楼房外有一大片施工工地,工地上有N栋待建的楼房.每天,这片工地上 ...
- 洛谷P2178 [NOI2015]品酒大会 后缀数组+单调栈
P2178 [NOI2015]品酒大会 题目链接 https://www.luogu.org/problemnew/show/P2178 题目描述 一年一度的"幻影阁夏日品酒大会" ...
- 洛谷 P2587 BZOJ 1034 [ZJOI2008]泡泡堂
题目描述 //不知道为什么BZOJ和洛谷都没有这幅图了,大牛们几年前的博客上都有这幅图的,把它贴上来吧 第XXXX届NOI期间,为了加强各省选手之间的交流,组委会决定组织一场省际电子竞技大赛,每一个省 ...
- 洛谷 P3307: bzoj 3202: [SDOI2013] 项链
题目传送门:洛谷P3307.这题在bzoj上是权限题. 题意简述: 这题分为两个部分: ① 有一些珠子,每个珠子可以看成一个无序三元组.三元组要满足三个数都在$1$到$m$之间,并且三个数互质,两个珠 ...
- 洛谷 4106 / bzoj 3614 [HEOI2014]逻辑翻译——思路+类似FWT
题目:https://www.luogu.org/problemnew/show/P4106 https://www.lydsy.com/JudgeOnline/problem.php?id=3614 ...
- 洛谷 P3332 BZOJ 3110 [ZJOI2013]K大数查询
题目链接 洛谷 bzoj 题解 整体二分 Code #include<bits/stdc++.h> #define LL long long #define RG register usi ...
随机推荐
- 转: 多版本并发控制(MVCC)在分布式系统中的应用 (from coolshell)
from: http://coolshell.cn/articles/6790.html 问题 最近项目中遇到了一个分布式系统的并发控制问题.该问题可以抽象为:某分布式系统由一个数据中心D和若干业务 ...
- visual c++ 2013进行MySQL编程(ODBC) -- (一) 套装安装
最近写了有些技术类文章了,只因为最近研究多了些东西,有一些项目用到了,所以自己记录一下,怕自己忘记,如果有哪位同学有自己的见解,可以邮件或者回复,技术类的探讨,不管对否,都是欢迎的. 操作之前,必须安 ...
- Java中HashMap遍历的两种方法(转)
第一种: Map map = new HashMap(); Iterator iter = map.entrySet().iterator(); while (iter.hasNext()) { Ma ...
- 如何将安卓手机中已安装的应用导出一个apk安装包
首先安装ES文件浏览器 然后应用里找到你要导出apk的应用,长按应用,备份应用, 然后在upbacks/apps文件夹里找吧 不得不说 ES浏览器还是屌啊
- 腾讯云centos,nginx安装
- windows下route命令详解(转载)
1.具体功能 该命令用于在本地IP路由表中显示和修改条目.使用不带参数的ROUTE可以显示帮助. 2.语法详解 route [-f] [-p] [co ...
- MonoTouch.Dialog简介
新建一个Single View Application项目 添加程序集 MonoTouch.Dialog.dll引用 删除 MainStoryboard.storyboard 添加空类Task.cs ...
- Android-Bundle的说明和用法
1.Bundle类的作用 Bundle类是一种数据载体,类似于Map,用于存放key-value名值对形式的值.相对于Map,它提供了各种常用类型的putXxx()/getXxx()方法, 如:put ...
- 02-合并frame
进入IJKMediaPlayer工程--->
- Xcode 6 IDE
本文转载至 http://www.cocoachina.com/ios/20140823/9442.html (via:苹果开发者中心) Xcode IDE 是 Apple 开发体验的核心.X ...