CDQ分治

要求可以计算前面的操作对后面询问的贡献

BZOJ1176

#include <cstdio>
#include <algorithm>
using namespace std;

  int tr[],p,ans[],s,q,opt,n,que;

  struct data1{
      int po,x,y,num,v,opt;
  }a[],tmp[];

  struct data{
      int num,po;
  }yt[];

  int mycomp(const data&a,const data&b){
      return(a.num<b.num);
  }

  int lowbit(int po){
      return(po&(-po));
  }

  int query(int po){
      int ret=;
      for (;po;po-=lowbit(po)) ret+=tr[po];
      return(ret);
  }

  void add(int po,int num){
      for (;po<=p;po+=lowbit(po)) tr[po]+=num;
  }

  void solve(int l,int r){
      if (l==r) return;
      int mid=(l+r)>>;
      solve(l,mid);solve(mid+,r);
      int po=l;
      for (int i=mid+;i<=r;i++)
        if (a[i].opt==){
          while (po<=mid&&a[po].x<=a[i].x){
            if (a[po].opt==) add(a[po].y,a[po].num);
          po++;
        }
        ans[a[i].po]+=query(a[i].y)*a[i].v;
      }
    for (int i=l;i<po;i++)
      if (a[i].opt==) add(a[i].y,-a[i].num);

    int po1=l,po2=mid+,po3=l;
    while (po1<=mid&&po2<=r){
      if (a[po1].x<a[po2].x) tmp[po3++]=a[po1++];else tmp[po3++]=a[po2++];
    }
    while (po1<=mid) tmp[po3++]=a[po1++];
    while (po2<=r) tmp[po3++]=a[po2++];
    for (int i=l;i<=r;i++) a[i]=tmp[i];
  }

  int main(){
      scanf("%d%d",&s,&q);
      while (scanf("%d",&opt),opt!=){
      if (opt==){
          n++;
          a[n].opt=;
          scanf("%d%d%d",&a[n].x,&yt[n].num,&a[n].num);
          yt[n].po=n;
      }else{
          int t1,t2,t3,t4;
        scanf("%d%d%d%d",&t1,&t2,&t3,&t4);
        que++;
        a[++n].po=que;a[n].v=;a[n].x=t3;yt[n].num=t4;yt[n].po=n;a[n].opt=;
        a[++n].po=que;a[n].v=-;a[n].x=t1-;yt[n].num=t4;yt[n].po=n;a[n].opt=;
        a[++n].po=que;a[n].v=-;a[n].x=t3;yt[n].num=t2-;yt[n].po=n;a[n].opt=;
        a[++n].po=que;a[n].v=;a[n].x=t1-;yt[n].num=t2-;yt[n].po=n;a[n].opt=;
      }
    }

    sort(yt+,yt+n+,mycomp);
    yt[].num=-1e9;
    p=;
    for (int i=;i<=n;i++){
      if (yt[i].num!=yt[i-].num) p++;
      a[yt[i].po].y=p;
    }

    solve(,n);

    for (int i=;i<=que;i++) printf("%d\n",ans[i]);
  }