noip 2016 D2T2 蚯蚓（队列）

题目描述

现在共有n只拥有长度[1~N]的蚯蚓，每一秒你会将最长的之一按固定比例p切成两条蚯蚓（向下取整），

同时其余蚯蚓的长度会增加q。请问m秒内每t秒被切断的蚯蚓被切断前的长度，以及m秒后，从大到小排

序，排名为t的整数倍的蚯蚓的长度。

思路

很容易可以想到将所有蚯蚓放入堆里，再模拟的方法，但时间复杂度不优。

可以观察到，如果没有增加长度这一条件，对于切断的蚯蚓，我们按切断的比例分成两对，在同一堆的

蚯蚓长度是单减的，所以每次选取就可以从两对的队首和还没切的蚯蚓里面选一个切。如果要增加长度

的话，也满足这个性质，因为他被切断的同时，别人长度都增加了。

关于如何维护长度，记录一个到目前需要增加的长度delt，再队列里面的长度为没增加的长度，每次将

一条蚯蚓取出，它现在的长度为x+delt。对于要切断的蚯蚓，放入时要减（delt+q），以达到别人会加，

他不加的目的。

细节

在模拟之前要先排一次序

在a[n+1]设置极小值，防止越界且方便分类

#include<bits/stdc++.h>
using namespace std;

const int maxn=100005;
const int maxm=7000005;
const int oo=-2100000000;
int n,m,q,u,v,t;
int ans[maxm+maxn];

template<class T>inline void read(T &x){
    x=0;char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) {x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
}

int a[maxn],k=1;
int ct1[maxm],ct2[maxm],h1=1,h2=1,t1,t2;

bool cmp(int x,int y){return x>y;}

int main(){
    read(n);read(m);read(q);read(u);read(v);read(t);
    double p=(double)u/v;
    for(int i=1;i<=n;i++) read(a[i]);
    sort(a+1,a+n+1,cmp);
    a[n+1]=oo;
    int delt=0,x;
    for(int i=1;i<=m;i++){
        if(a[k]>=ct1[h1]&&a[k]>=ct2[h2]) x=a[k++];
        else if(ct1[h1]>=a[k]&&ct1[h1]>=ct2[h2]) x=ct1[h1++];
        else if(ct2[h2]>=a[k]&&ct2[h2]>=ct1[h1]) x=ct2[h2++];
        x+=delt;
        int x1=floor((double)x*p),x2=x-x1;
        delt+=q;
        ct1[++t1]=x1-delt;
        ct2[++t2]=x2-delt;
        if(i%t==0) printf("%d ",x);
    }
    putchar(10);
    for(int i=k;i<=n;i++) ans[++ans[0]]=a[i];
    for(int i=h1;i<=t1;i++) ans[++ans[0]]=ct1[i];
    for(int i=h2;i<=t2;i++) ans[++ans[0]]=ct2[i];
    sort(ans+1,ans+ans[0]+1,cmp);
    for(int i=t;i<=ans[0];i+=t) printf("%d ",ans[i]+delt);
}

 然后你就会发现bzoj过不了会T，就是在取出各个队列的蚯蚓时会放入n+m个，排序就很难受。

其实队列本来就是有序的，所以只要照着之前的操作取一遍就好了，但是不好维护取的是哪个，所以三个队列写在一起。（好像也照着之前的操作打可以，不过这样也清晰一点吧）。

#include<bits/stdc++.h>
using namespace std;

const int maxn=100005;
const int maxm=7000005;
const int oo=-2100000000;
int n,m,q,u,v,t;
int ans[maxm+maxn];

template<class T>inline void read(T &x){
    x=0;char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) {x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
}

int Q[3][maxm],head[3],tail[3];

bool cmp(int x,int y){return x>y;}

int main(){
    read(n);read(m);read(q);read(u);read(v);read(t);
    double p=(double)u/v;
    for(int i=1;i<=n;i++) read(Q[0][i]);
    head[0]=head[1]=head[2]=1;
    tail[1]=tail[2]=0;tail[0]=n;
    sort(Q[0]+1,Q[0]+n+1,cmp);
    int delt=0,x,k=0;
    for(int i=1;i<=m;i++){
        x=oo;
        for(int j=0;j<3;j++)
         if(head[j]<=tail[j]&&x<Q[j][head[j]])
          x=Q[j][head[j]],k=j;
        head[k]++;
        x+=delt;
        int x1=floor((double)x*p),x2=x-x1;
        delt+=q;
        Q[1][++tail[1]]=x1-delt;
        Q[2][++tail[2]]=x2-delt;
        if(i%t==0) printf("%d ",x);
    }
    putchar(10);
    for(int i=1;i<=n+m;i++){
        x=oo;
        for(int j=0;j<3;j++)
         if(head[j]<=tail[j]&&x<Q[j][head[j]])
          x=Q[j][head[j]],k=j;
        head[k]++;
        if(!(i%t)) printf("%d ",x+delt);
    }
}