Luogu P2827 蚯蚓

https://www.luogu.org/problemnew/show/P2827

题目描述

本题中，我们将用符号 \lfloor c \rfloor⌊c⌋ 表示对 cc 向下取整，例如：\lfloor 3.0 \rfloor = \lfloor 3.1 \rfloor = \lfloor 3.9 \rfloor = 3⌊3.0⌋=⌊3.1⌋=⌊3.9⌋=3。

蛐蛐国最近蚯蚓成灾了！隔壁跳蚤国的跳蚤也拿蚯蚓们没办法，蛐蛐国王只好去请神刀手来帮他们消灭蚯蚓。

蛐蛐国里现在共有 nn 只蚯蚓（nn 为正整数）。每只蚯蚓拥有长度，我们设第 ii 只蚯蚓的长度为 a_iai​ (i=1,2,\dots,ni=1,2,…,n)，并保证所有的长度都是非负整数（即：可能存在长度为 00 的蚯蚓）。

每一秒，神刀手会在所有的蚯蚓中，准确地找到最长的那一只（如有多个则任选一个）将其切成两半。神刀手切开蚯蚓的位置由常数 pp（是满足 0 < p < 10<p<1 的有理数）决定，设这只蚯蚓长度为 xx，神刀手会将其切成两只长度分别为 \lfloor px \rfloor⌊px⌋ 和 x - \lfloor px \rfloorx−⌊px⌋ 的蚯蚓。特殊地，如果这两个数的其中一个等于 00，则这个长度为 00 的蚯蚓也会被保留。此外，除了刚刚产生的两只新蚯蚓，其余蚯蚓的长度都会增加 qq（是一个非负整常数）。

蛐蛐国王知道这样不是长久之计，因为蚯蚓不仅会越来越多，还会越来越长。蛐蛐国王决定求助于一位有着洪荒之力的神秘人物，但是救兵还需要 mm 秒才能到来……（mm 为非负整数）

蛐蛐国王希望知道这 mm 秒内的战况。具体来说，他希望知道：

• mm 秒内，每一秒被切断的蚯蚓被切断前的长度（有 mm 个数）；
• mm 秒后，所有蚯蚓的长度（有 n + mn+m 个数）。

蛐蛐国王当然知道怎么做啦！但是他想考考你……

输入输出格式

输入格式：

 

第一行包含六个整数 n,m,q,u,v,tn,m,q,u,v,t，其中：n,m,qn,m,q 的意义见【问题描述】；u,v,tu,v,t 均为正整数；你需要自己计算 p=u / vp=u/v（保证 0 < u < v0<u<v）；tt 是输出参数，其含义将会在【输出格式】中解释。

第二行包含 nn 个非负整数，为 a_1, a_2, \dots, a_na1​,a2​,…,an​，即初始时 nn 只蚯蚓的长度。

同一行中相邻的两个数之间，恰好用一个空格隔开。

保证 1 \leq n \leq 10^51≤n≤105，0 \leq m \leq 7 \times 10^60≤m≤7×106，0 < u < v \leq 10^90<u<v≤109，0 \leq q \leq 2000≤q≤200，1 \leq t \leq 711≤t≤71，0 \leq a_i \leq 10^80≤ai​≤108。

 

输出格式：

 

第一行输出 \left \lfloor \frac{m}{t} \right \rfloor⌊tm​⌋ 个整数，按时间顺序，依次输出第 tt 秒，第 2t2t 秒，第 3t3t 秒，……被切断蚯蚓（在被切断前）的长度。

第二行输出 \left \lfloor \frac{n+m}{t} \right \rfloor⌊tn+m​⌋ 个整数，输出 mm 秒后蚯蚓的长度；需要按从大到小的顺序，依次输出排名第 tt，第 2t2t，第 3t3t，……的长度。

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

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

 

输入输出样例

输入样例#1： 复制

3 7 1 1 3 1
3 3 2

输出样例#1： 复制

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

输入样例#2： 复制

3 7 1 1 3 2
3 3 2

输出样例#2： 复制

4 4 5
6 5 4 3 2

输入样例#3： 复制

3 7 1 1 3 9
3 3 2

输出样例#3： 复制

//空行
2

说明

【样例解释1】

在神刀手到来前：33只蚯蚓的长度为3,3,23,3,2。

11秒后：一只长度为33的蚯蚓被切成了两只长度分别为11和22的蚯蚓，其余蚯蚓的长度增加了11。最终44只蚯蚓的长度分别为(1,2),4,3(1,2),4,3。括号表示这个位置刚刚有一只蚯蚓被切断

22秒后：一只长度为44的蚯蚓被切成了11和33。55只蚯蚓的长度分别为：2,3,(1,3),42,3,(1,3),4。

3秒后：一只长度为44的蚯蚓被切断。66只蚯蚓的长度分别为：3,4,2,4,(1,3)3,4,2,4,(1,3)。

44秒后：一只长度为44的蚯蚓被切断。77只蚯蚓的长度分别为：4,(1,3),3,5,2,44,(1,3),3,5,2,4。

55秒后：一只长度为55的蚯蚓被切断。88只蚯蚓的长度分别为：5,2,4,4,(1,4),3,55,2,4,4,(1,4),3,5。

66秒后：一只长度为55的蚯蚓被切断。99只蚯蚓的长度分别为：(1,4),3,5,5,2,5,4,6(1,4),3,5,5,2,5,4,6。

77秒后：一只长度为66的蚯蚓被切断。1010只蚯蚓的长度分别为：2,5,4,6,6,3,6,5,(2,4)2,5,4,6,6,3,6,5,(2,4)。所以，77秒内被切断的蚯蚓的长度依次为3,4,4,4,5,5,63,4,4,4,5,5,6。77秒后，所有蚯蚓长度从大到小排序为6,6,6,5,5,4,4,3,2,26,6,6,5,5,4,4,3,2,2

【样例解释2】

这个数据中只有t=2t=2与上个数据不同。只需在每行都改为每两个数输出一个数即可。

虽然第一行最后有一个66没有被输出，但是第二行仍然要重新从第二个数再开始输出。

【样例解释3】

这个数据中只有t=9t=9与上个数据不同。

注意第一行没有数要输出，但也要输出一个空行。

【数据范围】

 

这题和合并果子差不多，可以先把那题做一做

https://www.luogu.org/problemnew/show/P1090

堆有80分，贼好写

开三个单调对面，1个是未经处理国，另两个是存切成的两个部分

可以证明原来比他大的切成两部分后的剩余部分仍然大小不变

记得开longlong

#include<cstdio>
#include<algorithm>
using namespace std;
int read()
{
    int ret=0; char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') 
        ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
    return ret; 
}
 
const int N=1e7+5;
int n,m,l,u,v,t;
long long x;
int l1,r1,l2,r2,l3,r3,a1[N],a2[N],a3[N];
 
bool cmp(int i,int j)
{
    return i>j;
}
  
int main()
{
    //freopen("earthworm.in","r",stdin);
    //freopen("earthworm.out","w",stdout);
    n=read(),m=read(),l=read(),
    u=read(),v=read(),t=read(); 
    for(int i=1;i<=n;i++) a1[i]=read();
    sort(a1+1,a1+n+1,cmp); 
    l1=l2=l3=1,r1=n,r2=r3=0;
    for(int i=1;i<=m;i++)
    {
        if(l1<=r1&&l2<=r2&&l3<=r3)
        {
            if(a1[l1]>=a2[l2]&&a1[l1]>=a3[l3]) x=a1[l1],l1++;
            else
            if(a2[l2]>=a1[l1]&&a2[l2]>=a3[l3]) x=a2[l2],l2++;
            else
            if(a3[l3]>=a1[l1]&&a3[l3]>=a2[l2]) x=a3[l3],l3++;
        }else
        if(l1<=r1&&l2<=r2)
        {
            if(a1[l1]>=a2[l2]) x=a1[l1],l1++;
                else x=a2[l2],l2++; 
        }else
        if(l1<=r1&&l3<=r3)
        {
            if(a1[l1]>=a3[l3]) x=a1[l1],l1++;
                else x=a3[l3],l3++;
        }else
        if(l2<=r2&&l3<=r3)
        {
            if(a2[l2]>=a3[l3]) x=a2[l2],l2++;
                else x=a3[l3],l3++;
        } else
        if(l1<=r1) x=a1[l1],l1++;
        else
        if(l2<=r2) x=a2[l2],l2++;
        else
        if(l3<=r3) x=a3[l3],l3++;
        x+=l*(i-1);
        a2[++r2]=x*u/v-l*i;
        a3[++r3]=x-x*u/v-l*i;
        if(i%t==0) printf("%lld ",x);
    }
    puts("");
    for(int i=1;i<=n+m;i++)
    {
        if(l1<=r1&&l2<=r2&&l3<=r3)
        {
            if(a1[l1]>=a2[l2]&&a1[l1]>=a3[l3]) x=a1[l1],l1++;
            else
            if(a2[l2]>=a1[l1]&&a2[l2]>=a3[l3]) x=a2[l2],l2++;
            else
            if(a3[l3]>=a1[l1]&&a3[l3]>=a2[l2]) x=a3[l3],l3++;
        }else
        if(l1<=r1&&l2<=r2)
        {
            if(a1[l1]>=a2[l2]) x=a1[l1],l1++;
                else x=a2[l2],l2++; 
        }else
        if(l1<=r1&&l3<=r3)
        {
            if(a1[l1]>=a3[l3]) x=a1[l1],l1++;
                else x=a3[l3],l3++;
        }else
        if(l2<=r2&&l3<=r3)
        {
            if(a2[l2]>=a3[l3]) x=a2[l2],l2++;
                else x=a3[l3],l3++;
        } else
        if(l1<=r1) x=a1[l1],l1++;
        else
        if(l2<=r2) x=a2[l2],l2++;
        else
        if(l3<=r3) x=a3[l3],l3++;
        x+=l*m;
        if(i%t==0) printf("%lld ",x);
    }
    //fclose(stdin); fclose(stdout);
    return 0;
}