自作セグメント木ライブラリ の使い方についての自分用のメモです． AtCoder Library にセグメント木はありますが， それができる前から使っていたものなので…

使用法

値の型を DAT, 更新演算の型を OP とする．

基本のセグメント木

作成

  auto st = make_seg_tree(unit_dat, add, init_vec);
  • unit_dat は，加法単位元
  • add には，加法の演算を行う関数を指定する． 関数ポインタ，クロージャ，関数オブジェクトが使える．
  • init_vec は初期ベクトル

初期ベクトル設定は，分けても良い:

  auto st = make_seg_tree(unit_dat, add);
  st.set_data(init_vec);

値の代入 (1点)

  st.set_single(i, x);
  # または
  st.rs(i) = x;

$i$ 番目の値として $x$ を設定する．

rs は，reference for substitution のつもり． これは，STSubst なるオブジェクトを作成して返す． STSubst オブジェクトには，代入演算子 = が再定義してあって，セグメント木の該当部分を更新するようになっている．

値の取得 (1点)

  st.at(i);

$i$ 番目の値を取得する．

値の取得 (範囲)

  DAT x = st.query(il, ir);

$il$ 以上 $ir$ 未満の値に add を適用した結果を返す．

遅延セグメント木

関数

3つの関数を指定する．

  • add は加法． $\verb!DAT! \times \verb!DAT! \to \verb!DAT!$
  • comp は更新演算の合成． $\verb!OP! \times \verb!OP! \to \verb!OP!$． comp(f, g) は，$f \circ g$ を表す．
  • apply は更新演算の適用．$\verb!OP! \times \verb!DAT! \to \verb!DAT!$．

当然，apply と add の順序交換ができる必要がある．具体的には:

  • $\text{apply}(f, \text{add}(x, y)) = \text{add}(\text{apply}(f, x), \text{apply}(f, y))$

作成

  auto st = make_seg_tree_lazy(unit_dat, unit_op, add, comp, apply, init_vec);

分けても良い:

  auto st = make_seg_tree_lazy(unit_dat, unit_op, add, comp, apply);
  st.set_data(init_vec);

値の代入 (1点)

  st.set_single(i, x);
  # または
  st.rs(i) = x;

$i$ 番目の値として $x$ を設定する．これは代入であって，OP とは関係無いことに注意．

値の更新 (範囲)

  st.update(il, ir, f);

半開区間 [il, ir) の各値 x を，apply(f, x) で更新する．

値の取得 (1点)

  st.at(i);

$i$ 番目の値を取得する．

値の取得 (範囲)

  DAT x = st.query(il, ir);

$il$ 以上 $ir$ 未満の値に add を適用した結果を返す．

設定例

1点更新，区間和取得

  auto st = make_seg_tree(0LL, plus<ll>());

1点更新，区間最小値取得

  auto st = make_seg_tree(LLONG_MAX, [](ll x, ll y) { return min(x, y); });

std::min は，オーバーロードがあるので，そのまま名前で渡すわけにはいかない．

区間加算，区間最小値取得

  auto mymin = [](ll x, ll y) -> ll { return min(x, y); };
  auto st = make_seg_tree_lazy(LLONG_MAX, 0LL, mymin, plus<ll>(), plus<ll>());
  • apply と add は自然に順序交換ができる．
  • オーバーフローしうる場合には， LLONG_MAX ではなく (ll)(1e18) などを用いる必要があるかもしれない．

区間代入，区間最小値取得

代入する値は optional<ll> で表す． nullptr が「何もしない」ことをあらわし，単位元となる．

  using OP = optional<ll>;
  auto mymin = [](ll x, ll y) -> ll { return min(x, y); };
  auto comp = [](OP f, OP g) -> OP { return f ? f : g; };
  auto appl = [](OP f, ll x) -> ll { return f ? *f : x; };
  auto st = make_seg_tree_lazy(LLONG_MAX, OP(), mymin, comp, appl);

区間代入，区間和取得

何も考えないと，add と apply の順序交換ができない． そこで， DAT として，long long の対をとり， 左に和を保持し，右には区間の長さを保持する． 右側は最初に設定したあとは変更されない． 加法演算としては，$\text{add}((v_1, l_1), (v_2, l_2)) = (v_1 + v_2, l_1 + l_2)$ を 用いる． 更新演算適用は，長さ $l$ の区間に一斉に $f$ を代入することになるので， 結果は $fl$ となる．ただし，$f$ が nullptr の時は，なにもしない．

  using DAT = pair<ll, ll>;
  using OP = optional<ll>;
  auto add = [](DAT x, DAT y) -> DAT { return DAT(x.first + y.first, x.second + y.second); };
  auto comp = [](OP f, OP g) -> OP { return f ? f : g; };
  auto appl = [](OP f, DAT x) -> DAT { return f ? DAT(*f * x.second, x.second) : x; };
  auto dat_init = vector<DAT>();
  REP(i, 0, init_vec.size()) dat_init.emplace_back(init_vec[i], 1);
  auto st = make_seg_tree_lazy(DAT(), OP(), add, comp, appl, dat_init);

区間加算，1点取得

上の区間加算・区間最小値取得で流用できる．Fp のような最小値を取れない場合には， 区間加算・区間和代入で代用しても良いが，区間の長さをとったりして面倒である． どうせモノイド計算は使わないので，次のように割り切る手もある． この場合，長さ2以上の区間に対する値は正しくない (が使わないので気にしない)．

  auto st = make_seg_tree_lazy(Fp(0), Fp(0), plus<Fp>(), plus<Fp>(), plus<Fp>());

なお，imos法ですむ場合もある．

区間代入，1点取得

上の「区間代入，区間和取得」を使っても問題無いが，どうせモノイド計算はしないので，以下でも良い．

  using DAT = ll
  using OP = optional<ll>;
  auto dummy_add = [](DAT x, DAT y) -> DAT { return 0; }
  auto comp = [](OP f, OP g) -> OP { return f ? f : g; };
  auto appl = [](OP f, DAT x) -> DAT { return f ? *f : x; };
  auto st = make_seg_tree_lazy(DAT(), OP(), dummy_add, comp, appl);

等差数列 (1次式) の代入・加算，区間和取得

  • DAT は，(和, 区間左端の添字，区間長)．
    • 単位元は，適当に (-1, -1, -1) とか決めておく．
    • モノイド加法は，$(v_1, i_1, l_1) + (v_2, i_2, l_2) = (v_1 + v_2, i_1, l_1 + l_2)$．
      • ただし，単位元を左・右に足すときは，特別扱いとする．
      • $i_1 + l_1 = i_2$ を仮定する．単位元を足すときを除いて，そういう演算しか出てこないはず．
  • OP は，$(a, b)$．
    • ただし，代入・加算する値が $v = ai + b$ ($v$ は値，$i$ は添字) であるとする．
    • 単位元
      • 加算のときには，$(0, 0)$
      • 代入のときには，適当に (INF, INF) とか決めておく．
    • 合成
      • 代入のときには，$f \circ g = f$．ただし，単位元は別扱い．
      • 加算のときには，$(a_1, b_1) \circ (a_2, b_2) = (a_1 + a_2, b_1 + b_2)$
    • apply
      • 下では，$w := 2(ai + b) + a(l - 1))l/2$ とする．
      • 代入: $\;\text{apply}((a, b), (v, i, l)) = (w, i, l)$
        • $(a, b)$ が単位元のときには別扱い．
      • 加算: $\;\text{apply}((a, b), (v, i, l)) = (w + v, i, l)$

二分探索

チェック関数は次のように書く:

  auto check = [&](DAT x) -> bool { ...; };

探索用に4つのメソッドが用意されている:

  • 範囲が広いほど check が成立しやすいときに使うメソッド:
    • st.binsearch_r_from(check, x)
      • check(st.query(x, y)) が成り立つ最小の y を返す．
      • そのような y が無いときには，init_vec.size() + 1 を返す．
    • st.binsearch_l_from(check, y)
      • check(st.query(x, y)) が成り立つ最大の x を返す．
      • そのような x が無いときには，-1 を返す．
  • 範囲が狭いほど check が成立しやすいときに使うメソッド:
    • st.binsearch_r_until(check, x)
      • check(st.query(x, y)) が成り立つ最大の y を返す．
      • そのような y が無いときには，x - 1 を返す．
    • st.binsearch_l_until(check, y)
      • check(st.query(x, y)) が成り立つ最小の x を返す．
      • そのような x が無いときには，y + 1 を返す．

デバッグ対応

vector<DAT> vec_view();   // 要素のベクトルを返す

セグメント木の vector

引数無しのコンストラクタが用意してあるので， 正の長さを持つ vector<SegTree<…» を宣言することもできる． その時に型を指定する必要があるが， make_seg_tree や make_seg_tree_lazy を使って一つインスタンスを 作った上で decltype を利用するのが現実的だと思う:

  using sg1_t = decltype(make_seg_tree_lazy(0, 0, plus<int>(), plus<int>(), plus<int>()));
  vector<sg1_t> sg1(10);

  auto sg4 = make_seg_tree(LLONG_MIN, [](ll x, ll y) -> ll { return max(x, y); });
  using sg4_t = decltype(sg4);
  vector<sg4_t> vec4(2, sg4);

ベクトルの要素は同じ型を持たなければならないので， 異なる加法関数などを用いる場合には，クロージャではうまくいかず， 関数ポインタを用いる必要があるかもしれない．

ノードと添字の対応

セグメント木自身を使用する際には不要だが，「セグメント木のような」構造を使用したいときに， 以下の関数が役に立つ可能性がある．

ノード $\to$ 添字の範囲

メソッド range_of_node(i) で，ノード i が表現する半開区間 [lo, hi) を取得できる．

  auto [lo, hi] = st.range_of_node(i);

このメソッドと同等のことを行う関数 segtree_range_of_node(ht, i) も定義されている． ht には，セグメント木の高さを指定する．init_vec のサイズを sz とすると，次のようになる:

  unsigned height = countr_zero(bit_ceil((unsigned)sz)));
  auto [lo, hi] = segtree_range_of_node(height, i);

添字の範囲 $\to$ ノード

メソッド nodes_for_range(lo, hi) で，半開区間 [lo, hi) をカバーするノードを vector で取得できる． 返り値中では，対応する半開区間が小さい順にノードが並んでいる．

  auto vec = st.nodes_for_range(lo, hi);

このメソッドと同等のことを行う関数 segtree_nodes_for_range(ht, lo, hi) も定義されている． ht には，セグメント木の高さを指定する．init_vec のサイズを sz とすると，次のようになる:

  unsigned height = countr_zero(bit_ceil((unsigned)sz)));
  auto vec = segtree_nodes_for_range(height, lo, hi);

コードへのリンク

GitHub へのリンク

keywords: segtree segment tree