3262: 陌上花开

来源:互联网 发布:友盟数据分析 编辑:程序博客网 时间:2024/06/02 22:57

题目链接

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

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

我的收获:CDQ神啊

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>using namespace std;#define M 200005int 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;}
原创粉丝点击