#include <stdio.h>
#include <stdlib.h>
#include <bits/stdc++.h>
#define MAX_N 500000
#define MAX_Q 100000
#define MAX_SUM_ST 1000000
#define MAX_VALUE 1000000000
#define ff first
#define ss second
#define pb push_back
#include "factories.h"
#define mp make_pair
typedef long long ll;
using namespace std;
ll n,a[1000001],b[1000001],d[1000001],distnce[1000001],euler[5000000],k=0,first[1000001],belong[25][1000001],sz[1000001],centroid,nearest[1000001],o=0,p[1000001];
bool vis[1000001];
vector<int>vec;
ll st[25][2000001];
vector<pair<ll,ll> >adj[1000001];
void dfs(ll node,ll par,ll dis){
distnce[node]=dis;
for(int i=0;i<adj[node].size();i++){
if(adj[node][i].ff!=par){
dfs(adj[node][i].ff,node,dis+adj[node][i].ss);
}
}
}
void build(ll node,ll par){
k++;
euler[k]=distnce[node];
first[node]=k;
for(int i=0;i<adj[node].size();i++){
if(adj[node][i].ff!=par){
build(adj[node][i].ff,node);
k++;
euler[k]=distnce[node];
}
}
}
void dfs1(ll node,ll par){
sz[node]=1;
for(int i=0;i<adj[node].size();i++){
if(adj[node][i].ff!=par && vis[adj[node][i].ff]==0){
dfs1(adj[node][i].ff,node);
sz[node]+=sz[adj[node][i].ff];
}
}
}
void solve(ll node,ll par,ll level,ll parent){
// cout << node << '\n';
o++;
dfs1(node,par);
ll cent=node;
while(1){
bool f=true;
for(int i=0;i<adj[cent].size();i++){
if(sz[adj[cent][i].ff]<sz[cent] && vis[adj[cent][i].ff]==0 && sz[adj[cent][i].ff]*2>sz[node]){
f=false;
cent=adj[cent][i].ff;
break;
}
}
if(f) break;
}
centroid=cent;
// cout << cent << ' ' << level << '\n';
p[cent]=parent;
vis[cent]=1;
for(int i=0;i<adj[cent].size();i++){
if(vis[adj[cent][i].ff]==0){
solve(adj[cent][i].ff,-1,level+1,cent);
}
}
}
void Init(int N, int A[], int B[], int D[]) {
n=(ll)N;
for(int i=0;i<(N-1);i++){
a[i]=(ll)A[i];
b[i]=(ll)B[i];
d[i]=(ll)D[i];
adj[a[i]].pb(mp(b[i],d[i]));
adj[b[i]].pb(mp(a[i],d[i]));
}
dfs(0,-1,0);
build(0,-1);
for(int i=1;i<=k;i++){
st[0][i]=euler[i];
}
int l=log2(k);
for(int i=1;i<=l;i++){
for(int j=1;j<=(k-(1<<i)+1);j++){
st[i][j]=min(st[i-1][j],st[i-1][j+(1<<(i-1))]);
}
}
solve(0,-1,1,-1);
for(int i=0;i<n;i++){
nearest[i]=1e18;
}
}
ll dist(int x,int y){
int xx=min(first[x],first[y]);
int yy=max(first[x],first[y]);
// cout << xx << ' ' << yy << '\n';
int z=log2(yy-xx+1);
ll mn=(min(st[z][xx],st[z][yy-(1<<z)+1]));
return distnce[x]+distnce[y]-2*mn;
}
void add(int node,int cent){
if(cent==-1) return;
nearest[cent]=min(nearest[cent],dist(node,cent));
// cout << node << ' ' << cent << ' ' << nearest[cent] << '\n';
vec.pb(cent);
add(node,p[cent]);
}
ll query(int node,int cent){
if(cent==-1) return 2e18;
return min(nearest[cent]+dist(node,cent),query(node,p[cent]));
}
long long Query(int S, int X[], int T, int Y[]) {
for(int i=0;i<S;i++){
add(X[i],X[i]);
}
ll ans=2e18;
for(int i=0;i<T;i++){
ans=min(ans,query(Y[i],Y[i]));
}
for(int i=0;i<vec.size();i++){
nearest[vec[i]]=1e18;
}
vec.clear();
return ans;
}
Compilation message
factories.cpp: In function 'void dfs(ll, ll, ll)':
factories.cpp:22:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<long long int, long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
22 | for(int i=0;i<adj[node].size();i++){
| ~^~~~~~~~~~~~~~~~~
factories.cpp: In function 'void build(ll, ll)':
factories.cpp:32:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<long long int, long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
32 | for(int i=0;i<adj[node].size();i++){
| ~^~~~~~~~~~~~~~~~~
factories.cpp: In function 'void dfs1(ll, ll)':
factories.cpp:42:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<long long int, long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
42 | for(int i=0;i<adj[node].size();i++){
| ~^~~~~~~~~~~~~~~~~
factories.cpp: In function 'void solve(ll, ll, ll, ll)':
factories.cpp:56:18: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<long long int, long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
56 | for(int i=0;i<adj[cent].size();i++){
| ~^~~~~~~~~~~~~~~~~
factories.cpp:69:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<long long int, long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
69 | for(int i=0;i<adj[cent].size();i++){
| ~^~~~~~~~~~~~~~~~~
factories.cpp: In function 'long long int Query(int, int*, int, int*)':
factories.cpp:127:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
127 | for(int i=0;i<vec.size();i++){
| ~^~~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
39 ms |
24428 KB |
Output is correct |
2 |
Correct |
950 ms |
38700 KB |
Output is correct |
3 |
Correct |
1248 ms |
38892 KB |
Output is correct |
4 |
Correct |
1236 ms |
39084 KB |
Output is correct |
5 |
Correct |
1438 ms |
39148 KB |
Output is correct |
6 |
Correct |
438 ms |
38636 KB |
Output is correct |
7 |
Correct |
1217 ms |
36824 KB |
Output is correct |
8 |
Correct |
1275 ms |
36588 KB |
Output is correct |
9 |
Correct |
1443 ms |
37096 KB |
Output is correct |
10 |
Correct |
428 ms |
36460 KB |
Output is correct |
11 |
Correct |
1217 ms |
36760 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
18 ms |
24172 KB |
Output is correct |
2 |
Correct |
3886 ms |
252212 KB |
Output is correct |
3 |
Correct |
5161 ms |
254920 KB |
Output is correct |
4 |
Correct |
1186 ms |
249544 KB |
Output is correct |
5 |
Correct |
7175 ms |
279704 KB |
Output is correct |
6 |
Correct |
5497 ms |
256128 KB |
Output is correct |
7 |
Correct |
4495 ms |
88308 KB |
Output is correct |
8 |
Correct |
878 ms |
88224 KB |
Output is correct |
9 |
Correct |
5574 ms |
92640 KB |
Output is correct |
10 |
Correct |
4379 ms |
89696 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
39 ms |
24428 KB |
Output is correct |
2 |
Correct |
950 ms |
38700 KB |
Output is correct |
3 |
Correct |
1248 ms |
38892 KB |
Output is correct |
4 |
Correct |
1236 ms |
39084 KB |
Output is correct |
5 |
Correct |
1438 ms |
39148 KB |
Output is correct |
6 |
Correct |
438 ms |
38636 KB |
Output is correct |
7 |
Correct |
1217 ms |
36824 KB |
Output is correct |
8 |
Correct |
1275 ms |
36588 KB |
Output is correct |
9 |
Correct |
1443 ms |
37096 KB |
Output is correct |
10 |
Correct |
428 ms |
36460 KB |
Output is correct |
11 |
Correct |
1217 ms |
36760 KB |
Output is correct |
12 |
Correct |
18 ms |
24172 KB |
Output is correct |
13 |
Correct |
3886 ms |
252212 KB |
Output is correct |
14 |
Correct |
5161 ms |
254920 KB |
Output is correct |
15 |
Correct |
1186 ms |
249544 KB |
Output is correct |
16 |
Correct |
7175 ms |
279704 KB |
Output is correct |
17 |
Correct |
5497 ms |
256128 KB |
Output is correct |
18 |
Correct |
4495 ms |
88308 KB |
Output is correct |
19 |
Correct |
878 ms |
88224 KB |
Output is correct |
20 |
Correct |
5574 ms |
92640 KB |
Output is correct |
21 |
Correct |
4379 ms |
89696 KB |
Output is correct |
22 |
Correct |
5842 ms |
254788 KB |
Output is correct |
23 |
Correct |
5913 ms |
259176 KB |
Output is correct |
24 |
Execution timed out |
8079 ms |
260252 KB |
Time limit exceeded |
25 |
Halted |
0 ms |
0 KB |
- |