#line 2 "/home/kowerkoint/workspace/CompPro-Make/library/KowerKoint/stl-expansion.hpp" #include template std::istream& operator>>(std::istream& is, std::pair& p) { is >> p.first >> p.second; return is; } template std::istream& operator>>(std::istream& is, std::array& a) { for (size_t i = 0; i < N; ++i) { is >> a[i]; } return is; } template std::istream& operator>>(std::istream& is, std::vector& v) { for (auto& e : v) is >> e; return is; } template std::ostream& operator<<(std::ostream& os, const std::pair& p) { os << p.first << " " << p.second; return os; } template std::ostream& operator<<(std::ostream& os, const std::array& a) { for (size_t i = 0; i < N; ++i) { os << a[i] << (i + 1 == a.size() ? "" : " "); } return os; } template std::ostream& operator<<(std::ostream& os, const std::vector& v) { for (size_t i = 0; i < v.size(); ++i) { os << v[i] << (i + 1 == v.size() ? "" : " "); } return os; } #line 3 "/home/kowerkoint/workspace/CompPro-Make/library/KowerKoint/base.hpp" using namespace std; #define REP(i, n) for(int i = 0; i < (int)(n); i++) #define FOR(i, a, b) for(ll i = a; i < (ll)(b); i++) #define ALL(a) (a).begin(),(a).end() #define RALL(a) (a).rbegin(),(a).rend() #define END(...) { print(__VA_ARGS__); return; } using VI = vector; using VVI = vector; using VVVI = vector; using ll = long long; using VL = vector; using VVL = vector; using VVVL = vector; using ull = unsigned long long; using VUL = vector; using VVUL = vector; using VVVUL = vector; using VD = vector; using VVD = vector; using VVVD = vector; using VS = vector; using VVS = vector; using VVVS = vector; using VC = vector; using VVC = vector; using VVVC = vector; using P = pair; using VP = vector

; using VVP = vector; using VVVP = vector; using LP = pair; using VLP = vector; using VVLP = vector; using VVVLP = vector; template using PQ = priority_queue; template using GPQ = priority_queue, greater>; constexpr int INF = 1001001001; constexpr ll LINF = 1001001001001001001ll; constexpr int DX[] = {1, 0, -1, 0}; constexpr int DY[] = {0, 1, 0, -1}; void print() { cout << '\n'; } template void print(const T &t) { cout << t << '\n'; } template void print(const Head &head, const Tail &... tail) { cout << head << ' '; print(tail...); } #ifdef DEBUG void dbg() { cerr << '\n'; } template void dbg(const T &t) { cerr << t << '\n'; } template void dbg(const Head &head, const Tail &... tail) { cerr << head << ' '; dbg(tail...); } #else template void dbg(const Args &... args) {} #endif template vector> split(typename vector::const_iterator begin, typename vector::const_iterator end, T val) { vector> res; vector cur; for(auto it = begin; it != end; it++) { if(*it == val) { res.push_back(cur); cur.clear(); } else cur.push_back(*it); } res.push_back(cur); return res; } vector split(typename string::const_iterator begin, typename string::const_iterator end, char val) { vector res; string cur = ""; for(auto it = begin; it != end; it++) { if(*it == val) { res.push_back(cur); cur.clear(); } else cur.push_back(*it); } res.push_back(cur); return res; } template< typename T1, typename T2 > inline bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); } template< typename T1, typename T2 > inline bool chmin(T1 &a, T2 b) { return a > b && (a = b, true); } template pair> compress(const vector &a) { int n = a.size(); vector x; REP(i, n) x.push_back(a[i]); sort(ALL(x)); x.erase(unique(ALL(x)), x.end()); VI res(n); REP(i, n) res[i] = lower_bound(ALL(x), a[i]) - x.begin(); return make_pair(res, x); } template auto rle(It begin, It end) { vector> res; if(begin == end) return res; auto pre = *begin; int num = 1; for(auto it = begin + 1; it != end; it++) { if(pre != *it) { res.emplace_back(pre, num); pre = *it; num = 1; } else num++; } res.emplace_back(pre, num); return res; } template vector> rle_sort(It begin, It end) { vector cloned(begin, end); sort(ALL(cloned)); auto e = rle(ALL(cloned)); sort(ALL(e), [](const auto& l, const auto& r) { return l.second < r.second; }); return e; } template pair, vector> factorial(int n) { vector res(n+1), rev(n+1); res[0] = 1; REP(i, n) res[i+1] = res[i] * (i+1); rev[n] = 1 / res[n]; for(int i = n; i > 0; i--) { rev[i-1] = rev[i] * i; } return make_pair(res, rev); } #line 2 "/home/kowerkoint/workspace/CompPro-Make/library/rainbou/cpp/lazy-segtree.hpp" /** * @file segtree.hpp * @brief 遅延伝搬セグメント木 */ #line 13 "/home/kowerkoint/workspace/CompPro-Make/library/rainbou/cpp/lazy-segtree.hpp" /** * @brief 遅延伝搬セグメント木のCRTP基底クラス * * @tparam S 値モノイドの型 * @tparam F 作用素モノイドの型 * @tparam ActualSegTree 派生クラス */ template class LazySegTreeBase { S op(const S& a, const S& b) const { return static_cast(*this).op(a, b); } S e() const { return static_cast(*this).e(); } S mapping(const F& f, const S& x, int l, int r) const { return static_cast(*this).mapping(f, x, l, r); } F composition(const F& f, const F& g) const { return static_cast(*this).composition(f, g); } F id() const { return static_cast(*this).id(); } int n, sz, height; std::vector data; std::vector lazy; void update(int k) { data[k] = op(data[2 * k], data[2 * k + 1]); } void apply_node(int k, int h, const F& f) { int l = (k << h) & (sz - 1); int r = l + (1 << h); data[k] = mapping(f, data[k], l, r); if(k < sz) lazy[k] = composition(f, lazy[k]); } void push(int k, int h) { apply_node(2 * k, h-1, lazy[k]); apply_node(2 * k + 1, h-1, lazy[k]); lazy[k] = id(); } class LazySegTreeReference { LazySegTreeBase& segtree; int k; public: LazySegTreeReference(LazySegTreeBase& segtree, int k) : segtree(segtree), k(k) {} LazySegTreeReference& operator=(const S& x) { segtree.set(k, x); return *this; } operator S() { return segtree.get(k); } }; protected: void construct_data() { sz = 1; height = 0; while (sz < n) { sz <<= 1; height++; } data.assign(sz * 2, e()); lazy.assign(sz * 2, id()); } void initialize(const std::vector& v) { for (int i = 0; i < n; i++) data[sz + i] = v[i]; for (int i = sz - 1; i > 0; i--) update(i); } public: // Warning: 継承先のコンストラクタでconstruct_data()を必ず呼び出す! LazySegTreeBase(int n = 0) : n(n) {} /** * @brief 指定された要素の値を返す * * @param k インデックス * @return S 値 */ S get(int k) { k += sz; for(int h = height; h > 0; h--) { push(k >> h, h); } return data[k]; } /** * @brief 指定された要素への参照を返す * * @param k * @return SegTreeReference 要素への参照 代入されるとset()が呼ばれる */ LazySegTreeReference operator[] (int k) { return LazySegTreeReference(*this, k); } /** * @brief 内容を出力する * * @tparam CharT 出力ストリームの文字型 * @tparam Traits 出力ストリームの文字型特性 * @param os 出力ストリーム * @param rhs セグメント木 * @return std::basic_ostream& 出力ストリーム */ template friend std::basic_ostream& operator<<(std::basic_ostream& os, LazySegTreeBase& rhs) { for(int i = 0; i < rhs.n; i++) { if(i != 0) os << CharT(' '); os << rhs[i]; } return os; } /** * @brief 指定された要素の値をxに更新する * * @param k インデックス * @param x 新しい値 */ void set(int k, const S& x) { k += sz; for(int h = height; h > 0; h--) { push(k >> h, h); } data[k] = x; while(k >>= 1) update(k); } /** * @brief [l, r)の区間の総積を返す * * @param l 半開区間の開始 * @param r 半開区間の終端 * @return S 総積 */ S prod(int l, int r) { l += sz; r += sz; for(int h = height; h > 0; h--) { if(((l >> h) << h) != l) push(l >> h, h); if(((r >> h) << h) != r) push(r >> h, h); } S left_prod = e(), right_prod = e(); while(l < r) { if(l & 1) left_prod = op(left_prod, data[l++]); if(r & 1) right_prod = op(data[--r], right_prod); l >>= 1; r >>= 1; } return op(left_prod, right_prod); } /** * @brief すべての要素の総積を返す * * @return S 総積 */ S all_prod() const { return data[1]; } /** * @brief 指定された要素の値にxを作用させる * * @param k インデックス * @param x 作用素 */ void apply(int k, const F& f) { k += sz; for(int h = height; h > 0; h--) { push(k >> h, h); } data[k] = mapping(f, data[k]); while(k >>= 1) update(k); } /** * @brief [l, r)の区間の値にxを作用させる * * @param l 半開区間の開始 * @param r 半開区間の終端 * @param f 作用素 */ void apply(int l, int r, const F& f) { if(l == r) return; l += sz; r += sz; for(int h = height; h > 0; h--) { if(((l >> h) << h) != l) push(l >> h, h); if(((r >> h) << h) != r) push(r >> h, h); } { int l2 = l, r2 = r; int h = 0; while(l < r) { if(l & 1) apply_node(l++, h, f); if(r & 1) apply_node(--r, h, f); l >>= 1; r >>= 1; h++; } l = l2; r = r2; } for(int h = 1; h <= height; h++) { if(((l >> h) << h) != l) update(l >> h); if(((r >> h) << h) != r) update((r - 1) >> h); } } /** * @brief (r = l or g(prod([l, r))) = true) and (r = n or g(prod([l, r+1))) = false)となるrを返す * gが単調なら、g(prod([l, r))) = trueとなる最大のr * * @tparam G * @param l 半開区間の開始 * @param g 判定関数 g(e) = true * @return int */ template int max_right(int l, G g) { assert(g(e())); if(l == n) return n; l += sz; for(int h = height; h > 0; h--) { push(l >> h, h); } int h = 0; while(l % 2 == 0) { l >>= 1; h++; } S sum = e(); while(g(op(sum, data[l]))) { sum = op(sum, data[l]); l++; while(l % 2 == 0) { l >>= 1; h++; } if(l == 1) return n; } while(l < sz) { push(l, h); if(!g(op(sum, data[l*2]))) { l = l*2; } else { sum = op(sum, data[l*2]); l = l*2+1; } h--; } return l - sz; } /** * @brief (l = 0 or g(prod([l, r))) = true) and (l = r or g(prod([l-1, r))) = false)となるlを返す * gが単調なら、g(prod([l, r))) = trueとなる最小のl * * @tparam G * @param r 半開区間の終端 * @param g 判定関数 g(e) = true * @return int */ template int min_left(int r, G g) { assert(g(e())); if (r == 0) return 0; r += sz; for(int h = height; h > 0; h--) { push(r >> h, h); } int h = 0; while(r % 2 == 0) { r >>= 1; h++; } S sum = e(); while(g(op(data[r-1], sum))) { sum = op(data[r-1], sum); r--; while(r % 2 == 0) { r >>= 1; h++; } if(r == 1) return 0; } while(r < sz) { push(r - 1, h); if(!g(op(data[r*2-1], sum))) r *= 2; else { sum = op(data[r*2-1], sum); r = r*2 - 1; } h--; } return r - sz; } }; /** * @brief ファンクタが静的な場合の遅延伝搬セグメント木の実装 * * @tparam S 値モノイドの型 * @tparam Op 値の積のファンクタ * @tparam E 積の単位元を返すファンクタ * @tparam F 作用素モノイドの型 * @tparam Mapping 作用素を適用するファンクタ 引数は(作用素, 値)または(作用素, 値, サイズ)(作用素, 値, 左の子, 右の子) * @tparam Composition 作用素の積のファンクタ * @tparam ID 作用素の単位元を返すファンクタ */ template class StaticLazySegTree : public LazySegTreeBase> { using BaseType = LazySegTreeBase>; inline static Op operator_object; inline static E identity_object; inline static Mapping mapping_object; inline static Composition lazy_operator_object; inline static ID lazy_identity_object; public: S op(const S& a, const S& b) const { return operator_object(a, b); } S e() const { return identity_object(); } S mapping(const F& f, const S& x, int l, int r) const { if constexpr(std::is_invocable_v) { return mapping_object(f, x, l, r); } else if constexpr(std::is_invocable_v) { return mapping_object(f, x, r - l); } else { return mapping_object(f, x); } } F composition(const F& f, const F& g) const { return lazy_operator_object(f, g); } F id() const { return lazy_identity_object(); } /** * @brief デフォルトコンストラクタ * */ StaticLazySegTree() = default; /** * @brief コンストラクタ * * @param n 要素数 */ explicit StaticLazySegTree(int n) : BaseType(n) { this->construct_data(); } /** * @brief コンストラクタ * * @param v 初期の要素 */ explicit StaticLazySegTree(const std::vector& v) : StaticLazySegTree(v.size()) { this->initialize(v); } }; /** * @brief コンストラクタで関数オブジェクトを与えることで積を定義する遅延伝搬セグメント木の実装 * テンプレート引数を省略することができ、ラムダ式が使える * * @tparam S モノイドの型 * @tparam Op 積の関数オブジェクトの型 * @tparam F 作用素モノイドの型 * @tparam Mapping 作用素を適用する関数オブジェクトの型 * @tparam Composition 作用素の積の関数オブジェクトの型 */ template class LazySegTree : public LazySegTreeBase> { using BaseType = LazySegTreeBase>; Op operator_object; S identity; Mapping mapping_object; Composition lazy_operator_object; F lazy_identity; public: S op(const S& a, const S& b) const { return operator_object(a, b); } S e() const { return identity; } S mapping(const F& f, const S& x, int l, int r) const { if constexpr(std::is_invocable_v) { return mapping_object(f, x, l, r); } else if constexpr(std::is_invocable_v) { return mapping_object(f, x, r - l); } else { return mapping_object(f, x); } } F composition(const F& f, const F& g) const { return lazy_operator_object(f, g); } F id() const { return lazy_identity; } /** * @brief デフォルトコンストラクタ */ LazySegTree() = default; /** * @brief コンストラクタ * * @param n 要素数 * @param op 積の関数オブジェクト * @param identity 単位元 * @param mapping 作用素を適用する関数オブジェクト * @param composition 作用素の積の関数オブジェクト * @param lazy_identity 作用素の単位元 */ explicit LazySegTree(int n, Op op, const S& identity, Mapping mapping, Composition composition, const F& lazy_identity) : BaseType(n), operator_object(std::move(op)), identity(identity), mapping_object(std::move(mapping)), lazy_operator_object(std::move(composition)), lazy_identity(lazy_identity) { this->construct_data(); } /** * @brief コンストラクタ * * @param v 初期の要素 * @param op 積の関数オブジェクト * @param identity 単位元 * @param mapping 作用素を適用する関数オブジェクト * @param composition 作用素の積の関数オブジェクト * @param lazy_identity 作用素の単位元 */ explicit LazySegTree(const std::vector& v, Op op, const S& identity, Mapping mapping, Composition composition, const F& lazy_identity) : LazySegTree(v.size(), std::move(op), identity, std::move(mapping), std::move(composition), lazy_identity) { this->initialize(v); } }; namespace lazy_segtree { template struct Max { const S operator() (const S& a, const S& b) const { return std::max(a, b); } }; template struct Min { const S operator() (const S& a, const S& b) const { return std::min(a, b); } }; template >* = nullptr> struct MaxLimit { constexpr S operator() () const { return std::numeric_limits::max(); } }; template >* = nullptr> struct MinLimit { constexpr S operator() () const { return std::numeric_limits::lowest(); } }; template struct AddWithSize { S operator() (const S& f, const S& x, int sz) const { return x + f * sz; } }; template struct Zero { S operator() () const { return S(0); } }; template struct Update { S operator() (const S& f, const S& x) const { return f == ID ? x : f; } }; template struct UpdateWithSize { S operator() (const S& f, const S& x, int sz) const { return f == ID ? x : f * sz; } }; template struct UpdateComposition { S operator() (const S& f, const S& g) const { return f == ID ? g : f; } }; } /** * @brief 区間最小値更新、区間最小値 * * @tparam S 型 */ template using RChminMinQ = StaticLazySegTree< S, lazy_segtree::Min, lazy_segtree::MaxLimit, S, lazy_segtree::Min, lazy_segtree::Min, lazy_segtree::MaxLimit >; /** * @brief 区間最大値更新、区間最大値 * * @tparam S 型 */ template using RChmaxMaxQ = StaticLazySegTree< S, lazy_segtree::Max, lazy_segtree::MinLimit, S, // F lazy_segtree::Max, lazy_segtree::Max, lazy_segtree::MinLimit >; /** * @brief 区間加算、区間最小値 * * @tparam S 型 */ template using RAddMinQ = StaticLazySegTree< S, // S lazy_segtree::Min, lazy_segtree::MaxLimit, S, std::plus, std::plus, lazy_segtree::Zero >; /** * @brief 区間加算、区間最大値 * * @tparam S 型 */ template using RAddMaxQ = StaticLazySegTree< S, lazy_segtree::Max, lazy_segtree::MinLimit, S, std::plus, std::plus, lazy_segtree::Zero >; /** * @brief 区間加算、区間和 * * @tparam S 型 */ template using RAddSumQ = StaticLazySegTree< S, std::plus, lazy_segtree::Zero, S, lazy_segtree::AddWithSize, std::plus, lazy_segtree::Zero >; /** * @brief 区間変更、区間最小値 * 注意: numeric_limits::max()での更新がないことをが要件 * * @tparam S 型 */ template using RUpdateMinQ = StaticLazySegTree< S, lazy_segtree::Min, lazy_segtree::MaxLimit, S, lazy_segtree::Update{}()>, lazy_segtree::UpdateComposition{}()>, lazy_segtree::MaxLimit >; /** * @brief 区間変更、区間最大値 * 注意: numeric_limits::max()での更新がないことをが要件 * * @tparam S 型 */ template using RUpdateMaxQ = StaticLazySegTree< S, lazy_segtree::Max, lazy_segtree::MinLimit, S, lazy_segtree::Update{}()>, lazy_segtree::UpdateComposition{}()>, lazy_segtree::MaxLimit >; /** * @brief 区間変更、区間和 * 注意: numeric_limits::max()での更新がないことをが要件 * * @tparam S 型 */ template using RUpdateSumQ = StaticLazySegTree< S, std::plus, lazy_segtree::Zero, S, lazy_segtree::UpdateWithSize{}()>, lazy_segtree::UpdateComposition{}()>, lazy_segtree::MaxLimit >; #line 3 "main.cpp" void solve(){ int n; cin >> n; VI a(n); cin >> a; using S = array, 4>; struct Op { S operator()(S a, S b) const { S res; REP(d, 4) REP(i, 4-d) { int j = i+d; res[i][j] = 0; if(d > 0) { chmax(res[i][j], res[i+1][j]); chmax(res[i][j], res[i][j-1]); } FOR(k, i, j+1) { chmax(res[i][j], a[i][k] + b[k][j]); } } return res; } }; struct E { S operator()() const { return {{0}}; } }; struct Mapping { S operator()(int f, S a, int len) const { if(f == -1) return a; S res; REP(i, 4) REP(j, 4) res[i][j] = i <= f && f <= j ? res[i][j] = len : res[i][j] = 0; return res; }; }; struct Composition { int operator()(int f, int g) const { if(f == -1) return g; return f; } }; struct Id { int operator()() const { return -1; } }; StaticLazySegTree seg(n); REP(i, n) seg.apply(i, i+1, a[i]-1); int q; cin >> q; while(q--) { int t, l, r; cin >> t >> l >> r; l--; if(t == 1) print(seg.prod(l, r)[0][3]); else { int x; cin >> x; seg.apply(l, r, x-1); } } } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(10); int t = 1; // cin >> t; for(int case_id = 1; case_id <= t; case_id++) solve(); return 0; }