【BZOJ3262】 陌上花开

Description

有n朵花，每朵花有三个属性：花形(s)、颜色(c)、气味(m)，又三个整数表示。现要对每朵花评级，一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美丽，当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb。显然，两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。

Input

第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值。

以下N行，每行三个整数si, ci, mi (1 <= si, ci, mi <= K)，表示第i朵花的属性

Output

包含N行，分别表示评级为0...N-1的每级花的数量。

Sample Input

10 3
3 3 3
2 3 3
2 3 1
3 1 1
3 1 2
1 3 1
1 1 2
1 2 2
1 3 2
1 2 1

Sample Output

3
1
3
0
1
0
1
0
0
1

HINT

1 <= N <= 100,000, 1 <= K <= 200,000

 

三维偏序集，黄学长告诉我们可以

1.树状数组+treap √

2.cdq分治

目前只会树套树的做法，然而我的树套树真心很弱

 #include<iostream>
 #include<cstdio>
 #include<cstdlib>
 #include<algorithm>
 #define M 5000005
 using namespace std;
 int n,m,sz,tmp;
 int root[],ans[],sum[];
 int size[M],ls[M],rs[M],v[M],w[M],rnd[M];
 struct node{int a,b,c;}a[];
 int lowbit(int x){return x&(-x);}
 bool cmp(node a,node b){
     if (a.a==b.a&&a.b==b.b) return (a.c<b.c);
     if (a.a==b.a) return (a.b<b.b);
     return (a.a<b.a);
 }
 void updata(int x){size[x]=size[ls[x]]+size[rs[x]]+w[k];}
 void lturn(int &k){int t=rs[k];rs[k]=ls[t];ls[t]=k;updata(k);updata(t);k=t;}
 void rturn(int &k){int t=ls[k];ls[k]=rs[t];rs[t]=k;updata(k);updata(t);k=t;}
 void ins(int &k,int num){
     if (!k){sz++;k=sz;rnd[k]=rand();v[k]=num;size[k]=w[k]=;return;}
     size[k]++;
     if (v[k]==num) {w[k]++;return;}
     else if (v[k]>num){ins(ls[k],num);if (rnd[ls[k]]<rnd[k]) rturn(k);}
     else {ins(rs[k],num);if (rnd[rs[k]]<rnd[k]) lturn(k);}
 }

 void getrank(int k,int num){
     if (!k) return;
     int l=ls[k],r=rs[k];
     if (v[k]==num) {tmp+=size[l]+w[k];return;}
     else if (v[k]>num) getrank(l,num);
     else if (v[k]<num){tmp+=size[l]+w[k];getrank(r,num);}
 }

 void ask(int x,int num){
     for (int i=x;i;i-=lowbit(i))
         getrank(root[i],num);
 }

 void insert(int x,int num){
     for (int i=x;i<=m;i+=lowbit(i))
         ins(root[i],num);
 }

 int main(){
     scanf("%d%d",&n,&m);
     for (int i=;i<=n;i++) scanf("%d%d%d",&a[i].a,&a[i].b,&a[i].c);
     sort(a+,a+n+,cmp);
     for (int i=;i<=n;i++){
         if (a[i].a==a[i+].a&&a[i].b==a[i+].b&&a[i].c==a[i+].c&&i!=n) sum[i+]=sum[i]+;
         else{
             tmp=;
             ask(a[i].b,a[i].c);
             ans[tmp]+=sum[i]+;
         }
         insert(a[i].b,a[i].c);
     }
     for (int i=;i<n;i++)printf("%d\n",ans[i]);
 }