[bzoj3262][cdq分治][树状数组]陌上花开

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

题解

cdq分治三维偏序裸题
很明显可以看出来三维就是x,y,z吧
首先对询问去重
第一维排序
第二维cdq
第三维树状数组统计答案并记录
最后根据记录的答案输出即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
int s[410000],zmax;
int lowbit(int x){return x&-x;}
void change(int x,int c){while(x<=zmax){s[x]+=c;x+=lowbit(x);}}
int findsum(int x){int ret=0;while(x>=1){ret+=s[x];x-=lowbit(x);}return ret;}
struct node
{
    int x,y,z,f,tmp;//表示和他相同的数目 
}a[210000],t[210000];int tot;
void ins(int x,int y,int z){tot++;a[tot].x=x;a[tot].y=y;a[tot].z=z;a[tot].tmp=1;a[tot].f=0;}
bool cmp(node n1,node n2)
{
    if(n1.x!=n2.x)return n1.x<n2.x;
    if(n1.y!=n2.y)return n1.y<n2.y;
    return n1.z<n2.z;
}
void cdq(int l,int r)
{
    if(l==r)return ;
    int mid=(l+r)/2;
    cdq(l,mid);cdq(mid+1,r);
    int i=l,j=mid+1,len=l;
    while(i<=mid && j<=r)
    {
        if(a[i].y<=a[j].y)change(a[i].z,a[i].tmp),t[len++]=a[i++];
        else a[j].f+=findsum(a[j].z),t[len++]=a[j++];
    }
    while(i<=mid)change(a[i].z,a[i].tmp),t[len++]=a[i++];
    while(j<=r)
    {
        a[j].f+=findsum(a[j].z);
        t[len++]=a[j++];
    }
    for(int i=l;i<=mid;i++)change(a[i].z,-a[i].tmp);
    for(int i=l;i<=r;i++)a[i]=t[i];
}
int n,ans[210000],N;
void pre()
{
    for(int i=1;i<=n;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        ins(x,y,z);
    }
    sort(a+1,a+1+n,cmp);
    N=n;
    n=1;
    for(int i=2;i<=N;i++)
    {
        if(a[i].x==a[n].x && a[i].y==a[n].y && a[i].z==a[n].z)a[n].tmp++;
        else a[++n]=a[i];
    }
}
int main()
{
    tot=0;
    scanf("%d%d",&n,&zmax);
    pre();
    cdq(1,n);
    for(int i=1;i<=n;i++)ans[a[i].f+a[i].tmp-1]+=a[i].tmp;
    for(int i=0;i<N;i++)printf("%d\n",ans[i]);
    return 0;
}