第五届CCPC中国大学生程序设计竞赛河南省赛-C

 

思路:离散化+主席树

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;
#define N 300003
#define mid (l+r>>1)
#define lc (tr[d].ld)
#define rc (tr[d].rd)
typedef long long ll;

int mx, k, val[N], bit[N], cnt;
int ro[N/3], tot, tid;
ll ans[N/3];
vector<int>V[N/3];

struct Tr{
	int ld, rd;
    int	c;
}tr[N*40];

int newTr() {
	tr[++tot].c = 0;
	tr[tot].ld = tr[tot].rd = -1;
	return tot;
}

void push(int d) {
	tr[d].c = tr[lc].c+tr[rc].c;
}

int build(int l, int r) {
	int d = newTr();
	if (l != r) {
		tr[d].ld = build(l, mid);
		tr[d].rd = build(mid+1, r);
		push(d);
	}
	return d;
}

int update(int l, int r, int pos, int ad, int pre) {
	int d = newTr();
	tr[d].c = tr[pre].c+ad;
	if (l != r) {
		if (pos <= mid) {
			tr[d].ld = update(l, mid, pos, ad, tr[pre].ld);
			tr[d].rd = tr[pre].rd;
		} else {
			tr[d].rd = update(mid+1, r, pos, ad, tr[pre].rd);
			tr[d].ld = tr[pre].ld;
		}
		push(d);
	}
	return d;
}

int query(int l, int r, int L, int R, int rt, int lt) {
	if (l == L && r == R) {
		return tr[rt].c-tr[lt].c;
	}
	if (R <= mid) return query(l, mid, L, R, tr[rt].ld, tr[lt].ld);
	else if (L > mid) return query(mid+1, r, L, R, tr[rt].rd, tr[lt].rd);
	else return query(l, mid, L, mid, tr[rt].ld, tr[lt].ld)+query(mid+1, r, mid+1, R, tr[rt].rd, tr[lt].rd);
}

ll dfs(int u) {
	int rk = lower_bound(bit+1, bit+cnt, val[u])-bit;
	int st = ++tid;
	ro[st] = update(1, cnt, rk, 1, ro[st-1]);
	int i, v;
	ans[u] = 0;
	for (i = 0;i < V[u].size();i++) {
		v = V[u][i];
		ans[u] += dfs(v);
	}
	int l, r;
	l = lower_bound(bit+1, bit+cnt, max(0,val[u]-k))-bit;
	r = lower_bound(bit+1, bit+cnt, min(mx, val[u]+k))-bit;
	ans[u] += query(1, cnt, l, r, ro[tid], ro[st-1]);
	return ans[u];
}

int main() {
	int i, j, n;
	while(~scanf("%d%d", &n, &k)) {
		cnt = 0;
		mx = 0;
		for (i = 1;i <= n;i++) {
			scanf("%d", val+i);
			mx = max(val[i], mx);
			V[i].clear();
		}
		for (i = 1;i <= n;i++) {
			bit[++cnt] = val[i];
			bit[++cnt] = max(0, val[i]-k);
			bit[++cnt] = min(mx, val[i]+k);
		}
		for (i = 2;i <= n;i++) {
			scanf("%d", &j);
			V[j].push_back(i);
		}
		sort(bit+1, bit+1+cnt);
		cnt = unique(bit+1, bit+1+cnt) - bit - 1;
		tot = tid = 0;
		ro[0] =	build(1, cnt);
		dfs(1);
		for (i = 1;i <= n;i++) printf("%lld\n", ans[i]);
	}
}

 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值