BZOJ3262 陌上花开

以前的树套树第一题，好像也是除了二维线段树和二维bit的唯一一道树套树。

三维偏序是很裸的cdq，注意排序时，一是cdq里面sort的时候不能说按x为第一关键字就只比较第一关键字，这可能导致y1处的修改在y2的询问后面(y1<=y2)，第二个是用重载的 < 比cmp快很多很多，把cdq里面的cmp去了之后快了1s。

//QWsin 
#include<cmath> 
#include<cstdio> 
#include<cstring> 
#include<iostream> 
#include<algorithm> 
using namespace std; 
const int maxn=100000+10; 
const int maxk=200000+10; 

inline int read() 
{ 
    int ret=0;char ch=getchar();     
    while(ch<'0'||ch>'9') ch=getchar(); 
    for(;ch>='0'&&ch<='9';ch=getchar()) ret=ret*10+ch-'0'; 
    return ret; 
} 

int n,k; 
struct Node{ 
    int x,y,z,t,id; 
    inline void input(int i){id=i;x=read();y=read();z=read();}   
    bool operator != (const Node &rhs)const{ 
        return  x!=rhs.x||y!=rhs.y||z!=rhs.z; 
    } 
    inline bool operator < (const Node &rhs)const{ 
        if(y!=rhs.y) return y<rhs.y; 
        if(x!=rhs.x) return x<rhs.x; 
        return z<rhs.z; 
    } 
}x[maxn],t[maxn]; 

int ans[maxn],cnt[maxn],num[maxn]; 

int cmp1(Node a,Node b){ 
    if(a.x!=b.x) return a.x<b.x; 
    if(a.y!=b.y) return a.y<b.y; 
    return a.z<b.z; 
} 
int id[maxn]; 

int C[maxk]; 

#define lowbit(i) ((i)&-(i)) 
inline void updata(int pos,int val){ 
    for(int i=pos;i<=k;i+=lowbit(i)) C[i]+=val;   
} 
inline int query(int pos){ 
    int ret=0; 
    for(int i=pos;i;i-=lowbit(i)) ret+=C[i]; 
    return ret;  
} 

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;i<=mid;++i) x[i].t=0; 
    for(int i=mid+1;i<=r;++i) x[i].t=1; 

    for(int i=l;i<=r;++i) t[i]=x[i]; 

    sort(t+l,t+r+1); 

    for(int i=l;i<=r;++i) 
        if(!t[i].t) updata(t[i].z,num[t[i].id]); 
        else ans[t[i].id]+=query(t[i].z); 

    for(int i=l;i<=r;++i)  
        if(!x[i].t) updata(x[i].z,-num[x[i].id]); 
} 

char s[20];int p; 
inline void output(int x){ 
    if(!x) {putchar('0');} 
    else{ 
        for(p=0;x;x/=10) s[++p]=x%10+'0'; 
        for(int i=p;i>=1;--i) putchar(s[i]); 
    } 
    putchar('\n'); 
} 

int main() 
{ 
    cin>>n>>k; 
    for(int i=1;i<=n;++i) t[i].input(i); 

    sort(t+1,t+n+1,cmp1); 

    int p=0; 
    for(int i=1;i<=n;++i)  
    { 
        if(i==1||t[i]!=t[i-1]) x[++p]=t[i],num[p]=1,x[p].id=p; 
        else ++num[p]; 
    } 
    solve(1,p); 

    for(int i=1;i<=n;++i) cnt[ans[i]+num[i]-1]+=num[i]; 

    for(int i=0;i<n;++i) output(cnt[i]); 
    return 0; 
}