差分

このページの2つのバージョン間の差分を表示します。

この比較画面へのリンク

両方とも前のリビジョン前のリビジョン
次のリビジョン
前のリビジョン
programming_algorithm:data_structure:sparse_table [2019/12/04] – [実装] ikatakosprogramming_algorithm:data_structure:sparse_table [2023/10/13] (現在) ikatakos
行 2: 行 2:
  
   * [[http://tookunn.hatenablog.com/entry/2016/07/13/211148|Sparse Tableを知ったので、忘れないように。 - tookunn’s diary]]   * [[http://tookunn.hatenablog.com/entry/2016/07/13/211148|Sparse Tableを知ったので、忘れないように。 - tookunn’s diary]]
 +
 +<note>2022/07/15 ブロック分割の項目を追加</note>
  
 ===== 概要 ===== ===== 概要 =====
行 13: 行 15:
   * 更新はできない   * 更新はできない
  
-足し算は冪等性を満たさないため、使えない。 +足し算は冪等性を満たさないため、使えない。\\ 
-(別に使ってもいいのだ、クエリ処理計算量 O(logN)O(logN) かかようになるし、 +データ持たせ方を工夫した DisjointSparseTable なら可能。(ただ、そもそも足し算なら累積和で区間クエリには答えられ 
-それならセグメント木と同じ上こちら更新も出来るのであまり優位性がい)+ 
 +  * [[programming_algorithm:data_structure:sparse_table:disjoint_sparse_table]] 
 + 
 +min, maxの他には、最小公倍数/最大公約数、bitwise-and/or ども冪等性を満たす。
  
-min, maxの他には、最小公倍数/最大公約数なども冪等性を満たす。 
  
 ===== データ構造 ===== ===== データ構造 =====
行 34: 行 38:
         2  |----------|                   3         2  |----------|                   3
         3  |----------------------|       1         3  |----------------------|       1
 +  i=1
 +      j=0     |-|                         6
 +        1     |----|                      3
 +        2     |----------|                3
 +        3     |----------------------|    1
 +  ...
  
 これを計算しておくと、例えば [2,8)[2,8) の最小値は、以下のように求められる。 これを計算しておくと、例えば [2,8)[2,8) の最小値は、以下のように求められる。
行 49: 行 59:
 ===== 実装 ===== ===== 実装 =====
  
-PythonならNumPyを使えば、tableの構築において k1k1 の計算結果を利用して kk を求める部分を高速化できる。 +以下の条件を満たす場合、PythonならNumPyを使えば、tableの構築において k1k1 の計算結果を利用して kk を求める部分を高速化できる。 
-(ただしminimum, maximumなど要素同士の演算を行う関数がNumPyに用意されているものに限る+ 
 +  * minimum, maximum, lcm, gcdなど2つの配列の要素同士の演算を行う関数がNumPyに用意されている 
 +  * 最小値の値そのものがわかればよく、最小値のindexは復元できなくてもよい 
 + 
 +(※未バリデーション
  
 <sxh python> <sxh python>
行 82: 行 96:
  
 </sxh> </sxh>
 +
 +
 +===== バリエーション =====
 +
 +==== ブロック分割 ====
 +
 +A0,...,AN1A0,...,AN1 をあらかじめ BB 個ごとのブロックに分け、ブロックに対しSparse Tableを構築する手法。
 +
 +=== 活用シーン ===
 +
 +主に以下の場合などに使われうる。
 +
 +  * 「[l,r)[l,r) の結果に ArAr を追加する」(区間を1伸ばす)計算量が、「[l1,r1)[l1,r1)[l2,r2)[l2,r2) の結果同士をマージする」計算量より小さい
 +  * 数列長 NN に比べ、クエリ数 QQ が小さい
 +  * 通常のSparse Tableで O(NlogN)O(NlogN) 個の情報を保存しておくにはメモリが足りない
 +
 +特に1つめが成り立つときに効果が大きいが、単純なSparse Tableに比べ処理が煩雑になり、
 +その分、定数倍もかさむので、劇的に改善することは期待しない方がよいかも。使いどころが難しい。
 +
 +=== 概要 ===
 +
 +「クエリ区間に完全に包括されるブロック区間」についてSparse Tableの前計算結果を利用し、はみ出た分は愚直に計算する。
 +
 +      0  1 ... B-1 | B ... 2B-1 | 2B ... | 3B ... | 4B ... | 5B ... | 6B ... |
 +   
 +  クエリ         |----------------------------------------------------------|
 +  
 +  愚直に計算     |---|
 +  table[1][2]        |---------------------------------------|
 +  table[2][2]                     |-----------------------------------|
 +  愚直に計算                                                          |-----|
 +  
 +  以上をマージしてクエリの答えとする
 +
 +=== 計算量 ===
 +
 +例えば求めたい演算処理が、「区間を1つだけ広げる計算量は O(x)O(x) だが、
 +結果のマージには O(y)O(y)x<yx<y)かかるような性質」を持っていると、
 +
 +  * 全てをSparse Table化する場合
 +    * 事前計算 O(NylogN)O(NylogN)
 +    * Q 回のクエリ O(Qy)
 +
 +  * ブロック分割する場合
 +    * 事前計算 O(Nx+NBylogNB)
 +    * Q 回のクエリ O(Q(Bx+y))
 +
 +クエリあたりの計算量を犠牲に、事前計算の計算量・メモリを減らすことができる。
 +
 +事前計算とクエリのバランスを取るブロック長 B はなかなか綺麗な形にはならないが、
 +いろいろ端折って計算すると B=NylogNQx を目安に設定すればよさそう。
 +
 +その場合、全体で O(Nx+NQxylogN+Qy) となり、
 +x<y,NQ の場合に、全てをSparse Table化するより高速になることがわかる。
 +
 +
 +==== ±1-RmQ ====
 +
 +上記のブロック分割の考え方にさらに制約を課して、省メモリ化・高速化する手法。
 +
 +=== 概要 ===
 +
 +A={A0,A1,...,AN1} の要素が、
 +必ず「左の値から1増えるか1減るかのどちらか」で変化していく配列の場合、
 +その区間最小値や最大値は、クエリ応答は O(1) のまま、
 +事前計算量と空間を O(NlogN) から O(N) に落とすことができる。\\
 +(結果のマージは O(1) とする)
 +
 +なかなか特殊な条件に見えるが、木の最小共通祖先だったり、たまに使える。
 +
 +(とはいえ競プロの文脈では、logのあるなしだけでTLEになるケースは稀で、実装の手間からこちらを選ぶのはよほどあとちょっとで通らない場合に限られる)
 +
 +  * [[https://www.slideshare.net/yumainoue965/lca-and-rmq|LCA and RMQ ~簡潔もあるよ!~]]
 +  * [[http://joisino.hatenablog.com/entry/2017/08/13/210000|前処理O(n)クエリO(1)のLCAと静的RMQ - ジョイジョイジョイ]]
 +
 +ブロック分割をベースの考え方として、
 +「はみ出た部分」の計算をクエリ毎に愚直に行うのではなく、
 +変化パターンが限られるのでパターン毎に事前計算しておくという発想。
 +
 +  B=3
 +  
 +      | 1  2  3 | 4  5  6 | 7  8  9 | 10 11 12 | 13 14 15 | 16 17 18 | 19 20 21 | 22 23 24 | 25 26 27 | 28 29 30 |
 +        |                                  |          |          |          |          |          |          |
 +  Ai    | 5  4  3 | 4  5  6 | 5  6  7 |  6  7  8 |  9  8  7 |  6  5  4 |  3  2  3 |  2  1  2 |  3  2  3 |  4  5  6 |
 +        |                                  |          |          |          |          |          |          |
 +  Min         3 | 4       | 5        6              7 |        4 |        |        |        |  4       |
 +  
 +  Query                            |----------------------------------------------------------------|
 +  ここはSparseTableで計算             |------------------------------------------------------|
 +  ここは個別に計算                 |-|                                                        |-----|
 +
 +== 個別に計算 ==
 +
 +ブロック内の変化は、今回1ブロック3要素なので、%%「++」「+-」「-+」「--」%%の4通りしかない。
 +
 +なので、全変化パターン・ブロック内の相対位置による区間 [l,r] ごとに、
 +最小値の左端要素からの差分配列を作っても、現実的な容量で間に合う。
 +
 +ブロックサイズを B とすると B22B1 の計算量とメモリがかかる。
 +
 +  先頭要素が0だったとして、変化パターン毎に最小値を計算
 +  
 +  [++]           [+-]           [-+]           [--]
 +  l\r 0  1  2    l\r 0  1  2    l\r 0  1  2    l\r 0  1  2
 +    0  0  0      0  0  0      0 -1 -1      0 -1 -2
 +        1          0        -1 -1        -1 -2
 +          2            0            0           -2
 +
 +例えば上の例で i=7,8,9 のブロックは「%%++%%」パターンであり、はみ出ている区間は [9,9]
 +ブロック内の相対的な位置になおすと [2,2] なので、
 +「%%++%% テーブル」から l=2,r=2 を参照すると 2 である。
 +ブロックの先頭要素 52 を足して、7 が区間の最小値とわかる。
 +
 +また、i=25,26,27 のブロックは「%%-+%%」パターンであり、区間 [0,1] のテーブル参照は 1 である。
 +ブロックの先頭要素 31 を足して、2 が区間の最小値である。
 +
 +=== 計算量 ===
 +
 +ブロックサイズを大きくしすぎると個別計算の計算量・メモリが増え、小さくしすぎるとSparseTableの計算量・メモリが増える。
 +
 +log2N2 が良いとされている。
 +
 +こうするとSparseTable, 個別計算のいずれも事前計算とメモリが O(N) に収まる。
 +
 +SparseTableの方は、ブロックの個数が 2Nlog2N となる。
 +この要素数でSparseTableを構築すると、2NlogNlog2NlogN=O(N) となる。
 +
 +(ちゃんとは理解してないけど、式変形すると 2N(logN2NlogNlogN) となり、対数の部分は極限取るとすぐ1になるので無視できる)
 +
 +また、個別計算の方も、logN222logN21=O(log2NN) となり、O(N) より小さくなる。
 +
 +=== 拡張性 ===
 +
 +必ずしも±1で遷移する場合に限定しなくても、ブロック内の変化パターン数が 2B 程度で収まる場合、同じ理由で活用できそう。具体的にあるかは知らんけど。
 +
  
  
programming_algorithm/data_structure/sparse_table.1575453699.txt.gz · 最終更新: by ikatakos
CC Attribution 4.0 International
Driven by DokuWiki Recent changes RSS feed Valid CSS Valid XHTML 1.0