This documentation is automatically generated by online-judge-tools/verification-helper
View the Project on GitHub plasmatic1/templates
#include "tree/lca_binary_lift.hpp"
#pragma once #include "../template.hpp" template <typename GraphType, typename EdgeType> struct LCABinaryLift { EdgeType E; int lg; vector<int> dep; vector<vector<int>> tb; void init(int N, GraphType &g, int rt) { lg = 31 - __builtin_clz(N) + 1; tb.assign(lg, vector<int>(N + 1, -1)); dep.assign(N + 1, 0); function<void(int, int)> dfs = [&] (int c, int p) { tb[0][c] = p; for (auto to : g[c]) { int v = E.v(to); if (v != p) { dep[v] = dep[c] + 1; dfs(v, c); } } }; dfs(rt, -1); for (auto i = 1; i < lg; i++) { for (auto j = 1; j <= N; j++) { int pp = tb[i - 1][j]; tb[i][j] = pp == -1 ? -1 : tb[i - 1][pp]; } } } int lca(int a, int b) { if (a == b) return a; if (dep[a] > dep[b]) swap(a, b); int delta = dep[b] - dep[a]; for (auto i = 0; i < lg; i++) if ((delta >> i) & 1) b = tb[i][b]; if (a == b) return a; for (auto i = lg-1; i >= 0; i--) { if (tb[i][a] != tb[i][b]) { a = tb[i][a]; b = tb[i][b]; } } return tb[0][a]; } #if __cplusplus == 201703L // CPP17 only things void bind(opt_ref<vector<int>> dep0, opt_ref<vector<vector<int>>> tb0) { if (dep0) dep.swap(*dep0); if (tb0) tb.swap(*tb0); } #endif };
#line 2 "template.hpp" #include <bits/stdc++.h> #define DEBUG 1 using namespace std; // Defines #define fs first #define sn second #define pb push_back #define eb emplace_back #define mpr make_pair #define mtp make_tuple #define all(x) (x).begin(), (x).end() // Basic type definitions #if __cplusplus == 201703L // CPP17 only things template <typename T> using opt_ref = optional<reference_wrapper<T>>; // for some templates #endif using ll = long long; using ull = unsigned long long; using ld = long double; using pii = pair<int, int>; using pll = pair<long long, long long>; #ifdef __GNUG__ // PBDS order statistic tree #include <ext/pb_ds/assoc_container.hpp> // Common file #include <ext/pb_ds/tree_policy.hpp> using namespace __gnu_pbds; template <typename T, class comp = less<T>> using os_tree = tree<T, null_type, comp, rb_tree_tag, tree_order_statistics_node_update>; template <typename K, typename V, class comp = less<K>> using treemap = tree<K, V, comp, rb_tree_tag, tree_order_statistics_node_update>; // HashSet #include <ext/pb_ds/assoc_container.hpp> template <typename T, class Hash> using hashset = gp_hash_table<T, null_type, Hash>; template <typename K, typename V, class Hash> using hashmap = gp_hash_table<K, V, Hash>; const ll RANDOM = chrono::high_resolution_clock::now().time_since_epoch().count(); struct chash { ll operator()(ll x) const { return x ^ RANDOM; } }; #endif // More utilities int SZ(string &v) { return v.length(); } template <typename C> int SZ(C &v) { return v.size(); } template <typename C> void UNIQUE(vector<C> &v) { sort(v.begin(), v.end()); v.resize(unique(v.begin(), v.end()) - v.begin()); } template <typename T, typename U> void maxa(T &a, U b) { a = max(a, b); } template <typename T, typename U> void mina(T &a, U b) { a = min(a, b); } const ll INF = 0x3f3f3f3f, LLINF = 0x3f3f3f3f3f3f3f3f; #line 3 "tree/lca_binary_lift.hpp" template <typename GraphType, typename EdgeType> struct LCABinaryLift { EdgeType E; int lg; vector<int> dep; vector<vector<int>> tb; void init(int N, GraphType &g, int rt) { lg = 31 - __builtin_clz(N) + 1; tb.assign(lg, vector<int>(N + 1, -1)); dep.assign(N + 1, 0); function<void(int, int)> dfs = [&] (int c, int p) { tb[0][c] = p; for (auto to : g[c]) { int v = E.v(to); if (v != p) { dep[v] = dep[c] + 1; dfs(v, c); } } }; dfs(rt, -1); for (auto i = 1; i < lg; i++) { for (auto j = 1; j <= N; j++) { int pp = tb[i - 1][j]; tb[i][j] = pp == -1 ? -1 : tb[i - 1][pp]; } } } int lca(int a, int b) { if (a == b) return a; if (dep[a] > dep[b]) swap(a, b); int delta = dep[b] - dep[a]; for (auto i = 0; i < lg; i++) if ((delta >> i) & 1) b = tb[i][b]; if (a == b) return a; for (auto i = lg-1; i >= 0; i--) { if (tb[i][a] != tb[i][b]) { a = tb[i][a]; b = tb[i][b]; } } return tb[0][a]; } #if __cplusplus == 201703L // CPP17 only things void bind(opt_ref<vector<int>> dep0, opt_ref<vector<vector<int>>> tb0) { if (dep0) dep.swap(*dep0); if (tb0) tb.swap(*tb0); } #endif };