Mokia(三维偏序)P4390

提到cdq,就不得不提这道该死的,挨千刀的题目了。

极简题面：

给定一个二维平面，在ti时刻会在（xi，yi）放一个点，会在tj时刻查询一个方框里面的点的数量

看道题就是二维线段树乱搞啊，这么水？？？

数据范围劝退警告

单是一维都快有点吃不消了...1e6*1e6的数组？几个GB？？？

。。。

于是，伟大的CDQ分治出场了。

题面其实可以这样翻译：

按时插入点，询问小于（x，y）且时间也小于当前点的点的个数

这不就是CDQ的事吗？比模板题还要裸。。。

但是可能要差分一下（二维差分）因为统计的是点与00组成的大矩形，所以要剪去两个矩形，再加上一个小矩形，所以要统计四个点的偏序

总结一下，就是cdq。

第一维时间，第二维x，第三维y

一定要离线做

于是开始了愉快的CDQ

#include<bits/stdc++.h>
using namespace std;
const int maxn=;

struct node
{
    int time,x,y,val,id;
}e[maxn];
int m,cnt,t[maxn<<],a[maxn],ans[maxn];
inline int lowbit(int x)
{
    return x & - x ;
}
void add(int x,int y)
{
    for(;x<=m;x+=lowbit(x))
    {
        t[x]+=y;
    }
}
int ask(int x)
{
    int res=;
    for(;x;x-=lowbit(x))
    {
        res+=t[x];
    }
    return res;
}
bool cmp2(node a,node b)
{
    if(a.x!=b.x)return a.x<b.x;
    if(a.y!=b.y)return a.y<b.y;
    //else        return a.time<b.time;
}
bool cmp(node a,node b)
{
    return a.time<b.time;
}
void cdq(int l,int r)
{
    if(l==r)return;
    int mid=l+r>>;
    cdq(l,mid);
    cdq(mid+,r);
    sort(e+l,e++r,cmp2);
    for(int i=l;i<=r;i++)
    {
        if(e[i].x<=mid&&e[i].id==)
        add(e[i].y,e[i].val);
        else e[i].val+=ask(e[i].y);
    }
    for(int i=l;i<=r;i++)
    {
        if(e[i].x<=mid&&e[i].id==)
        add(e[i].y,-e[i].val);
    }
}
int read()
{
    int f=,x=;char s=getchar();
    while(s>''||s<''){if(s=='-')f=-;s=getchar();}
    while(s<=''&&s>=''){x=x*+s-'';s=getchar();}
    return x*f;
}
int main()
{
    read();
    m=read();
    int flag=read();
    while(flag!=)
    {
        if(flag==)
        {
            int x=read()+,y=read()+,val=read();
            e[++cnt]=(node){cnt,x,y,val,};
        }
        else
        {
            int x1=read(),yl=read(),x2=read()+,y2=read()+;
            e[++cnt]=(node){cnt,x1,yl,,};//数据结构体化
            e[++cnt]=(node){cnt,x2,y2,,};
            e[++cnt]=(node){cnt,x2,yl,,};
            e[++cnt]=(node){cnt,x1,y2,,};
        }
        flag=read();
    }
    cdq(,cnt);然后硬cdq就行了
    sort(e+,e+cnt+,cmp);
    for(int i=;i<=cnt;++i)
    {
        if(e[i].id==)
        {
            printf("%d\n",e[i].val+e[i+].val-e[i+].val-e[i+].val);
            i+=;
        }
    }
    return ;
}

（完）