[20250226] BOJ / P2 / 트리의 외심 / 권혁준

khj20006/202502/26 BOJ P2 점프.md

@@ -0,0 +1,62 @@
+```cpp
+
+#include <iostream>
+#include <vector>
+#include <algorithm>
+#include <cmath>
+using namespace std;
+using ll = long long;
+
+int D[32]{1};
+
+int sol(int d, int s, int e);
+int ans(int x, int y);
+
+int main()
+{
+	cin.tie(0)->sync_with_stdio(0);
+
+	for (int i = 1; i < 32; i++) D[i] = D[i - 1] + i;
+	int T, x, y;
+	for (cin >> T; T--;) {
+		cin >> x >> y;
+		cout << ans(x, y) << '\n';
+	}
+
+}
+
+int J(int x) {
+	if (x <= 1) return x;
+	int lg = log2(x);
+	if (x == (1 << (lg + 1)) - 1) return lg + 1;
+	if (x == (1 << (lg + 1)) - 2) return lg * 2;
+	return lg + J(x - (1 << lg) + 1);
+}
+
+int sol(int d, int s, int e) {
+	if (s == e) return J(s + (1 << d) - 1);
+	if (e == (1 << d)) {
+		int res = max(d * 2, d + 1);
+		if (s == (1 << d) - 1) return res;
+		return max(res, d + ans(s, e - 2));
+	}
+	return d + ans(s, e);
+}
+
+int ans(int x, int y) {
+	if (y <= 5) {
+		int arr[6] = { 0,1,2,2,3,4 };
+		int mx = 0;
+		for (int i = x; i <= y; i++) mx = max(mx, arr[i]);
+		return mx;
+	}
+	int depX = log2(x);
+	int depY = log2(y);
+	int idxX = x - (1 << depX) + 1;
+	int idxY = y - (1 << depY) + 1;
+	if (depY - depX > 1) return max(D[depY - 1], sol(depY, 1, idxY));
+	if (depY - depX == 1) return max(sol(depY, 1, idxY), sol(depX, idxX, (1 << depX)));
+	return sol(depX, idxX, idxY);
+}
+
+```

khj20006/202502/26 BOJ P2 트리의 외심.md

@@ -0,0 +1,99 @@
+```cpp
+
+#include <iostream>
+#include <vector>
+#include <algorithm>
+using namespace std;
+using ll = long long;
+
+int N, Q;
+vector<vector<int>> V(100001);
+int par[100001][17]{}, dep[100001]{};
+
+void dfs(int n, int p, int d) {
+	par[n][0] = p, dep[n] = d;
+	for (int i : V[n]) if (i != p) dfs(i, n, d + 1);
+}
+
+void makeTable() {
+	for (int k = 1; k < 17; k++) for (int i = 1; i <= N; i++) par[i][k] = par[par[i][k - 1]][k - 1];
+}
+
+pair<int, int> solve(int a, int b) {
+	int diff = abs(dep[a] - dep[b]);
+	int ma = 0, mb = 0;
+	int oa = a, ob = b;
+	for (int i = 0; i < 17; i++) if (diff & (1 << i)) {
+		if (dep[a] > dep[b]) a = par[a][i], ma += (1 << i);
+		else b = par[b][i], mb += (1 << i);
+	}
+
+	while (par[a][0] != par[b][0]) {
+		for (int i = 16; i >= 0; i--) if (par[a][i] != par[b][i]) {
+			a = par[a][i], ma += (1 << i);
+			b = par[b][i], mb += (1 << i);
+			break;
+		}
+	}
+	if (a != b) {
+		a = par[a][0], ma++;
+		b = par[b][0], mb++;
+	}
+
+	int d = ma + mb;
+	if (ma == mb) return { a, d };
+	if ((ma + mb) % 2) return { -1,d };
+
+	int from = ma > mb ? oa : ob, to = a, len = (ma + mb) / 2;
+	for (int i = 0; i < 17; i++) if (len & (1 << i)) from = par[from][i];
+	return { from, d };
+
+}
+
+int main()
+{
+	cin.tie(0)->sync_with_stdio(0);
+
+	cin >> N;
+	for (int i = 1, a, b; i < N; i++) {
+		cin >> a >> b;
+		V[a].push_back(b);
+		V[b].push_back(a);
+	}
+
+	dfs(1, 0, 0);
+	makeTable();
+
+	for (cin >> Q; Q--;) {
+		int a, b, c;
+		cin >> a >> b >> c;
+		int res = -1;
+
+		auto[m1, _x] = solve(a, b);
+		if (m1 != -1) {
+			auto[_1, l1] = solve(a, m1);
+			auto[_2, l2] = solve(b, m1);
+			auto[_3, l3] = solve(c, m1);
+			if (l1 == l2 && l2 == l3) res = m1;
+		}
+		auto[m2, _y] = solve(b, c);
+		if (m2 != -1) {
+			auto[_1, l1] = solve(a, m2);
+			auto[_2, l2] = solve(b, m2);
+			auto[_3, l3] = solve(c, m2);
+			if (l1 == l2 && l2 == l3) res = m2;
+		}
+		auto[m3, _z] = solve(a, c);
+		if (m3 != -1) {
+			auto[_1, l1] = solve(a, m3);
+			auto[_2, l2] = solve(b, m3);
+			auto[_3, l3] = solve(c, m3);
+			if (l1 == l2 && l2 == l3) res = m3;
+		}
+
+		cout << res << '\n';
+	}
+
+}
+
+```