【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[..,..]of longint;
     t,w:array[..,..]of longint;
     a:array[..]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)>>];
  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,]>w[x,] then exit(-maxlongint);
  exit(q[x,t[x,]]);
 end;

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

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