树套树【bzoj3262】陌上花开

/*

【bzoj3262】陌上花开

2014年6月19日1,2430

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的每级花的数量。*/

//树套树  cdq分治不会

//卡时函数大视野竟然不认，给我报了n遍RE，气死我了。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#define M 5000005
using namespace std;
struct shu
{
 int x,y,z;
}a[100005];
int n,m,sum[100005],ans[100005],sz,root[200005],w[M],v[M],sz1[M],lr[M],rr[M],rou[M];
int tmp;
bool cmp(shu a1,shu a2)
{
 if(a1.x==a2.x&&a1.y==a2.y)
   return a1.z<a2.z;
 if(a1.x==a2.x)
   return a1.y<a2.y;
 return a1.x<a2.x;
}
void update(int a1)
{
 sz1[a1]=sz1[lr[a1]]+sz1[rr[a1]]+w[a1];
}
void zuo(int &a1)
{
 int t=lr[a1];
 lr[a1]=rr[t];
 rr[t]=a1;
 update(a1);
 update(t);
 a1=t;
}
void you(int &a1)
{
 int t=rr[a1];
 rr[a1]=lr[t];
 lr[t]=a1;
 update(a1);
 update(t);
 a1=t;
}
void jiaa(int &a1,int a2)
{
 if(a1==0)
   {
    sz++;
    a1=sz;
    v[a1]=a2;
    sz1[a1]=1;
    w[a1]=1;
    rou[a1]=rand();
    return;
   } 
 sz1[a1]++;
 if(v[a1]==a2)
   {
    w[a1]++;
    return;
   }
 if(v[a1]>a2)
   {
    jiaa(lr[a1],a2);
    if(rou[lr[a1]]<rou[a1])
      zuo(a1);
   }
 else
   {
    jiaa(rr[a1],a2);
    if(rou[rr[a1]]<rou[a1])
      you(a1);
   }
}
void jia(int a1,int a2)
{
 for(;a1<=m;a1+=(a1&(-a1)))
   jiaa(root[a1],a2);
}
void askk(int a1,int a2)
{
 if(a1==0)
   return;
 if(v[a1]==a2)
    {
     tmp+=w[a1]+sz1[lr[a1]];
     return;
    }
 if(v[a1]>a2)
   askk(lr[a1],a2);
 else
   {
    tmp+=w[a1]+sz1[lr[a1]];
    askk(rr[a1],a2);
   }
 return;
}
void ask(int a1,int a2)
{
 for(;a1;a1-=(a1&(-a1)))
   askk(root[a1],a2);
}
int main()
{
 scanf("%d%d",&n,&m);
 for(int i=1;i<=n;i++)
   scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
 sort(a+1,a+n+1,cmp);
 for(int i=1;i<=n;i++)
  {
   if(a[i].y==a[i+1].y&&a[i].z==a[i+1].z&&a[i].x==a[i+1].x&&i!=n)
     sum[i+1]+=sum[i]+1;
   else
     {
      tmp=0;
      ask(a[i].y,a[i].z);
      ans[tmp]+=sum[i]+1;
  }
   jia(a[i].y,a[i].z);
     }
 for(int i=0;i<n;i++)
   printf("%d\n",ans[i]);
    return 0;
}