边权

P3038 [USACO11DEC] Grass Planting G
树链剖分解决的是点权问题,那么遇到边权问题该怎么办?
其实只要把边权问题改为点权问题就行了,这不是废话

首先我们考虑一条边对应两个点,因为父节点可能对应多个儿子,但儿子节点只有一位父亲(如果是颗树的话),所以我们考虑让儿子节点对应边权。
但这样结果会wa

原因是lca(最近公共祖先)的点也被算了进来,那我们该这么办。
我们考虑如何去掉lca点,根据树链剖分模板,我们知道,在最后一步深度浅的那个点就是lca点
因此只要在最后一步入手即可

我们又知道,最后一步必然在同一条重链上,并且在线段树上是一段连续的区间,因此lca+1就是其儿子节点,这样就成功避免了计算lca

因此只要在最后一步深度浅的那个点的编号+1,就完美解决
是不是很简单呢

点击查看代码
#include<bits/stdc++.h>
using namespace std;
const int N = 1e5 + 10;
int n, m;
vector<int> g[N];
int dep[N], siz[N], son[N], fa[N];
void dfs1(int u, int father) {
	fa[u] = father; dep[u] = dep[father] + 1;
	siz[u] = 1;
	for (auto v : g[u]) {
		if (v == father) continue;
		dfs1(v, u);
		siz[u] += siz[v];
		if (siz[son[u]] < siz[v]) son[u] = v;
	}
}
int top[N], id[N], nw[N], cnt;
void dfs2(int u, int t) {
	top[u] = t; id[u] = ++cnt;
	if (!son[u]) return;
	dfs2(son[u], t);
	for (auto v : g[u]) {
		if (v == fa[u] || v == son[u]) continue;
		dfs2(v, v);
	}
}

struct tree {
	int l, r;
	int w, lz;
}t[N*4];
#define lc k<<1
#define rc k<<1|1
void updata(int k) {
	t[k].w = t[lc].w + t[rc].w;
}
void build(int k, int l, int r) {
	t[k] = { l,r };
	if (l == r) return;
	int mid = l + r >> 1;
	build(lc, l, mid); build(rc, mid + 1, r);
	updata(k);
}
void updown(int k) {
	if (t[k].lz) {
		t[lc].w += (t[lc].r - t[lc].l + 1) * t[k].lz;
		t[rc].w += (t[rc].r - t[rc].l + 1) * t[k].lz;
		t[lc].lz += t[k].lz;
		t[rc].lz += t[k].lz;
		t[k].lz = 0;
	}
}
int query(int k, int l, int r) {
	if (t[k].l >= l && t[k].r <= r) return t[k].w;
	updown(k);
	int mid = t[k].l + t[k].r >> 1;
	int ans = 0;
	if (l <= mid) ans = query(lc, l, r);
	if (r > mid) ans += query(rc, l, r);
	return ans;
}
int query_path(int u, int v) {
	int res = 0;
	while (top[v] != top[u]) {
		if (dep[top[u]] < dep[top[v]]) swap(u, v);
		res += query(1, id[top[u]], id[u]);
		u =fa[top[u]];
	}
	if (dep[u] < dep[v]) swap(v, u);
	res += query(1, id[v]+1, id[u]);//避免多加
	return res;
}
void modify(int k, int l, int r, int z) {
	if (t[k].l >= l && t[k].r <= r) {
		t[k].w += (t[k].r - t[k].l + 1) * z;
		t[k].lz += z;
		return;
	}
	updown(k);
	int mid = t[k].l + t[k].r >> 1;
	if (l <= mid) modify(lc, l, r, z);
	if (r > mid) modify(rc, l, r, z);
	updata(k);
}
void modify_path(int u, int v,int z) {
	while (top[v] != top[u]) {
		if (dep[top[u]] < dep[top[v]]) swap(u, v);
		modify(1, id[top[u]], id[u],z);
		u = fa[top[u]];
	}
	if (dep[u] < dep[v]) swap(v, u);
	modify(1, id[v]+1, id[u],z);//避免多改
}
int main() {
	ios::sync_with_stdio(0), cin.tie(0),cout.tie(0);
	cin >> n >> m;
	for (int i = 1; i < n; i++) {
		int x, y;
		cin >> x >> y;
		g[x].push_back(y);
		g[y].push_back(x);
	}
	dfs1(1,0);
	dfs2(1, 1);
	build(1, 1, n);
	while (m--) {
		char c;
		cin >> c;
		int u, v;
		cin >> u >> v;
		if (c == 'P') {
			modify_path(u, v,1);
		}
		else {
			cout << query_path(u, v)<<'\n';
		}
	}
	return 0;
}
posted @ 2023-09-25 10:30  不o凡  阅读(42)  评论(0)    收藏  举报