ac-library-rb 0.5.0

data/document_ja/priority_queue.md

@@ -0,0 +1,103 @@
+# PriorityQueue
+
+優先度付きキューです。
+デフォルトでは降順で要素を保持します。
+
+**エイリアス** `HeapQueue`
+
+## 特異メソッド
+
+### new(array = []) -> PriorityQueue
+
+```ruby
+PriorityQueue.new
+pq = PriorityQueue.new([1, -1, 100])
+pq.pop # => 100
+pq.pop # => 1
+pq.pop # => -1
+```
+
+引数に与えた配列から優先度付きキューを構築します。
+
+**計算量** `O(n)` (`n` は配列の要素数)
+
+### new(array = []) {|x, y| ... } -> PriorityQueue
+
+```ruby
+PriorityQueue.new([1, -1, 100]) {|x, y| x < y }
+pq.pop # => -1
+pq.pop # => 1
+pq.pop # => 100
+```
+
+引数に与えた配列から優先度付きキューを構築します。
+ブロックで比較関数を渡すことができます。比較関数が渡されると、それを使って要素の優先度が計算されます。
+
+比較関数で比較できるなら、任意のオブジェクトをキューの要素にできます。
+
+**計算量** `O(n)` (`n` は配列の要素数)
+
+## インスタンスメソッド
+
+### push(item)
+
+```ruby
+pq.push(10)
+```
+
+要素を追加します。
+
+**エイリアス** `<<`, `append`
+
+**計算量** `O(log n)` (`n` はキューの要素数)
+
+### pop -> 最も優先度の高い要素 | nil
+
+```ruby
+pq.pop # => 100
+```
+
+最も優先度の高い要素をキューから取り除き、それを返します。キューが空のときは `nil` を返します。
+
+**計算量** `O(log n)` (`n` はキューの要素数)
+
+### get -> 最も優先度の高い要素 | nil
+
+```ruby
+pq.get # => 10
+```
+
+最も優先度の高い要素を**取り除くことなく**、その値を返します。
+
+キューが空のときは `nil` を返します。
+
+**エイリアス** `top`
+
+**計算量** `O(1)`
+
+### empty? -> bool
+
+```ruby
+pq.empty? # => false
+```
+
+キューの中身が空なら `true`、そうでないなら `false` を返します。
+
+**計算量** `O(1)`
+
+### heap -> Array(キューの要素)
+
+```ruby
+pq.heap
+```
+
+キューが内部で持っている二分ヒープを返します。
+
+## Verified
+
+- [Aizu Online Judge ALDS1_9_C Priority Queue](https://onlinejudge.u-aizu.ac.jp/problems/ALDS1_9_C)
+  - [code](https://onlinejudge.u-aizu.ac.jp/solutions/problem/ALDS1_9_C/review/4835681/qsako6/Ruby)
+
+## 参考
+
+- [cpython/heapq.py at master - python/cpython](https://github.com/python/cpython/blob/master/Lib/heapq.py)

data/document_ja/scc.md

@@ -0,0 +1,65 @@
+# SCC
+
+Strongly Connected Components
+
+有向グラフを強連結成分分解します。
+
+## 特異メソッド
+
+### new(n = 0) -> SCC
+
+```ruby
+graph = SCC.new(6)
+```
+
+`n` 頂点 `0` 辺の有向グラフを作ります。
+
+頂点番号は、0-based indexです。
+
+**制約** `0 ≦ n`
+
+**計算量** `O(n)`
+
+## インスタンスメソッド
+
+### add_edge(from, to)
+
+```ruby
+graph.add_edge(1, 4)
+```
+
+頂点 `from` から頂点 `to` へ有向辺を足します。
+
+**制約** `0 ≦ from < n`, `0 ≦ to < n`
+
+**計算量** `ならしO(1)`
+
+### scc -> Array[Array[Integer]]
+
+```ruby
+graph.scc
+```
+
+強連結成分分解して、頂点のグループの配列を返します。
+
+強連結成分分解は、お互いが行き来できる頂点を同じグループとなるように分解します。
+
+また、グループの順番は、トポロジカルソートとなっており、頂点uから頂点vに一方的に到達できるとき、頂点uに属するグループは頂点vに属するグループよりも先頭に来ます。
+
+**計算量** `O(n + m)` ※ `m` は辺数です。
+
+## Verified
+
+[ALPC: G \- SCC](https://atcoder.jp/contests/practice2/tasks/practice2_g)
+- [2538ms 2020/9/8](https://atcoder.jp/contests/practice2/submissions/16569175)
+
+# 参考リンク
+
+- 当ライブラリ
+  - [当ライブラリの実装コード scc\.rb](https://github.com/universato/ac-library-rb/blob/master/lib/scc.rb)
+  - [当ライブラリの実装コード scc_test\.rb](https://github.com/universato/ac-library-rb/blob/master/test/scc_test.rb)
+- 本家ライブラリ
+  - [本家ライブラリの実装コード scc\.hpp](https://github.com/atcoder/ac-library/blob/master/atcoder/scc.hpp)
+  - [本家ライブラリのドキュメント scc\.md](https://github.com/atcoder/ac-library/blob/master/document_ja/scc.md)
+- その他
+  - [強連結成分分解の意味とアルゴリズム \| 高校数学の美しい物語](https://mathtrain.jp/kyorenketsu)

data/document_ja/segtree.md

@@ -0,0 +1,145 @@
+# Segtree
+
+セグメント木です。
+
+## 特異メソッド
+
+### new(arg, e, &op) -> Segtree
+
+```rb
+seg = Segtree.new(arg, e) { |x, y| ... }
+```
+
+第1引数は、`Integer`または`Array`です。
+
+- 第1引数が`Integer`の`n`のとき、長さ`n`・初期値`e`のセグメント木を作ります。
+- 第1引数が長さ`n`の`Array`の`a`のとき、`a`をもとにセグメント木を作ります。
+
+第2引数は単位元`e`で、ブロックで二項演算`op(x, y)`を定義することで、モノイドを定義する必要があります。
+
+**計算量** `O(n)`
+
+### モノイドの設定コード例
+
+```ruby
+n   = 10**5
+inf = (1 << 60) - 1
+
+Segtree.new(n, 0) { |x, y| x.gcd y }       # gcd
+Segtree.new(n, 1) { |x, y| x.lcm y }       # lcm
+Segtree.new(n, -inf) { |x, y| [x, y].max } # max
+Segtree.new(n,  inf) { |x, y| [x, y].min } # min
+Segtree.new(n, 0) { |x, y| x | y }         # or
+Segtree.new(n, 1) { |x, y| x * y }         # prod
+Segtree.new(n, 0) { |x, y| x + y }         # sum
+```
+
+## インスタンスメソッド
+
+### set(pos, x)
+
+```rb
+seg.set(pos, x)
+```
+
+`a[pos]` に `x` を代入します。
+
+**計算量**
+
+- `O(logn)`
+
+### get(pos)
+
+```rb
+seg.get(pos)
+```
+
+`a[pos]` を返します。
+
+**計算量**
+
+- `O(1)`
+
+### prod(l, r)
+
+```rb
+seg.prod(l, r)
+```
+
+`op(a[l], ..., a[r - 1])` を返します。
+
+**制約**
+
+- `0 ≦ l ≦ r ≦ n`
+
+**計算量**
+
+- `O(logn)`
+
+### all_prod
+
+```rb
+seg.all_prod
+```
+
+`op(a[0], ..., a[n - 1])` を返します。
+
+**計算量**
+
+- `O(1)`
+
+### max_right(l, &f) -> Integer
+
+```ruby
+seg.max_right(l, &f)
+```
+
+Segtree上で二分探索をします。
+
+### min_left(r, &f) -> Integer
+
+```ruby
+seg.min_left(r, &f)
+```
+
+Segtree上で二分探索をします。
+
+## Verified
+
+- [ALPC: J \- Segment Tree](https://atcoder.jp/contests/practice2/tasks/practice2_j)
+
+- [F \- Range Xor Query](https://atcoder.jp/contests/abc185/tasks/abc185_f)
+  - xorのセグメントツリーの基本的な典型問題です。FenwickTree(BIT)をxorに改造するだけでも解けます。
+  - [ACコード(1538ms)](https://atcoder.jp/contests/abc185/submissions/18746817): 通常のSegtree解。
+  - [ACコード(821ms)](https://atcoder.jp/contests/abc185/submissions/18769200): FenwickTree(BIT)のxor改造版。
+
+## 参考リンク
+
+- 当ライブラリ
+  - [当ライブラリの実装コード segtree.rb](https://github.com/universato/ac-library-rb/blob/master/lib/segtree.rb)
+  - [当ライブラリのテストコード segtree.rb](https://github.com/universato/ac-library-rb/blob/master/test/segtree_test.rb)
+    - テストコードも具体的な使い方として役に立つかもしれまん。
+- 本家
+  - [本家ライブラリのドキュメント segtree.md(GitHub)](https://github.com/atcoder/ac-library/blob/master/document_ja/segtree.md)
+  - [本家のドキュメント appendix.md(GitHub)](https://github.com/atcoder/ac-library/blob/master/document_ja/appendix.md)
+  - [本家ライブラリの実装コード segtree.hpp(GitHub)](https://github.com/atcoder/ac-library/blob/master/atcoder/segtree.hpp)
+  - [本家ライブラリのテストコード segtree_test.cpp(GitHub)](https://github.com/atcoder/ac-library/blob/master/test/unittest/segtree_test.cpp)
+  - [本家ライブラリのサンプルコード segtree_practice.cpp(GitHub)](https://github.com/atcoder/ac-library/blob/master/test/example/segtree_practice.cpp)
+- セグメントツリーについて
+  - [セグメント木をソラで書きたいあなたに \- hogecoder](https://tsutaj.hatenablog.com/entry/2017/03/29/204841)
+
+## 本家ライブラリとの違い等
+
+基本的に、本家の実装に合わせています。
+
+内部実装に関しても、変数`@d`の0番目の要素には必ず単位元`@e`が入ります。
+
+### 変数名の違い
+
+Rubyには`p`メソッドがあるので、引数`p`について、`p`ではなくpositionの`pos`を変数名として使いました。
+
+同様に、本家の変数`size`を、わかりやすさから`leaf_size`としています。
+
+### `ceil_pow2`ではなく、`bit_length`
+
+本家C++ライブラリは独自定義の`internal::ceil_pow2`関数を用いてますが、本ライブラリではRuby既存のメソッド`Integer#bit_length`を用いています。そちらの方が計測した結果、高速でした。

data/document_ja/string.md

@@ -0,0 +1,105 @@
+# String
+
+文字列`s`の`a`番目から`b-1`番目の要素の`substring`を、`s[a..b)`と表記します。
+
+## suffix_array
+
+```ruby
+(1) suffix_array(s)
+(2) suffix_array(a)
+(3) suffix_array(a, upper)
+```
+
+長さ`n`の文字列`s`に対し、Suffix Arrayとして長さ`n`の`Integer`の配列を返します。
+
+Suffix Array `sa`は`0 ... n`の順列であって、各`i`について`s[sa[i]..n) < s[sa[i+1]..n)`を満たします。
+
+このSuffix Arrayの概念は一般の列にも適用でき、`<=>`演算子で比較可能な要素の配列`a`に対しても動作します。
+
+1. `suffix_array(s)`
+
+内部で`s.bytes`によってバイト列に変換します。
+
+**制約**
+
+- `s`は文字コード255以下の文字のみからなる文字列
+
+**計算量**
+
+`O(|s|)`
+
+2. `suffix_array(a)`
+
+内部で、いわゆる座標圧縮を行います。要素は大小関係を保ったまま非負整数に変換されます。
+
+**制約**
+
+- `a`の要素は互いに`<=>`演算子で比較可能
+
+**計算量**
+
+要素の比較が定数時間で行えるとして、`O(|a| log |a|)`
+
+3. `suffix_array(a, upper)`
+
+**制約**
+
+- `a`は`Integer`の配列
+
+- `a`の全ての要素`x`について`0≦x≦upper`
+
+**計算量**
+
+`O(|a| + upper)`
+
+## lcp_array
+
+```ruby
+lcp_array(s, sa)
+```
+
+長さ`n`の文字列`s`のLCP Arrayとして、長さ`n-1`の`Integer`の配列を返します。`i`番目の要素は`s[sa[i]..n)`と`s[sa[i+1]..n)`のLCP ( Longest Common Prefix ) の長さです。
+
+こちらも`suffix_array`と同様一般の列に対して動作します。
+
+**制約**
+
+- `sa`は`s`のSuffix Array
+
+- (`s`が文字列の場合) `s`は文字コード255以下の文字のみからなる
+
+**計算量**
+
+`O(|s|)`
+
+## z_algorithm
+
+```ruby
+z_algorithm(s)
+```
+
+入力された配列の長さを`n`として、長さ`n`の`Integer`の配列を返します。`i`番目の要素は`s[0..n)`と`s[i..n)`のLCP ( Longest Common Prefix ) の長さです。
+
+**制約**
+
+- `s`の要素は互いに`==`演算子で比較可能
+
+**計算量**
+
+`O(|s|)`
+
+<hr>
+
+## Verified
+
+- suffix_array, lcp_array
+  - [I - Number of Substrings](https://atcoder.jp/contests/practice2/tasks/practice2_i)
+    - [AC 1362 ms](https://atcoder.jp/contests/practice2/submissions/17194669)
+
+- z_algorithm
+  - [ABC135 F - Strings of Eternity](https://atcoder.jp/contests/abc135/tasks/abc135_f)
+    - [AC 1076 ms](https://atcoder.jp/contests/abc135/submissions/16656965)
+
+## 参考
+
+- [本家ACLのドキュメント String](https://atcoder.github.io/ac-library/master/document_ja/string.html)

data/document_ja/two_sat.md

@@ -0,0 +1,87 @@
+# 2-SAT
+
+2-SAT を解きます。変数 x_0, x_1, ..., x_N-1 に関して、
+
+(x_i = f) ∨ (x_j = g)
+
+というクローズを足し、これをすべて満たす変数の割当があるかを解きます。
+
+**エイリアス**
+
+- `TwoSAT`
+- `TwoSat`
+
+<hr>
+
+## 特異メソッド
+
+### new(n = 0) -> TwoSAT
+
+```ruby
+ts = TwoSAT.new(n)
+```
+
+n 変数の 2-SAT を作ります。
+
+**計算量**
+
+- `O(n)`
+
+<hr>
+
+## インスタンスメソッド
+
+### add_clause(i, f, j, g) -> nil
+
+`i: Integer`, `f: bool`, `j: Integer`, `g: bool`
+
+```ruby
+ts.add_clause(i, true, j, false)
+```
+
+(x_i = f) ∨ (x_j = g)　というクローズを足します。
+
+**計算量**
+
+- ならし `O(1)`
+
+### satisfiable? -> bool
+
+```ruby
+ts.satisfiable?
+```
+
+条件を満たす割当が存在するかどうかを判定します。割当が存在するならば `true`、そうでないなら `false` を返します。
+
+**エイリアス**
+
+- `satisfiable?`
+- `satisfiable`
+
+**計算量**
+
+- 足した制約の個数を `m` として `O(n + m)`
+
+### answer -> Array[bool]
+
+```ruby
+ts.answer
+```
+
+最後に呼んだ `satisfiable?` のクローズを満たす割当を返します。
+`satisfiable?` を呼ぶ前や、`satisfiable?` で割当が存在しなかったときにこの関数を呼ぶと、長さ n の意味のない配列を返します。
+
+**計算量**
+
+- `O(n)`
+
+<hr>
+
+## Verified
+
+- [H - Two SAT](https://atcoder.jp/contests/practice2/tasks/practice2_h)
+  - [163 ms](https://atcoder.jp/contests/practice2/submissions/16655036)
+
+## 参考
+
+[本家 ACL のドキュメント 2-SAT](https://atcoder.github.io/ac-library/master/document_ja/twosat.html)