cogs 577 蝗灾 CDQ分治

第一道CDQ，抄了下helenkeller的代码，感觉和归并排序差不多。。。

因为左半边的修改肯定在右半边的询问之前，所以就不用管时间的限制了，可以直接x轴排序树状数组处理y轴。。。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<cmath>
 using namespace std;
 int w,n;
 const int maxn=;
 inline int read()
 {
     int p=;char c=getchar();
     while(c<''||c>'')c=getchar();
     while(c>=''&&c<='')p=p*+c-'',c=getchar();
     return p;
 }
 struct node
 {
     int x1,y1,x2,y2,op,k,num;
     int ans;
 }p[maxn],cc[maxn<<];
 int c[];
 bool cmp(const node&a,const node&b)
 {
     if(a.x1==b.x1)return a.op<b.op;
     return a.x1<b.x1;
 }
 int tot=;
 void add(int x,int y)
 {
     for(int i=x;i<=w;i+=(i&(-i)))
     {
         c[i]+=y;
     }
 }
 int sum(int x)
 {
     int q=;
     for(int i=x;i;i-=(i&(-i)))
     {
         q+=c[i];
     }
     return q;
 }
 void solve(int l,int r)
 {
     if(l==r)return ;
     int mid=(l+r)>>;int cnt=;
     for(int i=l;i<=mid;i++)
     {
         if(p[i].op==)cc[cnt++]=p[i];
     }
     for(int i=mid+;i<=r;i++)
     {
         if(p[i].op==)
         {
             cc[cnt++]=p[i];
             cc[cnt++]=p[i];
             cc[cnt-].x1--;
             cc[cnt-].x1=p[i].x2;
             cc[cnt-].op=;
         }
     }
     sort(cc,cc+cnt,cmp);
     for(int i=;i<cnt;i++)
     {
         if(cc[i].op==)
         {
             add(cc[i].y1,cc[i].k);
         }
         else if(cc[i].op==)
         {
             p[cc[i].num].ans-=sum(cc[i].y2)-sum(cc[i].y1-);
         }
         else
         {
             p[cc[i].num].ans+=sum(cc[i].y2)-sum(cc[i].y1-);
         }
     }
     for(int i=;i<cnt;i++)
     {
         if(cc[i].op==)
         {
             add(cc[i].y1,-cc[i].k);
         }
     }
     solve(l,mid);
     solve(mid+,r);
     return ;
 }
 int main()
 {
     freopen("locust.in","r",stdin);
     freopen("locust.out","w",stdout);
     scanf("%d%d",&w,&n);
     int flag,a,b,c,d;
     for(int i=;i<=n;i++)
     {
         flag=read();
         if(flag==)
         {
             a=read();b=read();c=read();
             p[i].op=;
             p[i].x1=a;
             p[i].y1=b;
             p[i].k=c;
         }
         else
         {
             a=read();b=read();c=read();d=read();
             p[i].op=;
             p[i].x1=min(a,c);p[i].x2=max(a,c);
             p[i].y1=min(b,d);p[i].y2=max(b,d);
         }
         p[i].num=i;
     }
     solve(,n);
     for(int i=;i<=n;i++)
     {
         if(p[i].op==)printf("%d\n",p[i].ans);
     }
     return ;
 }