結果
問題 | No.1675 Strange Minimum Query |
ユーザー | iiljj |
提出日時 | 2021-09-10 21:51:09 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 618 ms / 2,000 ms |
コード長 | 26,360 bytes |
コンパイル時間 | 2,330 ms |
コンパイル使用メモリ | 155,792 KB |
実行使用メモリ | 48,000 KB |
最終ジャッジ日時 | 2024-06-11 23:31:13 |
合計ジャッジ時間 | 15,621 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 68 ms
12,032 KB |
testcase_04 | AC | 77 ms
25,600 KB |
testcase_05 | AC | 27 ms
22,912 KB |
testcase_06 | AC | 77 ms
12,672 KB |
testcase_07 | AC | 86 ms
26,752 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 4 ms
5,376 KB |
testcase_10 | AC | 310 ms
44,128 KB |
testcase_11 | AC | 119 ms
42,328 KB |
testcase_12 | AC | 350 ms
44,672 KB |
testcase_13 | AC | 247 ms
47,940 KB |
testcase_14 | AC | 457 ms
47,968 KB |
testcase_15 | AC | 409 ms
48,000 KB |
testcase_16 | AC | 2 ms
5,376 KB |
testcase_17 | AC | 83 ms
13,952 KB |
testcase_18 | AC | 111 ms
14,464 KB |
testcase_19 | AC | 211 ms
43,904 KB |
testcase_20 | AC | 359 ms
45,312 KB |
testcase_21 | AC | 344 ms
45,420 KB |
testcase_22 | AC | 439 ms
46,488 KB |
testcase_23 | AC | 258 ms
16,768 KB |
testcase_24 | AC | 292 ms
26,112 KB |
testcase_25 | AC | 181 ms
15,488 KB |
testcase_26 | AC | 125 ms
23,808 KB |
testcase_27 | AC | 313 ms
45,196 KB |
testcase_28 | AC | 193 ms
43,244 KB |
testcase_29 | AC | 167 ms
43,264 KB |
testcase_30 | AC | 129 ms
23,680 KB |
testcase_31 | AC | 341 ms
26,880 KB |
testcase_32 | AC | 541 ms
47,928 KB |
testcase_33 | AC | 549 ms
47,924 KB |
testcase_34 | AC | 537 ms
48,000 KB |
testcase_35 | AC | 609 ms
47,976 KB |
testcase_36 | AC | 618 ms
47,872 KB |
ソースコード
/* #region Head */ // #include <bits/stdc++.h> #include <algorithm> #include <array> #include <bitset> #include <cassert> // assert.h #include <cmath> // math.h #include <cstring> #include <ctime> #include <deque> #include <fstream> #include <functional> #include <iomanip> #include <iostream> #include <list> #include <map> #include <memory> #include <numeric> #include <queue> #include <random> #include <set> #include <sstream> #include <stack> #include <string> #include <unordered_map> #include <unordered_set> #include <vector> using namespace std; using ll = long long; using ull = unsigned long long; using ld = long double; using pll = pair<ll, ll>; template <class T> using vc = vector<T>; template <class T> using vvc = vc<vc<T>>; using vll = vc<ll>; using vvll = vvc<ll>; using vld = vc<ld>; using vvld = vvc<ld>; using vs = vc<string>; using vvs = vvc<string>; template <class T, class U> using um = unordered_map<T, U>; template <class T> using pq = priority_queue<T>; template <class T> using pqa = priority_queue<T, vc<T>, greater<T>>; template <class T> using us = unordered_set<T>; #define TREP(T, i, m, n) for (T i = (m), i##_len = (T)(n); i < i##_len; ++(i)) #define TREPM(T, i, m, n) for (T i = (m), i##_max = (T)(n); i <= i##_max; ++(i)) #define TREPR(T, i, m, n) for (T i = (m), i##_min = (T)(n); i >= i##_min; --(i)) #define TREPD(T, i, m, n, d) for (T i = (m), i##_len = (T)(n); i < i##_len; i += (d)) #define TREPMD(T, i, m, n, d) for (T i = (m), i##_max = (T)(n); i <= i##_max; i += (d)) #define REP(i, m, n) for (ll i = (m), i##_len = (ll)(n); i < i##_len; ++(i)) #define REPM(i, m, n) for (ll i = (m), i##_max = (ll)(n); i <= i##_max; ++(i)) #define REPR(i, m, n) for (ll i = (m), i##_min = (ll)(n); i >= i##_min; --(i)) #define REPD(i, m, n, d) for (ll i = (m), i##_len = (ll)(n); i < i##_len; i += (d)) #define REPMD(i, m, n, d) for (ll i = (m), i##_max = (ll)(n); i <= i##_max; i += (d)) #define REPI(itr, ds) for (auto itr = ds.begin(); itr != ds.end(); itr++) #define REPIR(itr, ds) for (auto itr = ds.rbegin(); itr != ds.rend(); itr++) #define ALL(x) begin(x), end(x) #define SIZE(x) ((ll)(x).size()) #define PERM(c) \ sort(ALL(c)); \ for (bool c##p = 1; c##p; c##p = next_permutation(ALL(c))) #define UNIQ(v) v.erase(unique(ALL(v)), v.end()); #define CEIL(a, b) (((a) + (b)-1) / (b)) #define endl '\n' constexpr ll INF = 1'010'000'000'000'000'017LL; constexpr int IINF = 1'000'000'007LL; constexpr ll MOD = 1'000'000'007LL; // 1e9 + 7 // constexpr ll MOD = 998244353; constexpr ld EPS = 1e-12; constexpr ld PI = 3.14159265358979323846; template <typename T> istream &operator>>(istream &is, vc<T> &vec) { // vector 入力 for (T &x : vec) is >> x; return is; } template <typename T> ostream &operator<<(ostream &os, const vc<T> &vec) { // vector 出力 (for dump) os << "{"; REP(i, 0, SIZE(vec)) os << vec[i] << (i == i_len - 1 ? "" : ", "); os << "}"; return os; } template <typename T> ostream &operator>>(ostream &os, const vc<T> &vec) { // vector 出力 (inline) REP(i, 0, SIZE(vec)) os << vec[i] << (i == i_len - 1 ? "\n" : " "); return os; } template <typename T, size_t _Nm> istream &operator>>(istream &is, array<T, _Nm> &arr) { // array 入力 REP(i, 0, SIZE(arr)) is >> arr[i]; return is; } template <typename T, size_t _Nm> ostream &operator<<(ostream &os, const array<T, _Nm> &arr) { // array 出力 (for dump) os << "{"; REP(i, 0, SIZE(arr)) os << arr[i] << (i == i_len - 1 ? "" : ", "); os << "}"; return os; } template <typename T, typename U> istream &operator>>(istream &is, pair<T, U> &pair_var) { // pair 入力 is >> pair_var.first >> pair_var.second; return is; } template <typename T, typename U> ostream &operator<<(ostream &os, const pair<T, U> &pair_var) { // pair 出力 os << "(" << pair_var.first << ", " << pair_var.second << ")"; return os; } // map, um, set, us 出力 template <class T> ostream &out_iter(ostream &os, const T &map_var) { os << "{"; REPI(itr, map_var) { os << *itr; auto itrcp = itr; if (++itrcp != map_var.end()) os << ", "; } return os << "}"; } template <typename T, typename U> ostream &operator<<(ostream &os, const map<T, U> &map_var) { return out_iter(os, map_var); } template <typename T, typename U> ostream &operator<<(ostream &os, const um<T, U> &map_var) { os << "{"; REPI(itr, map_var) { auto [key, value] = *itr; os << "(" << key << ", " << value << ")"; auto itrcp = itr; if (++itrcp != map_var.end()) os << ", "; } os << "}"; return os; } template <typename T> ostream &operator<<(ostream &os, const set<T> &set_var) { return out_iter(os, set_var); } template <typename T> ostream &operator<<(ostream &os, const us<T> &set_var) { return out_iter(os, set_var); } template <typename T> ostream &operator<<(ostream &os, const pq<T> &pq_var) { pq<T> pq_cp(pq_var); os << "{"; if (!pq_cp.empty()) { os << pq_cp.top(), pq_cp.pop(); while (!pq_cp.empty()) os << ", " << pq_cp.top(), pq_cp.pop(); } return os << "}"; } // tuple 出力 template <size_t N = 0, bool end_line = false, typename... Args> ostream &operator<<(ostream &os, tuple<Args...> &a) { if constexpr (N < std::tuple_size_v<tuple<Args...>>) { os << get<N>(a); if constexpr (N + 1 < std::tuple_size_v<tuple<Args...>>) { os << ' '; } else if constexpr (end_line) { os << '\n'; } return operator<<<N + 1, end_line>(os, a); } return os; } template <typename... Args> void print_tuple(tuple<Args...> &a) { operator<<<0, true>(cout, a); } void pprint() { cout << endl; } template <class Head, class... Tail> void pprint(Head &&head, Tail &&...tail) { cout << head; if (sizeof...(Tail) > 0) cout << ' '; pprint(move(tail)...); } // dump #define DUMPOUT cerr void dump_func() { DUMPOUT << endl; } template <class Head, class... Tail> void dump_func(Head &&head, Tail &&...tail) { DUMPOUT << head; if (sizeof...(Tail) > 0) DUMPOUT << ", "; dump_func(move(tail)...); } // chmax (更新「される」かもしれない値が前) template <typename T, typename U, typename Comp = less<>> bool chmax(T &xmax, const U &x, Comp comp = {}) { if (comp(xmax, x)) { xmax = x; return true; } return false; } // chmin (更新「される」かもしれない値が前) template <typename T, typename U, typename Comp = less<>> bool chmin(T &xmin, const U &x, Comp comp = {}) { if (comp(x, xmin)) { xmin = x; return true; } return false; } // ローカル用 #ifndef ONLINE_JUDGE #define DEBUG_ #endif #ifdef DEBUG_ #define DEB #define dump(...) \ DUMPOUT << " " << string(#__VA_ARGS__) << ": " \ << "[" << to_string(__LINE__) << ":" << __FUNCTION__ << "]" << endl \ << " ", \ dump_func(__VA_ARGS__) #else #define DEB if (false) #define dump(...) #endif #define VAR(type, ...) \ type __VA_ARGS__; \ cin >> __VA_ARGS__; template <typename T> istream &operator,(istream &is, T &rhs) { return is >> rhs; } template <typename T> ostream &operator,(ostream &os, const T &rhs) { return os << ' ' << rhs; } struct AtCoderInitialize { static constexpr int IOS_PREC = 15; static constexpr bool AUTOFLUSH = false; AtCoderInitialize() { ios_base::sync_with_stdio(false), cin.tie(nullptr), cout.tie(nullptr); cout << fixed << setprecision(IOS_PREC); if (AUTOFLUSH) cout << unitbuf; } } ATCODER_INITIALIZE; void Yn(bool p) { cout << (p ? "Yes" : "No") << endl; } void YN(bool p) { cout << (p ? "YES" : "NO") << endl; } /* #endregion */ // #include <atcoder/all> // using namespace atcoder; /* #region SegmentTreeBeats */ template <class T> class SegmentTreeBeats { T inf, zero; size_t nn; // 外から見たサイズ size_t length; // 葉ノードの個数(2べきになるようにパディング) vc<T> node_max_first; // ノードに対応する区間の最大値 vc<T> node_max_second; // ノードに対応する区間の二番目の最大値,区間内が最大値のみのときは -inf vc<int> count_max_first; // ノードに対応する区間の最大値の個数 vc<T> node_min_first; // ノードに対応する区間の最小値 vc<T> node_min_second; // ノードに対応する区間の二番目の最小値,区間内が最大値のみのときは inf vc<int> count_min_first; // ノードに対応する区間の最小値の個数 vc<T> node_sum; // ノードに対応する区間の総和 vc<T> lazy_add, lazy_update; // ノードに対応する遅延配列 vc<pair<int, int>> range; // ノードに対応する区間の左端・右端(半開区間) stack<int> down, up; // クエリを処理するためのスタック.インスタンスを使い回すためにここに置く. // 木のノード k に対する区間 chmin 更新に応じて,区間最大値を更新する(区間和も更新) // 区間に含まれる x より大きい要素は,最大値の要素のみである前提 // 更新前後の差分(負数)を返す inline T internal_chmax(const int k, const T x) noexcept { T delta = (x - node_max_first[k]) * count_max_first[k]; node_sum[k] += delta; // これから減る分を反映 if (node_max_first[k] == node_min_first[k]) // 区間内の値が全て同じ値である(つまり x より大きい)とき node_max_first[k] = node_min_first[k] = x; // 新しい最大値・最小値が x になる else if (node_max_first[k] == node_min_second[k]) // 区間内の値が 2 種類である(大きい方は x より大きい)とき node_max_first[k] = node_min_second[k] = x; // 2種類のうち大きい方の値が x になる else // 区間内の値が 3 種類以上あるとき node_max_first[k] = x; // 最大値のみが x になる // 必要に応じて,遅延配列にも chmin をかける if (lazy_update[k] != inf && x < lazy_update[k]) { delta += (x - lazy_update[k]); lazy_update[k] = x; } // assert(lazy_add[k] == zero); // 加算については先に伝播処理が終わっている前提 return delta; } // 木のノード k に対する区間 chmax 更新に応じて,区間最小値を更新する(区間和も更新) // 区間に含まれる x 未満の要素は,最小値の要素のみである前提 inline T internal_chmin(const int k, const T x) noexcept { T delta = (x - node_min_first[k]) * count_min_first[k]; node_sum[k] += delta; // これから増える分を反映 if (node_max_first[k] == node_min_first[k]) // 区間内の値が全て同じ値(つまり x 未満)のとき node_max_first[k] = node_min_first[k] = x; // 新しい最大値・最小値が x になる else if (node_max_second[k] == node_min_first[k]) // 区間内の値が 2 種類(小さい方は x 未満)のとき node_min_first[k] = node_max_second[k] = x; // 2種類のうち小さい方の値が x になる else // 区間内の値が 3 種類以上あるとき node_min_first[k] = x; // 最小値のみが x になる // 必要に応じて,遅延配列にも chmax をかける if (lazy_update[k] != inf && lazy_update[k] < x) { delta += (x - lazy_update[k]); lazy_update[k] = x; } // assert(lazy_add[k] == zero); // 加算については先に伝播処理が終わっている前提 return delta; } // 木のノード k に対する加算処理(と区間和の更新) inline void internal_add(const int k, const T x) noexcept { node_max_first[k] += x; if (node_max_second[k] != -inf) node_max_second[k] += x; node_min_first[k] += x; if (node_min_second[k] != inf) node_min_second[k] += x; node_sum[k] += x * (range[k].second - range[k].first); (lazy_update[k] != inf ? lazy_update[k] : lazy_add[k]) += x; } // 木のノード k の値の更新(と区間和の更新) inline T internal_update(const int k, const T x) noexcept { int cnt = (range[k].second - range[k].first); T delta = x * cnt - node_sum[k]; node_max_first[k] = x; node_max_second[k] = -inf; node_min_first[k] = x; node_min_second[k] = inf; count_max_first[k] = count_min_first[k] = cnt; node_sum[k] = x * cnt; lazy_update[k] = x; lazy_add[k] = zero; return delta; } // 木のノード k に関して,子ノードへの伝播処理を行う inline void propagate(const int k) noexcept { if ((int)length - 1 <= k) return; if (lazy_update[k] != inf) { internal_update((k << 1) | 0, lazy_update[k]); internal_update((k << 1) | 1, lazy_update[k]); lazy_update[k] = inf; return; } if (lazy_add[k] != zero) { internal_add((k << 1) | 0, lazy_add[k]); internal_add((k << 1) | 1, lazy_add[k]); lazy_add[k] = zero; } if (node_max_first[k] < node_max_first[(k << 1) | 0]) { internal_chmax((k << 1) | 0, node_max_first[k]); } if (node_min_first[(k << 1) | 0] < node_min_first[k]) { internal_chmin((k << 1) | 0, node_min_first[k]); } if (node_max_first[k] < node_max_first[(k << 1) | 1]) { internal_chmax((k << 1) | 1, node_max_first[k]); } if (node_min_first[(k << 1) | 1] < node_min_first[k]) { internal_chmin((k << 1) | 1, node_min_first[k]); } } // ノード k について,子ノードの情報をマージする inline void merge(const int k) noexcept { // 和の更新 node_sum[k] = node_sum[(k << 1) | 0] + node_sum[(k << 1) | 1]; // 最大値の更新 if (node_max_first[(k << 1) | 0] < node_max_first[(k << 1) | 1]) { // 右側の方が大きい node_max_first[k] = node_max_first[(k << 1) | 1]; count_max_first[k] = count_max_first[(k << 1) | 1]; node_max_second[k] = max(node_max_first[(k << 1) | 0], node_max_second[(k << 1) | 1]); } else if (node_max_first[(k << 1) | 0] > node_max_first[(k << 1) | 1]) { // 左側の方が大きい node_max_first[k] = node_max_first[(k << 1) | 0]; count_max_first[k] = count_max_first[(k << 1) | 0]; node_max_second[k] = max(node_max_second[(k << 1) | 0], node_max_first[(k << 1) | 1]); } else { // 大きさが同じ node_max_first[k] = node_max_first[(k << 1) | 0]; count_max_first[k] = count_max_first[(k << 1) | 0] + count_max_first[(k << 1) | 1]; node_max_second[k] = max(node_max_second[(k << 1) | 0], node_max_second[(k << 1) | 1]); } // 最小値の更新 if (node_min_first[(k << 1) | 0] < node_min_first[(k << 1) | 1]) { // 左側の方が小さい node_min_first[k] = node_min_first[(k << 1) | 0]; count_min_first[k] = count_min_first[(k << 1) | 0]; node_min_second[k] = min(node_min_second[(k << 1) | 0], node_min_first[(k << 1) | 1]); } else if (node_min_first[(k << 1) | 0] > node_min_first[(k << 1) | 1]) { // 右側の方が小さい node_min_first[k] = node_min_first[(k << 1) | 1]; count_min_first[k] = count_min_first[(k << 1) | 1]; node_min_second[k] = min(node_min_first[(k << 1) | 0], node_min_second[(k << 1) | 1]); } else { // 小ささが同じ node_min_first[k] = node_min_first[(k << 1) | 0]; count_min_first[k] = count_min_first[(k << 1) | 0] + count_min_first[(k << 1) | 1]; node_min_second[k] = min(node_min_second[(k << 1) | 0], node_min_second[(k << 1) | 1]); } } // up スタックに記録された木のノードに順番にマージ処理をかける // up スタックからは木の葉側から順にノードが取り出される前提 inline void up_merge() noexcept { while (!up.empty()) { merge(up.top()); up.pop(); } } // 木のノード k に伝播処理をかけて,子ノードを処理対象に追加する // 伝播が必要な場合のみ呼ばれるので,子ノードは必ず存在する前提 inline void down_propagate(const int k) noexcept { propagate(k); down.push((k << 1) | 0); down.push((k << 1) | 1); } public: // コンストラクタ SegmentTreeBeats(const int n_, const T inf = INF, const T zero = 0) noexcept : inf(inf), zero(zero) { init(n_); } // コンストラクタ SegmentTreeBeats(const vector<T> &a, const T inf = INF, const T zero = 0) noexcept : inf(inf), zero(zero) { build(a); } // 指定要素数で初期化 void init(const size_t n_) noexcept { nn = n_; length = 1; while (length <= n_) length <<= 1; node_max_first.resize(2 * length); node_max_second.resize(2 * length); count_max_first.resize(2 * length); node_min_first.resize(2 * length); node_min_second.resize(2 * length); count_min_first.resize(2 * length); node_sum.assign(2 * length, zero); lazy_add.assign(2 * length, zero); lazy_update.resize(2 * length, inf); range.resize(2 * length); REP(i, 0, length) range[i + length] = {i, i + 1}; REPR(i, length - 1, 0) range[i] = {range[(i << 1) + 0].first, range[(i << 1) + 1].second}; // 実際に要素がある区間 REP(i, 0, n_) { node_max_first[length + i] = node_min_first[length + i] = zero; node_max_second[length + i] = -inf; node_min_second[length + i] = inf; count_max_first[length + i] = count_min_first[length + i] = 1; } // 2冪にするためにパディングした区間 REP(i, n_, length) { node_max_first[length + i] = node_max_second[length + i] = -inf; node_min_first[length + i] = node_min_second[length + i] = inf; count_max_first[length + i] = count_min_first[length + i] = 0; } // 上へ順にマージ REPR(i, length - 1, 1) merge(i); } // 指定ベクトルで初期化 void build(const vector<T> &a) noexcept { init(a.size()); // 実際に要素がある区間だけ更新 REP(i, 0, nn) node_max_first[length + i] = node_min_first[length + i] = node_sum[length + i] = a[i]; // 上へ順にマージ REPR(i, length - 1, 1) merge(i); } // 区間 chmin,半開区間 // x より大きい要素を x に置き換える inline T range_chmin(const int a, const int b, const T x) noexcept { T delta = zero; down.push(1); while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていないか,更新が不要な場合(x 以下の値しか区間に含まれない場合)は飛ばす if (b <= range[k].first || range[k].second <= a || node_max_first[k] <= x) continue; if (a <= range[k].first && range[k].second <= b && node_max_second[k] < x) { // ノード k の区間がクエリ区間に内包されていて,最大値以外は x よりも小さい場合 delta += internal_chmax(k, x); } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); up.push(k); } } up_merge(); return delta; } // 区間 chmax,半開区間 // x 未満の要素を x に置き換える inline T range_chmax(const int a, const int b, const T x) noexcept { T delta = zero; down.push(1); while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていないか,更新が不要な場合(x 以上の値しか区間に含まれない場合)は飛ばす if (b <= range[k].first || range[k].second <= a || x <= node_min_first[k]) continue; if (a <= range[k].first && range[k].second <= b && x < node_min_second[k]) { // ノード k の区間がクエリ区間に内包されていて,最小値以外は x よりも大きい場合 delta += internal_chmin(k, x); } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); up.push(k); } } up_merge(); return delta; } // 区間 add,半開区間 inline void range_add(const int a, const int b, const T x) noexcept { down.push(1); while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていない場合は飛ばす if (b <= range[k].first || range[k].second <= a) continue; if (a <= range[k].first && range[k].second <= b) { // ノード k の区間がクエリ区間に内包されている場合 internal_add(k, x); } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); up.push(k); } } up_merge(); } // 区間 update,半開区間 inline void range_update(const int a, const int b, const T x) noexcept { T delta = zero; down.push(1); while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていない場合は飛ばす if (b <= range[k].first || range[k].second <= a) continue; if (a <= range[k].first && range[k].second <= b) { // ノード k の区間がクエリ区間に内包されている場合 delta += internal_update(k, x); } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); up.push(k); } } up_merge(); return delta; } // 区間最大値取得,半開区間 inline T get_max(const int a, const int b) noexcept { down.push(1); T v = -inf; while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていない場合は飛ばす if (b <= range[k].first || range[k].second <= a) continue; if (a <= range[k].first && range[k].second <= b) { // ノード k の区間がクエリ区間に内包されている場合 v = max(v, node_max_first[k]); } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); } } return v; } // 区間最小値取得,半開区間 inline T get_min(const int a, const int b) noexcept { down.push(1); T v = inf; while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていない場合は飛ばす if (b <= range[k].first || range[k].second <= a) continue; if (a <= range[k].first && range[k].second <= b) { // ノード k の区間がクエリ区間に内包されている場合 v = min(v, node_min_first[k]); } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); } } return v; } // 区間和取得,半開区間 inline T get_sum(const int a, const int b) noexcept { down.push(1); T v = zero; while (!down.empty()) { int k = down.top(); down.pop(); // 区間が被っていない場合は飛ばす if (b <= range[k].first || range[k].second <= a) continue; if (a <= range[k].first && range[k].second <= b) { // ノード k の区間がクエリ区間に内包されている場合 v += node_sum[k]; } else { // 区間をさらに分割して処理する必要がある場合 down_propagate(k); } } return v; } }; /* #endregion */ // Problem void solve() { VAR(ll, n, q); vll l(q), r(q), b(q); REP(i, 0, q) { cin >> l[i], r[i], b[i]; --l[i]; } vll data(n, 1); SegmentTreeBeats<ll> seg(data); REP(i, 0, q) { seg.range_chmax(l[i], r[i], b[i]); } REP(i, 0, q) { if (seg.get_min(l[i], r[i]) != b[i]) { pprint(-1); return; } } vll ret(n); REP(i, 0, n) ret[i] = seg.get_min(i, i + 1); cout >> ret; } // entry point int main() { solve(); return 0; }