NOIP2016 DAY2 T2蚯蚓
传送门
Description
本题中,我们将用符号[c]表示对c向下取整,例如:[3.0」= [3.1」=[3.9」=3。蛐蛐国最近蚯蚓成灾了!隔壁跳
蚤国的跳蚤也拿蚯蚓们没办法,蛐蛐国王只好去请神刀手来帮他们消灭蚯蚓。蛐蛐国里现在共有n只蚯蚓(n为正整
数)。每只蚯蚓拥有长度,我们设第i只蚯蚓的长度为a_i(i=1,2,...,n),并保证所有的长度都是非负整数(即:可
能存在长度为0的蚯蚓)。每一秒,神刀手会在所有的蚯蚓中,准确地找到最长的那一只(如有多个则任选一个)
将其切成两半。神刀手切开蚯蚓的位置由常数p(是满足0<p<1的有理数)决定,设这只蚯蚓长度为x,神刀手会将其
切成两只长度分别为[px]和x-[px]的蚯蚓。特殊地,如果这两个数的其中一个等于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个非负整数,为ai,a2,...,an,即初始时n只蚯蚓的长度。
同一行中相邻的两个数之间,恰好用一个空格隔开。
保证1<=n<=10^5,0<m<7*10^6,0<u<v<10^9,0<=q<=200,1<t<71,0<ai<10^8。
Output
第一行输出[m/t]个整数,按时间顺序,依次输出第t秒,第2t秒,第3t秒……被切断蚯蚓(在被切断前)的长度。
第二行输出[(n+m)/t]个整数,输出m秒后蚯蚓的长度;需要按从大到小的顺序
依次输出排名第t,第2t,第3t……的长度。
同一行中相邻的两个数之间,恰好用一个空格隔开。即使某一行没有任何数需要 输出,你也应输出一个空行。
请阅读样例来更好地理解这个格式。
Sample Input
3 7 1 1 3 1
3 3 2
3 3 2
Sample Output
3 4 4 4 5 5 6
6 6 6 5 5 4 4 3 2 2
6 6 6 5 5 4 4 3 2 2
感觉挺莫名奇妙的.....
在洛谷上AC了这题,然后去下载官方数据, 居然超时??
也许是我这台机子太慢了.....................
正解:
第一遍看完题目时就觉得用优先队列,但想想.... CCF会这么容易给你满分吗?? 做梦吧。。。
优先队列65分。
正解是用三个普通的队列做的。
一开始将所有蚯蚓读入,放进一个队列并且排序,从大到小。
然后拿出最长的一条,进行切割,将 x*q 放在一个队列,将x - x*q放在另一个队列。
经过证明,这两个队列都是单调的。
利用这一点,我们只需每次从这三个队列里面挑出一个最大的蚯蚓进行切割然后放进那两个队列就行了。
下面贴代码,有问题留言。
#include<cstdio> #include<algorithm> #define N 20000900 using namespace std; int que[3][N],head[3],tail[3]; int cmp(const int &a,const int &b){ return a > b; } int read(){ int x = 0; char ch = getchar(); while(ch < '0' || ch > '9')ch = getchar(); while(ch >= '0' && ch <= '9'){ x = x * 10 + ch -'0'; ch = getchar(); } return x; } int main(){ int n = read(),m = read(),q = read(),u = read(),v = read(),t = read(); double p = 1.0*u/v; for(int i = 1; i <= n; i++)que[0][i] = read(); sort(que[0]+1,que[0]+n+1,cmp); int first = 1,len = 0; head[0] = head[1] = head[2] = 1; tail[0] = n,tail[1] = tail[2] = 0; for(int i = 1; i <= m; i++){ long long mx = -99999999999,pos = 0; for(int j = 0; j <= 2; j++) if(head[j] <= tail[j] && que[j][head[j]] > mx){ mx = que[j][head[j]]; pos = j; } head[pos]++; mx += len; if(i % t == 0){ if(first)first = 0; else printf(" "); printf("%d",mx); } int x1 = (int)(mx*u/v); int x2 = mx - x1; que[1][++tail[1]] = x1 - len-q; que[2][++tail[2]] = x2 - len-q; len += q; } printf("\n"); first = 1; for(int i = 1; i <= n + m; i++){ long long mx = -9999999999999,pos = 0; //mx要设小一点,不然会爆 for(int j = 0; j <= 2; j++) if(head[j] <= tail[j] && que[j][head[j]] > mx){ mx = que[j][head[j]]; pos = j; } head[pos]++; mx += len; if(i % t == 0){ if(first)first = 0; else printf(" "); printf("%lld",mx); } } printf("\n"); return 0; }