#include "factories.h"
#include<iostream>
#include<stdio.h>
#include<vector>
using namespace std ;
#define MAXN 500007
#define LOG 20
int n ;
vector < pair < int , int > > v[ MAXN ] ;
long long dist[ MAXN ][ LOG + 1 ] ;
int LCA[ MAXN ][ LOG + 1 ] ;
int lvl[ MAXN ] ;
int subtree[ MAXN ] ;
int lst[ MAXN ] ;
int used[ MAXN ] ;
int tot = 0 ;
long long ans[ MAXN ] ;
int w[ MAXN ] ;
void dfs ( int vertex , int prv , int ori ) {
if ( used[ vertex ] == 1 ) { return ; }
tot ++ ;
int i ;
int sz = v[ vertex ].size ( ) ;
subtree[ vertex ] = 1 ;
for ( i = 0 ; i < sz ; i ++ ) {
if ( v[ vertex ][ i ].first == prv ) { continue ; }
if ( used[ v[ vertex ][ i ].first ] != 0 ) { continue ; }
dist[ v[ vertex ][ i ].first ][ w[ ori ] ] = dist[ vertex ][ w[ ori ] ] + v[ vertex ][ i ].second ;
dfs ( v[ vertex ][ i ].first , vertex , ori ) ;
subtree[ vertex ] += subtree[ v[ vertex ][ i ].first ] ;
}
}
int find_centroid ( int vertex , int prv ) {
int i ;
int sz = v[ vertex ].size ( ) ;
for ( i = 0 ; i < sz ; i ++ ) {
if ( v[ vertex ][ i ].first == prv ) { continue ; }
if ( used[ v[ vertex ][ i ].first ] == 0 && subtree[ v[ vertex ][ i ].first ] > ( tot / 2 ) ) {
return find_centroid ( v[ vertex ][ i ].first , vertex ) ;
}
}
return vertex ;
}
void decompose ( int vertex , int prv ) {
tot = 0 ;
if ( prv == -1 ) {
dfs ( vertex , -1 , 0 ) ;
}
else {
dfs ( vertex , -1 , prv ) ;
}
int cen = find_centroid ( vertex , -1 ) ;
lst[ cen ] = prv ;
used[ cen ] = 1 ;
if ( prv != -1 ) {
w[ cen ] = w[ prv ] + 1 ;
}
else {
w[ cen ] = 1 ;
}
int i , sz ;
sz = v[ cen ].size ( ) ;
for ( i = 0 ; i < sz ; i ++ ) {
if ( used[ v[ cen ][ i ].first ] == 0 ) {
dist[ v[ cen ][ i ].first ][ w[ cen ] ] = v[ cen ][ i ].second ;
decompose ( v[ cen ][ i ].first , cen ) ;
}
}
}
void update ( int vertex , int x ) {
while ( vertex != -1 ) {
long long cur = dist[ x ][ w[ vertex ] ] ;
if ( ans[ vertex ] == -1 ) { ans[ vertex ] = cur ; }
if ( ans[ vertex ] > cur ) { ans[ vertex ] = cur ; }
vertex = lst[ vertex ] ;
}
}
void reset ( int vertex ) {
while ( vertex != -1 ) {
ans[ vertex ] = -1 ;
vertex = lst[ vertex ] ;
}
}
long long query ( int vertex , int x ) {
long long ret = -1 ;
while ( vertex != -1 ) {
if ( ans[ vertex ] > -1 ) {
long long cur = dist[ x ][ w[ vertex ] ] + ans[ vertex ] ;
if ( ret == -1 ) { ret = cur ; }
if ( ret > cur ) { ret = cur ; }
}
vertex = lst[ vertex ] ;
}
return ret ;
}
void Init(int N, int A[], int B[], int D[]) {
n = N ;
int i ;
for ( i = 0 ; i < N - 1 ; i ++ ) {
A[ i ] ++ ;
B[ i ] ++ ;
v[ A[ i ] ].push_back ( make_pair ( B[ i ] , D[ i ] ) ) ;
v[ B[ i ] ].push_back ( make_pair ( A[ i ] , D[ i ] ) ) ;
}
decompose ( 1 , -1 ) ;
for ( i = 1 ; i <= n ; i ++ ) {
ans[ i ] = -1 ;
}
}
long long Query(int S, int X[], int T, int Y[]) {
int i ;
for ( i = 0 ; i < S ; i ++ ) {
X[ i ] ++ ;
}
for ( i = 0 ; i < T ; i ++ ) {
Y[ i ] ++ ;
}
for ( i = 0 ; i < S ; i ++ ) {
update ( X[ i ] , X[ i ] ) ;
}
long long ret = -1 ;
for ( i = 0 ; i < T ; i ++ ) {
long long cur = query ( Y[ i ] , Y[ i ] ) ;
if ( cur > -1 ) {
if ( ret == -1 ) { ret = cur ; }
if ( ret > cur ) { ret = cur ; }
}
}
for ( i = 0 ; i < S ; i ++ ) {
reset ( X[ i ] ) ;
}
return ret ;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
13 ms |
172724 KB |
Output is correct |
2 |
Correct |
489 ms |
172988 KB |
Output is correct |
3 |
Correct |
416 ms |
172988 KB |
Output is correct |
4 |
Correct |
439 ms |
172988 KB |
Output is correct |
5 |
Correct |
443 ms |
173076 KB |
Output is correct |
6 |
Correct |
323 ms |
173024 KB |
Output is correct |
7 |
Correct |
469 ms |
172988 KB |
Output is correct |
8 |
Correct |
503 ms |
172988 KB |
Output is correct |
9 |
Correct |
546 ms |
173080 KB |
Output is correct |
10 |
Correct |
356 ms |
173024 KB |
Output is correct |
11 |
Correct |
479 ms |
172988 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
172724 KB |
Output is correct |
2 |
Correct |
2129 ms |
191468 KB |
Output is correct |
3 |
Correct |
2719 ms |
195360 KB |
Output is correct |
4 |
Correct |
1219 ms |
192436 KB |
Output is correct |
5 |
Correct |
4389 ms |
221740 KB |
Output is correct |
6 |
Correct |
3463 ms |
194828 KB |
Output is correct |
7 |
Correct |
1109 ms |
177136 KB |
Output is correct |
8 |
Correct |
573 ms |
176956 KB |
Output is correct |
9 |
Correct |
1419 ms |
179832 KB |
Output is correct |
10 |
Correct |
999 ms |
177016 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
3016 ms |
191468 KB |
Output is correct |
2 |
Correct |
2316 ms |
191468 KB |
Output is correct |
3 |
Correct |
3943 ms |
194136 KB |
Output is correct |
4 |
Correct |
3469 ms |
194856 KB |
Output is correct |
5 |
Correct |
4073 ms |
194680 KB |
Output is correct |
6 |
Correct |
4943 ms |
214512 KB |
Output is correct |
7 |
Correct |
1376 ms |
192436 KB |
Output is correct |
8 |
Correct |
4243 ms |
194224 KB |
Output is correct |
9 |
Correct |
3599 ms |
193560 KB |
Output is correct |
10 |
Correct |
3319 ms |
194272 KB |
Output is correct |
11 |
Correct |
1413 ms |
180428 KB |
Output is correct |
12 |
Correct |
643 ms |
176956 KB |
Output is correct |
13 |
Correct |
966 ms |
176420 KB |
Output is correct |
14 |
Correct |
996 ms |
176420 KB |
Output is correct |
15 |
Correct |
1089 ms |
177020 KB |
Output is correct |
16 |
Correct |
1123 ms |
177004 KB |
Output is correct |