#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pll;
const int inf = 1e9;
struct SegmentTree{
vector<int> tree;
void init(int n) {
int sz = 1 << __lg(n-1) + 2;
tree.assign(sz, inf);
}
void update(int node, int s, int e, int l, int r, int val) {
if (l > e || s > r) return;
if (l <= s && e <= r) {
tree[node] = min(tree[node], val);
return;
}
int mid = s + e >> 1;
update(node*2, s, mid, l, r, val);
update(node*2+1, mid+1, e, l, r, val);
}
int query(int node, int s, int e, int tar) {
if (s == e) return tree[node];
int mid = s + e >> 1;
if (tar <= mid) return min(tree[node], query(node*2, s, mid, tar));
return min(tree[node], query(node*2+1, mid+1, e, tar));
}
};
struct Node{
int l, r, val;
};
struct DynamicSegmentTree{
vector<Node> tree;
void init() {
tree.clear();
tree.push_back({-1, -1, inf});
}
void update(int node, int s, int e, int tar, int val, bool f = false) {
tree[node].val = min(tree[node].val, val);
if (s == e) {
return;
}
int mid = s + e >> 1;
if (tar <= mid) {
if (tree[node].l == -1) {
tree[node].l = tree.size();
tree.push_back({-1, -1, inf});
}
update(tree[node].l, s, mid, tar, val, f);
}
else {
if (tree[node].r == -1) {
tree[node].r = tree.size();
tree.push_back({-1, -1, inf});
}
update(tree[node].r, mid+1, e, tar, val, f);
}
}
int query(int node, int s, int e, int l, int r, bool f = false) {
if (l > e || s > r || node == -1) return inf;
if (l <= s && e <= r) return tree[node].val;
int mid = s + e >> 1;
return min(query(tree[node].l, s, mid, l, r, f), query(tree[node].r, mid+1, e, l, r, f));
}
};
SegmentTree seg;
int N, R, A, B;
vector<int> adj[300010];
vector<int> g[300010];
int down[300010], par[300010];
int dep[300010];
int in[300010], out[300010];
int pv;
int sz[300010];
int top[300010];
int chk[300010];
int getSize(int x, int p) {
sz[x] = 1;
for (auto &y : adj[x]) {
if (chk[y] || y == p) continue;
sz[x] += getSize(y, x);
}
return sz[x];
}
int getCent(int x, int p, int cap) {
for (auto &y : adj[x]) {
if (chk[y] || y == p) continue;
if (sz[y] * 2 > cap)
return getCent(y, x, cap);
}
return x;
}
set<int> ps[300010];
int cent[300010][20], dist[300010][20];
int sub[300010][20];
void dfs3(int x, int p, int lv, int S) {
sub[x][lv] = S;
for (auto &y : adj[x]) {
if (y == p || chk[y]) continue;
cent[y][lv] = cent[x][lv];
dist[y][lv] = dist[x][lv] + 1;
dfs3(y, x, lv, S);
}
}
void dnc(int x, int lv) {
x = getCent(x, -1, getSize(x, -1));
int t = x;
while (t && !chk[t]) {
ps[x].insert(t);
t = par[t];
}
sub[x][lv] = x;
cent[x][lv] = x;
dist[x][lv] = 0;
for (auto &y : adj[x]) {
if (chk[y]) continue;
cent[y][lv] = x;
dist[y][lv] = 1;
dfs3(y, x, lv, y);
}
chk[x] = 1;
for (auto &y : adj[x]) {
if (chk[y]) continue;
dnc(y, lv+1);
}
}
void dfs(int x, int p) {
sz[x] = 1;
down[x] = 0;
for (auto &y : adj[x]) {
if (y == p) continue;
par[y] = x;
dep[y] = dep[x] + 1;
g[x].push_back(y);
dfs(y, x);
down[x] = max(down[x], down[y] + 1);
if (sz[g[x][0]] < sz[y]) swap(g[x][0], g[x].back());
sz[x] += sz[y];
}
}
void dfs2(int x) {
in[x] = ++ pv;
for (auto &y : g[x]) {
top[y] = y == g[x][0] ? top[x] : y;
dfs2(y);
}
out[x] = pv;
}
void HLDupdate(int x) {
int ori = x;
int cnt = B;
while (x) {
int tp = top[x];
if (dep[x] - dep[tp] + 1 < cnt) {
seg.update(1, 1, N, in[tp], in[x], ori);
cnt -= dep[x] - dep[tp] + 1;
}
else {
seg.update(1, 1, N, in[x]-cnt+1, in[x], ori);
break;
}
x = par[tp];
}
}
int dp[300010];
int dp2[300010];
vector<int> vtx[300010];
vector<pair<int,int>> poss[300010], imp[300010];
DynamicSegmentTree seg1[300010], seg2[300010];
void CentUpdate(int x) {
for (int lv=19; lv>=0; lv--) {
if (!cent[x][lv]) continue;
if (cent[x][lv] == x) continue;
int c = cent[x][lv];
int val;
if (poss[x].empty()) val = imp[x][0].first;
else val = poss[x][0].first;
if (sub[x][lv] == par[c]) {
if (ps[c].find(x) != ps[c].end()) {
if (poss[x].empty()) {
if (ps[c].find(imp[x][0].second) != ps[c].end())
val = imp[x][1].first;
}
else if (poss[x].size() == 1) {
if (ps[c].find(poss[x][0].second) != ps[c].end())
val = imp[x][0].first;
}
else {
if (ps[c].find(poss[x][0].second) != ps[c].end())
val = poss[x][1].first;
}
}
seg1[c].update(0, 1, N, dist[x][lv], val);
}
else {
seg1[c].update(0, 1, N, dist[x][lv], val);
seg2[c].update(0, 1, N, dist[x][lv], val);
}
}
}
int Centquery(int x) {
int ret = N + 1;
for (int lv=19; lv>=0; lv--) {
if (!cent[x][lv]) continue;
int d = B - dist[x][lv];
if (d < 0) continue;
int c = cent[x][lv];
int val;
if (poss[c].empty()) val = imp[c][0].first;
else val = poss[c][0].first;
if (x == c) {
ret = min(ret, val);
}
else {
if (poss[c].empty()) {
if (sub[x][lv] == imp[c][0].second)
val = imp[c][1].first;
}
else if (poss[c].size() == 1) {
if (sub[x][lv] == poss[c][0].second)
val = imp[c][0].first;
}
else {
if (sub[x][lv] == poss[c][0].second)
val = poss[c][1].first;
}
ret = min(ret, val);
}
if (d == 0) continue;
if (sub[x][lv] == par[c]) {
ret = min(ret, seg2[c].query(0, 1, N, 1, d));
}
else {
ret = min(ret, seg1[c].query(0, 1, N, 1, d));
}
}
assert(ret != N+1);
return ret;
}
//DynamicSegmentTree final[300010];
int simulate(int x, int p, int dist) {
if (dist > B) return N + 1;
int ret;
int d = 0;
for (auto &y : adj[x]) {
if (y == p || y == par[x]) continue;
d = max(d, down[y] + 1);
}
if (d < A) {
ret = x;
for (auto &y : adj[x]) {
if (y == p || y == par[x]) continue;
ret = min(ret, dp[y]);
}
}
else {
ret = 0;
for (auto &y : adj[x]) {
if (y == p || y == par[x]) continue;
if (down[y] >= A-1)
ret = max(ret, dp[y]);
}
}
for (auto &y : adj[x]) {
if (y == p) continue;
ret = min(ret, simulate(y, x, dist + 1));
}
return ret;
}
int main() {
ios_base :: sync_with_stdio(false); cin.tie(NULL);
cin >> N >> R >> A >> B;
if (A <= B) {
cout << "1\n";
return 0;
}
for (int i=0; i<N-1; i++) {
int u, v;
cin >> u >> v;
adj[u].push_back(v);
adj[v].push_back(u);
}
dfs(R, 0);
top[R] = R;
dfs2(R);
seg.init(N);
int mxdep = 0;
for (int i=1; i<=N; i++) {
mxdep = max(mxdep, dep[i]);
vtx[dep[i]].push_back(i);
HLDupdate(i);
imp[i].push_back({i, -1});
seg1[i].init();
seg2[i].init();
}
dnc(1, 0);
for (int i=mxdep; i>=1; i--) {
for (auto &x : vtx[i]) {
sort(imp[x].begin(), imp[x].end());
sort(poss[x].begin(), poss[x].end());
reverse(poss[x].begin(), poss[x].end());
CentUpdate(x);
}
if (i+A-1 <= mxdep) {
for (auto &x : vtx[i+A-1]) {
int ret1 = Centquery(x);
int ret2 = simulate(x, -1, 0);
assert(ret1 >= ret2);
dp2[x] = ret2;
}
}
for (auto &x : vtx[i]) {
if (down[x] < A-1) {
dp[x] = min(par[x], seg.query(1, 1, N, in[x]));
imp[par[x]].push_back({dp[x], x});
}
else {
for (auto &y : vtx[i+A-1]) {
if (in[x] <= in[y] && in[y] <= out[x])
dp[x] = max(dp[x], dp2[y]);
}
poss[par[x]].push_back({dp[x], x});
}
}
}
if (down[R] < A) {
int ans = R;
for (auto &x : vtx[1]) {
ans = min(ans, dp[x]);
}
cout << ans;
}
else {
int ans = 0;
for (auto &x : vtx[1]) {
if (down[x] >= A-1)
ans = max(ans, dp[x]);
}
cout << ans;
}
}
Compilation message
Main.cpp: In member function 'void SegmentTree::init(int)':
Main.cpp:12:33: warning: suggest parentheses around '+' inside '<<' [-Wparentheses]
12 | int sz = 1 << __lg(n-1) + 2;
| ~~~~~~~~~~^~~
Main.cpp: In member function 'void SegmentTree::update(int, int, int, int, int, int)':
Main.cpp:21:21: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
21 | int mid = s + e >> 1;
| ~~^~~
Main.cpp: In member function 'int SegmentTree::query(int, int, int, int)':
Main.cpp:27:21: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
27 | int mid = s + e >> 1;
| ~~^~~
Main.cpp: In member function 'void DynamicSegmentTree::update(int, int, int, int, int, bool)':
Main.cpp:48:21: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
48 | int mid = s + e >> 1;
| ~~^~~
Main.cpp: In member function 'int DynamicSegmentTree::query(int, int, int, int, int, bool)':
Main.cpp:67:21: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
67 | int mid = s + e >> 1;
| ~~^~~
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
13 ms |
67936 KB |
Output is correct |
| 2 |
Correct |
10 ms |
67932 KB |
Output is correct |
| 3 |
Correct |
9 ms |
67932 KB |
Output is correct |
| 4 |
Correct |
10 ms |
67928 KB |
Output is correct |
| 5 |
Correct |
12 ms |
68084 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Runtime error |
1664 ms |
728500 KB |
Execution killed with signal 6 |
| 2 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
12 ms |
80476 KB |
Output is correct |
| 2 |
Correct |
11 ms |
80472 KB |
Output is correct |
| 3 |
Correct |
13 ms |
80392 KB |
Output is correct |
| 4 |
Correct |
13 ms |
80476 KB |
Output is correct |
| 5 |
Correct |
14 ms |
80384 KB |
Output is correct |
| 6 |
Correct |
12 ms |
80344 KB |
Output is correct |
| 7 |
Correct |
15 ms |
80472 KB |
Output is correct |
| 8 |
Correct |
11 ms |
80508 KB |
Output is correct |
| 9 |
Correct |
12 ms |
67932 KB |
Output is correct |
| 10 |
Correct |
14 ms |
80476 KB |
Output is correct |
| 11 |
Runtime error |
62 ms |
163204 KB |
Execution killed with signal 6 |
| 12 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
12 ms |
80476 KB |
Output is correct |
| 2 |
Correct |
11 ms |
80472 KB |
Output is correct |
| 3 |
Correct |
13 ms |
80392 KB |
Output is correct |
| 4 |
Correct |
13 ms |
80476 KB |
Output is correct |
| 5 |
Correct |
14 ms |
80384 KB |
Output is correct |
| 6 |
Correct |
12 ms |
80344 KB |
Output is correct |
| 7 |
Correct |
15 ms |
80472 KB |
Output is correct |
| 8 |
Correct |
11 ms |
80508 KB |
Output is correct |
| 9 |
Correct |
12 ms |
67932 KB |
Output is correct |
| 10 |
Correct |
14 ms |
80476 KB |
Output is correct |
| 11 |
Runtime error |
62 ms |
163204 KB |
Execution killed with signal 6 |
| 12 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
1172 ms |
157264 KB |
Output is correct |
| 2 |
Correct |
974 ms |
163156 KB |
Output is correct |
| 3 |
Correct |
828 ms |
158572 KB |
Output is correct |
| 4 |
Correct |
787 ms |
163196 KB |
Output is correct |
| 5 |
Correct |
1057 ms |
158556 KB |
Output is correct |
| 6 |
Correct |
945 ms |
163152 KB |
Output is correct |
| 7 |
Runtime error |
251 ms |
266228 KB |
Execution killed with signal 6 |
| 8 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
12 ms |
80476 KB |
Output is correct |
| 2 |
Correct |
11 ms |
80472 KB |
Output is correct |
| 3 |
Correct |
13 ms |
80392 KB |
Output is correct |
| 4 |
Correct |
13 ms |
80476 KB |
Output is correct |
| 5 |
Correct |
14 ms |
80384 KB |
Output is correct |
| 6 |
Correct |
12 ms |
80344 KB |
Output is correct |
| 7 |
Correct |
15 ms |
80472 KB |
Output is correct |
| 8 |
Correct |
11 ms |
80508 KB |
Output is correct |
| 9 |
Correct |
12 ms |
67932 KB |
Output is correct |
| 10 |
Correct |
14 ms |
80476 KB |
Output is correct |
| 11 |
Runtime error |
62 ms |
163204 KB |
Execution killed with signal 6 |
| 12 |
Halted |
0 ms |
0 KB |
- |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
13 ms |
67936 KB |
Output is correct |
| 2 |
Correct |
10 ms |
67932 KB |
Output is correct |
| 3 |
Correct |
9 ms |
67932 KB |
Output is correct |
| 4 |
Correct |
10 ms |
67928 KB |
Output is correct |
| 5 |
Correct |
12 ms |
68084 KB |
Output is correct |
| 6 |
Runtime error |
1664 ms |
728500 KB |
Execution killed with signal 6 |
| 7 |
Halted |
0 ms |
0 KB |
- |
