#952917 (C++17) No.2697 Range LIS Query

提出ソース
結果

問題	No.2697 Range LIS Query
コンテスト
ユーザー	KowerKoint2010
提出日時	2024-02-19 22:37:54
言語	C++17 (gcc 13.3.0 + boost 1.89.0)
結果	AC
実行時間	508 ms / 10,000 ms
コード長	24,654 bytes
記録記録タグの例: 初AC ショートコード純ショートコード純主流ショートコード最速実行時間
コンパイル時間	2,530 ms
コンパイル使用メモリ	211,336 KB
最終ジャッジ日時	2025-02-19 17:22:11
ジャッジサーバーID （参考情報）	judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン	結果
sample	AC * 3
other	AC * 15
権限があれば一括ダウンロードができます
ソースコード

raw source code
#line 2 "/home/kowerkoint/workspace/CompPro-Make/library/KowerKoint/stl-expansion.hpp"
#include <bits/stdc++.h>

template <typename T1, typename T2>
std::istream& operator>>(std::istream& is, std::pair<T1, T2>& p) {
    is >> p.first >> p.second;
    return is;
}
template <typename T, size_t N>
std::istream& operator>>(std::istream& is, std::array<T, N>& a) {
    for (size_t i = 0; i < N; ++i) {
        is >> a[i];
    }
    return is;
}
template <typename T>
std::istream& operator>>(std::istream& is, std::vector<T>& v) {
    for (auto& e : v) is >> e;
    return is;
}
template <typename T1, typename T2>
std::ostream& operator<<(std::ostream& os, const std::pair<T1, T2>& p) {
    os << p.first << " " << p.second;
    return os;
}
template <typename T, size_t N>
std::ostream& operator<<(std::ostream& os, const std::array<T, N>& a) {
    for (size_t i = 0; i < N; ++i) {
        os << a[i] << (i + 1 == a.size() ? "" : " ");
    }
    return os;
}
template <typename T>
std::ostream& operator<<(std::ostream& os, const std::vector<T>& 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<int>;
using VVI = vector<VI>;
using VVVI = vector<VVI>;
using ll = long long;
using VL = vector<ll>;
using VVL = vector<VL>;
using VVVL = vector<VVL>;
using ull = unsigned long long;
using VUL = vector<ull>;
using VVUL = vector<VUL>;
using VVVUL = vector<VVUL>;
using VD = vector<double>;
using VVD = vector<VD>;
using VVVD = vector<VVD>;
using VS = vector<string>;
using VVS = vector<VS>;
using VVVS = vector<VVS>;
using VC = vector<char>;
using VVC = vector<VC>;
using VVVC = vector<VVC>;
using P = pair<int, int>;
using VP = vector<P>;
using VVP = vector<VP>;
using VVVP = vector<VVP>;
using LP = pair<ll, ll>;
using VLP = vector<LP>;
using VVLP = vector<VLP>;
using VVVLP = vector<VVLP>;

template <typename T>
using PQ = priority_queue<T>;
template <typename T>
using GPQ = priority_queue<T, vector<T>, greater<T>>;

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<typename T>
void print(const T &t) { cout << t << '\n'; }
template<typename Head, typename... Tail>
void print(const Head &head, const Tail &... tail) {
    cout << head << ' ';
    print(tail...);
}

#ifdef DEBUG
void dbg() { cerr << '\n'; }
template<typename T>
void dbg(const T &t) { cerr << t << '\n'; }
template<typename Head, typename... Tail>
void dbg(const Head &head, const Tail &... tail) {
    cerr << head << ' ';
    dbg(tail...);
}
#else
template<typename... Args>
void dbg(const Args &... args) {}
#endif

template<typename T>
vector<vector<T>> split(typename vector<T>::const_iterator begin, typename vector<T>::const_iterator end, T val) {
    vector<vector<T>> res;
    vector<T> 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<string> split(typename string::const_iterator begin, typename string::const_iterator end, char val) {
    vector<string> 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 <typename T>
pair<VI, vector<T>> compress(const vector<T> &a) {
    int n = a.size();
    vector<T> 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 <typename It>
auto rle(It begin, It end) {
    vector<pair<typename It::value_type, int>> 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 <typename It>
vector<pair<typename It::value_type, int>> rle_sort(It begin, It end) {
    vector<typename It::value_type> 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 <typename T>
pair<vector<T>, vector<T>> factorial(int n) {
    vector<T> 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 <typename S, typename F, typename ActualLazySegTree>
class LazySegTreeBase {
    S op(const S& a, const S& b) const { return static_cast<const ActualLazySegTree&>(*this).op(a, b); }
    S e() const { return static_cast<const ActualLazySegTree&>(*this).e(); }
    S mapping(const F& f, const S& x, int l, int r) const { return static_cast<const ActualLazySegTree&>(*this).mapping(f, x, l, r); }
    F composition(const F& f, const F& g) const { return static_cast<const ActualLazySegTree&>(*this).composition(f, g); }
    F id() const { return static_cast<const ActualLazySegTree&>(*this).id(); }

    int n, sz, height;
    std::vector<S> data;
    std::vector<F> 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<S>& 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<CharT, Traits>& 出力ストリーム 
     */
    template <class CharT, class Traits>
    friend std::basic_ostream<CharT, Traits>& operator<<(std::basic_ostream<CharT, Traits>& 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 <typename G>
    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 <typename G>
    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 <typename S, typename Op, typename E, typename F, typename Mapping, typename Composition, typename ID>
class StaticLazySegTree : public LazySegTreeBase<S, F, StaticLazySegTree<S, Op, E, F, Mapping, Composition, ID>> {
    using BaseType = LazySegTreeBase<S, F, StaticLazySegTree<S, Op, E, F, Mapping, Composition, ID>>;

    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<Mapping, F, S, int, int>) {
            return mapping_object(f, x, l, r);
        } else if constexpr(std::is_invocable_v<Mapping, F, S, int>) {
            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<S>& v) : StaticLazySegTree(v.size()) {
        this->initialize(v);
    }
};

/**
 * @brief コンストラクタで関数オブジェクトを与えることで積を定義する遅延伝搬セグメント木の実装
 * テンプレート引数を省略することができ、ラムダ式が使える
 * 
 * @tparam S モノイドの型
 * @tparam Op 積の関数オブジェクトの型
 * @tparam F 作用素モノイドの型
 * @tparam Mapping 作用素を適用する関数オブジェクトの型
 * @tparam Composition 作用素の積の関数オブジェクトの型
 */
template <typename S, typename Op, typename F, typename Mapping, typename Composition>
class LazySegTree : public LazySegTreeBase<S, F, LazySegTree<S, Op, F, Mapping, Composition>> {
    using BaseType = LazySegTreeBase<S, F, LazySegTree<S, Op, F, Mapping, Composition>>;

    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<Mapping, F, S, int, int>) {
            return mapping_object(f, x, l, r);
        } else if constexpr(std::is_invocable_v<Mapping, F, S, int>) {
            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<S>& 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 <typename S>
    struct Max {
        const S operator() (const S& a, const S& b) const { return std::max(a, b); }
    };
    template <typename S>
    struct Min {
        const S operator() (const S& a, const S& b) const { return std::min(a, b); }
    };
    template <typename S, std::enable_if_t<std::is_scalar_v<S>>* = nullptr>
    struct MaxLimit {
        constexpr S operator() () const { return std::numeric_limits<S>::max(); }
    };
    template <typename S, std::enable_if_t<std::is_scalar_v<S>>* = nullptr>
    struct MinLimit {
        constexpr S operator() () const { return std::numeric_limits<S>::lowest(); }
    };
    template <typename S>
    struct AddWithSize {
        S operator() (const S& f, const S& x, int sz) const { return x + f * sz; }
    };
    template <typename S>
    struct Zero {
        S operator() () const { return S(0); }
    };
    template <typename S, S ID>
    struct Update {
        S operator() (const S& f, const S& x) const { return f == ID ? x : f; }
    };
    template <typename S, S ID>
    struct UpdateWithSize {
        S operator() (const S& f, const S& x, int sz) const { return f == ID ? x : f * sz; }
    };
    template <typename S, S ID>
    struct UpdateComposition {
        S operator() (const S& f, const S& g) const { return f == ID ? g : f; }
    };
}

/**
 * @brief 区間最小値更新、区間最小値
 * 
 * @tparam S 型
 */
template <typename S>
using RChminMinQ = StaticLazySegTree<
    S,
    lazy_segtree::Min<S>,
    lazy_segtree::MaxLimit<S>,
    S,
    lazy_segtree::Min<S>,
    lazy_segtree::Min<S>,
    lazy_segtree::MaxLimit<S>
>;
/**
 * @brief 区間最大値更新、区間最大値
 * 
 * @tparam S 型
 */
template <typename S>
using RChmaxMaxQ = StaticLazySegTree<
    S,
    lazy_segtree::Max<S>,
    lazy_segtree::MinLimit<S>,
    S, // F
    lazy_segtree::Max<S>,
    lazy_segtree::Max<S>,
    lazy_segtree::MinLimit<S>
>;
/**
 * @brief 区間加算、区間最小値
 * 
 * @tparam S 型
 */
template <typename S>
using RAddMinQ = StaticLazySegTree<
    S, // S
    lazy_segtree::Min<S>,
    lazy_segtree::MaxLimit<S>,
    S,
    std::plus<S>,
    std::plus<S>,
    lazy_segtree::Zero<S>
>;
/**
 * @brief 区間加算、区間最大値
 * 
 * @tparam S 型
 */
template <typename S>
using RAddMaxQ = StaticLazySegTree<
    S,
    lazy_segtree::Max<S>,
    lazy_segtree::MinLimit<S>,
    S,
    std::plus<S>,
    std::plus<S>,
    lazy_segtree::Zero<S>
>;
/**
 * @brief 区間加算、区間和
 * 
 * @tparam S 型
 */
template <typename S>
using RAddSumQ = StaticLazySegTree<
    S,
    std::plus<S>,
    lazy_segtree::Zero<S>,
    S,
    lazy_segtree::AddWithSize<S>,
    std::plus<S>,
    lazy_segtree::Zero<S>
>;
/**
 * @brief 区間変更、区間最小値
 * 注意: numeric_limits<S>::max()での更新がないことをが要件
 * 
 * @tparam S 型
 */
template <typename S>
using RUpdateMinQ = StaticLazySegTree<
    S,
    lazy_segtree::Min<S>,
    lazy_segtree::MaxLimit<S>,
    S,
    lazy_segtree::Update<S, lazy_segtree::MaxLimit<S>{}()>,
    lazy_segtree::UpdateComposition<S, lazy_segtree::MaxLimit<S>{}()>,
    lazy_segtree::MaxLimit<S>
>;
/**
 * @brief 区間変更、区間最大値
 * 注意: numeric_limits<S>::max()での更新がないことをが要件
 * 
 * @tparam S 型
 */
template <typename S>
using RUpdateMaxQ = StaticLazySegTree<
    S,
    lazy_segtree::Max<S>,
    lazy_segtree::MinLimit<S>,
    S,
    lazy_segtree::Update<S, lazy_segtree::MaxLimit<S>{}()>,
    lazy_segtree::UpdateComposition<S, lazy_segtree::MaxLimit<S>{}()>,
    lazy_segtree::MaxLimit<S>
>;
/**
 * @brief 区間変更、区間和
 * 注意: numeric_limits<S>::max()での更新がないことをが要件
 * 
 * @tparam S 型
 */
template <typename S>
using RUpdateSumQ = StaticLazySegTree<
    S,
    std::plus<S>,
    lazy_segtree::Zero<S>,
    S,
    lazy_segtree::UpdateWithSize<S, lazy_segtree::MaxLimit<S>{}()>,
    lazy_segtree::UpdateComposition<S, lazy_segtree::MaxLimit<S>{}()>,
    lazy_segtree::MaxLimit<S>
>;
#line 3 "main.cpp"

void solve(){
    int n; cin >> n;
    VI a(n); cin >> a;
    using S = array<array<int, 4>, 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<S, Op, E, int, Mapping, Composition, Id> 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;
}
yukicoder

結果

ソースコード
