content/graph/virtualTree.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22

// needs dfs in- and out- time and lca function
vector<int> in, out;

void virtualTree(vector<int> ind) { // indices of used nodes
	sort(all(ind), [&](int x, int y) { return in[x] < in[y]; });
	for (int i = 1, n = sz(ind); i < n; i++) {
		ind.push_back(lca(ind[i - 1], ind[i]));
	}
	sort(all(ind), [&](int x, int y) { return in[x] < in[y]; });
	ind.erase(unique(all(ind)), ind.end());

	int n = sz(ind);
	vector<vector<int>> tree(n);
	vector<int> st = {0};
	for (int i = 1; i < n; i++) {
		while (in[ind[i]] >= out[ind[st.back()]]) st.pop_back();
		tree[st.back()].push_back(i);
		st.push_back(i);
	}
	// virtual directed tree with n nodes, original indices in ind
	// weights can be calculated, e.g. with binary lifting
}