[bzoj1176]Mokia[CDQ分治]

啃了一天论文，发现CDQ分治的原理其实很简单，大概就是这样的一类分治：将左右区间按一定规律排序后分开处理，递归到底时直接计算答案，对于一个区间，按照第二关键字split成两个区间，先处理左区间，之后因为整个区间是有序的，就可以根据左区间来推算右区间的答案，最后递归处理右区间即可。拿此题做比方，先把全区间按照x坐标排序，然后自左向右用前一半（按时间排序）的修改来更新后一半的查询，之后将整个区间按照时间分成两个部分，递归处理。归纳起来就是split->left->push->right->merge->return。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <cmath>

using namespace std;

struct node { int    v,x,y,op,d,pos; }a[210000],temp[210000];

int    s,w,cnt;
int    L[2100000],Ans[11000];

bool    operator<(const node temp1,const node temp2)
{
    if(temp1.x!=temp2.x)return temp1.x<temp2.x;
    if(temp1.y!=temp2.y)return temp1.y<temp2.y;
    return temp1.op<temp2.op;
}

void    Add(const int x,const int d) { for(int i=x;i<=w;i+=i&-i)L[i]+=d; return ; }
int    Query(const int x) { int t=0; for(int i=x;i;i-=i&-i) t+=L[i]; return t; }

void    CDQ(const int l,const int r)
{
    if(l==r)return ;

    int    i,mid=l+((r-l)>>1);

    for(i=l;i<=r;++i)
    {
        if(a[i].v<=mid && a[i].op==0)Add(a[i].y,a[i].d);
        if(a[i].v>mid  && a[i].op==1)Ans[a[i].pos]+=Query(a[i].y)*a[i].d;
    }
    for(i=l;i<=r;++i)
        if(a[i].v<=mid && a[i].op==0)Add(a[i].y,-a[i].d);

    int l1=l,l2=mid+1;
    for(i=l;i<=r;++i)
    {
        if(a[i].v<=mid)temp[l1++]=a[i];
        else    temp[l2++]=a[i];
    }
    for(i=l;i<=r;++i)a[i]=temp[i];
    CDQ(l,mid),CDQ(mid+1,r);
    return ;
}


int main()
{
    int    i,op,x,y,z,xa,xb,ya,yb;

    scanf("%d%d",&s,&w);
    while(1)
    {
        scanf("%d",&op);
        if(op==1)
        {
            scanf("%d%d%d",&x,&y,&z);
            a[++cnt]=(node){cnt,x,y,0,z,0};
        }
        else if(op==2)
        {
            scanf("%d%d%d%d",&xa,&ya,&xb,&yb);
            a[++cnt]=(node){cnt,xb  ,yb  ,1,+1,++Ans[0]};
            a[++cnt]=(node){cnt,xb  ,ya-1,1,-1,Ans[0]};
            a[++cnt]=(node){cnt,xa-1,yb  ,1,-1,Ans[0]};
            a[++cnt]=(node){cnt,xa-1,ya-1,1,+1,Ans[0]};
        }
        else    break;
    }

    sort(a+1,a+cnt+1);

    CDQ(1,cnt);

    for(i=1;i<=Ans[0];++i)
        printf("%d
",Ans[i]);

    return 0;
}