BZOJ 4756 [Usaco2017 Jan]Promotion Counting 線段樹合併

題目大意：給出一棵有根樹，每個點有一個權值，問對於每一個點的子樹中比這個點大的點的個數。

首先權值很大需要離散化。（2017-7-13UPD:不用離散化也可反正也是動態開點
對於每一個點建一個權值線段樹。遍歷每一棵子樹，然後將子樹的線段樹合併到自己的線段樹上，在查詢之後將自己插入。時間複雜度O(nlogn)

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 100005
using namespace std;
struct Node {
    Node *ch[2];
    int l,r,sum;
    Node(int _l,int _r):l(_l),r(_r) {
        sum=0;
        ch[0]=ch[1]=NULL;
    }
    Node() {}
    void maintain() {
        sum=0;
        if(ch[0]) sum+=ch[0]->sum;
        if(ch[1]) sum+=ch[1]->sum;
        return ;
    }
    void* operator new (size_t) {
        static Node *mempool,*C;
        if(mempool==C) mempool=(C=new Node[1<<20])+(1<<20);
        return C++;
    }
}*root[N];
struct Edge {
    int to,nxt;
    Edge() {}
    Edge(int _to,int _nxt):to(_to),nxt(_nxt) {}
}e[N];
int n,m,tot=-1,fir[N],a[N],b[N],ans[N];
void Add_Edge(int x,int y) {
    e[++tot]=Edge(y,fir[x]), fir[x]=tot;
    return ;
}
void Merge(Node*& x,Node*& y) {
    if(!y) return ;
    if(!x) {
        x=y;
        return ;
    }
    x->sum+=y->sum;
    Merge(x->ch[0],y->ch[0]), Merge(x->ch[1],y->ch[1]);
    return ;
}
void Insert(Node*& o,int v,int L,int R) {
    if(!o) o=new Node(L,R);
    if(o->l==o->r) {
        o->sum++;
        return ;
    }
    int mid=L+R>>1;
    if(v<=mid) Insert(o->ch[0],v,L,mid);
    else Insert(o->ch[1],v,mid+1,R);
    o->maintain();
    return ;
}
int Query(Node* o,int v) {
    if(!o) return 0;
    if(v==o->r) return o->sum;
    int mid=o->l+o->r>>1;
    if(v<=mid) return Query(o->ch[0],v);
    return (o->ch[0] ? o->ch[0]->sum : 0)+Query(o->ch[1],v);
}
void dfs(int x) {
    for(int i=fir[x];~i;i=e[i].nxt) {
        dfs(e[i].to);
        Merge(root[x],root[e[i].to]);
    }
    ans[x]=Query(root[x],n-a[x]);
    Insert(root[x],n-a[x]+1,1,n);
    return ;
}
int main() {
    memset(fir,-1,sizeof fir);
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",a+i), b[i]=a[i];
    sort(b+1,b+n+1);
    m=unique(b+1,b+n+1)-b-1;
    for(int i=1;i<=n;i++) a[i]=lower_bound(b+1,b+1+m,a[i])-b;
    for(int i=2;i<=n;i++) {
        int x;
        scanf("%d",&x);
        Add_Edge(x,i);
    }
    dfs(1);
    for(int i=1;i<=n;i++) printf("%d\n",ans[i]);
    return 0;
}