結果
問題 | No.2634 Tree Distance 3 |
ユーザー | ecottea |
提出日時 | 2024-02-18 03:55:24 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 19,256 bytes |
コンパイル時間 | 7,165 ms |
コンパイル使用メモリ | 334,936 KB |
実行使用メモリ | 149,284 KB |
最終ジャッジ日時 | 2024-09-29 00:08:17 |
合計ジャッジ時間 | 115,552 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | TLE | - |
testcase_01 | TLE | - |
testcase_02 | TLE | - |
testcase_03 | TLE | - |
testcase_04 | TLE | - |
testcase_05 | AC | 2,266 ms
114,480 KB |
testcase_06 | AC | 2,232 ms
114,052 KB |
testcase_07 | AC | 2,314 ms
114,676 KB |
testcase_08 | AC | 1,776 ms
111,364 KB |
testcase_09 | AC | 1,803 ms
111,772 KB |
testcase_10 | AC | 1,765 ms
111,296 KB |
testcase_11 | AC | 1,687 ms
111,072 KB |
testcase_12 | AC | 1,643 ms
111,248 KB |
testcase_13 | AC | 1,171 ms
111,748 KB |
testcase_14 | AC | 1,112 ms
111,160 KB |
testcase_15 | AC | 1,128 ms
111,824 KB |
testcase_16 | AC | 1,618 ms
123,364 KB |
testcase_17 | AC | 874 ms
74,332 KB |
testcase_18 | AC | 1,161 ms
95,296 KB |
testcase_19 | AC | 1,368 ms
107,612 KB |
testcase_20 | AC | 373 ms
36,404 KB |
testcase_21 | AC | 1,568 ms
135,592 KB |
testcase_22 | AC | 1,605 ms
135,680 KB |
testcase_23 | AC | 1,680 ms
149,284 KB |
testcase_24 | AC | 2,043 ms
115,724 KB |
testcase_25 | AC | 2,055 ms
115,436 KB |
testcase_26 | AC | 2,005 ms
118,588 KB |
testcase_27 | AC | 1,499 ms
117,032 KB |
testcase_28 | AC | 1,519 ms
117,928 KB |
testcase_29 | AC | 1,516 ms
112,692 KB |
testcase_30 | AC | 2,233 ms
112,800 KB |
testcase_31 | AC | 2,211 ms
114,384 KB |
testcase_32 | AC | 2,176 ms
113,336 KB |
testcase_33 | AC | 259 ms
77,428 KB |
testcase_34 | AC | 49 ms
18,460 KB |
testcase_35 | AC | 149 ms
47,292 KB |
testcase_36 | AC | 85 ms
28,388 KB |
testcase_37 | AC | 176 ms
54,552 KB |
testcase_38 | AC | 4 ms
6,816 KB |
testcase_39 | AC | 4 ms
6,820 KB |
testcase_40 | AC | 4 ms
6,816 KB |
testcase_41 | AC | 3 ms
6,820 KB |
testcase_42 | AC | 3 ms
6,816 KB |
testcase_43 | AC | 741 ms
38,296 KB |
testcase_44 | AC | 473 ms
26,460 KB |
testcase_45 | AC | 2,356 ms
95,076 KB |
testcase_46 | AC | 1,458 ms
63,420 KB |
testcase_47 | AC | 2,711 ms
106,168 KB |
testcase_48 | AC | 2,894 ms
111,848 KB |
testcase_49 | AC | 2,857 ms
112,108 KB |
testcase_50 | AC | 2,914 ms
111,996 KB |
testcase_51 | AC | 2,875 ms
111,928 KB |
testcase_52 | AC | 2,875 ms
111,736 KB |
testcase_53 | AC | 8 ms
6,820 KB |
testcase_54 | AC | 8 ms
6,816 KB |
testcase_55 | AC | 7 ms
6,820 KB |
testcase_56 | AC | 7 ms
6,816 KB |
testcase_57 | AC | 8 ms
6,816 KB |
testcase_58 | AC | 1 ms
6,816 KB |
testcase_59 | AC | 2 ms
6,820 KB |
testcase_60 | AC | 2,859 ms
111,772 KB |
testcase_61 | AC | 2,916 ms
111,804 KB |
testcase_62 | AC | 2,920 ms
111,812 KB |
testcase_63 | AC | 1,285 ms
107,900 KB |
testcase_64 | AC | 791 ms
71,640 KB |
testcase_65 | AC | 1,210 ms
101,808 KB |
testcase_66 | AC | 300 ms
32,772 KB |
testcase_67 | AC | 243 ms
27,796 KB |
testcase_68 | AC | 982 ms
111,204 KB |
testcase_69 | AC | 976 ms
111,336 KB |
testcase_70 | AC | 998 ms
111,260 KB |
ソースコード
// いろいろ高速化 #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #ifndef HIDDEN_IN_VS // 折りたたみ用 // 警告の抑制 #define _CRT_SECURE_NO_WARNINGS // ライブラリの読み込み #include <bits/stdc++.h> 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<int, int>; using pll = pair<ll, ll>; using pil = pair<int, ll>; using pli = pair<ll, int>; using vi = vector<int>; using vvi = vector<vi>; using vvvi = vector<vvi>; using vvvvi = vector<vvvi>; using vl = vector<ll>; using vvl = vector<vl>; using vvvl = vector<vvl>; using vvvvl = vector<vvvl>; using vb = vector<bool>; using vvb = vector<vb>; using vvvb = vector<vvb>; using vc = vector<char>; using vvc = vector<vc>; using vvvc = vector<vvc>; using vd = vector<double>; using vvd = vector<vd>; using vvvd = vector<vvd>; template <class T> using priority_queue_rev = priority_queue<T, vector<T>, greater<T>>; 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 <class T> inline ll powi(T n, int k) { ll v = 1; rep(i, k) v *= n; return v; } template <class T> inline bool chmax(T& M, const T& x) { if (M < x) { M = x; return true; } return false; } // 最大値を更新(更新されたら true を返す) template <class T> inline bool chmin(T& m, const T& x) { if (m > x) { m = x; return true; } return false; } // 最小値を更新(更新されたら true を返す) template <class T> inline T get(T set, int i) { return (set >> i) & T(1); } // 演算子オーバーロード template <class T, class U> inline istream& operator>>(istream& is, pair<T, U>& p) { is >> p.first >> p.second; return is; } template <class T> inline istream& operator>>(istream& is, vector<T>& v) { repea(x, v) is >> x; return is; } template <class T> inline vector<T>& operator--(vector<T>& v) { repea(x, v) --x; return v; } template <class T> inline vector<T>& operator++(vector<T>& v) { repea(x, v) ++x; return v; } #endif // 折りたたみ用 #if __has_include(<atcoder/all>) #include <atcoder/all> 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<mint>; using vvm = vector<vm>; using vvvm = vector<vvm>; using vvvvm = vector<vvvm>; #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<T>(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<T>::max(), -1} を返す) */ template <class T> 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<vector<T>> acc; vvi id; public: // 配列 a[0..n) で初期化する.min_flag = true[false] のときは最小値[最大値] を求める. Sparse_table_indexed(const vector<T>& a) : n(sz(a)), m(msb(n) + 1), acc(m, vector<T>(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<int>; SEG seg; void dfs(const Graph& g, int rt) { int time = 0; function<void(int, int)> 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<int> 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<int> 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 <class T, T(*merge)(T, T, int), T(*leaf)(int), T(*apply)(T, int, int)> vector<T> rerooting(const Graph& g) { // 参考 : https://atcoder.jp/contests/abc222/editorial/2749 int n = sz(g); vector<T> res(n); // sub[s][i] : 頂点 s と接続する i 番目の頂点を t としたとき, // s-t 間の辺を切断し,t を根と見たときの答え vector<vector<T>> sub(n); rep(s, n) sub[s].resize(sz(g[s])); // 大きさ 1 の木に対する例外処理 if (n == 1) return vector<T>{ leaf(0) }; // p-s 間の辺を切断し,s を根と見たときの答えを計算する. // p : 0 を根としたときの s の親 // si : s が p に接続する何番目の頂点か function<void(int, int, int)> 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<void(int, int, const T&)> dfs2 = [&](int s, int p, const T& val) { // K : 根 s から出る辺の数 int K = sz(g[s]); // ds[i] : 根 s から出る i 番目の辺だけを s に接続したときの答え vector<T> 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<T> 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<T> 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<T> solve_by_rerooting(const Graph& g, vector<vector<T>>* sub = nullptr) { return rerooting<T, merge, leaf, apply>(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<T_hut> height_of_undirected_tree(Graph& g) { return rerooting<T_hut, merge_hut, leaf_hut, apply_hut>(g); } int main() { // input_from_file("input.txt"); // output_to_file("output.txt"); int n; cin >> n; vi a(n); cin >> a; auto g = read_Graph(n); mt19937_64 mt((int)time(NULL)); uniform_int_distribution<int> rnd(0, n - 1); Auxiliary_tree G(g, rnd(mt)); vector<pii> ai(n); rep(i, n) ai[i] = { a[i], i }; sort(all(ai)); vi is(n); vi as(n); rep(i, n) tie(as[i], is[i]) = ai[i]; vi res(n, 0); function<void(int, int, int, int)> rf = [&](int l, int r, int lb, int ub) { dump("----", l, r, lb, ub, "----"); if (r - l <= 16) { repi(i1, l, r - 1) repi(i2, i1 + 1, r - 1) { int d = ET.dist(is[i1], is[i2]); if (as[i1] <= as[i2]) chmax(res[is[i1]], d); if (as[i1] >= as[i2]) chmax(res[is[i2]], d); } return; } int mb = as[(l + r) / 2]; int m1 = (int)distance(as.begin(), lower_bound(as.begin() + l, as.begin() + r, mb)); int m2 = (int)distance(as.begin(), upper_bound(as.begin() + l, as.begin() + r, mb)); int m; if (m1 - l < r - m2) { m = m2; mb++; } else { m = m1; } dump(mb, m); vi vs(is.begin() + l, is.begin() + r); dump(vs); 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] = a[id[i2]] >= mb; } auto hs = height_of_undirected_tree(g2); dump(hs); if (l < m && m < r) { rep(i2, n2) { if (col[i2] == 0) chmax(res[id[i2]], hs[i2]); } dump(res); rf(l, m, lb, mb); rf(m, r, mb, ub); } else { rep(i2, n2) { if (col[i2] != -1) chmax(res[id[i2]], hs[i2]); } } }; rf(0, n, 0, as.back() + 1); rep(s, n) cout << res[s] << " \n"[s == n - 1]; } // 5286 ms