【例题记录】树链剖分·线段树·边权点权化

这篇博客主要讨论了一种具有n个节点的有根树,每条边都有初始权值为0,并进行m次操作,每次操作选择两个节点u和v,使路径上所有边权值加1。文章通过构建 segtree 数据结构实现了快速更新和查询边权值的功能,对于每次操作后,可以高效地计算出每条边的权值。
摘要由CSDN通过智能技术生成

有一棵以 1 为根节点,具有 n 个节点的有根树,树上的每条边都具有边权,初始时所有边权均为 0。共有 m 次操作,每次操作将选定两个不同的节点 u,v,使从 u 到 v 的简单路径上的所有边权值均 +1。求每次操作结束后,每条边的权值。

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5 + 5;
int top[N], son[N], fa[N], siz[N], dep[N], dfn[N], rnk[N], dfncnt, tmp[N], w[N];
int head[N], cnt, n, m;
struct node {
	int v, ne;
}e[N<<1];
void add(int u, int v) {
	e[++cnt].v = v, e[cnt].ne = head[u], head[u] = cnt;
}
struct e_node {
	int u, v;
}id[N];
struct segtree {
	int sum[N << 2], add[N << 2];
	void build(int root, int left, int right) {
		if (left == right) {
			sum[root] = 0; add[root] = 0; return;
		}
		int mid = left + right >> 1;
		build(root << 1, left, mid);
		build(root << 1 | 1, mid + 1, right);
		sum[root] = sum[root << 1] + sum[root << 1 | 1];
	}
	void pushdown(int root, int left, int right) {
		if (!add[root]) return;
		int mid = left + right >> 1;
		sum[root << 1] += (mid - left + 1) * add[root];
		add[root << 1] += add[root];
		sum[root << 1 | 1] += (right - mid) * add[root];
		add[root << 1 | 1] += add[root];
		add[root] = 0;
	}
	void update(int root, int left, int right, int l, int r, int val) {
		if (l <= left && right <= r) {
			sum[root] += (right - left + 1) * val;
			add[root] += val; return;
		}
		pushdown(root, left, right);
		int mid = left + right >> 1;
		if (l <= mid) update(root << 1, left, mid, l, r, val);
		if (r > mid) update(root << 1 | 1, mid + 1, right, l, r, val);
		sum[root] = sum[root << 1] + sum[root << 1 | 1];
	}
	int query(int root, int left, int right, int id) {
		if (left == right && id == left) {
			return sum[root];
		}
		pushdown(root, left, right);
		int mid = left + right >> 1;
		if (id <= mid) return query(root<<1, left, mid, id);
		else return query(root<<1|1, mid + 1, right, id);
	}
}tr;
void dfs1(int now) {
	siz[now] = 1;
	for (int i = head[now]; i; i = e[i].ne) {
		if (dep[e[i].v]) continue;
		dep[e[i].v] = dep[now] + 1;
		fa[e[i].v] = now;
		dfs1(e[i].v);
		tmp[e[i].v] = 0; //用子节点映射边权值
		siz[now] += siz[e[i].v];
		if (!son[now] || siz[son[now]] < siz[e[i].v]) son[now] = e[i].v;
	}
}
void dfs2(int now, int tp) {
	top[now] = tp;
	dfn[now] = ++dfncnt;
	w[dfncnt] = tmp[now]; //映射当前dfn对应的原标号的对应边的值
	rnk[dfncnt] = now; //映射当前dfn对应的原标号
	if (son[now]) dfs2(son[now], tp);
	for (int i = head[now]; i; i = e[i].ne) {
		if (e[i].v == son[now] || e[i].v == fa[now]) continue;
		dfs2(e[i].v, e[i].v);
	}
}
void update(int x, int y) {
	int fx = top[x], fy = top[y];
	while (fx != fy) {
		if (dep[fx] >= dep[fy]) tr.update(1, 1, n, dfn[fx], dfn[x], 1), x = fa[fx];
		else tr.update(1, 1, n, dfn[fy], dfn[y], 1), y = fa[fy];
		fx = top[x], fy = top[y];
	}
	if (dep[x] > dep[y]) swap(x, y);
	if (x != y) tr.update(1, 1, n, dfn[x] + 1, dfn[y], 1);
}
int main() {
	cin >> n;
	for (int i = 1; i <= n - 1; i++) {
		int u, v; cin >> u >> v;
		add(u, v), add(v, u);
		id[i].u = u, id[i].v = v;
	}
	dep[1] = 1, dfs1(1), dfs2(1, 1);
	tr.build(1, 1, n);
	cin >> m;
	while (m--) {
		int u, v; cin >> u >> v;
		update(u, v);
		for (int i = 1; i <= n - 1; i++) {
			int po = dep[id[i].u] > dep[id[i].v] ? id[i].u : id[i].v;
			cout << tr.query(1, 1, n, dfn[po]) << " " << endl;
		}
		cout << endl;
	}
	return 0;
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Ctrl AC

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值