bzoj3262陌上花开 cdq分治 三维偏序

cdq分治的入门题目，学习一下这个很实用的玩意儿。
大概思想就是把在线变为离线，去掉时间的限制，然后每次对于询问区间分治处理，块内之间的影响可以递归处理，然后剩下的就是块与块之间的影响，这个其实就很简单了，随便搞搞就好，反正是静态查询，双指针啊啥的随便上。。
这题的话第一维排序，第二维cdq，第三维树状数组。
直接做就好啦。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=2e5+5;
int m;
typedef long long ll;
struct node
{
    int a,b,c,cnt,ans;
    bool operator < (const node &x) const 
    {
        if (x.a!=a)return a<x.a;
        if (x.b!=b)return b<x.b;
        return c<x.c;
    }
}a[N],tmp[N];
bool cmp(const node &x,const node &y)
{
    if (x.b!=y.b)return x.b<y.b;
    if (x.c!=y.c)return x.c<y.c;
    return x.a<y.a;
}
int n,K,tot,c[N+10],tim[N],cnt[N],T;
inline int lowbit(int x)
{
    return x&(-x);
}
inline void add(int x,int val)
{
    while (x<=K)
    {
        if (tim[x]!=T)c[x]=0;
        tim[x]=T;
        c[x]+=val;
        x+=lowbit(x);
    }
}
inline int query(int x)
{
    int ret=0;
    while (x>0)
    {
        if (tim[x]==T)
        ret+=c[x];
        x-=lowbit(x);
    }
    return ret;
}
inline void cdq(int l,int r)
{
    if (l==r)
    {
        a[l].ans+=a[l].cnt-1;
        return;
    }
    int mid=(l+r)>>1;
    int l1,l2;
    l1=l,l2=mid+1;
    fo(i,l,r)
    {
        if (a[i].a<=mid)tmp[l1++]=a[i];
        else tmp[l2++]=a[i];
    }
    fo(i,l,r)a[i]=tmp[i];
    cdq(l,mid);
    int j=l;T++;
    fo(i,mid+1,r)
    {
        for(;j<=mid&&a[j].b<=a[i].b;j++)add(a[j].c,a[j].cnt);
        a[i].ans+=query(a[i].c);
    }
    cdq(mid+1,r);
    l1=l,l2=mid+1;
    fo(i,l,r)
    if (l1<=mid&&(l2>r||cmp(a[l1],a[l2])))tmp[i]=a[l1++];
    else tmp[i]=a[l2++];
    fo(i,l,r)a[i]=tmp[i];
} 
int main()
{
    scanf("%d%d",&n,&K);
    fo(i,1,n)
    {
        scanf("%d%d%d",&a[i].a,&a[i].b,&a[i].c);
        a[i].cnt++;
    }
    sort(a+1,a+1+n);
    fo(i,1,n)
    {
        if (i==1||a[i-1]<a[i])a[++tot]=a[i];
        else a[tot].cnt++;
    }
    fo(i,1,tot)a[i].a=i;
    sort(a+1,a+1+tot,cmp);
    cdq(1,tot);
    fo(i,1,tot)cnt[a[i].ans]+=a[i].cnt;
    fo(i,0,n-1)printf("%d\n",cnt[i]);
}