luogu3810 陌上花开 (cdq分治)

求三维偏序

设三维为a,b,c。先对a排序，这样i的偏序就只能<i。

然而排序的时候需要三个维度都判断一遍，最后还要去重，不然会出现实际应该记答案的数出现在它后面的情况。

（排序用的函数里不要写类似于<=之类的东西啊..会出奇奇怪怪的问题的(RE)）

然后分治来做，我们在做区间[l,r]的时候，先去做[l,m]和[m+1,r]

之后左区间[l,m],右区间[m+1,r]都已经按照b排好序了，而且左右两区间内部的答案已经统计过了，所以现在只要考虑左区间中满足（右区间的数）的数量就好了。

那么就也把[l,r]按照b排好序，在排的时候再用一个权值树状数组维护c，

也就是，如果这个点是左区间的点，就把它的c值对应的树状数组中+=这个点的重复数（刚才去重了）

　　　　如果这个点是右区间的点，就询问树状数组中<=它的c值的数量，然后加到这个点的答案里。

而且每次做的时候树状数组都要清空，但不能用memset来清，复杂度有问题。（一直迷信memset的速度，结果一查告诉我也就比循环清快一倍？？）

所以只要把刚才加过的再减回去就可以了。

复杂度$O(n*log_2n*log_2k)$

也可以cdq套cdq，然后不用树状数组，复杂度是一样的。

这样一直套下去，k维的话复杂度也就是$O(n*log^{k-1}_2n)$啦

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<cmath>
#include<ctime>
#define LL long long int
#define inf 0x3f3f3f3f
#define lowbit(x) ((x)&(-(x)))
using namespace std;
const int maxn=,maxk=;

LL rd(){
   LL x=;char c=getchar();int neg=;
   while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
   while(c>=''&&c<='') x=x*+c-'',c=getchar();
   return x*neg;
}

struct Node{
    int a,b,c,i;
}inp[maxn],num[maxn],tmp[maxn];
int iniN,N,K;
int tr[maxk],cnt[maxn],siz[maxn],ans[maxn];

inline bool cmp(Node a,Node b){return a.a==b.a?(a.b==b.b?a.c<b.c:a.b<b.b):a.a<b.a;}

inline void add(int x,int y){
    while(x&&x<=K) tr[x]+=y,x+=lowbit(x);
}
inline int query(int x){
    int re=;while(x) re+=tr[x],x-=lowbit(x);return re;
}

void cdq(int l,int r){
    int m=l+r>>,p=l,q=m+,t=;
    if(l>=r) return;
    cdq(l,m);cdq(m+,r);
    while(p<=m&&q<=r){
        if(num[p].b<=num[q].b){
            tmp[++t]=num[p];add(num[p].c,siz[num[p].i]);p++;
        }else{
            tmp[++t]=num[q];cnt[num[q].i]+=query(num[q].c);q++;
        }
    }    while(q<=r){
        tmp[++t]=num[q];cnt[num[q].i]+=query(num[q].c);q++;
    }for(int i=l;i<p;i++) add(num[i].c,-siz[num[i].i]);
    while(p<=m) tmp[++t]=num[p++];

    memcpy(num+l,tmp+,sizeof(Node)*t);
}

int main(){
    int i,j,k;
    iniN=N=rd();K=rd();
    for(i=;i<=N;i++){
        int a=rd(),b=rd(),c=rd();
        inp[i].a=a;inp[i].b=b;inp[i].c=c;num[i].i=i;
    }
    sort(inp+,inp+N+,cmp);//printf("ll");
    for(i=,j=;i<=N;i++){
        if(inp[i].a==inp[i-].a&&inp[i].b==inp[i-].b&&inp[i].c==inp[i-].c) inp[i].i=j,cnt[j]++,siz[j]++;
        else{
            inp[i].i=++j;num[j]=inp[i];siz[j]=;
        }
    }N=j;
    cdq(,N);
    for(i=;i<=N;i++) ans[cnt[i]]+=siz[i];
    for(i=;i<iniN;i++) printf("%d\n",ans[i]);
    return ;
}