cdq分治入门--BZOJ3262: 陌上花开

n<=100000个人，每个人三个属性Ai,Bi,Ci，一个人i的等级为Ai>=Aj,Bi>=Bj,Ci>=Cj的人数，求每个等级有多少人。

裸的三维偏序。按照常规思路，一维排序，一维归并，一维利用单调性或用树状数组维护，这里选择后者。

先按Ai排序，然后在分治过程中，solve(l,mid)，solve(mid+1,r)，然后考虑(l,mid)对(mid+1,r)答案的贡献，先把这两部分分别按B排序，然后两个指针一起扫，扫的过程中，把C的值作下标丢进树状数组，查询时相当于树状数组前缀和。

思路很好，可样例过不了。

原因：在计算过程中默认后面对前面是没有贡献的，但按Ai排序后，相同的Ai值可能引起后面对前面的贡献。

方法：保证后面对前面无贡献，一开始就先按A再按B最后按C排序即可。

那还有重复的呢？去重呗！

 #include<stdio.h>
 #include<string.h>
 #include<stdlib.h>
 #include<algorithm>
 //#include<iostream>
 using namespace std;

 int n,m;
 #define maxn 200011
 int ord[maxn],tmpord[maxn];
 struct BIT
 {
     int a[maxn];
     BIT() {memset(a,,sizeof(a));}
     void add(int x,int v) {for (;x<=n;x+=x&-x) a[x]+=v;}
     int query(int x) {int ans=;for (;x;x-=x&-x) ans+=a[x];return ans;}
 }t;

 struct Point
 {
     int x,y,z,cnt;
     bool operator < (const Point &b) const
     {return x<b.x || (x==b.x && y<b.y) || (x==b.x && y==b.y && z<b.z);}
 }q[maxn],p[maxn];
 int ans[maxn],ansnum[maxn];
 void solve(int L,int R)
 {
     if (L==R) {ans[L]+=p[L].cnt-;ord[L]=L;return;}
     const int mid=(L+R)>>;
     solve(L,mid);
     solve(mid+,R);
     int i=L,j=mid+,k=L;
     while (i<=mid && j<=R)
     {
         if (p[ord[i]].y<=p[ord[j]].y)
         {
             tmpord[k++]=ord[i];
             t.add(p[ord[i]].z,p[ord[i]].cnt);
             i++;
         }
         else
         {
             tmpord[k++]=ord[j];
             ans[ord[j]]+=t.query(p[ord[j]].z);
             j++;
         }
     }
     for (;j<=R;j++) ans[ord[j]]+=t.query(p[ord[j]].z),tmpord[k++]=ord[j];
     for (int ii=L;ii<i;ii++) t.add(p[ord[ii]].z,-p[ord[ii]].cnt);
     for (;i<=mid;i++) tmpord[k++]=ord[i];
     for (int x=L;x<=R;x++) ord[x]=tmpord[x];
 }

 int lisan[maxn];
 int main()
 {
     scanf("%d%d",&n,&m);
     for (int i=;i<=n;i++)
     {
         scanf("%d%d%d",&q[i].x,&q[i].y,&q[i].z);
         lisan[i]=q[i].z;
     }
     sort(lisan+,lisan++n);
     for (int i=;i<=n;i++) q[i].z=lower_bound(lisan+,lisan++n,q[i].z)-lisan;
     sort(q+,q++n);
     int tot=;q[].x=-0x3f3f3f3f;
     for (int i=;i<=n;i++)
     {
         if (q[i-].x==q[i].x && q[i-].y==q[i].y && q[i-].z==q[i].z) p[tot].cnt++;
         else p[++tot]=q[i],p[tot].cnt=;
     }
     solve(,tot);
     for (int i=;i<=tot;i++) ansnum[ans[i]]+=p[i].cnt;
     for (int i=;i<n;i++) printf("%d\n",ansnum[i]);
     return ;
 }