【BZOJ 3262】 3262: 陌上花开 （CDQ分治）

3262: 陌上花开

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

Source

树套树 CDQ分治

【分析】

　　三维偏序。。经典CDQ。。

　　如此文艺的名字。。

　　先把属性完全相同的合并，然后三个关键字一起排序，可以避免后影响前的问题。

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 using namespace std;
 #define Maxn 100010
 #define Maxk 200010

 struct node
 {
     int a,b,c,sm,ans;
 }t[Maxn];

 int mx=;

 int mymax(int x,int y) {return x>y?x:y;}

 bool cmp(node x,node y)
 {
     if(x.a!=y.a) return x.a<y.a;
     if(x.b!=y.b) return x.b<y.b;
     return x.c<y.c;
 }

 bool cmp2(int x,int y) {return (t[x].b==t[y].b)?(x<y):(t[x].b<t[y].b);}

 int c[Maxk],num[Maxn];
 void add(int x,int y)
 {
     for(int i=x;i<=mx;i+=i&(-i)) c[i]+=y;
 }

 int query(int x)
 {
     int ans=;
     for(int i=x;i>=;i-=i&(-i))
         ans+=c[i];
     return ans;
 }

 int nw[Maxn];
 void ffind(int l,int r)
 {
     if(l==r) return;
     int mid=(l+r)>>;
     nw[]=;
     for(int i=l;i<=r;i++) nw[++nw[]]=i;
     sort(nw+,nw++nw[],cmp2);
     for(int i=;i<=nw[];i++)
     {
         if(nw[i]<=mid) add(t[nw[i]].c,t[nw[i]].sm);
         else
         {
             t[nw[i]].ans+=query(t[nw[i]].c);
         }
     }
     for(int i=;i<=nw[];i++) if(nw[i]<=mid) add(t[nw[i]].c,-t[nw[i]].sm);
     ffind(l,mid);ffind(mid+,r);
 }

 int main()
 {
     int n,k;
     scanf("%d%d",&n,&k);
     for(int i=;i<=n;i++)
     {
         scanf("%d%d%d",&t[i].a,&t[i].b,&t[i].c);
         mx=mymax(mx,t[i].c);
         t[i].sm=;
     }
     sort(t+,t++n,cmp);
     int len=;
     for(int i=;i<=n;i++)
     {
         if(t[i].a!=t[len].a||t[i].b!=t[len].b||t[i].c!=t[len].c) t[++len]=t[i];
         else t[len].sm++;
     }
     for(int i=;i<=len;i++) t[i].ans=;
     memset(c,,sizeof(c));
     ffind(,len);
     memset(num,,sizeof(num));
     // for(int i=1;i<=n;i++) printf("%d\n",t[i].ans);
     for(int i=;i<=len;i++) num[t[i].ans+t[i].sm-]+=t[i].sm;
     for(int i=;i<n;i++) printf("%d\n",num[i]);
     return ;
 }

2017-01-19 09:17:11