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

CDQ分治。先按照x排序。然后我们用归并排序的方法排序y。同时把z放进树状数组

因为x排序好了。所以归并的时候保证左半边的x小于等于右半边。那么我们只需要在放入右半边数的时候记录有几个z比它小就可以了

【果然状态不行各种傻逼错调了好久】

#include<cstdio>#include<string>#include<cstring>#include<algorithm>using namespace std;struct number{     int x,y,z;     int p;     int s;}a[100001],b[100001],c[100001];int rank[100001];int sx[100001];inline bool cmpx(number x,number y){     if(x.x<y.x)          return true;     if(x.x==y.x)     {          if(x.y<y.y)               return true;          if(x.y==y.y)               if(x.z<y.z)                    return true;     }     return false;}int tr[200001];int k;inline int lowbit(int x){     return x&(-x);}inline void add(int x,int xx){     int i;     for(i=x;i<=k;i+=lowbit(i))          tr[i]+=xx;}inline int ask(int x){     int i;     int ans=0;     for(i=x;i>=1;i-=lowbit(i))          ans+=tr[i];     return ans;}inline void solve(int l,int r){ int mid=(l+r)/2; if(l!=r) {      solve(l,mid);      solve(mid+1,r);      int i1=l,i2=mid+1;      int p=0;      //memset(tr,0,sizeof(tr));      while(i1<=mid&&i2<=r)      {           if(a[i1].y<=a[i2].y)           {                p++;                c[p]=a[i1];                add(c[p].z,c[p].s);                i1++;           }           else           {                p++;                c[p]=a[i2];                rank[c[p].p]+=ask(c[p].z);                i2++;           }      }      if(i1<=mid)      {           int i;               for(i=i1;i<=mid;i++)               {                    p++;                    c[p]=a[i];                    add(c[p].z,c[p].s);               }      }      else if(i2<=r)      {           int i;               for(i=i2;i<=r;i++)               {                    p++;                    c[p]=a[i];                    rank[c[p].p]+=ask(c[p].z);               }      }      int i;      for(i=l;i<=mid;i++)           add(a[i].z,-a[i].s);      for(i=l;i<=r;i++)           a[i]=c[i-l+1]; }}int main(){//     freopen("test.in","r",stdin);//     freopen("test.out","w",stdout);     int n;     scanf("%d%d",&n,&k);     int i;     for(i=1;i<=n;i++)     {          scanf("%d%d%d",&b[i].x,&b[i].y,&b[i].z);          b[i].p=i;     }     sort(b+1,b+1+n,cmpx);     int t=0;     for(i=1;i<=n;i++)     {          if(b[i].x!=b[i-1].x||b[i].y!=b[i-1].y||b[i].z!=b[i-1].z)          {               t++;               a[t]=b[i];               a[t].s=1;               a[t].p=t;          }          else               a[t].s++;     }     for(i=1;i<=t;i++)          b[i]=a[i];     solve(1,t);     //for(i=1;i<=t;i++)       //   printf("%d\n",rank[i]);     for(i=1;i<=t;i++)          sx[rank[i]+b[i].s]+=b[i].s;     for(i=1;i<=n;i++)          printf("%d\n",sx[i]);     return 0;}