#include <bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define dl double
#define st first
#define nd second
#define II pair <int, int>
#include "factories.h"
using namespace std;
const ll inf = 7 + 1e18;
int n, logn, tsize;
vector <int> sz, pa, h;
vector <ll> f, ans;
vector <vector <int>> p;
vector <vector <II>> adj;
///--------------------------LCA-------------------------------
void dfs_lca(int u) {
for (auto z : adj[u]) {
int v = z.st;
ll w = z.nd;
if (h[v])
continue;
h[v] = h[u] + 1;
f[v] = f[u] + w;
p[0][v] = u;
dfs_lca(v);
}
}
int lca(int u, int v) {
if (h[u] < h[v])
swap(u, v);
int d = h[u] - h[v];
for (int si = d; si > 0; si ^= si & -si) {
int i = __builtin_ctz(si & -si);
u = p[i][u];
}
if (u == v)
return u;
for (int i = logn - 1; i >= 0; i --)
if (p[i][u] != p[i][v]) {
u = p[i][u];
v = p[i][v];
}
return p[0][u];
}
ll dis(int u, int v) {
return f[u] + f[v] - f[lca(u, v)] * 2;
}
///-------------------------CENTROID-----------------------------
void dfs(int u, int p) {
tsize ++;
sz[u] = 1;
for (auto z : adj[u]) {
int v = z.st;
if (v == p || pa[v])
continue;
dfs(v, u);
sz[u] += sz[v];
}
}
int centroid(int u, int p) {
for (auto z : adj[u]) {
int v = z.st;
if (v == p || pa[v])
continue;
if (sz[v] * 2 > tsize)
return centroid(v, u);
}
return u;
}
void build(int u, int p) {
tsize = 0;
dfs(u, 0);
int c = centroid(u, 0);
if (p == 0)
pa[c] = c;
else
pa[c] = p;
for (auto z : adj[c]) {
int v = z.st;
if (pa[v])
continue;
build(v, c);
}
}
void Init(int N, int A[], int B[], int D[]) {
n = N;
adj.resize(n + 1);
for (int i = 0; i < n - 1; i ++) {
A[i] ++; B[i] ++;
adj[A[i]].push_back({B[i], D[i]});
adj[B[i]].push_back({A[i], D[i]});
}
for (logn = 1; (1 << logn) <= n; logn ++);
h.assign(n + 1, 0);
p.assign(logn, vector <int> (n + 1, 0));
f.assign(n + 1, 0);
h[1] = 1;
dfs_lca(1);
for (int i = 1; i < logn; i ++)
for (int u = 1; u <= n; u ++)
p[i][u] = p[i - 1][p[i - 1][u]];
sz.assign(n + 1, 0);
pa.assign(n + 1, 0);
build(1, 0);
ans.assign(n + 1, inf);
}
long long Query(int S, int X[], int T, int Y[]) {
vector <int> del;
for (int i = 0; i < S; i ++) {
int u = X[i] + 1;
ans[u] = 0;
del.push_back(u);
for (int v = u; v != pa[v]; ) {
v = pa[v];
ans[v] = min(ans[v], dis(v, u));
del.push_back(v);
}
}
ll res = inf;
for (int i = 0; i < T; i ++) {
int u = Y[i] + 1;
res = min(res, ans[u]);
for (int v = u; v != pa[v]; ) {
v = pa[v];
res = min(res, dis(v, u) + ans[v]);
}
}
for (int u : del)
ans[u] = inf;
return res;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
16728 KB |
Output is correct |
2 |
Correct |
350 ms |
21704 KB |
Output is correct |
3 |
Correct |
664 ms |
21584 KB |
Output is correct |
4 |
Correct |
631 ms |
22064 KB |
Output is correct |
5 |
Correct |
633 ms |
22080 KB |
Output is correct |
6 |
Correct |
150 ms |
21584 KB |
Output is correct |
7 |
Correct |
640 ms |
21692 KB |
Output is correct |
8 |
Correct |
652 ms |
21584 KB |
Output is correct |
9 |
Correct |
606 ms |
22080 KB |
Output is correct |
10 |
Correct |
144 ms |
21584 KB |
Output is correct |
11 |
Correct |
637 ms |
21852 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
16728 KB |
Output is correct |
2 |
Correct |
2074 ms |
104908 KB |
Output is correct |
3 |
Correct |
4824 ms |
125644 KB |
Output is correct |
4 |
Correct |
471 ms |
124096 KB |
Output is correct |
5 |
Correct |
5884 ms |
147912 KB |
Output is correct |
6 |
Correct |
5277 ms |
126540 KB |
Output is correct |
7 |
Correct |
2188 ms |
51960 KB |
Output is correct |
8 |
Correct |
250 ms |
51804 KB |
Output is correct |
9 |
Correct |
1598 ms |
54376 KB |
Output is correct |
10 |
Correct |
2168 ms |
52172 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
9 ms |
16728 KB |
Output is correct |
2 |
Correct |
350 ms |
21704 KB |
Output is correct |
3 |
Correct |
664 ms |
21584 KB |
Output is correct |
4 |
Correct |
631 ms |
22064 KB |
Output is correct |
5 |
Correct |
633 ms |
22080 KB |
Output is correct |
6 |
Correct |
150 ms |
21584 KB |
Output is correct |
7 |
Correct |
640 ms |
21692 KB |
Output is correct |
8 |
Correct |
652 ms |
21584 KB |
Output is correct |
9 |
Correct |
606 ms |
22080 KB |
Output is correct |
10 |
Correct |
144 ms |
21584 KB |
Output is correct |
11 |
Correct |
637 ms |
21852 KB |
Output is correct |
12 |
Correct |
3 ms |
16728 KB |
Output is correct |
13 |
Correct |
2074 ms |
104908 KB |
Output is correct |
14 |
Correct |
4824 ms |
125644 KB |
Output is correct |
15 |
Correct |
471 ms |
124096 KB |
Output is correct |
16 |
Correct |
5884 ms |
147912 KB |
Output is correct |
17 |
Correct |
5277 ms |
126540 KB |
Output is correct |
18 |
Correct |
2188 ms |
51960 KB |
Output is correct |
19 |
Correct |
250 ms |
51804 KB |
Output is correct |
20 |
Correct |
1598 ms |
54376 KB |
Output is correct |
21 |
Correct |
2168 ms |
52172 KB |
Output is correct |
22 |
Correct |
3130 ms |
132448 KB |
Output is correct |
23 |
Correct |
3277 ms |
136620 KB |
Output is correct |
24 |
Execution timed out |
8071 ms |
139140 KB |
Time limit exceeded |
25 |
Halted |
0 ms |
0 KB |
- |