BZOJ 4721 [Noip2016]蚯蚓 单调队列

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

Sample Output

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

HINT

传送门

堆可以拿70，80左右的分数……

满分其实不是很复杂，

我们假设有两条蚯蚓长度为x和y，且x>y。

如果说x先分裂，y第二个分裂，

那么x分裂后是px，(1-p)x，y变成y+q

y分裂后是p(y+q)，(1-p)(y+q)，同时px变为px+q，(1-p)x变为(1-p)x+q

我们发现，既然0<=p<1， 那么：

px+q>py+pq

(1-p)x+q>(1-p)y+q(1-p)

也就是说假如初始的时候x>y，那么x分裂的蚯蚓比y分裂的蚯蚓还是要长……

假如x>y，那么x+q>y+q，所以x和y都不处理的时候x还是比y大。

所以，我们建立3个队列分别存储原来的，分成px的，分成(1-p)x的蚯蚓，

那么分裂后存放入后两个队列，

我们可以知道一定是满足单调递减的了。

每次从3个队列的头取max的即可。

那么+q怎么处理呢？我们对每个蚯蚓加入时间点Time[u]，

原来的蚯蚓都是Time[u]=0，

然后第i个时间点分出来两个新蚯蚓x,y，则Time[x]=Time[y]=i。

这样子的话，具体计算可以是now*q-Time[u]*q。now是当前时间点。

此题有坑……long long的问题

计算分出来的长度的时候，我们都这样写：len*x/y。

但是！len*x是完全可以爆int的

所以这个地方需要long long。

但是假如所有数组等等都开long long，是会TLE的！！

……坑啊= =

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int read(){
    int f=1,x=0;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int 
    N=100005,
    M=7000005;
int n,m,Q,u,v,T,num;
int h[3],t[3];
int q[3][M];
int LEN[N+M],Time[N+M];
int ANS1[M],ANS2[N+M];
bool cmp(int x,int y){return LEN[x]>LEN[y];}
int count(int x,int i){return LEN[x]+(i-1-Time[x])*Q;}

int MAX(int ttime){
    int tt=-1,tt1=-1,tmp;
    for (int i=0;i<3;i++)
        if (h[i] && h[i]<=t[i]){
             tmp=count(q[i][h[i]],ttime);
             if (tt1<tmp) tt1=tmp,tt=i;
        }
    return tt;
}
void solve2(){
    int len=0,Ttime=0,tmp;
    while (Ttime++<=n+m){
        tmp=MAX(m+1);
        if (tmp<0) break;
        if (!(Ttime%T)) ANS2[++len]=count(q[tmp][h[tmp]],m+1);
        h[tmp]++;
    }
    for (int i=1;i<len;i++) printf("%d ",ANS2[i]);
    printf("%d",ANS2[len]);
}
int main(){
    n=read(),m=read(),Q=read(),u=read(),v=read(),T=read();

    t[0]=n,t[1]=t[2]=0;
    h[0]=1,h[1]=h[2]=1;
    for (int i=1;i<=n;i++) LEN[i]=read(),q[0][i]=i;
    sort(q[0]+1,q[0]+1+n,cmp);

    int x,x1,x2,tmp,len=0;num=n;
    for (int i=1;i<=m;i++){
        x=MAX(i);
        tmp=count(q[x][h[x]],i);
        if (!(i%T)) ANS1[++len]=tmp;

        x1=(ll)tmp*u/v,x2=tmp-x1;

        LEN[q[x][h[x]]]=x1,Time[q[x][h[x]]]=i;
        q[1][++t[1]]=q[x][h[x]],q[2][++t[2]]=++num;
        LEN[num]=x2,Time[num]=i;

        h[x]++;
    }
    if (T<=m){
        for (int i=1;i<len;i++) printf("%d ",ANS1[i]);
        printf("%d",ANS1[len]);
    }
    putchar('\n');
    solve2();
    return 0;
}