Submission #958504

#TimeUsernameProblemLanguageResultExecution timeMemory
958504JwFXoizFactories (JOI14_factories)C++14
15 / 100
8068 ms158056 KiB
#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<long long, 2>> adj[MXN]; long long dist[MXN], dep[MXN]; long long p[LOG][MXN]; long long in[MXN], out[MXN], tim = -1; void dfs(int a) { in[a] = ++tim; for (array<long long, 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[]) { for (int i = 0; i < n; i++) { adj[i].clear(); dist[i] = dep[i] = 0; for (int j = 0; j < LOG; j++) p[j][i] = 0; in[i] = out[i] = 0; } 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<long long, 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]); long long res = inf; if (max(S, T) >= 1) { 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 timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...