【队列】【P2827】【NOIP2016D2T3】蚯蚓
Description
本题中,我们将用符号 $\lfloor c \rfloor$ 表示对 $c$ 向下取整,例如:$\lfloor 3.0 \rfloor = \lfloor 3.1 \rfloor = \lfloor 3.9 \rfloor = 3$。
蛐蛐国最近蚯蚓成灾了!隔壁跳蚤国的跳蚤也拿蚯蚓们没办法,蛐蛐国王只好去请神刀手来帮他们消灭蚯蚓。
蛐蛐国里现在共有 $n$ 只蚯蚓($n$ 为正整数)。每只蚯蚓拥有长度,我们设第 $i$ 只蚯蚓的长度为 $a_i$ ($i=1,2,\dots,n$),并保证所有的长度都是非负整数(即:可能存在长度为 $0$ 的蚯蚓)。
每一秒,神刀手会在所有的蚯蚓中,准确地找到最长的那一只(如有多个则任选一个)将其切成两半。神刀手切开蚯蚓的位置由常数 $p$(是满足 $0 < p < 1$ 的有理数)决定,设这只蚯蚓长度为 $x$,神刀手会将其切成两只长度分别为 $\lfloor px \rfloor$ 和 $x - \lfloor px \rfloor$ 的蚯蚓。特殊地,如果这两个数的其中一个等于 $0$,则这个长度为 $0$ 的蚯蚓也会被保留。此外,除了刚刚产生的两只新蚯蚓,其余蚯蚓的长度都会增加 $q$(是一个非负整常数)。
蛐蛐国王知道这样不是长久之计,因为蚯蚓不仅会越来越多,还会越来越长。蛐蛐国王决定求助于一位有着洪荒之力的神秘人物,但是救兵还需要 $m$ 秒才能到来……($m$ 为非负整数)
蛐蛐国王希望知道这 $m$ 秒内的战况。具体来说,他希望知道:
- $m$ 秒内,每一秒被切断的蚯蚓被切断前的长度(有 $m$ 个数);
- $m$ 秒后,所有蚯蚓的长度(有 $n + m$ 个数)。
蛐蛐国王当然知道怎么做啦!但是他想考考你……
## Input
第一行包含六个整数 $n,m,q,u,v,t$,其中:$n,m,q$ 的意义见【问题描述】;$u,v,t$ 均为正整数;你需要自己计算 $p=u / v$(保证 $0 < u < v$);$t$ 是输出参数,其含义将会在【输出格式】中解释。
第二行包含 $n$ 个非负整数,为 $a_1, a_2, \dots, a_n$,即初始时 $n$ 只蚯蚓的长度。
同一行中相邻的两个数之间,恰好用一个空格隔开。
## Output
第一行输出 $\left \lfloor \frac{m}{t} \right \rfloor$ 个整数,按时间顺序,依次输出第 $t$ 秒,第 $2t$ 秒,第 $3t$ 秒,……被切断蚯蚓(在被切断前)的长度。
第二行输出 $\left \lfloor \frac{n+m}{t} \right \rfloor$ 个整数,输出 $m$ 秒后蚯蚓的长度;需要按从大到小的顺序,依次输出排名第 $t$,第 $2t$,第 $3t$,……的长度。
同一行中相邻的两个数之间,恰好用一个空格隔开。即使某一行没有任何数需要输出,你也应输出一个空行。
请阅读样例来更好地理解这个格式。
Sample Input
3 7 1 1 3 1
3 3 2
Sample Output
3 4 4 4 5 5 6
6 6 6 5 5 4 4 3 2 2
Hint
保证 $1 \leq n \leq 10^5$,$0 \leq m \leq 7 \times 10^6$,$0 < u < v \leq 10^9$,$0 \leq q \leq 200$,$1 \leq t \leq 71$,$0 \leq a_i \leq 10^8$。
Solution
看看题意每次取出一堆数中最大的进行操作,貌似是个堆,直接操作的话显然复杂度爆炸。考虑优化
对于所有的蚯蚓每秒都是在增长的。所以不妨将被切的蚯蚓看做短了q。每次输出的时候加上总共伸长的长度即可。使用堆维护他们的原长度。压入堆的时候将长度-q。这样时间复杂度达到了\(IO(~mlong(m+n)~)\)期望得分在70-95之间。
考虑整个序列是满足一定单调性的。
对于一条蚯蚓\(x\),我们将切下来的第一块叫做\(x_1\),另一块叫做\(x_2\)。考虑有两条蚯蚓\(a,b\),首先将\(a\)切成两节,在\(t\)秒后将\(b\)切成两节,那么\(Len_{a_1}=p*a*t,Len_{b_1}=p*b*t\)。因为\(a \geq b\)所以\(Len_{a_1}~\geq~Len_{b_1}\)。类似的可以证明另一半之间的大小关系。
通过简单的数学归纳我们可以得到,将条蚯蚓分成的两部分组成序列,序列都是具有单调性的。所以可以直接采用三个队列,分别是排序后的原长度,切得第一部分,切得第二部分,每次取出队头进行比较。就可以完成了。因为所有元素都入队出队1次。所以复杂度是\(O(n+m)\),考虑到排序与m近似同阶,可以认为复杂度就是\(O(n+m)\)
Code
#include<cmath>
#include<queue>
#include<cstdio>
#include<algorithm>
#define rg register
#define ci const int
#define cl const long long int
typedef long long int ll;
namespace IO {
char buf[50];
}
template<typename T>
inline void qr(T &x) {
char ch=getchar(),lst=' ';
while(ch>'9'||ch<'0') lst=ch,ch=getchar();
while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
if (lst=='-') x=-x;
}
template<typename T>
inline void write(T x,const char aft,const bool pt) {
if(x<0) {putchar('-');x=-x;}
int top=0;
do {
IO::buf[++top]=x%10+'0';
x/=10;
} while(x);
while(top) putchar(IO::buf[top--]);
if(pt) putchar(aft);
}
template <typename T>
inline T mmax(const T a,const T b) {if(a>b) return a;return b;}
template <typename T>
inline T mmin(const T a,const T b) {if(a<b) return a;return b;}
template <typename T>
inline T mabs(const T a) {if(a<0) return -a;return a;}
template <typename T>
inline void mswap(T &a,T &b) {T temp=a;a=b;b=temp;}
const int maxn = 100010;
int n,m,q,u,v,t,ad;
int MU[maxn];
std::queue<double>beg,big,sml;
inline bool cmp(const int &_a,const int &_b) {
return _a>_b;
}
int main() {
qr(n);qr(m);qr(q);qr(u);qr(v);qr(t);
for(rg int i=1;i<=n;++i) qr(MU[i]);
std::sort(MU+1,MU+1+n,cmp);
for(rg int i=1;i<=n;++i) beg.push(MU[i]);
for(rg int i=1;i<=m;++i) {
rg int maxx=-0x7fffffff,p=0;
if(!beg.empty()) {
if(beg.front()>maxx) {
maxx=beg.front();p=1;
}
}
if(!big.empty()) {
if(big.front()>maxx) {
maxx=big.front();p=2;
}
}
if(!sml.empty()) {
if(sml.front()>maxx) {
maxx=sml.front();p=3;
}
}
switch(p) {
case 1:
beg.pop();break;
case 2:
big.pop();break;
case 3:
sml.pop();break;
}
maxx+=ad;
if(!(i%t)) write(maxx,' ',true);
int _a=floor(1.0*(maxx)*u/v),_b=maxx-_a;
_a-=ad;_b-=ad;
big.push(_a-q);sml.push(_b-q);
ad+=q;
}
putchar('\n');
rg int ss=n+m;
for(rg int i=1;i<=ss;++i) {
rg int maxx=-0x7fffffff,p=0;
if(!beg.empty()) {
if(beg.front()>maxx) {
maxx=beg.front();p=1;
}
}
if(!big.empty()) {
if(big.front()>maxx) {
maxx=big.front();p=2;
}
}
if(!sml.empty()) {
if(sml.front()>maxx) {
maxx=sml.front();p=3;
}
}
switch(p) {
case 1:
beg.pop();break;
case 2:
big.pop();break;
case 3:
sml.pop();break;
}
if(!(i%t)) {
write(maxx+ad,' ',true);
}
}
putchar('\n');
return 0;
}
Summary
1、等效思想的运用:全部增长等于单个缩短
2、局部特殊性质的发现:对于每种蚯蚓都隐含单调性
【队列】【P2827】【NOIP2016D2T3】蚯蚓的更多相关文章
- 【题解】洛谷P2827 [NOIP2016TG] 蚯蚓(优先队列)
题目来源:洛谷P2827 思路 阅读理解题 一开始以为是裸的优先队列而已 但是发现维护一个切开并且其他的要分别加上一个值很不方便 而且如果直接用优先队列会TLE3到4个点 自测85分 所以我们需要发现 ...
- 洛谷 P2827 蚯蚓 题解
每日一题 day32 打卡 Analysis 我们可以想一下,对于每一秒除了被切的哪一个所有的蚯蚓都增长Q米,我们来维护3个队列,队列1表示最开始的蚯蚓,队列2表示每一次被切的蚯蚓被分开的较长的那一部 ...
- [noip2016]蚯蚓<单调队列+模拟>
题目链接:https://vijos.org/p/2007 题目链接:https://www.luogu.org/problem/show?pid=2827#sub 说实话当两个网站给出AC后,我很感 ...
- 【luoguP2827】 蚯蚓
题目描述 本题中,我们将用符号\rfloor⌊c⌋ 表示对 cc 向下取整,例如:\lfloor 3.0 \rfloor = \lfloor 3.1 \rfloor = \lfloor 3.9 \rf ...
- C++头文件居然可以这么打!!!! 长见识了!!!
返回主页 longdie 这人,生于天,立于地,为的就是顶天立地. 未来的答案早已被宇宙计算好了,人类自出现,答案就在那里,人类灭亡了,答案也在那里,,但是人活着,不就是为了看看未来发生了什么吗?如果 ...
- 【BZOJ】4721: [Noip2016]蚯蚓 / 【洛谷】P2827 蚯蚓(单调队列)
Description 本题中,我们将用符号[c]表示对c向下取整,例如:[3.0」= [3.1」=[3.9」=3.蛐蛐国最近蚯蚓成灾了!隔壁跳 蚤国的跳蚤也拿蚯蚓们没办法,蛐蛐国王只好去请神刀手来帮 ...
- 【Luogu】P2827蚯蚓(堆转队列)
按照国际惯例先发题目链接 woc从4月就开始做这sb题.最开始30分升到65分不管了,直到最近几天升到85分,再到今天AC.激动的心情自然是那些一遍就A或者一小时以内就A的神犇难以想象的. 下面说说 ...
- Luogu P2827 [NOIp2016提高组]蚯蚓 | 神奇的队列
题目链接 80分思路: 弄一个优先队列,不停地模拟,切蚯蚓时就将最长的那一条出队,然后一分为二入队,简单模拟即可.还要弄一个标记,表示从开始到当前时间每一条蚯蚓应该加上的长度,操作时就加上,入队时就减 ...
- 洛谷P2827 [NOIP2016 提高组] 蚯蚓 (二叉堆/队列)
容易想到的是用二叉堆来解决,切断一条蚯蚓,其他的都要加上一个值,不妨用一个表示偏移量的delta. 1.取出最大的x,x+=delta: 2.算出切断后的两个新长度,都减去delta和q: 3.del ...
随机推荐
- Siki_Unity_2-10_数据结构与算法
Unity 2-10 数据结构与算法 任务1-1:数据结构简介 数据结构:数据存储的结构,数据之间的关系 数据结构分类: 集合:同属于一个集合 线性结构:数据元素存在一对一的关系 树形结构:数据元素存 ...
- VMware实现控制台功能(VMware Remote Console)
说明: 刚开始一脸懵逼,google了一些资料,发现基本没有能快速落地的,自己做完后梳理了一下发上来供大家参考. 如果帮到你了,请点赞评论关注,以资鼓励,多谢~ 实现VMware控制台功能主要有两种方 ...
- 传入中文参数-->服务器_转码的方法
如果要传入 中文参数到 服务器 使用lr_convert_string_encoding() LR_ENC_SYSTEM_LOCALE , 转为 ...
- Python入门(4)
一.while循环 有时候,你可能需要计算机来帮重复做一件事,这时就需要循环. while condition: statements (else: statements ) 当condition条件 ...
- 关于mysql开元数据库的几个随想
现在已经是凌晨了,昨天晚上写了我人生中的第一篇笔记,觉得没什么可写的,写了一个多小时都没写出什么,现在突然想写点东西了,这是一个比较有趣的问题,前两个月换了新工作,记得当初面试这份工作的时候面试到第三 ...
- 京东2018秋招c++岗 神奇数
题意大概是: 一个数比如242,把所有数字分成两组,而且两组的和相等,那么这个数就是神奇数,此时242,能够分成{2,2}和{4},所以242是神奇数. 题目要求输入n和m求[n,m]区间内神奇数的个 ...
- POJ 2208 Pyramids(求四面体体积)
Description Recently in Farland, a country in Asia, a famous scientist Mr. Log Archeo has discovered ...
- Windows Phone编程回顾
前言 已有一年多没有碰WP相关的开发了. 近期经常看博客园的文章, 发现开发WP应用的同学很多, 其中博问频道关于"WPF", "C#", "WP8& ...
- "Hello world!"团队第八次会议
Scrum会议 今天是我们"Hello world!"团队第八次召开会议,博客内容是: 1.会议时间 2.会议成员 3.会议地点 4.会议内容 5.todo list 6.会议照片 ...
- 欢迎来怼--第二十一次Scrum会议
一.小组信息 队名:欢迎来怼 小组成员 队长:田继平 成员:李圆圆,葛美义,王伟东,姜珊,邵朔,冉华 小组照片 二.开会信息 时间:2017/11/2 17:05~17:15,总计10min. 地点: ...