HDU - 5618 Jam's problem again 三维偏序 CDQ分治+树状数组

题目链接点这里

弱不会写树套树，，只能学·CDQ分治了、、、

关于CDQ分治，，推荐这么几篇博客。。

(1)【教程】简易CDQ分治教程&学习笔记

(2)【教程】CDQ套CDQ——四维偏序问题

(3) 关于偏序问题的一些总结

一篇篇按顺序往下看，，就懂了，，

然后，，如果用分治或者树套树（虽然我不会树套树），，每加一维都会多一个log，，所以到5维偏序的时候，，时间复杂度已经到 n (logn)^4了，，已经接近n^2的暴力复杂度了，，所以一般在5维及以上偏序问题，，分块bitset暴力的效果还是不错的

#include<algorithm>
#include<iostream>
#include<stdio.h>
#include<cstring>
#include<cmath>
using namespace std;
#define mem(x,y) memset(x,y,sizeof(x))
#define FIN freopen("input.txt","r",stdin)
#define fuck(x) cout<<x<<endl
const double eps=1e-7;
const int MX=111111;
#define INF 0x3f3f3f3f
#define INFLL 0x3f3f3f3f3f3f3f3f
typedef long long LL;
typedef pair<LL,LL> PLL;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
int n,k,block;
struct QQ
{
    int x,y,z,id;
    bool operator <(const QQ a)const
    {
        if(x!=a.x) return x<a.x;
        if(y!=a.y) return y<a.y;
        return z<a.z;
    }
    bool operator ==(const QQ a)const
    {
        return x==a.x&&y==a.y&&z==a.z;
    }
} Q[MX],tmp[MX];
bool cmp(QQ a,QQ b)
{
    return a.x<b.x;
}
int ans[MX];
int tree[MX*11];
int lowbit(int x)
{
    return x&-x;
}
int tree_sum(int x)
{
    int ans=0;
    while(x)
    {
        ans+=tree[x];
        x-=lowbit(x);
    }
    return ans;
}
void tree_add(int x,int d)
{
    while(x<=k)
    {
        tree[x]+=d;
        x+=lowbit(x);
    }
}
void CDQ(int l,int r)
{
    if(l==r)return ;
    int m=(l+r)>>1;
    CDQ(l,m);
    CDQ(m+1,r);
    int q=l,p=m+1,cnt=l;
    while(q<=m||p<=r)
    {
        if(p>r||(q<=m&&Q[q].y<=Q[p].y))
        {
            tree_add(Q[q].z,1);
            tmp[cnt++]=Q[q];
            q++;
        }
        else
        {
            ans[Q[p].id]+=tree_sum(Q[p].z);
            tmp[cnt++]=Q[p];
            p++;
        }
    }
    for(int i=l; i<=m; i++) tree_add(Q[i].z,-1);
    for(int i=l; i<=r; i++) Q[i]=tmp[i];
}
int arr[MX];
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        cin>>n;
        k=0;
        for(int i=0; i<n; i++)
        {
            scanf("%d%d%d",&Q[i].x,&Q[i].y,&Q[i].z);
            k=max(k,Q[i].z);
            Q[i].id=i;
        }
        sort(Q,Q+n);
        mem(ans,0);
        CDQ(0,n-1);
        for(int i=n-2; i>=0; i--)if(Q[i]==Q[i+1]) ans[Q[i].id]=ans[Q[i+1].id];
        for(int i=0; i<n; i++)printf("%d\n",ans[i]);
    }
    return 0;
}