[BZOJ]3262: 陌上花开 CDQ分治+树状数组

Description
有n朵花，每朵花有三个属性：花形(s)、颜色(c)、气味(m)，又三个整数表示。现要对每朵花评级，一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美丽，当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb。显然，两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。

题解：

三维偏序，CDQ分治裸题，大概就是先对第一维排序，忽略其影响，然后CDQ分治搞第二维，第三维用权值树状数组搞，注意对两朵完全相同的花的处理。写归并似乎比快排快很多？

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=100010;
const int Maxk=200010;
const int inf=2147483647;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int n,k,cnt=0,Ans[Maxn];
struct flower{int a,b,c,id,tot,ans;}A[Maxn],q[Maxn],temp[Maxn];
bool cmp(flower x,flower y){return((x.b==y.b)?x.c<y.c:x.b<y.b);}
bool operator < (flower x,flower y)
{
    if(x.a!=y.a)return x.a<y.a;
    if(x.b!=y.b)return x.b<y.b;
    return x.c<y.c;
}
bool operator == (flower x,flower y){return(x.a==y.a&&x.b==y.b&x.c==y.c);}
int s[Maxk];
void add(int x,int y){for(;x<=k;x+=(x&(-x)))s[x]+=y;}
int getsum(int x){int re=0;for(;x;x-=(x&(-x)))re+=s[x];return re;}
void solve(int l,int r)
{
    if(l==r)return;
    int mid=l+r>>1;
    solve(l,mid);solve(mid+1,r);
    int i=l,j=mid+1,len=0;
    while(i<=mid&&j<=r)
    {
        if(q[i].b<=q[j].b)add(q[i].c,q[i].tot),temp[++len]=q[i++];
        else q[j].ans+=getsum(q[j].c),temp[++len]=q[j++];
    }
    int t=i;
    while(i<=mid)temp[++len]=q[i++];
    while(j<=r)q[j].ans+=getsum(q[j].c),temp[++len]=q[j++];
    for(int p=l;p<t;p++)add(q[p].c,-q[p].tot);
    for(int p=1;p<=len;p++)q[p+l-1]=temp[p];
}
int main()
{
    n=read(),k=read();
    for(int i=1;i<=n;i++)A[i].a=read(),A[i].b=read(),A[i].c=read();
    sort(A+1,A+1+n);
    A[0].a=A[0].b=A[0].c=0;
    for(int i=1;i<=n;i++)
    {
        if(A[i]==A[i-1])q[cnt].tot++;
        else q[++cnt]=A[i],q[cnt].tot=1,q[cnt].id=cnt;
    }
    solve(1,cnt);
    for(int i=1;i<=cnt;i++)Ans[q[i].ans+q[i].tot-1]+=q[i].tot;
    for(int i=0;i<n;i++)printf("%d\n",Ans[i]);
}