【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分治

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

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstdlib>
 4 #include<algorithm>
 5 #define M 5000005
 6 using namespace std;
 7 int n,m,sz,tmp;
 8 int root[200005],ans[100005],sum[100005];
 9 int size[M],ls[M],rs[M],v[M],w[M],rnd[M];
10 struct node{int a,b,c;}a[100005];
11 int lowbit(int x){return x&(-x);}
12 bool cmp(node a,node b){
13     if (a.a==b.a&&a.b==b.b) return (a.c<b.c);
14     if (a.a==b.a) return (a.b<b.b);
15     return (a.a<b.a);
16 }
17 void updata(int x){size[x]=size[ls[x]]+size[rs[x]]+w[k];}
18 void lturn(int &k){int t=rs[k];rs[k]=ls[t];ls[t]=k;updata(k);updata(t);k=t;}
19 void rturn(int &k){int t=ls[k];ls[k]=rs[t];rs[t]=k;updata(k);updata(t);k=t;}
20 void ins(int &k,int num){
21     if (!k){sz++;k=sz;rnd[k]=rand();v[k]=num;size[k]=w[k]=1;return;}
22     size[k]++;
23     if (v[k]==num) {w[k]++;return;}
24     else if (v[k]>num){ins(ls[k],num);if (rnd[ls[k]]<rnd[k]) rturn(k);}
25     else {ins(rs[k],num);if (rnd[rs[k]]<rnd[k]) lturn(k);}
26 }
27 
28 void getrank(int k,int num){
29     if (!k) return;
30     int l=ls[k],r=rs[k];
31     if (v[k]==num) {tmp+=size[l]+w[k];return;}
32     else if (v[k]>num) getrank(l,num);
33     else if (v[k]<num){tmp+=size[l]+w[k];getrank(r,num);}
34 }
35 
36 void ask(int x,int num){
37     for (int i=x;i;i-=lowbit(i))
38         getrank(root[i],num);
39 }
40 
41 void insert(int x,int num){
42     for (int i=x;i<=m;i+=lowbit(i))
43         ins(root[i],num);
44 }
45 
46 int main(){
47     scanf("%d%d",&n,&m);
48     for (int i=1;i<=n;i++) scanf("%d%d%d",&a[i].a,&a[i].b,&a[i].c);
49     sort(a+1,a+n+1,cmp);
50     for (int i=1;i<=n;i++){
51         if (a[i].a==a[i+1].a&&a[i].b==a[i+1].b&&a[i].c==a[i+1].c&&i!=n) sum[i+1]=sum[i]+1;
52         else{
53             tmp=0;
54             ask(a[i].b,a[i].c);
55             ans[tmp]+=sum[i]+1;
56         }
57         insert(a[i].b,a[i].c);
58     }
59     for (int i=0;i<n;i++)printf("%d
",ans[i]);
60 }