【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[..]of extended;
     t:array[..]of longint;
     c,q:array[..]of record
                              l,r,op,t:longint;
                              a:int64;
                              s:extended;
                             end;
     b:array[..]of extended;
     head,vet,next,ans,flag:array[..]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]:=; s2[i]:=; 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:=; i:=x;
  while i> do
  begin
   if t[i]=time then query:=query+s1[i]*(x+)-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= then ans[c[i].a]:=l;
   exit;
  end;
  mid:=(l+r)>>;
  pre:=x; now:=x; inc(time);
  for i:=x to y do
   if c[i].op= 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(,c[i].a);
      add(c[i].r+,-c[i].a);
     end
      else
      begin
       add(c[i].l,c[i].a);
       add(c[i].r+,-c[i].a);
      end;
     flag[i]:=;
     inc(now);
    end
     else flag[i]:=;
   end;
  for i:=x to y do
   if c[i].op= then
   begin
    u:=c[i].a; tmp:=;
    e:=head[u];
    while e<> do
    begin
     v:=vet[e];
     tmp:=tmp+query(v)-query(v-);
     if tmp>=c[i].s then break;
     e:=next[e];
    end;
    if tmp>=c[i].s then
    begin
     inc(now); flag[i]:=;
    end
     else
     begin
      flag[i]:=; c[i].s:=c[i].s-tmp;
     end;
   end;

  for i:=x to y do
   if flag[i]= 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-,l,mid);
  solve(pre,y,mid+,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:= to m do
  begin
   read(x);
   add(x,i);
  end;
  for i:= to n do read(b[i]);
  read(k);
  for i:= to k do
  begin
   read(c[i].l,c[i].r,c[i].a);
   c[i].op:=; c[i].t:=i;
  end;
  for i:= to n do
  begin
   c[i+k].op:=;
   c[i+k].a:=i; c[i+k].s:=b[i];
  end;

  solve(,k+n,,k+);
  for i:= to n do
   if ans[i]<=k then writeln(ans[i])
    else writeln('NIE');

 end.