[NOIP2016]蚯蚓 D2 T2 队列

[NOIP2016]蚯蚓 D2 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

Sample Output

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

题解：开三个队列，p1装的是原来的n个蚯蚓（排好序），p2，p3分别表示每一次切割后分成u/v和(v-u)/v的段，每次取3个队列中队首最大的元素进行切割，容易发现在每次切割后p2,p3始终具有单调性。每次对所有元素加P，只需要将now+=p，计算时在统一加上now就可以了。

考场上用的STL的queue，结果只拿了60多分，被人家用手写堆和读入优化80多分完虐了~

所以必须用手写队列+读入优化，否则还是会TLE

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m,q,u,v,t,now;
int a[100010];
struct queue
{
    int l,r,v[8000010];
    void push(int x)
    {
        v[++r]=x;
    }
    int front()
    {
        return v[l];
    }
    void pop()
    {
        l++;
    }
    int empty()
    {
        return l>r?1:0;
    }
}p[3];
int getbig()
{
    for(int i=0;i<=2;i++)
    {
        if(p[i].empty()) continue;
        if(i!=0&&!p[0].empty()&&p[0].front()>p[i].front()) continue;
        if(i!=1&&!p[1].empty()&&p[1].front()>p[i].front()) continue;
        if(i!=2&&!p[2].empty()&&p[2].front()>p[i].front()) continue;
        int x=p[i].front(),p[i].pop();
        return x;
    }
}
int readin()
{
    int ret=0;    char gc;
    while(gc<'0'||gc>'9')    gc=getchar();
    while(gc>='0'&&gc<='9')    ret=ret*10+gc-'0',gc=getchar();
    return ret;
}
int main()
{
    int i,x;
    n=readin(),m=readin(),q=readin(),u=readin(),v=readin(),t=readin();
    for(i=1;i<=n;i++)    a[i]=readin(); 
    sort(a+1,a+n+1);
    for(i=n;i>=1;i--) p[0].push(a[i]);
    p[0].l=p[1].l=p[2].l=1;
    for(i=1;i<=m;i++)
    {
        x=getbig();
        if(i%t==0)
        {
            if(i!=t)
            printf(" ");
            printf("%d",x+now);
        }
        p[1].push((long long)(x+now)*u/v-now-q);
        p[2].push((x+now)-(long long)(x+now)*u/v-now-q);
        now+=q;
    }
    printf("
");
    i=0;
    while(!p[0].empty()||!p[1].empty()||!p[2].empty())
    {
        i++;
        x=getbig();
        if(i%t==0)
        {
            if(i!=t) printf(" ");
            printf("%d",x+now);
        }
    }
    return 0;
}