// いろいろ高速化 #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #ifndef HIDDEN_IN_VS // 折りたたみ用 // 警告の抑制 #define _CRT_SECURE_NO_WARNINGS // ライブラリの読み込み #include using namespace std; // 型名の短縮 using ll = long long; using ull = unsigned long long; // -2^63 ~ 2^63 = 9 * 10^18(int は -2^31 ~ 2^31 = 2 * 10^9) using pii = pair; using pll = pair; using pil = pair; using pli = pair; using vi = vector; using vvi = vector; using vvvi = vector; using vvvvi = vector; using vl = vector; using vvl = vector; using vvvl = vector; using vvvvl = vector; using vb = vector; using vvb = vector; using vvvb = vector; using vc = vector; using vvc = vector; using vvvc = vector; using vd = vector; using vvd = vector; using vvvd = vector; template using priority_queue_rev = priority_queue, greater>; using Graph = vvi; // 定数の定義 const double PI = acos(-1); const vi DX = { 1, 0, -1, 0 }; // 4 近傍(下,右,上,左) const vi DY = { 0, 1, 0, -1 }; int INF = 1001001001; ll INFL = 4004004003104004004LL; // (int)INFL = 1010931620; // 入出力高速化 struct fast_io { fast_io() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(18); } } fastIOtmp; // 汎用マクロの定義 #define all(a) (a).begin(), (a).end() #define sz(x) ((int)(x).size()) #define lbpos(a, x) (int)distance((a).begin(), std::lower_bound(all(a), x)) #define ubpos(a, x) (int)distance((a).begin(), std::upper_bound(all(a), x)) #define Yes(b) {cout << ((b) ? "Yes\n" : "No\n");} #define YES(b) {cout << ((b) ? "YES\n" : "NO\n");} #define rep(i, n) for(int i = 0, i##_len = int(n); i < i##_len; ++i) // 0 から n-1 まで昇順 #define repi(i, s, t) for(int i = int(s), i##_end = int(t); i <= i##_end; ++i) // s から t まで昇順 #define repir(i, s, t) for(int i = int(s), i##_end = int(t); i >= i##_end; --i) // s から t まで降順 #define repe(v, a) for(const auto& v : (a)) // a の全要素(変更不可能) #define repea(v, a) for(auto& v : (a)) // a の全要素(変更可能) #define repb(set, d) for(int set = 0, set##_ub = 1 << int(d); set < set##_ub; ++set) // d ビット全探索(昇順) #define repis(i, set) for(int i = lsb(set), bset##i = set; i >= 0; bset##i -= 1 << i, i = lsb(bset##i)) // set の全要素(昇順) #define repp(a) sort(all(a)); for(bool a##_perm = true; a##_perm; a##_perm = next_permutation(all(a))) // a の順列全て(昇順) #define smod(n, m) ((((n) % (m)) + (m)) % (m)) // 非負mod #define uniq(a) {sort(all(a)); (a).erase(unique(all(a)), (a).end());} // 重複除去 #define EXIT(a) {cout << (a) << endl; exit(0);} // 強制終了 #define inQ(x, y, u, l, d, r) ((u) <= (x) && (l) <= (y) && (x) < (d) && (y) < (r)) // 矩形内判定 // 汎用関数の定義 template inline ll powi(T n, int k) { ll v = 1; rep(i, k) v *= n; return v; } template inline bool chmax(T& M, const T& x) { if (M < x) { M = x; return true; } return false; } // 最大値を更新(更新されたら true を返す) template inline bool chmin(T& m, const T& x) { if (m > x) { m = x; return true; } return false; } // 最小値を更新(更新されたら true を返す) template inline T get(T set, int i) { return (set >> i) & T(1); } // 演算子オーバーロード template inline istream& operator>>(istream& is, pair& p) { is >> p.first >> p.second; return is; } template inline istream& operator>>(istream& is, vector& v) { repea(x, v) is >> x; return is; } template inline vector& operator--(vector& v) { repea(x, v) --x; return v; } template inline vector& operator++(vector& v) { repea(x, v) ++x; return v; } #endif // 折りたたみ用 #if __has_include() #include using namespace atcoder; #ifdef _MSC_VER #include "localACL.hpp" #endif //using mint = modint1000000007; using mint = modint998244353; //using mint = modint; // mint::set_mod(m); namespace atcoder { inline istream& operator>>(istream& is, mint& x) { ll x_; is >> x_; x = x_; return is; } inline ostream& operator<<(ostream& os, const mint& x) { os << x.val(); return os; } } using vm = vector; using vvm = vector; using vvvm = vector; using vvvvm = vector; #endif #ifdef _MSC_VER // 手元環境(Visual Studio) #include "local.hpp" #else // 提出用(gcc) inline int popcount(int n) { return __builtin_popcount(n); } inline int popcount(ll n) { return __builtin_popcountll(n); } inline int lsb(int n) { return n != 0 ? __builtin_ctz(n) : -1; } inline int lsb(ll n) { return n != 0 ? __builtin_ctzll(n) : -1; } inline int msb(int n) { return n != 0 ? (31 - __builtin_clz(n)) : -1; } inline int msb(ll n) { return n != 0 ? (63 - __builtin_clzll(n)) : -1; } #define dump(...) #define dumpel(v) #define dump_list(v) #define dump_mat(v) #define input_from_file(f) #define output_to_file(f) #define Assert(b) { if (!(b)) while (1) cout << "OLE"; } #endif //【グラフの入力】O(n + m) /* * (始点, 終点) の組からなる入力を受け取り,n 頂点 m 辺のグラフを構築して返す. * * n : グラフの頂点の数 * m : グラフの辺の数(省略すれば n-1) * undirected : 無向グラフか(省略すれば true) * one_indexed : 入力が 1-indexed か(省略すれば true) */ Graph read_Graph(int n, int m = -1, bool undirected = true, bool one_indexed = true) { // verify : https://atcoder.jp/contests/tessoku-book/tasks/tessoku_book_bi Graph g(n); if (m == -1) m = n - 1; rep(i, m) { int a, b; cin >> a >> b; if (one_indexed) { --a; --b; } g[a].push_back(b); if (undirected && a != b) g[b].push_back(a); } return g; } //【Sparse Table(最小値と位置)】(の改変) /* * Sparse_table_indexed(vS a, bool min_flag = true) : O(n log n) * 配列 a[0..n) で初期化する.min_flag = true[false] のときは最小値[最大値] を求める. * * pTi get(int l, int r) : O(1) * min a[l..r) とそれを与える位置の組を返す.(空なら {numeric_limits::max(), -1} を返す) */ template class Sparse_table_indexed { // 参考 : https://tookunn.hatenablog.com/entry/2016/07/13/211148 int n, m; // acc[j][i] : min a[i..i+2^j) vector> acc; vvi id; public: // 配列 a[0..n) で初期化する.min_flag = true[false] のときは最小値[最大値] を求める. Sparse_table_indexed(const vector& a) : n(sz(a)), m(msb(n) + 1), acc(m, vector(n)), id(m, vi(n)) { // verify : https://atcoder.jp/contests/agc026/tasks/agc026_d rep(i, n) { acc[0][i] = a[i]; id[0][i] = i; } repi(j, 1, m - 1) { int d = 1 << (j - 1); rep(i, n - d) { if (acc[j - 1][i] < acc[j - 1][i + d]) { acc[j][i] = acc[j - 1][i]; id[j][i] = id[j - 1][i]; } else { acc[j][i] = acc[j - 1][i + d]; id[j][i] = id[j - 1][i + d]; } } } } Sparse_table_indexed() : n(0), m(0) {} // min a[l..r) を返す. int get(int l, int r) const { // verify : https://atcoder.jp/contests/agc026/tasks/agc026_d int j = msb(r - l); int d = 1 << j; if (acc[j][l] < acc[j][r - d]) { return id[j][l]; } else { return id[j][r - d]; } } #ifdef _MSC_VER friend ostream& operator<<(ostream& os, const Sparse_table_indexed& st) { rep(j, st.m) { rep(i, st.n) os << st.acc[j][i] << " "; os << "\n"; } return os; } #endif }; //【オイラーツアー】(の改変) /* * Euler_tour(Graph g, int rt) : O(n) * rt を根とする根付き木 g で初期化する. * * int lca(int s, int t) : O(log n) * 頂点 s, t の最小共通祖先を返す. * * int dist(int s, int t) : O(log n) * 頂点 s, t 間の距離を返す. * * int jump(int s, int t, int i) : O(log n) * 頂点 s から t までのパスの i 番目(0-indexed)の頂点を返す(なければ -1) * * sort_by_DFS_order(vi& vs) : O(log |vs|) * 頂点集合 vs を DFS 昇順にソートする. */ struct Euler_tour { int n; // in[s] : rt からの DFS で最初に頂点 s を訪れた時刻(根なら 0) // out[s] : rt からの DFS で最後に頂点 s から離れた時刻(根なら 2n-1) // pos[t] : rt からの DFS で時刻 t に居た頂点の番号(長さ 2n-1) // dep[s] : 頂点 s の深さ vi in, out, pos, dep; // seg[t] : 時刻 t に居た頂点の (深さ, 番号) using SEG = Sparse_table_indexed; SEG seg; void dfs(const Graph& g, int rt) { int time = 0; function rf = [&](int s, int p) { // s を最初に訪れた in[s] = time; pos[time] = s; time++; repe(t, g[s]) { if (t == p) continue; dep[t] = dep[s] + 1; rf(t, s); pos[time] = s; time++; } // s から最後に離れる out[s] = time; }; // 根から順に探索する. rf(rt, -1); } public: // rt を根とする根付き木 g で初期化する. Euler_tour(const Graph& g, int rt) : n(sz(g)), in(n), out(n), pos(2 * n - 1), dep(n) { // verify : https://judge.yosupo.jp/problem/lca dfs(g, rt); vector ini(2 * n - 1); rep(t, 2 * n - 1) ini[t] = dep[pos[t]]; seg = SEG(ini); } Euler_tour() : n(0) {} // 頂点 s, t の最小共通祖先を返す. int lca(int s, int t) const { // verify : https://judge.yosupo.jp/problem/lca // 初めて s または t に訪れたとき int l = min(in[s], in[t]); // 最後に s または t から離れたとき int r = max(out[s], out[t]); // その途中で訪れたことのある最も浅い頂点が最小共通祖先 return pos[seg.get(l, r)]; } // 頂点 s, t 間の距離を返す. int dist(int s, int t) const { // verify : https://yukicoder.me/problems/no/2337 int p = lca(s, t); // 根からの距離(深さ)の和を求め,ダブっている分を引く. return dep[s] + dep[t] - 2 * dep[p]; } // 頂点集合 vs を DFS 昇順にソートする. void sort_by_DFS_order(vi& vs) { // verify : https://yukicoder.me/problems/no/2588 sort(all(vs), [&](int s, int t) { return in[s] < in[t]; }); } }; Euler_tour ET; //【木の座標圧縮】(の改変) /* * Auxiliary_tree(Graph g, int rt) : O(n) * rt を根とする根付き木 g で初期化する. * * Graph create(vi vs, vi& id) : O(k (log k + log n)) (k = |vs|) * 頂点集合 vs とそれらの LCA からなる座標圧縮された木 gc(根は 0)を構築して返す. * gc[i] は g[id[i]] と対応する. * * 利用:【オイラーツアー】 */ int id[200000]; struct Auxiliary_tree { // 参考 : https://tjkendev.github.io/procon-library/python/graph/auxiliary_tree.html public: Auxiliary_tree(const Graph& g, int rt) { // verify : https://mojacoder.app/users/Tonegawac/problems/lca_tree ET = Euler_tour(g, rt); } // 頂点集合 vs とそれらの LCA からなる座標圧縮された木 gc(根は 0)を構築して返す. // gc[i] は g[id[i]] と対応する. Graph create(vi vs) { // verify : https://mojacoder.app/users/Tonegawac/problems/lca_tree int K = sz(vs); // 頂点集合 vs をオイラーツアーの行きがけ順にソートする. ET.sort_by_DFS_order(vs); // 行きがけ順で隣り合う 2 頂点の LCA は必要なので頂点集合に追加する. rep(k, K - 1) vs.emplace_back(ET.lca(vs[k], vs[k + 1])); // LCA も含めた頂点集合 vs をオイラーツアーの行きがけ順にソートし重複を除去する. ET.sort_by_DFS_order(vs); auto it = unique(all(vs)); vs.erase(it, vs.end()); K = sz(vs); Graph gc(K); stack stk; rep(si, K) { // v = vs[si] とし,スタックトップが v の先祖になるまで走査済の頂点をポップする. while (!stk.empty() && ET.out[vs[stk.top()]] < ET.in[vs[si]]) stk.pop(); // v に先祖が居ればそれは直近の先祖であるから辺で繋ぐ. if (!stk.empty()) { gc[stk.top()].push_back(si); gc[si].push_back(stk.top()); } stk.push(si); } rep(k, K) id[k] = vs[k]; return gc; } }; //【全方位木 DP】O(n) /* * 与えられた木 g に対し,各 s∈[0..n) について, * g の頂点 s を根と見たときの問題の答えを格納したリストを返す. * また必要なら各 s∈[0..n) と s に隣接する各頂点 t(j 番目)について, * s-t 間の辺を切断し t を根と見たときの問題の答えを sub[s][j] に格納する. * * T merge(T x, T y, int s) : * 根 s のみを共有する部分木 2 つに対する答えがそれぞれ x, y のとき, * これらをマージした部分木について同じく s を根と見たときの答えを返す. * * T leaf(int s) : * 木 g の葉 s のみからなる部分木について,s を根と見たときの答えを返す. * * T apply(T x, int p, int s) : * 頂点 s を根とする部分木の暫定の答えが x のとき, * 辺 p→s を追加して p を根と見たときの答えを返す. */ template vector rerooting(const Graph& g, vector>* sub = nullptr) { // 参考 : https://atcoder.jp/contests/abc222/editorial/2749 int n = sz(g); vector res(n); // sub[s][i] : 頂点 s と接続する i 番目の頂点を t としたとき, // s-t 間の辺を切断し,t を根と見たときの答え if (sub == nullptr) sub = new vector>; sub->resize(n); rep(s, n) (*sub)[s] = vector(sz(g[s])); // 大きさ 1 の木に対する例外処理 if (n == 1) return vector{ leaf(0) }; // p-s 間の辺を切断し,s を根と見たときの答えを計算する. // p : 0 を根としたときの s の親 // si : s が p に接続する何番目の頂点か function dfs1 = [&](int s, int p, int si) { // is_leaf : s が葉か bool is_leaf = true; rep(ti, sz(g[s])) { int t = g[s][ti]; if (t == p) continue; // s-t 間の辺を切断し,t を根と見たときの答えを計算する. dfs1(t, s, ti); // 先の部分木に対して辺 s→t を接続した場合の答えを得る. T val = apply((*sub)[s][ti], s, t); // それを暫定の答えとマージして自身の答えを計算していく. if (p != -1) { if (is_leaf) (*sub)[p][si] = move(val); else (*sub)[p][si] = merge((*sub)[p][si], val, s); } is_leaf = false; } // s が葉の場合は専用の答えを代入しておく. if (is_leaf && p != -1) (*sub)[p][si] = leaf(s); }; dfs1(0, -1, -1); // s を根と見たときの答えを計算する. // p : 0 を根としたときの s の親 // val : s-p 間の辺を切断し,p を根と見たときの答え function dfs2 = [&](int s, int p, const T& val) { // K : 根 s から出る辺の数 int K = sz(g[s]); // ds[i] : 根 s から出る i 番目の辺だけを s に接続したときの答え vector ds(K); rep(ti, K) { const auto& t = g[s][ti]; if (t == p) { (*sub)[s][ti] = val; ds[ti] = apply(val, s, p); continue; } // s-t 間の辺を切断し,t を根と見たときの答えは計算し終えているので, // その部分木に対して辺 s→t を接続し s を根と見た場合の答えを得る. ds[ti] = apply((*sub)[s][ti], s, t); } // acc_l[i] : 根 s の [0..i] 番目の辺を s に接続したときの答え vector acc_l(K); acc_l[0] = ds[0]; repi(i, 1, K - 1) acc_l[i] = merge(acc_l[i - 1], ds[i], s); // acc_r[i] : 根 s の [i..K) 番目の辺を s に接続したときの答え vector acc_r(K); acc_r[K - 1] = ds[K - 1]; repir(i, K - 2, 0) acc_r[i] = merge(acc_r[i + 1], ds[i], s); // 根 s から出る全ての辺を s に接続したときの答えが求めるものである. res[s] = acc_l[K - 1]; rep(ti, K) { const auto& t = g[s][ti]; if (t == p) continue; // 根 s に辺 s→t 以外の全ての辺を接続したときの答え, // すなわち,辺 t-s を切断し,s を根と見たときの答えを再帰関数に渡す. if (K == 1) dfs2(t, s, leaf(s)); else if (ti == 0) dfs2(t, s, acc_r[1]); else if (ti == K - 1) dfs2(t, s, acc_l[K - 2]); else dfs2(t, s, merge(acc_l[ti - 1], acc_r[ti + 1], s)); } }; dfs2(0, -1, T()); // 後ろ 1 つの引数はダミー return res; /* 雛形 using T = int; T merge(T x, T y, int s) { return max(x, y); } T leaf(int s) { return 0; } T apply(T x, int p, int s) { return x + 1; } vector solve_by_rerooting(const Graph& g, vector>* sub = nullptr) { return rerooting(g, sub); } */ }; //【木の高さ】O(n)(の改変) /* * 与えられた重み付き木 g に対し,各 s∈[0..n) について * 頂点 s を根にしたときの高さ(最も遠い葉までのコスト)を格納したリストを返す. * * 利用:【全方位木 DP】 */ int col[200000]; using T_hut = int; T_hut merge_hut(T_hut x, T_hut y, int s) { return max(x, y); } T_hut leaf_hut(int s) { return col[s] == 1 ? 0 : -INF; } T_hut apply_hut(T_hut x, int p, int s) { x += ET.dist(id[p], id[s]); if (col[p] == 1) chmax(x, 0); return x; } vector height_of_undirected_tree(Graph& g) { return rerooting(g); } //【木の高さ】O(n)(の改変) /* * 与えられた重み付き木 g に対し,各 s∈[0..n) について * 頂点 s を根にしたときの高さ(最も遠い葉までのコスト)を格納したリストを返す. * * 利用:【全方位木 DP】 */ using T_hut2 = pii; T_hut2 merge_hut2(T_hut2 x, T_hut2 y, int s) { return { max(x.first, y.first), max(x.second, y.second) }; } T_hut2 leaf_hut2(int s) { return col[s] == 1 ? make_pair(-INF, 0) : make_pair(0, -INF); } T_hut2 apply_hut2(T_hut2 x, int p, int s) { int d = ET.dist(id[p], id[s]); x.first += d; x.second += d; if (col[p] == 1) chmax(x.second, 0); else chmax(x.first, 0); return x; } vector height_of_undirected_tree2(Graph& g) { return rerooting(g); } //【座標圧縮】O(n log n) /* * a[0..n) を座標圧縮した結果を a_cp[0..n) に格納し,その値域の大きさを返す. * また xs[j] に圧縮された座標 j に対応する元の座標を格納する. * * a に重複する要素がなければ,a_cp[i] は a[i] が昇順で何番目かを表し, * xs[j] は昇順で j 番目の要素が何かを表す. */ template int coordinate_compression(const vector& a, vi& a_cp, vector* xs = nullptr) { // verify : https://atcoder.jp/contests/tessoku-book/tasks/tessoku_book_o int n = sz(a); if (xs == nullptr) xs = new vector; // *xs : a の x 座標のユニークな昇順列 *xs = a; uniq(*xs); // a[i] が xs において何番目かを求める. a_cp.resize(n); rep(i, n) a_cp[i] = lbpos(*xs, a[i]); return sz(*xs); } int main() { input_from_file("input.txt"); output_to_file("output.txt"); int n; cin >> n; vi a0(n); cin >> a0; auto g = read_Graph(n); mt19937_64 mt((int)time(NULL)); uniform_int_distribution rnd(0, n - 1); Auxiliary_tree G(g, rnd(mt)); vi a; int A = coordinate_compression(a0, a); int B = msb(A); // dump(a); vector ai(n); rep(i, n) ai[i] = { a[i], i }; sort(all(ai)); vi is(n); vi a_sorted(n); rep(i, n) tie(a_sorted[i], is[i]) = ai[i]; vi res(n, -INF); auto rf2 = [&](int b, int l, int r, int lb, int ub) { vi vs(is.begin() + l, is.begin() + r); auto g2 = G.create(move(vs)); int n2 = sz(g2); rep(i2, n2) { if (a[id[i2]] < lb || ub <= a[id[i2]]) { col[i2] = -1; continue; } col[i2] = get(a[id[i2]], b); } auto hs = height_of_undirected_tree2(g2); rep(i2, n2) { if (col[i2] == 0) chmax(res[id[i2]], hs[i2].first); if (col[i2] != -1) chmax(res[id[i2]], hs[i2].second); } }; function rf = [&](int b, int l, int r, int lb, int ub) { // dump("---", b, l, r, lb, ub); vi vs(is.begin() + l, is.begin() + r); // dump(vs); auto g2 = G.create(move(vs)); int n2 = sz(g2); // dump(id); rep(i2, n2) { if (a[id[i2]] < lb || ub <= a[id[i2]]) { col[i2] = -1; continue; } col[i2] = get(a[id[i2]], b); } // dump(col); auto hs = height_of_undirected_tree(g2); // dump(hs); rep(i2, n2) { if (col[i2] == 0) chmax(res[id[i2]], hs[i2]); } // dump(res); int mb = (lb + ub) >> 1; if (b > 1) { auto it = lower_bound(a_sorted.begin() + l, a_sorted.begin() + r, mb); int m = (int)distance(a_sorted.begin(), it); if (l < m) rf(b - 1, l, m, lb, mb); if (m < r) rf(b - 1, m, r, mb, ub); } else { auto it = lower_bound(a_sorted.begin() + l, a_sorted.begin() + r, (lb + ub) / 2); int m = (int)distance(a_sorted.begin(), it); if (l < m) rf2(b - 1, l, m, lb, mb); if (m < r) rf2(b - 1, m, r, mb, ub); } }; rf(B, 0, n, 0, 1 << (B + 1)); rep(s, n) cout << res[s] << " \n"[s == n - 1]; } // 5286 ms