【BZOJ 3262三维偏序】陌上花开

【BZOJ 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分治维护三维偏序,第一维sort即可，第二维分治即可，第三维套树状数组求解；

　　　!：对于细节处理要注意！在分治时原先的树状数组记录的由于y的情况不同会有误！

　　　所以要多一个数组记录。

}

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<ctime>
#include<cmath>
#include<map>
#include<set>
#define MAXX 1000010
#define lowbit(k) k&-k
using namespace std;
struct f{
  int x,y,z,num,ans;
}e[MAXX*2];
int fa[MAXX*2],tree[MAXX*2],now,n,k,anss[MAXX*2];
//                                   shuzhuangshuzu
int sum(int p){
  int num=0;
  for(int k=p;k;k-=lowbit(k))
    if(fa[k]==now)num+=tree[k];
  return num;
}
void add(int p,int num){
  for(int k=p;k<MAXX;k+=lowbit(k)){
    if(fa[k]!=now)tree[k]=0;//yi jing ji suan guo
    fa[k]=now;tree[k]+=num;
  }
}
//                                   sort!!!!!!
bool Comp(const f & a,const f & b){
  if(a.y==b.y)return a.z<a.z;
  return a.y<b.y;
}
bool comp(const f & a,const f & b){
  if(a.x==b.x){
    if(a.y==b.y)return a.z<b.z;
    return a.y<b.y;
  }
  return a.x<b.x;
}
//                                     CDQ   fenzhi!!
void CDQ(int l,int r){
  if(l>=r)return;
  int mid=(l+r)>>1;
  CDQ(l,mid),CDQ(mid+1,r);
  sort(e+l,e+mid+1,Comp),sort(e+mid+1,e+r+1,Comp);now++;int i,j;
  for(j=mid+1,i=l;j<=r;++j){
    for(;e[i].y<=e[j].y&&i<=mid;++i)add(e[i].z,e[i].num);
    e[j].ans+=sum(e[j].z);
   }
}
int main(){
  scanf("%d%d",&n,&k);int nn=1;
  for(int i=1;i<=n;++i)scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].z),e[i].num++;
  sort(e+1,e+n+1,comp);
  for(int i=2;i<=n;++i)
    if(e[i].x!=e[nn].x||e[i].y!=e[nn].y||e[i].z!=e[nn].z)e[++nn]=e[i];
    else e[nn].num++;
  CDQ(1,nn);
  for(int i=1;i<=nn;++i)anss[e[i].ans+e[i].num-1]+=e[i].num;
  for(int i=0;i<n;++i)printf("%d
",anss[i]);
  return 0;
}

 

 si