【NOIP 2016 提高组】蚯蚓

传送门

---

problem

本题中，我们将用符号 $\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$ 个数）。

数据范围：$n\le 10^5$，$m\le 7\times 10^6$，$t\le 71$，$a_i\le 10^8$，$q\le 200$。

---

solution

有一个比较显然的做法是用优先队列模拟这一过程。比较不好处理的是每条蚯蚓的增长，我们可以反着考虑，除了被切开的两条蚯蚓长度增加 $q$，可以看成被切开的两条蚯蚓长度减少 $q$，所以我们维护一个累计增加量，在被切开的两条蚯蚓上动手脚就可以了。

这样复杂度是 $O(m\log n)$ 的，可以拿到大约 $80$ 分。

然后 $100$ 分的做法就是注意到题目中本身的单调性，即先被切掉的蚯蚓分成的蚯蚓一定比后被切掉的蚯蚓分成的对应蚯蚓大。

证明如下：
有两只蚯蚓 $a,b$ 且原长度 $l_a>l_b$。
假设在某一时刻 $a$ 被切割，此时 $a$ 被分成了 $p{l}_a$ 和 $l_a(1-p)$ 两部分。
假设 $t$ 秒后 $b$ 被切割，此时 $b$ 的长度是 $l_b+tq$，分成的两条蚯蚓就是 $l_{b1}=p(l_b+tq)$，$b_{b2}=(l_b+tq)(1-p)$。而此时 $l_{a1}=p{l}_a+tq$，$l_{a2}=l_a(1-p)+tq$。容易得到 $l_{a1}>l_{b1}$，$l_{a_2}>l_{b2}$。
证毕。

因此我们维护三个队列，分别记录未被切割的蚯蚓，被分成 $\lfloor px\rfloor$ 段的蚯蚓和被分成 $x-\lfloor px\rfloor$ 段的蚯蚓，在里面找最长的即可。

时间复杂度 $O(m)$。

---

code

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,q,u,v,t,a[100005];
queue<int>Q[4];
int main(){
	scanf("%d%d%d%d%d%d",&n,&m,&q,&u,&v,&t);
	for(int i=1;i<=n;++i)  scanf("%d",&a[i]);
	sort(a+1,a+n+1);
	for(int i=n;i>=1;--i)  Q[1].push(a[i]);
	int del=0,tot=0;
	for(int i=1;i<=m;++i){
		int Max=-1000000000,p=0;
		if(!Q[1].empty()&&Q[1].front()>Max)  Max=Q[1].front(),p=1;
		if(!Q[2].empty()&&Q[2].front()>Max)  Max=Q[2].front(),p=2;
		if(!Q[3].empty()&&Q[3].front()>Max)  Max=Q[3].front(),p=3;
		Q[p].pop(),Max+=del,del+=q;
		if(i%t==0)  printf("%d ",Max);
		int L=1ll*Max*u/v,R=Max-L;
		L-=del,R-=del;
		Q[2].push(L),Q[3].push(R);
	}puts("");
	for(int i=1;i<=n+m;++i){
		int Max=-1000000000,p=0;
		if(!Q[1].empty()&&Q[1].front()>Max)  Max=Q[1].front(),p=1;
		if(!Q[2].empty()&&Q[2].front()>Max)  Max=Q[2].front(),p=2;
		if(!Q[3].empty()&&Q[3].front()>Max)  Max=Q[3].front(),p=3;
		Q[p].pop();
		if(i%t==0)  printf("%d ",Max+del);
	}
	return 0;
}