三维偏序 模板

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

ll read(){
    ll x=0,f=1;char c=getchar();
    while(c<'0' || c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0' && c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

struct Node{
    int a,b,c;
    int ans,cnt;
    bool operator < (const Node &ano) const{
        if(a!=ano.a) return a<ano.a;
        else if(b!=ano.b) return b<ano.b;
        else return c<ano.c;
    }
    bool operator <=(const Node &ano) const{
        if(a<=ano.a && b<=ano.b && c<=ano.c) return 1;
    }
    bool operator ==(const Node &ano) const{
        return a==ano.a && b==ano.b && c==ano.c;
    }
} a[100100],put[100100];
int n,k;
int ans[100100];

#define lowbit(x) (x&-x)
struct Bit{
    int a[200200];

    void add(int pos,int val){
        //cout<<"add "<<pos<<' '<<val<<endl;
        while(pos<=200000){
            a[pos]+=val;
            pos+=lowbit(pos);
        }
    }

    int ask(int pos,int ret=0){
        //cout<<"ask "<<pos<<' ';
        while(pos){
            ret+=a[pos];
            pos-=lowbit(pos);
        }
        //cout<<ret<<endl;
        return ret;
    }
} BIT;

void solve(int l,int r){
    //cout<<l<<' '<<r<<endl;
    if(l>=r) return;
    int md=(l+r)>>1;
    solve(l,md),solve(md+1,r);
    int pos=l,top=l;
    for(int i=md+1;i<=r;i++){
        while(a[i].b>=a[pos].b && pos<=md){BIT.add(a[pos].c,a[pos].cnt);put[top++]=a[pos++];}
        a[i].ans+=BIT.ask(a[i].c);put[top++]=a[i];
    }
    for(int i=l;i<pos;i++)
        BIT.add(a[i].c,-a[i].cnt);
    while(pos<=md) put[top++]=a[pos++];
    for(int i=l;i<=r;i++)
        a[i]=put[i];
}

int main(){
#ifdef LZT
    freopen("in","r",stdin);
#endif
    n=read(),k=read();
    for(int i=1;i<=n;i++)
        a[i].a=read(),a[i].b=read(),a[i].c=read(),a[i].cnt=1;
    sort(a+1,a+n+1);
    int cnt=1;
    for(int i=2;i<=n;i++){
        if(a[cnt]==a[i]) a[cnt].cnt++; else a[++cnt]=a[i];
    }
    //for(int i=1;i<=cnt;i++)
    //    cout<<a[i].a<<' '<<a[i].b<<' '<<a[i].c<<' '<<a[i].cnt<<endl;
    solve(1,cnt);
    for(int i=1;i<=cnt;i++)
        ans[a[i].ans+a[i].cnt-1]+=a[i].cnt;
    for(int i=0;i<n;i++)
        printf("%d
",ans[i]);
    return 0;
}