結果
問題 |
No.2337 Equidistant
|
ユーザー |
![]() |
提出日時 | 2023-06-02 22:38:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 313 ms / 4,000 ms |
コード長 | 2,372 bytes |
コンパイル時間 | 2,062 ms |
コンパイル使用メモリ | 199,232 KB |
最終ジャッジ日時 | 2025-02-13 19:33:19 |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 28 |
ソースコード
#include <bits/stdc++.h> using ll = long long; using namespace std; struct Graph { static constexpr int N = 18; // n <= 2^N struct Vertex { int n, d, k, p[N]; }; struct Edge { int i, n; }; Graph(int n, int m) : v(n, { -1, -1, -1, {} }), e(m), n(n), m(0) {} void add_edge(int i, int j) { e[m] = { j, v[i].n }; v[i].n = m; m++; } int dfs(int i, int p, int d) { v[i].d = d; v[i].p[0] = p; int k = 1; for (int j = v[i].n; j >= 0; j = e[j].n) { Edge &o = e[j]; if (o.i == p) continue; k += dfs(o.i, i, d + 1); } return v[i].k = k; } void doubling() { dfs(0, 0, 0); for (int k = 1; k < N; k++) { for (int i = 0; i < n; i++) { v[i].p[k] = v[v[i].p[k - 1]].p[k - 1]; } } } int f() { int s, t; cin >> s >> t; s--; t--; int d = v[s].d - v[t].d; if (d < 0) d *= -1, swap(s, t); int i = s, j = t; if (d > 0) { for (int k = N - 1; k >= 0; k--) { if (d & 1 << k) i = v[i].p[k]; } } if (i != j) { for (int k = N - 1; k >= 0; k--) { if (v[i].p[k] == v[j].p[k]) continue; i = v[i].p[k]; j = v[j].p[k]; } i = v[i].p[0]; } int l = v[s].d + v[t].d - v[i].d * 2; if (l % 2 != 0) return 0; l /= 2; l--; if (l > 0) { for (int k = N - 1; k >= 0; k--) { if (l & 1 << k) s = v[s].p[k]; } } int u = v[s].p[0]; int k = v[u].k - v[s].k; if (d == 0) { for (int k = N - 1; k >= 0; k--) { if (l & 1 << k) t = v[t].p[k]; } return n - v[s].k - v[t].k; } return k; } vector<Vertex> v; vector<Edge> e; int n, m; }; int main() { ios::sync_with_stdio(false); cin.tie(0); int n, q; cin >> n >> q; Graph g(n, (n - 1) * 2); for (int i = 0; i < n - 1; i++) { int a, b; cin >> a >> b; a--; b--; g.add_edge(a, b); g.add_edge(b, a); } g.doubling(); while (q--) { cout << g.f() << '\n'; } return 0; }