【BZOJ2527】Meteors（整体二分）

题意：

Byteotian Interstellar Union有N个成员国。现在它发现了一颗新的星球，这颗星球的轨道被分为M份（第M份和第1份相邻），第i份上有第Ai个国家的太空站。

这个星球经常会下陨石雨。BIU已经预测了接下来K场陨石雨的情况。
BIU的第i个成员国希望能够收集Pi单位的陨石样本。你的任务是判断对于每个国家，它需要在第几次陨石雨之后，才能收集足够的陨石。
输入：
第一行是两个数N,M。
第二行有M个数，第i个数Oi表示第i段轨道上有第Oi个国家的太空站。
第三行有N个数，第i个数Pi表示第i个国家希望收集的陨石数量。
第四行有一个数K，表示BIU预测了接下来的K场陨石雨。
接下来K行，每行有三个数Li,Ri,Ai，表示第K场陨石雨的发生地点在从Li顺时针到Ri的区间中（如果Li<=Ri，就是Li,Li+1,...,Ri，否则就是Ri,Ri+1,...,m-1,m,1,...,Li），向区间中的每个太空站提供Ai单位的陨石样本。
输出：
N行。第i行的数Wi表示第i个国家在第Wi波陨石雨之后能够收集到足够的陨石样本。如果到第K波结束后仍然收集不到，输出NIE。

数据范围： 1<=n,m,k<=3*10^5 1<=Pi<=10^9 1<=Ai<10^9

 

思路：RYZ作业

整体二分很显然，卡了5-6天，线段树标记持久化都过不去，树状数组暴力清空也T，实在是没办法了……

var s1,s2:array[1..400000]of extended;
    t:array[1..400000]of longint;
    c,q:array[1..710000]of record
                             l,r,op,t:longint;
                             a:int64;
                             s:extended;
                            end;
    b:array[1..710000]of extended;
    head,vet,next,ans,flag:array[1..710000]of longint;
    n,m,i,k,tot,x,time:longint;
 
procedure add(x:longint;y:extended);
var i:longint;
begin
 i:=x;
 while i<=m do
 begin
  if t[i]<>time then
  begin
   s1[i]:=0; s2[i]:=0; t[i]:=time;
  end;
  s1[i]:=s1[i]+y;
  s2[i]:=s2[i]+y*x;
  i:=i+i and (-i);
 end;
end;
 
function query(x:longint):extended;
var i:longint;
begin
 query:=0; i:=x;
 while i>0 do
 begin
  if t[i]=time then query:=query+s1[i]*(x+1)-s2[i];
  i:=i-i and (-i);
 end;
end;
 
procedure solve(x,y,l,r:longint);
var i,j,mid,pre,now,u,e,v:longint;
    tmp:extended;
begin
 if l>r then exit;
 if l=r then
 begin
  for i:=x to y do
   if c[i].op=2 then ans[c[i].a]:=l;
  exit;
 end;
 mid:=(l+r)>>1;
 pre:=x; now:=x; inc(time);
 for i:=x to y do
  if c[i].op=1 then
  begin
   if c[i].t<=mid then
   begin
    if c[i].l>c[i].r then
    begin
     add(c[i].l,c[i].a);
     add(1,c[i].a);
     add(c[i].r+1,-c[i].a);
    end
     else
     begin
      add(c[i].l,c[i].a);
      add(c[i].r+1,-c[i].a);
     end;
    flag[i]:=0;
    inc(now);
   end
    else flag[i]:=1;
  end;
 for i:=x to y do
  if c[i].op=2 then
  begin
   u:=c[i].a; tmp:=0;
   e:=head[u];
   while e<>0 do
   begin
    v:=vet[e];
    tmp:=tmp+query(v)-query(v-1);
    if tmp>=c[i].s then break;
    e:=next[e];
   end;
   if tmp>=c[i].s then
   begin
    inc(now); flag[i]:=0;
   end
    else
    begin
     flag[i]:=1; c[i].s:=c[i].s-tmp;
    end;
  end;
  
 for i:=x to y do
  if flag[i]=1 then
  begin
   q[now]:=c[i]; inc(now);
  end
   else
   begin
    q[pre]:=c[i]; inc(pre);
   end;
 for i:=x to y do c[i]:=q[i];
 solve(x,pre-1,l,mid);
 solve(pre,y,mid+1,r);
end;
 
procedure add(a,b:longint);
begin
 inc(tot);
 next[tot]:=head[a];
 vet[tot]:=b;
 head[a]:=tot;
end;
 
begin
 
 read(n,m);
 for i:=1 to m do
 begin
  read(x);
  add(x,i);
 end;
 for i:=1 to n do read(b[i]);
 read(k);
 for i:=1 to k do
 begin
  read(c[i].l,c[i].r,c[i].a);
  c[i].op:=1; c[i].t:=i;
 end;
 for i:=1 to n do
 begin
  c[i+k].op:=2;
  c[i+k].a:=i; c[i+k].s:=b[i];
 end;
 
 solve(1,k+n,1,k+1);
 for i:=1 to n do
  if ans[i]<=k then writeln(ans[i])
   else writeln('NIE');
 
end.