bzoj4721[Noip2016]蚯蚓

Posted 2020-08-24 GXZlegend

题目描述

本题中，我们将用符号[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个数)。蛐蛐国王当然知道怎么做啦！但是他想考考你......

输入

第一行包含六个整数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。

输出

第一行输出[m/t]个整数，按时间顺序，依次输出第t秒，第2t秒，第3t秒……被切断蚯蚓（在被切断前）的长度。

第二行输出[(n+m)/t]个整数，输出m秒后蚯蚓的长度；需要按从大到小的顺序

依次输出排名第t，第2t，第3t……的长度。

同一行中相邻的两个数之间，恰好用一个空格隔开。即使某一行没有任何数需要 输出，你也应输出一个空行。

请阅读样例来更好地理解这个格式。

样例输入

3 7 1 1 3 1
3 3 2

样例输出

3 4 4 4 5 5 6
6 6 6 5 5 4 4 3 2 2

---

题解

队列题

由于先切的两段一定大于等于后切的对应的两段，

于是可以用两个队列分别存储[px]段和x-[px]段，还有一个存储原始长度。

每次把所有没切的长度+q，相当于把切了的长度-q，然后查找时再一起加上。

注意除法需要long long。

还有，在bzoj上，可能会PE无数次，只能按照样例格式输出，不能多任何的空格和换行。

 1 #include <cstdio>
 2 #include <algorithm>
 3 #define inf 0x7fffffff
 4 using namespace std;
 5 int q[3][8000001] , l[3] , r[3] = {-1 , -1 , -1} , a[100001];
 6 bool cmp(int a , int b)
 7 {
 8     return a > b;
 9 }
10 int getmax()
11 {
12     int k = -inf , ans = 0;
13     if(l[0] <= r[0] && q[0][l[0]] > k) k = q[0][l[0]] , ans = 0;
14     if(l[1] <= r[1] && q[1][l[1]] > k) k = q[1][l[1]] , ans = 1;
15     if(l[2] <= r[2] && q[2][l[2]] > k) k = q[2][l[2]] , ans = 2;
16     return ans;
17 }
18 int main()
19 {
20     int n , m , g , u , v , t , i , p , c;
21     scanf("%d%d%d%d%d%d" , &n , &m , &g , &u , &v , &t);
22     for(i = 0 ; i < n ; i ++ )
23         scanf("%d" , &a[i]);
24     sort(a , a + n , cmp);
25     for(i = 0 ; i < n ; i ++ )
26         q[0][++r[0]] = a[i];
27     for(i = 1 ; i <= m ; i ++ )
28     {
29         p = getmax();
30         c = q[p][l[p]++] + g * (i - 1);
31         if(i % t == 0)
32         {
33             if(i != t)
34                 printf(" ");
35             printf("%d" , c);
36         }
37         q[1][++r[1]] = (int)((long long)c * u / v) - g * i;
38         q[2][++r[2]] = c - (int)((long long)c * u / v) - g * i;
39     }
40     printf("\n");
41     for(i = 1 ; i <= m + n ; i ++ )
42     {
43         p = getmax();
44         c = q[p][l[p]++] + g * m;
45         if(i % t == 0)
46         {
47             if(i != t)
48                 printf(" ");
49             printf("%d" , c);
50         }
51     }
52     printf("\n");
53     return 0;
54 }