BZOJ_4756_[Usaco2017 Jan]Promotion Counting_树状数组

BZOJ_4756_[Usaco2017 Jan]Promotion Counting_树状数组

Description

n只奶牛构成了一个树形的公司,每个奶牛有一个能力值pi,1号奶牛为树根。
问对于每个奶牛来说,它的子树中有几个能力值比它大的。

Input

n,表示有几只奶牛 n<=100000
接下来n行为1-n号奶牛的能力值pi
接下来n-1行为2-n号奶牛的经理(树中的父亲)

Output

共n行,每行输出奶牛i的下属中有几个能力值比i大

Sample Input

5
804289384
846930887
681692778
714636916
957747794
1
1
2
3

Sample Output

2
0
1
0
0
 

这道题的思路比较巧妙。
我们对整棵树进行DFS时每个点有两次处理的机会。
一次是刚刚遍历到,一次是子树的回溯。
分别求出两次能力值大于这个点的点的个数,用第二次的减去第一次的就是这个点的答案。
然后离散用树状数组处理一下。
 
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 100050
struct A {
    int num,id,v;
}a[N];
bool cmp1(const A &x,const A &y){return x.num>y.num; }
bool cmp2(const A &x,const A &y){return x.id<y.id; }
int n,head[N],to[N<<1],nxt[N<<1],cnt,ans[N],c[N];
inline void add(int u,int v) {
    to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;
}
void fix(int x,int v) {
    for(;x<=n;x+=x&-x) c[x]+=v;
}
int inq(int x) {
    int re=0;
    for(;x;x-=x&-x) re+=c[x];
    return re;
}
void dfs(int x,int y) {
    int i;
    fix(a[x].v,1);
    int tmp=inq(a[x].v-1);
    for(i=head[x];i;i=nxt[i]) {
        if(to[i]!=y) {
            dfs(to[i],x);
        }
    }
    ans[x]=inq(a[x].v-1)-tmp;
}
int main() {
    //freopen("tt.in","r",stdin);
    scanf("%d",&n);
    int i,x,j;
    for(i=1;i<=n;i++) scanf("%d",&a[i].num),a[i].id=i;
    sort(a+1,a+n+1,cmp1); a[0].num=-245345;
    for(i=1,j=0;i<=n;i++) { if(a[i].num!=a[i-1].num) j++; a[i].v=j; }
    sort(a+1,a+n+1,cmp2);
    for(i=2;i<=n;i++) {
        scanf("%d",&x);
        add(i,x); add(x,i);
    }
    dfs(1,0);
    for(i=1;i<=n;i++) printf("%d\n",ans[i]);
}

 

转载于:https://www.cnblogs.com/suika/p/8734488.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值