BZOJ-3262-陌上花开

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

第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值。
以下N行，每行三个整数si, ci, mi (1 <= si, ci, mi <= K)，表示第i朵花的属性

包含N行，分别表示评级为0…N-1的每级花的数量。

CDQ分治经典题目，求解三维偏序问题

CDQ分治与普通分治法的区别在于：普通分治法分为两个区间互不影响，CDQ分治左区间队有区间产生影响，有区间队左区间无影响，需额外计算左区间队右区间的影响。
CDQ分治
参考题解

#include<bits/stdc++.h>
#define lowbit(x) x&(-x)
#define de(x) cout<<#x<<" = "<<endl;
using namespace std;
typedef long long ll;
const ll maxn=2e5+10;

ll Tree[maxn];//树状数组
ll ans[maxn];

struct node{
    ll x,y,z,cnt,ans;
}a[maxn];


bool cmpx(node a,node b)
{
	if(a.x<b.x) return 1;
	if(a.x>b.x) return 0;
	if(a.y<b.y) return 1;
	if(a.y>b.y) return 0;
	if(a.z<b.z) return 1;
	return 0;
}

bool cmpy(node a,node b)
{
	if(a.y<b.y) return 1;
	if(a.y>b.y) return 0;
	if(a.z<b.z) return 1;
	if(a.z>b.z) return 0;
	if(a.x<b.x) return 1;
	return 0;
}

void update(ll x,ll val){
    while(x<maxn){
        Tree[x]+=val;
        x+=lowbit(x);
    }
}

ll query(ll x){
    ll sum=0;
    while(x>0){
        sum+=Tree[x];
        x-=lowbit(x);
    }
    return sum;
}

void CDQ(ll l,ll r){
    if(l==r){
        a[l].ans+=a[l].cnt-1;//自己不算
        return;
    }
    ll mid=(l+r)>>1;
    CDQ(l,mid);
    CDQ(mid+1,r);
    sort(a+l,a+mid+1,cmpy);
    sort(a+mid+1,a+r+1,cmpy);//分别按y排序
    ll j=l;
    for(ll i=mid+1;i<=r;i++){
        while(j<=mid&&a[j].y<=a[i].y){
            update(a[j].z,a[j].cnt);
            j++;
        }
        a[i].ans+=query(a[i].z);
    }
    for(ll i=l;i<j;i++){
        update(a[i].z,-a[i].cnt);//清空树状数组
    }
}

int main(){
    ll n,k;
    scanf("%lld%lld",&n,&k);
    for(ll i=1;i<=n;i++){
        scanf("%lld%lld%lld",&a[i].x,&a[i].y,&a[i].z);
        a[i].ans=0;
        a[i].cnt=1;
    }
    sort(a+1,a+1+n,cmpx);
    ll tot=0;
    for(ll i=1;i<=n;i++){
        if(i!=1&&a[i].x==a[i-1].x&&a[i].y==a[i-1].y&&a[i].z==a[i-1].z)
            a[tot].cnt++;
        else
            a[++tot]=a[i];
    }
    CDQ(1,tot);
    for(ll i=1;i<=tot;i++){
        ans[a[i].ans]+=a[i].cnt;
    }
    for(ll i=0;i<n;i++){
        printf("%lld\n",ans[i]);
    }
}

提交不了代码，样例过了就算过了吧