3262: 陌上花开

题目链接

题目大意：给定一堆花，每个花有三个属性，定义一朵花比另一朵花美丽当期仅当三个值都大于等于另一朵花 定义花的评级为没有它美丽的花的数量 求评级为0~N-1的花的数量

题解：三维偏序，x排序，yCDQ分治，z树状数组

我的收获：CDQ神啊

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

#define M 200005

int n,k,cnt;
int c[M],num[M];

struct flo{int x,y,z,v,ans;}a[M],p[M];

bool cmpa(flo a,flo b){return a.x<b.x||(a.x==b.x&&a.y<b.y)||(a.x==b.x&&a.y==b.y&&a.z<b.z);}
bool cmpp(flo a,flo b){return a.y<b.y||(a.y==b.y&&a.z<b.z);}

void add(int x,int v){for(int i=x;i<=k;i+=i&(-i)) c[i]+=v;}
int query(int x){int ret=0;for(int i=x;i>0;i-=i&(-i)) ret+=c[i];return ret;}

void solve(int l,int r){//一维排序,二维CDQ,三维树状数组 
    if(l==r){p[l].ans+=p[l].v-1;return ;}
    int mid=(l+r)>>1,i=l,j=mid+1;
    solve(l,mid);solve(mid+1,r);
    sort(p+l,p+mid+1,cmpp);sort(p+mid+1,p+r+1,cmpp);
    for(;j<=r;j++){
        for(;i<=mid&&p[i].y<=p[j].y;i++) add(p[i].z,p[i].v);//x有序，z树状数组维护，只考虑y 
        p[j].ans+=query(p[j].z);
    }
    for(j=l;j<i;j++) add(p[j].z,-p[j].v);//删掉 
}

void work()
{
    solve(1,cnt);
    for(int i=1;i<=cnt;i++) num[p[i].ans]+=p[i].v;
    for(int i=1;i<=n;i++) printf("%d\n",num[i]);
}

void init()
{
    scanf("%d%d",&n,&k);
    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,cmpa);
    for(int i=1,w=1;i<=n;i++,w++)
        if(a[i].x!=a[i+1].x||a[i].y!=a[i+1].y||a[i].z!=a[i+1].z)//把重复的合并,它们互相比对方美丽 
            p[++cnt]=a[i],p[cnt].v=w,w=0,p[cnt].ans=1;
}

int main()
{
    init();
    work();
    return 0;
}