bzoj 3262 陌上花开

本质是一个三维偏序，一位排序后cdq分治，一维在子函数里排序，一维用树状数组维护。

把三维相等的合并到一个里面。

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #define N 100005
 using namespace std;
 int n,k,m;
 struct node
 {
     int x,y,z,num,op;
     int yuan;
     bool operator == (const node &b)
     {
         return x==b.x&&y==b.y&&z==b.z;
     }
 }a[N],b[N];
 int shan[N];
 int ans[N];
 bool cmp1(node aa,node bb)
 {
     if(aa.num==)return ;
     if(bb.num==)return ;
     if(aa.x!=bb.x)return aa.x<bb.x;
     if(aa.y!=bb.y)return aa.y<bb.y;
     if(aa.z!=bb.z)return aa.z<bb.z;

 }
 bool cmp2(node aa,node bb)
 {
     if(aa.y!=bb.y)return aa.y<bb.y;
     return aa.op<bb.op;
 }
 int c[];
 void add(int x,int y)
 {
     for(int i=x;i<=k;i+=(i&(-i)))c[i]+=y;
 }
 int qur(int x)
 {
     int now=;
     for(int i=x;i;i-=(i&(-i)))now+=c[i];
     return now;
 }
 void solve(int l,int r)
 {
     if(l==r)return ;
     int mid=(l+r)>>;
     solve(l,mid);solve(mid+,r);
     int cnt=;
     for(int i=l;i<=mid;i++)
     {
         b[++cnt]=a[i];
         b[cnt].op=;
     }
     for(int i=mid+;i<=r;i++)
     {
         b[++cnt]=a[i];
         b[cnt].op=;
     }
     sort(b+,b+cnt+,cmp2);
     for(int i=;i<=cnt;i++)
     {
         if(b[i].op==)add(b[i].z,b[i].num);
         else ans[b[i].yuan]+=qur(b[i].z);
     }
     for(int i=;i<=cnt;i++)if(b[i].op==)add(b[i].z,-b[i].num);
 }
 int zui[N];
 int main()
 {
     scanf("%d%d",&n,&k);m=n;
     for(int i=;i<=n;i++)
     {
         scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
         a[i].num=;
     }
     sort(a+,a+n+,cmp1);
     for(int i=;i<=n;i++)
     {
         if(a[i]==a[i-])
         {
             a[shan[i-]].num++;
             shan[i]=shan[i-];
             a[i].num--;
         }
         else shan[i]=i;
     }
     sort(a+,a+n+,cmp1);
     while(a[n].num==)n--;
     for(int i=;i<=n;i++)a[i].yuan=i;
     solve(,n);
     for(int i=;i<=n;i++)zui[ans[a[i].yuan]+a[i].num-]+=a[i].num;
     for(int i=;i<=m-;i++)printf("%d\n",zui[i]);
     return ;
 }