#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
using namespace __gnu_pbds;
template<typename T> using Tree = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
typedef long long int ll;
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
#define fastio ios_base::sync_with_stdio(false); cin.tie(NULL)
#define pb push_back
#define endl '\n'
#define sz(a) (int)a.size()
#define setbits(x) __builtin_popcountll(x)
#define ff first
#define ss second
#define conts continue
#define ceil2(x,y) ((x+y-1)/(y))
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
#define yes cout << "Yes" << endl
#define no cout << "No" << endl
#define rep(i,n) for(int i = 0; i < n; ++i)
#define rep1(i,n) for(int i = 1; i <= n; ++i)
#define rev(i,s,e) for(int i = s; i >= e; --i)
#define trav(i,a) for(auto &i : a)
template<typename T>
void amin(T &a, T b) {
a = min(a,b);
}
template<typename T>
void amax(T &a, T b) {
a = max(a,b);
}
#ifdef LOCAL
#include "debug.h"
#else
#define debug(x) 42
#endif
/*
refs:
edi
https://codeforces.com/blog/entry/65042?#comment-491880
*/
const int MOD = 1e9 + 7;
const int N = 2e5 + 5;
const int inf1 = int(1e9) + 5;
const ll inf2 = ll(1e18) + 5;
vector<ll> adj[N];
vector<ll> a(N);
pll diam = {-1,-1};
void dfs1(ll u, ll p, ll d){
pll px = {d,u};
amax(diam,px);
trav(v,adj[u]){
if(v == p) conts;
dfs1(v,u,d+1);
}
}
ll dis[N][2];
ll deepest[N][2];
void dfs2(ll u, ll p, ll t){
trav(v,adj[u]){
if(v == p) conts;
dis[v][t] = dis[u][t]+1;
dfs2(v,u,t);
amax(deepest[u][t],deepest[v][t]+1);
}
}
vector<ll> ans(N);
vector<ll> mx_without(N);
vector<ll> first_occ(N,inf2);
vector<pll> active(N);
ll siz = 0;
void dfs3(ll u, ll p, ll d, ll t){
ll mx = deepest[u][t];
ll cnt = upper_bound(active.begin(),active.begin()+siz,make_pair(d-mx,-1ll))-active.begin();
if(dis[u][t] >= dis[u][t^1]){
ans[u] = cnt;
}
mx = 0;
trav(v,adj[u]){
if(v == p) conts;
mx_without[v] = mx;
amax(mx,deepest[v][t]+1);
}
mx = 0;
reverse(all(adj[u]));
trav(v,adj[u]){
if(v == p) conts;
amax(mx_without[v],mx);
amax(mx,deepest[v][t]+1);
}
reverse(all(adj[u]));
trav(v,adj[u]){
if(v == p) conts;
vector<array<ll,4>> history;
mx = mx_without[v];
ll new_siz = upper_bound(active.begin(),active.begin()+siz,make_pair(d-mx,-1ll))-active.begin();
ll old_siz = siz;
history.pb({3,siz,-1,-1});
siz = new_siz;
ll x = a[u];
if(first_occ[x] >= siz){
history.pb({1,siz,active[siz].ff,active[siz].ss});
history.pb({2,x,first_occ[x],-1});
history.pb({3,siz,-1,-1});
ll w = active[siz].ss;
history.pb({2,a[w],first_occ[a[w]],-1});
first_occ[a[w]] = inf2;
active[siz] = {d,u};
first_occ[x] = siz;
siz++;
}
dfs3(v,u,d+1,t);
reverse(all(history));
trav(ar,history){
auto [t,ind,val1,val2] = ar;
if(t == 1){
active[ind] = {val1,val2};
}
else if(t == 2){
first_occ[ind] = val1;
}
else{
siz = ind;
}
}
}
}
void solve(int test_case)
{
ll n,m; cin >> n >> m;
rep1(i,n-1){
ll u,v; cin >> u >> v;
adj[u].pb(v), adj[v].pb(u);
}
rep1(i,n) cin >> a[i];
dfs1(1,-1,0);
ll s = diam.ss;
diam = {-1,-1};
dfs1(s,-1,0);
ll t = diam.ss;
dfs2(s,-1,0);
dfs2(t,-1,1);
dfs3(s,-1,0,0);
dfs3(t,-1,0,1);
rep1(i,n){
cout << ans[i] << endl;
}
}
int main()
{
fastio;
int t = 1;
// cin >> t;
rep1(i, t) {
solve(i);
}
return 0;
}
Compilation message
joi2019_ho_t5.cpp: In function 'void dfs3(ll, ll, ll, ll)':
joi2019_ho_t5.cpp:126:12: warning: unused variable 'old_siz' [-Wunused-variable]
126 | ll old_siz = siz;
| ^~~~~~~
| # |
결과 |
실행 시간 |
메모리 |
Grader output |
| 1 |
Correct |
5 ms |
16476 KB |
Output is correct |
| 2 |
Correct |
7 ms |
16732 KB |
Output is correct |
| 3 |
Correct |
8 ms |
16820 KB |
Output is correct |
| 4 |
Incorrect |
7 ms |
16988 KB |
Output isn't correct |
| 5 |
Halted |
0 ms |
0 KB |
- |
| # |
결과 |
실행 시간 |
메모리 |
Grader output |
| 1 |
Correct |
124 ms |
24908 KB |
Output is correct |
| 2 |
Correct |
194 ms |
54924 KB |
Output is correct |
| 3 |
Correct |
33 ms |
22364 KB |
Output is correct |
| 4 |
Correct |
303 ms |
30536 KB |
Output is correct |
| 5 |
Correct |
439 ms |
82256 KB |
Output is correct |
| 6 |
Correct |
424 ms |
56148 KB |
Output is correct |
| 7 |
Correct |
297 ms |
30216 KB |
Output is correct |
| 8 |
Correct |
392 ms |
35464 KB |
Output is correct |
| 9 |
Correct |
374 ms |
33836 KB |
Output is correct |
| 10 |
Correct |
352 ms |
33796 KB |
Output is correct |
| 11 |
Correct |
216 ms |
30664 KB |
Output is correct |
| 12 |
Correct |
478 ms |
63228 KB |
Output is correct |
| 13 |
Correct |
414 ms |
55588 KB |
Output is correct |
| 14 |
Correct |
479 ms |
54256 KB |
Output is correct |
| 15 |
Correct |
167 ms |
28472 KB |
Output is correct |
| 16 |
Correct |
394 ms |
65148 KB |
Output is correct |
| 17 |
Correct |
420 ms |
56268 KB |
Output is correct |
| # |
결과 |
실행 시간 |
메모리 |
Grader output |
| 1 |
Correct |
263 ms |
28628 KB |
Output is correct |
| 2 |
Correct |
702 ms |
122704 KB |
Output is correct |
| 3 |
Correct |
44 ms |
26772 KB |
Output is correct |
| 4 |
Correct |
417 ms |
30872 KB |
Output is correct |
| 5 |
Correct |
660 ms |
127312 KB |
Output is correct |
| 6 |
Correct |
650 ms |
78428 KB |
Output is correct |
| 7 |
Correct |
435 ms |
30892 KB |
Output is correct |
| 8 |
Correct |
465 ms |
48756 KB |
Output is correct |
| 9 |
Correct |
467 ms |
42372 KB |
Output is correct |
| 10 |
Correct |
432 ms |
37460 KB |
Output is correct |
| 11 |
Correct |
334 ms |
30408 KB |
Output is correct |
| 12 |
Correct |
701 ms |
108168 KB |
Output is correct |
| 13 |
Correct |
560 ms |
71292 KB |
Output is correct |
| 14 |
Correct |
567 ms |
72636 KB |
Output is correct |
| 15 |
Correct |
191 ms |
29128 KB |
Output is correct |
| 16 |
Correct |
493 ms |
95224 KB |
Output is correct |
| 17 |
Correct |
550 ms |
78556 KB |
Output is correct |
| # |
결과 |
실행 시간 |
메모리 |
Grader output |
| 1 |
Correct |
5 ms |
16476 KB |
Output is correct |
| 2 |
Correct |
7 ms |
16732 KB |
Output is correct |
| 3 |
Correct |
8 ms |
16820 KB |
Output is correct |
| 4 |
Incorrect |
7 ms |
16988 KB |
Output isn't correct |
| 5 |
Halted |
0 ms |
0 KB |
- |
