#include <bits/stdc++.h>
#include "factories.h"
using namespace std;
#define inf 0x3F3F3F3F3F3F3F3F
const int MXN = 5e5 + 5;
const int LOG = 21;
const int B = 600;
int n;
vector<array<int, 2>> adj[MXN];
long long dist[MXN], dep[MXN];
int p[LOG][MXN];
int in[MXN], out[MXN], tim = -1;
void dfs(int a)
{
in[a] = ++tim;
for (array<int, 2> &v : adj[a])
{
if (v[0] == p[0][a]) continue;
dep[v[0]] = dep[a] + v[1];
p[0][v[0]] = a;
dfs(v[0]);
}
out[a] = tim;
}
void Init(int N, int A[], int B[], int D[])
{
n = N;
for (int i = 0; i < n - 1; i++)
{
adj[A[i]].push_back({B[i], D[i]});
adj[B[i]].push_back({A[i], D[i]});
}
p[0][0] = 0;
dfs(0);
for (int i = 1; i < LOG; i++) for (int j = 0; j < n; j++) p[i][j] = p[i - 1][p[i - 1][j]];
}
void bfs(vector<int> &s)
{
for (int i = 0; i < n; i++) dist[i] = inf;
priority_queue<array<long long, 2>, vector<array<long long, 2>>, greater<array<long long, 2>>> pq;
for (int x : s)
{
dist[x] = 0;
pq.push({0, x});
}
while (!pq.empty())
{
array<long long, 2> f = pq.top();
pq.pop();
if (f[0] > dist[f[1]]) continue;
for (array<int, 2> &v : adj[f[1]])
{
if (dist[v[0]] > f[0] + v[1])
{
dist[v[0]] = f[0] + v[1];
pq.push({dist[v[0]], v[0]});
}
}
}
}
int anc(int u, int v)
{
return in[u] <= in[v] && out[v] <= out[u];
}
int LCA(int u, int v)
{
if (anc(u, v)) return u;
if (anc(v, u)) return v;
for (int i = LOG - 1; i >= 0; i--)
{
if (anc(p[i][u], v)) continue;
u = p[i][u];
}
return p[0][u];
}
long long DIST(int u, int v)
{
return dep[u] + dep[v] - 2 * dep[LCA(u, v)];
}
long long Query(int S, int X[], int T, int Y[])
{
vector<int> v[2];
for (int i = 0; i < S; i++) v[0].push_back(X[i]);
for (int i = 0; i < T; i++) v[1].push_back(Y[i]);
if (v[0].size() < v[1].size()) swap(v[0], v[1]);
long long res = inf;
if (v[0].size() >= B)
{
bfs(v[0]);
for (int x : v[1]) res = min(res, dist[x]);
}
else
{
for (int x : v[0])
{
for (int y : v[1])
{
res = min(res, DIST(x, y));
}
}
}
return res;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
128 ms |
60228 KB |
Output is correct |
2 |
Correct |
2340 ms |
70232 KB |
Output is correct |
3 |
Execution timed out |
8004 ms |
69712 KB |
Time limit exceeded |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
10 ms |
59996 KB |
Output is correct |
2 |
Correct |
2245 ms |
121344 KB |
Output is correct |
3 |
Correct |
4569 ms |
123452 KB |
Output is correct |
4 |
Correct |
2428 ms |
120956 KB |
Output is correct |
5 |
Correct |
1170 ms |
145660 KB |
Output is correct |
6 |
Correct |
3920 ms |
124320 KB |
Output is correct |
7 |
Correct |
4322 ms |
97712 KB |
Output is correct |
8 |
Correct |
1824 ms |
96132 KB |
Output is correct |
9 |
Correct |
1191 ms |
98932 KB |
Output is correct |
10 |
Correct |
2447 ms |
98324 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
128 ms |
60228 KB |
Output is correct |
2 |
Correct |
2340 ms |
70232 KB |
Output is correct |
3 |
Execution timed out |
8004 ms |
69712 KB |
Time limit exceeded |
4 |
Halted |
0 ms |
0 KB |
- |