BZOJ 3262 陌上花开 (三维偏序CDQ+树状数组)

题目大意：

题面传送门

三维偏序裸题

首先，把三元组关于$a_{i}$排序

然后开始$CDQ$分治，回溯后按$b_{i}$排序

现在要处理左侧对右侧的影响了，显然现在左侧三元组的$a_{i}$都小于等于右侧

而$c_{i}$这一维需要用权值树状数组维护

归并排序时，已知左侧右侧两个指针分别是$i,j$

如果$b_{i} \leq bj$，将左侧已经遍历过的三元组的$c_{i}$推入树状数组，然后$i++$

如果$b_{i}>bj$，那么右侧能取到的贡献就是树状数组内$\leq c_{j}$的三元组数量，然后$j++$

可三元组有重复的啊！

如果我们不去重，假设有两个相同的三元组，靠前的三元组会得不到后面的贡献

还需要去重，并额外记录相同三元组数量

这道题的判定条件三个元是小于等于

发现我们只关于$a_{i}$排序也不行啊

在归并过程中会发现有一些$b_{i}$或者是$c_{i}$明明比较大，却得不到后面较小的贡献

这是由于小于等于的情况也合法，而我们必须让较大的在后面，才能保证后面能取到前面的贡献

所以排序过程中第二维还要按$b_{i}$排序，第三维按$c_{i}$

 #include <vector>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define N1 100100
 #define ll long long
 #define dd double
 #define inf 0x3f3f3f3f3f3f3f3fll
 using namespace std;

 int gint()
 {
     int ret=,fh=;char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
     return ret*fh;
 }
 int n,nn,K;
 struct node{
 int id,num,a,b,c;
 friend bool operator == (const node &s1,const node &s2)
 {return (s1.a==s2.a)&&(s1.b==s2.b)&&(s1.c==s2.c);}
 friend bool operator < (const node &s1,const node &s2)
 {
     if(s1.a!=s2.a) return s1.a<s2.a;
     if(s1.b!=s2.b) return s1.b<s2.b;
     return s1.c<s2.c;
 }
 }p[N1],t[N1],tmp[N1];

 struct b_{i}T{
 int s[N1<<];
 void update(int x,int w){for(int i=x;i<=K;i+=(i&(-i))) s[i]+=w;}
 int query(int x){int ans=; for(int i=x;i>;i-=(i&(-i))) ans+=s[i]; return ans;}
 }b;
 int f[N1],hs[N1],que[N1],tl;
 void CDQ(int L,int R)
 {
     if(R-L<=) return;
     int M=(L+R)>>;
     CDQ(L,M); CDQ(M,R);
     int i=L,j=M,cnt=;
     while(i<M&&j<R)
     {
         if(t[i].b<=t[j].b){
             b.update(t[i].c,t[i].num);
             tmp[++cnt]=t[i]; i++; que[++tl]=cnt;
         }else{
             f[t[j].id]+=b.query(t[j].c);
             tmp[++cnt]=t[j]; j++;
         }
     }
     while(i<M){tmp[++cnt]=t[i]; i++;}
     while(j<R){f[t[j].id]+=b.query(t[j].c); tmp[++cnt]=t[j]; j++;}
     while(tl) i=que[tl--],b.update(tmp[i].c,-tmp[i].num);
     for(i=L;i<R;i++) t[i]=tmp[i-L+];
 }

 int main()
 {
     scanf("%d%d",&n,&K);
     int i;
     for(i=;i<=n;i++) p[i].a=gint(),p[i].b=gint(),p[i].c=gint();
     sort(p+,p+n+);
     for(i=;i<=n;i++)
         if(!(p[i]==p[i-])) t[++nn]=p[i],t[nn].num=;
         else t[nn].num++;
     for(i=;i<=nn;i++) t[i].id=i;
     CDQ(,nn+);
     for(i=;i<=nn;i++) hs[f[t[i].id]+t[i].num-]+=t[i].num;
     for(i=;i<n;i++) printf("%d\n",hs[i]);
     return ;
 }