【NOIP2016】蚯蚓（单调队列）

题意：

思路：

我们发现，
对于任意两次切割i和j，i<j，
在进行完第j次切割后，第i次切割的u/v部分一定大于等于第j次切割的u/v部分，
第i次的1-u/v部分也一定大于等于第j次的1-u/v部分
证明很显然，
假设在第i次切割的时候，切割的蚯蚓长度为x，
第j次的时候为y+(i-j)q，
那么有x>=y，px+(i-j)q显然大于等于p(y+(i-j)q)
那么我们维护三个队列，
刚开始的时候把所有蚯蚓排序后推到第一个队列里，
然后每次取三个队头里最长的，
把他弹出来，再把切开的两半分别推入第二和第三个队列即可

From http://blog.csdn.net/neither_nor/article/details/53300255

线性复杂度还是被卡一个点 也是醉了

var q:array[1..3,1..8000000]of longint;
    t,w:array[1..3,0..0]of longint;
    a:array[1..200000]of longint;
    n,m,q1,u,v,t1,i:longint;
    now,s1,s2:longint;

procedure swap(var x,y:longint);
var t:longint;
begin
 t:=x; x:=y; y:=t;
end;

procedure qsort(l,r:longint);
var i,j,mid:longint;
begin
 i:=l; j:=r; mid:=a[(l+r)>>1];
 repeat
  while mid<a[i] do inc(i);
  while mid>a[j] do dec(j);
  if i<=j then
  begin
   swap(a[i],a[j]);
   inc(i); dec(j);
  end;
 until i>j;
 if l<j then qsort(l,j);
 if i<r then qsort(i,r);
end;

function h(x:longint):longint;
begin
 if t[x,0]>w[x,0] then exit(-maxlongint);
 exit(q[x,t[x,0]]);
end;

begin
 assign(input,'bzoj4721.in'); reset(input);
 assign(output,'bzoj4721.out'); rewrite(output);
 readln(n,m,q1,u,v,t1);
 for i:=1 to n do read(a[i]);
 qsort(1,n);
 for i:=1 to n do q[1,i]:=a[i];
 t[1,0]:=1; t[2,0]:=1; t[3,0]:=1;
 w[1,0]:=n; w[2,0]:=0; w[3,0]:=0;
 for i:=1 to m do
 begin
  if (h(1)<>-maxlongint)and(h(1)>=h(2))and(h(1)>=h(3)) then
  begin
   now:=h(1); inc(t[1,0]);
  end
   else if (h(2)<>-maxlongint)and(h(2)>=h(1))and(h(2)>=h(3)) then
   begin
    now:=h(2); inc(t[2,0]);
   end
    else
    begin
     now:=h(3); inc(t[3,0]);
    end;
  now:=now+(i-1)*q1;
  if i mod t1=0 then
  begin
   write(now);
   if i+t1<=m then write(' ');
  end;
  s1:=int64(now)*u div v;
  s2:=int64(now)-s1;
  inc(w[2,0]); q[2,w[2,0]]:=s1-i*q1;
  inc(w[3,0]); q[3,w[3,0]]:=s2-i*q1;
 end;
 writeln;
 for i:=1 to n+m do
 begin
  if (h(1)<>-maxlongint)and(h(1)>=h(2))and(h(1)>=h(3)) then
  begin
   now:=h(1); inc(t[1,0]);
  end
   else if (h(2)<>-maxlongint)and(h(2)>=h(1))and(h(2)>=h(3)) then
   begin
    now:=h(2); inc(t[2,0]);
   end
    else
    begin
     now:=h(3); inc(t[3,0]);
    end;
  now:=now+m*q1;
  if i mod t1=0 then
  begin
   write(now);
   if i+t1<=n+m then write(' ');
  end;
 end;



 close(input);
 close(output);
end.