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=200005;
inline int read()
{
    int p=0;char c=getchar();
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')p=p*10+c-'0',c=getchar();
    return p;
}
struct node
{
    int x1,y1,x2,y2,op,k,num;
    int ans;
}p[maxn],cc[maxn<<1];
int c[500005];
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=0;
void add(int x,int y)
{
    for(int i=x;i<=w;i+=(i&(-i)))
    {
        c[i]+=y;
    }
}
int sum(int x)
{
    int q=0;
    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)>>1;int cnt=0;
    for(int i=l;i<=mid;i++)
    {
        if(p[i].op==1)cc[cnt++]=p[i];
    }
    for(int i=mid+1;i<=r;i++)
    {
        if(p[i].op==2)
        {
            cc[cnt++]=p[i];
            cc[cnt++]=p[i];
            cc[cnt-2].x1--;
            cc[cnt-1].x1=p[i].x2;
            cc[cnt-1].op=3;
        }
    }
    sort(cc,cc+cnt,cmp);
    for(int i=0;i<cnt;i++)
    {
        if(cc[i].op==1)
        {
            add(cc[i].y1,cc[i].k);
        }
        else if(cc[i].op==2)
        {
            p[cc[i].num].ans-=sum(cc[i].y2)-sum(cc[i].y1-1);
        }
        else 
        {
            p[cc[i].num].ans+=sum(cc[i].y2)-sum(cc[i].y1-1);
        }
    }
    for(int i=0;i<cnt;i++)
    {
        if(cc[i].op==1)
        {
            add(cc[i].y1,-cc[i].k);
        }
    }
    solve(l,mid);
    solve(mid+1,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=1;i<=n;i++)
    {
        flag=read();
        if(flag==1)
        {
            a=read();b=read();c=read();
            p[i].op=1;
            p[i].x1=a;
            p[i].y1=b;
            p[i].k=c;
        }
        else 
        {
            a=read();b=read();c=read();d=read();
            p[i].op=2;
            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(1,n);
    for(int i=1;i<=n;i++)
    {
        if(p[i].op==2)printf("%d
",p[i].ans);
    }
    return 0;
}