#include "factories.h"
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")
const int nx=5e5+5, kx=19;
ll lvl[nx], pa[nx][kx], ds[nx], c[nx], sz[nx], dp[nx];
vector<pair<ll, ll>> d[nx];
bool used[nx];
vector<ll> rem;
void dfs(int u, int p)
{
lvl[u]=lvl[p]+1;
pa[u][0]=p;
for (int i=1; i<kx; i++) pa[u][i]=pa[pa[u][i-1]][i-1];
for (auto [v, w]:d[u]) if (v!=p) ds[v]=ds[u]+w, dfs(v, u);
}
int lca(int u, int v)
{
if (lvl[u]>lvl[v]) swap(u, v);
for (int i=kx-1; i>=0; i--) if (lvl[pa[v][i]]>=lvl[u]) v=pa[v][i];
if (u==v) return v;
for (int i=kx-1; i>=0; i--) if (pa[u][i]!=pa[v][i]) u=pa[u][i], v=pa[v][i];
return pa[u][0];
}
ll distance(int u, int v)
{
return ds[u]+ds[v]-2*ds[lca(u, v)];
}
int dfssz(int u, int p)
{
sz[u]=1;
for (auto [v, w]:d[u]) if (v!=p&&!used[v]) sz[u]+=dfssz(v, u);
return sz[u];
}
int findcentroid(int u, int p, int rtsz)
{
for (auto [v, w]:d[u]) if (v!=p&&!used[v]&&2*sz[v]>rtsz) return findcentroid(v, u, rtsz);
return u;
}
void decomposition(int u, int p)
{
u=findcentroid(u, u, dfssz(u, u));
used[u]=1;
c[u]=p;
for (auto [v, w]:d[u]) if (v!=p&&!used[v]) decomposition(v, u);
}
void update(int u)
{
int rt=u;
while (u!=-1) dp[u]=min(dp[u], distance(u, rt)), rem.push_back(u), u=c[u];
}
ll query(int u)
{
ll tmp=LLONG_MAX, rt=u;
while (u!=-1) tmp=(dp[u]!=LLONG_MAX)?min(tmp, distance(u, rt)+dp[u]):tmp, u=c[u];
return tmp;
}
void Init(int N, int A[], int B[], int D[]) {
for (int i=0; i<N-1; i++) d[A[i]].push_back({B[i], D[i]}), d[B[i]].push_back({A[i], D[i]});
for (int i=0; i<N; i++) dp[i]=LLONG_MAX;
dfs(0, 0);
decomposition(0, -1);
//for (int i=0; i<N; i++) cout<<i<<' '<<c[i]<<'\n';
}
long long Query(int S, int X[], int T, int Y[]) {
ll res=LLONG_MAX;
for (int i=0; i<S; i++) update(X[i]);
for (int i=0; i<T; i++) res=min(res, query(Y[i]));
for (auto x:rem) dp[x]=LLONG_MAX;
rem.clear();
return res;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
17 ms |
39256 KB |
Output is correct |
2 |
Correct |
553 ms |
44012 KB |
Output is correct |
3 |
Correct |
1219 ms |
43888 KB |
Output is correct |
4 |
Correct |
992 ms |
44368 KB |
Output is correct |
5 |
Correct |
1600 ms |
44624 KB |
Output is correct |
6 |
Correct |
221 ms |
43872 KB |
Output is correct |
7 |
Correct |
1246 ms |
43996 KB |
Output is correct |
8 |
Correct |
1300 ms |
43972 KB |
Output is correct |
9 |
Correct |
1614 ms |
44492 KB |
Output is correct |
10 |
Correct |
220 ms |
43852 KB |
Output is correct |
11 |
Correct |
1251 ms |
43892 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
7 ms |
39256 KB |
Output is correct |
2 |
Correct |
1705 ms |
155468 KB |
Output is correct |
3 |
Correct |
4864 ms |
160628 KB |
Output is correct |
4 |
Correct |
655 ms |
153400 KB |
Output is correct |
5 |
Correct |
7615 ms |
189880 KB |
Output is correct |
6 |
Correct |
5182 ms |
179188 KB |
Output is correct |
7 |
Correct |
3043 ms |
87496 KB |
Output is correct |
8 |
Correct |
372 ms |
86472 KB |
Output is correct |
9 |
Correct |
4607 ms |
91300 KB |
Output is correct |
10 |
Correct |
3403 ms |
88052 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
17 ms |
39256 KB |
Output is correct |
2 |
Correct |
553 ms |
44012 KB |
Output is correct |
3 |
Correct |
1219 ms |
43888 KB |
Output is correct |
4 |
Correct |
992 ms |
44368 KB |
Output is correct |
5 |
Correct |
1600 ms |
44624 KB |
Output is correct |
6 |
Correct |
221 ms |
43872 KB |
Output is correct |
7 |
Correct |
1246 ms |
43996 KB |
Output is correct |
8 |
Correct |
1300 ms |
43972 KB |
Output is correct |
9 |
Correct |
1614 ms |
44492 KB |
Output is correct |
10 |
Correct |
220 ms |
43852 KB |
Output is correct |
11 |
Correct |
1251 ms |
43892 KB |
Output is correct |
12 |
Correct |
7 ms |
39256 KB |
Output is correct |
13 |
Correct |
1705 ms |
155468 KB |
Output is correct |
14 |
Correct |
4864 ms |
160628 KB |
Output is correct |
15 |
Correct |
655 ms |
153400 KB |
Output is correct |
16 |
Correct |
7615 ms |
189880 KB |
Output is correct |
17 |
Correct |
5182 ms |
179188 KB |
Output is correct |
18 |
Correct |
3043 ms |
87496 KB |
Output is correct |
19 |
Correct |
372 ms |
86472 KB |
Output is correct |
20 |
Correct |
4607 ms |
91300 KB |
Output is correct |
21 |
Correct |
3403 ms |
88052 KB |
Output is correct |
22 |
Correct |
3043 ms |
183404 KB |
Output is correct |
23 |
Correct |
2525 ms |
188016 KB |
Output is correct |
24 |
Execution timed out |
8095 ms |
192188 KB |
Time limit exceeded |
25 |
Halted |
0 ms |
0 KB |
- |