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排序即可。
那还有重复的呢？去重呗！
1 #include<stdio.h> 2 #include<string.h> 3 #include<stdlib.h> 4 #include<algorithm> 5 //#include<iostream> 6 using namespace std; 7 8 int n,m; 9 #define maxn 200011 10 int ord[maxn],tmpord[maxn]; 11 struct BIT 12 { 13 int a[maxn]; 14 BIT() {memset(a,0,sizeof(a));} 15 void add(int x,int v) {for (;x<=n;x+=x&-x) a[x]+=v;} 16 int query(int x) {int ans=0;for (;x;x-=x&-x) ans+=a[x];return ans;} 17 }t; 18 19 struct Point 20 { 21 int x,y,z,cnt; 22 bool operator < (const Point &b) const 23 {return x<b.x || (x==b.x && y<b.y) || (x==b.x && y==b.y && z<b.z);} 24 }q[maxn],p[maxn]; 25 int ans[maxn],ansnum[maxn]; 26 void solve(int L,int R) 27 { 28 if (L==R) {ans[L]+=p[L].cnt-1;ord[L]=L;return;} 29 const int mid=(L+R)>>1; 30 solve(L,mid); 31 solve(mid+1,R); 32 int i=L,j=mid+1,k=L; 33 while (i<=mid && j<=R) 34 { 35 if (p[ord[i]].y<=p[ord[j]].y) 36 { 37 tmpord[k++]=ord[i]; 38 t.add(p[ord[i]].z,p[ord[i]].cnt); 39 i++; 40 } 41 else 42 { 43 tmpord[k++]=ord[j]; 44 ans[ord[j]]+=t.query(p[ord[j]].z); 45 j++; 46 } 47 } 48 for (;j<=R;j++) ans[ord[j]]+=t.query(p[ord[j]].z),tmpord[k++]=ord[j]; 49 for (int ii=L;ii<i;ii++) t.add(p[ord[ii]].z,-p[ord[ii]].cnt); 50 for (;i<=mid;i++) tmpord[k++]=ord[i]; 51 for (int x=L;x<=R;x++) ord[x]=tmpord[x]; 52 } 53 54 int lisan[maxn]; 55 int main() 56 { 57 scanf("%d%d",&n,&m); 58 for (int i=1;i<=n;i++) 59 { 60 scanf("%d%d%d",&q[i].x,&q[i].y,&q[i].z); 61 lisan[i]=q[i].z; 62 } 63 sort(lisan+1,lisan+1+n); 64 for (int i=1;i<=n;i++) q[i].z=lower_bound(lisan+1,lisan+1+n,q[i].z)-lisan; 65 sort(q+1,q+1+n); 66 int tot=0;q[0].x=-0x3f3f3f3f; 67 for (int i=1;i<=n;i++) 68 { 69 if (q[i-1].x==q[i].x && q[i-1].y==q[i].y && q[i-1].z==q[i].z) p[tot].cnt++; 70 else p[++tot]=q[i],p[tot].cnt=1; 71 } 72 solve(1,tot); 73 for (int i=1;i<=tot;i++) ansnum[ans[i]]+=p[i].cnt; 74 for (int i=0;i<n;i++) printf("%d\n",ansnum[i]); 75 return 0; 76 }
