hdu_5618_Jam's problem again(cdq分治+lowbit)

题目链接：hdu_5618_Jam's problem again

题意：

给你n个点，每个点有一个坐标(x,y,z)，找出有ans个点，3个坐标都比该点小，这个点的level就为ans,然后让你输出所有点的ans.

题解：

对于第一维，直接排序，后面两维的处理可以用线段树套lowbit，但空间用的很大，这里满足运用cdq分治的条件，所有用cdq分治套lowbit，常数很小。

注意相同点的处理。

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;i++)
using namespace std;

const int N=1e5+7;
int t,n,sum[N];

struct node
{
    int x,y,z,cnt,id;
    bool operator < (const node &b)const
    {
        if(x==b.x&&y==b.y)return z<b.z;
        if(x==b.x)return y<b.y;
        return x<b.x;
    }
    bool operator == (const node &b)const{return x==b.x&&y==b.y&&z==b.z;}
}a[N],b[N];

bool cmp(node a,node b){return a.id<b.id;}

inline void add(int x,int c){while(x<N)sum[x]+=c,x+=x&-x;}
inline  int ask(int x){int an=0;while(x)an+=sum[x],x-=x&-x;return an;}

void solve(int l,int r)
{
    if(l==r)return;
    int mid=l+r>>1;
    solve(l,mid),solve(mid+1,r);
    for(int i=l,j=l,k=mid+1;i<=r;i++)
    {
        if((a[j].y<=a[k].y||r<k)&&j<=mid)b[i]=a[j++],add(b[i].z,1);
        else b[i]=a[k++],b[i].cnt+=ask(b[i].z);
    }//当左半边的y大于当前右半边的y时，先将这个点更新了，这样就保证了前面的y都是比这个点小
    F(i,l,mid)add(a[i].z,-1),a[i]=b[i];
    F(i,mid+1,r)a[i]=b[i];
}

int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        F(i,1,n)
        {
            scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
            a[i].id=i,a[i].cnt=0;
        }
        sort(a+1,a+1+n);
        for(int i=n-1;i>0;i--)if(a[i]==a[i+1])a[i].cnt=a[i+1].cnt+1;
        solve(1,n);
        sort(a+1,a+1+n,cmp);
        F(i,1,n)printf("%d
",a[i].cnt);
    }
    return 0;
}