メモ帳がわり

個人的なメモを残します。主に競プロ

【AOJ】RSQ and RUQ

AOJの問題

区間和の取得は普通にできるが、更新は工夫が必要。
二分木は区間和を保持する。
また、遅延配列はINFで初期化しておこう。これはフラグとして利用する。(更新後の値として与えられる可能性があるため、0はフラグの値として使えない。)
更新すべき区間が子をもつ場合、(l-r)*xで更新する必要がある。これは、この区間を根とした部分木の葉は全てxに更新されるため区間の和は、その要素の数かけるxになるため。
子は、親が内包する葉の数の半分の葉をもつため、親の半分の値に更新される。

実装

#include <bits/stdc++.h>
using namespace std;

long long INF = (1LL<<31)-1;

struct RMQ_RUQ_LazySegTree {
    int n; // 葉の数
    vector<long long> node; // 完全2分木
    vector<long long> lazy; // 遅延配列

    // 初期化
    RMQ_RUQ_LazySegTree(vector<long long> v) {
        int sz = v.size();
        n = 1;
        while (n < sz) n *= 2; // 与えられた数列の項数以上の2^n個、葉を作る
        node.resize(n*2-1,  0); // 0で初期化
        lazy.resize(n*2-1, INF); // 遅延配列はINFで初期化

        for (int i = 0; i < sz; i++) node[n-1+i] = v[i]; // 葉を初期化
        // 下から順に葉以外のnodeを初期化
        for (int i = n-2; i >= 0; i--) node[i] = node[i*2+1] + node[i*2+2];
    }

    void eval(int k) {
        if (lazy[k] != INF) { // lazy[k]がINFでない場合、伝播&更新の反映
            node[k] = lazy[k];
            if (k < n-1) { // 子に伝搬
                lazy[k*2+1] = lazy[k]/2;
                lazy[k*2+2] = lazy[k]/2;
            }
            lazy[k] = INF;
        }
    }

    // クエリ処理
    // [s, t)の和を取得
    // 再帰的に探索するために呼び出す側を別の関数に
    void sum(int s, int t, long long &ans) { sum_query(s, t, 0, n, 0, ans); }

    void sum_query(int s, int t, int l, int r, int n, long long &ans) {
        eval(n); // 普通のセグメント木と違うのはここだけ
        if (r <= s || t <= l) return; // 範囲外なら終了
        // [s, t)が[l, r)を内包しているとき
        else if (s <= l && t >= r) ans += node[n];
        else {
            // (r+l)/2は区間の中心, 区間の中心を左端にするか、右端にするかで分岐する
            sum_query(s, t, l, (r+l)/2, n*2+1, ans); // 左下の子を探索
            sum_query(s, t, (r+l)/2, r, n*2+2, ans);  // 右下の子を探索
        }
    }

    // [s, t)の全ての値をxに変更
    // 再帰的に更新するために呼び出す側を別の関数に
    void update(int s, int t, long long x) { update_query(s, t, 0, n, 0, x); }

    void update_query(int s, int t, int l, int r, int n, long long x) {
        eval(n);
        if (r <= s || t <= l) return; // 範囲外なら終了
        // [s, t)が[l, r)を内包しているとき
        else if (s <= l && t >= r) {
            lazy[n] = (r - l) * x; // 子にもつ葉の数 * x に変更する
            eval(n);
        } else {
            // (r+l)/2は区間の中心, 区間の中心を左端にするか、右端にするかで分岐する
            update_query(s, t, l, (r+l)/2, n*2+1, x); // 左下の子を更新
            update_query(s, t, (r+l)/2, r, n*2+2, x);  // 右下の子を更新
            node[n] = node[n*2+1] + node[n*2+2];
        }
    }

    // デバック用
    void output() {
        for (int i = 0; i < n*2-1; i++) cout << node[i] << " ";
    }

};

int main() {
    int n, q;
    cin >> n >> q;
    vector<long long> v(n,  0); // 0で初期化
    auto segtree = RMQ_RUQ_LazySegTree(v);

    for (int i = 0; i < q; i++) {
        int com;
        cin >> com;
        if (com == 0) {
            int s, t, x;
            cin >> s >> t >> x;
            segtree.update(s, t+1, x); // 半開区間に直して呼び出し
        } else {
            int s, t;
            cin >> s >> t;
            long long ans = 0;
            segtree.sum(s, t+1, ans); // 半開区間に直して呼び出し
            cout << ans << endl;
        }
    }
}