ac-library-rb 0.5.0

data/document_en/string.md

@@ -0,0 +1,106 @@
+# String
+
+The substring of the `a`-th to `b-1`-th element of the string `s` is written as `s[a..b)`.
+
+
+## suffix_array
+
+```ruby
+(1) suffix_array(s)
+(2) suffix_array(a)
+(3) suffix_array(a, upper)
+```
+For a string `s` of length `n`, returns an array of `Integer` of length `n` as a Suffix Array.
+
+Suffix Array `sa` is a permutation of `0 ... ... n`, and for each `i`, it is a permutation of `s[sa[i]. .n) < s[sa[i+1]. .n)` is satisfied.
+
+This concept of Suffix Array can also be applied to general columns and works for `a`, an array of elements comparable with `<=>` operator. 1.
+
+### 1. `suffix_array(s)`
+
+Convert internally to a sequence of bytes by `s.bytes`.
+
+**Constraints**
+
+- `s` is a string consisting only of characters with character code 255 or less.
+
+**Complexity**
+
+- `O(|s|)`
+
+### 2. `suffix_array(a)`
+
+Inside, so-called coordinate compression is performed. Elements are converted to non-negative integers keeping their size relation.
+
+**Constraints**
+
+- The elements of `a` are comparable to each other with the `<=>` operator.
+
+**Complexity**
+
+- Assuming elements can be compared in constant time, `O(|a| log |a|)`
+
+### 3. `suffix_array(a, upper)`
+
+**Constraints**
+
+- `a` is an array of `Integer`.
+
+- For all elements `x` of `a`, `0 ≤ x ≤ upper`.
+
+**Complexity**
+
+- `O(|a| + upper)`
+
+## lcp_array
+
+```ruby
+lcp_array(s, sa)
+```
+
+It returns an array of `Integer` of length `n-1` as an LCP Array of length `n` for strings `s`. The `i` th element is the length of LCP ( Longest Common Prefix ) for `s[sa[i]. .n)` and `s[sa[i+1]. .n)`.
+
+This also works for general columns as well as `suffix_array`.
+
+**Constraints**
+
+- `sa` is a Suffix Array of `s`.
+
+- If `s` is a string, `s` consists only of characters with character code 255 or less.
+
+**Complexity**
+
+- `O(|s|)`
+
+## z_algorithm
+
+```ruby
+z_algorithm(s)
+```
+
+It returns an array of `Integer` of length `n`, where `n` is the length of the input array. The `i` th element is the length of LCP ( Longest Common Prefix ) of `s[0..n)` and `s[i..n)`.
+
+**Constraints**
+
+- The elements of `s` can be compared with each other with the `==` operator
+
+
+**Complexity**
+
+- `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)
+
+## Links
+
+- AtCoder Library
+  - [Document String(HTML)](https://atcoder.github.io/ac-library/master/document_ja/string.html)

data/document_en/two_sat.md

@@ -0,0 +1,91 @@
+# 2-SAT
+
+It solves 2-SAT.
+
+For variables x_0, x_1, ..., x_N-1 and clauses with form
+
+(x_i = f) ∨ (x_j = g)
+
+it decides whether there is a truth assignment that satisfies all clauses.
+
+
+**Alias**
+
+`TwoSAT`, `TwoSat`
+
+<hr>
+
+## Class Methods
+
+### new(n = 0) -> TwoSAT
+
+```ruby
+ts = TwoSAT.new(n)
+```
+
+It creates a 2-SAT of `n` variables and `0` clauses.
+
+**Complexity**
+
+O(n)
+
+<hr>
+
+## Instance Methods
+
+### add_clause(i, f, j, g) -> nil
+
+`i: Integer`, `f: bool`, `j: Integer`, `g: bool`
+
+```ruby
+ts.add_clause(i, true, j, false)
+```
+
+It adds a clause (x_i = f) ∨ (x_j = g).
+
+**Constraints**
+
+- `0 <= i < n`
+- `0 <= j < n`
+
+**Complexity**
+
+- `O(1)` amortized
+
+### satisfiable? -> bool
+
+```ruby
+ts.satisfiable?
+```
+
+If there is a truth assignment that satisfies all clauses, it returns `true`. Otherwise, it returns `false`.
+
+**Alias**
+
+`satisfiable?`, `satisfiable`
+
+**Complexity**
+
+`O(n + m)`, where `m` is the number of added clauses.
+
+### answer -> Array(bool)
+
+```ruby
+ts.answer
+```
+
+It returns a truth assignment that `satisfie` all clauses of the last call of satisfiable. If we call it before calling `satisfiable` or when the last call of `satisfiable` returns `false`, it returns the vector of length `n` with undefined elements.
+
+**Complexity**
+
+- `O(n)`
+
+## Verified
+
+- [H - Two SAT](https://atcoder.jp/contests/practice2/tasks/practice2_h)
+  - [163 ms](https://atcoder.jp/contests/practice2/submissions/16655036 )
+
+## Links
+
+- AtCoder Library(ACL)
+  - [Documetn twosat.html](https://atcoder.github.io/ac-library/master/document_ja/twosat.html)

data/document_en/union_find.md

@@ -0,0 +1,3 @@
+See [DSU document (dsu.md)](https://github.com/universato/ac-library-rb/blob/main/document_en/dsu.md).
+
+DSU is an alias of UnionFind.

data/document_ja/convolution.md

@@ -0,0 +1,64 @@
+# Convolution
+
+畳み込みを行います。長さ`N`の`Integer`の配列`a[0],a[1],...,a[N-1]`と長さ`M`の`Integer`の配列`b[0],b[1],...,b[M-1]`から畳み込みを計算し、長さ`N+M-1`の`Integer`の配列`c`として返します。
+
+## convolution
+
+```ruby
+(1) conv = Convolution.new(m = 998244353)
+(2) conv = Convolution.new(m, primitive_root)
+```
+
+畳み込みを`mod m`で計算するためのオブジェクトを作成します。
+
+`primitive_root`は法`m`に対する原始根である必要があります。与えられなかった場合は、内部で最小の原始根を計算します。
+
+**制約**
+
+- `2≦m`
+
+- `m`は素数
+
+- ((2)のみ) `primitive_root`は法`m`に対する原始根
+
+**計算量**
+
+1. `O(sqrt m)`
+
+2. `O(log m)`
+
+### 使用方法
+
+実際の計算は、上で作成したオブジェクト`conv`を用いて次のように行います。
+
+```ruby
+c = conv.convolution(a, b)
+```
+
+`a`と`b`の少なくとも一方が空配列の場合は空配列を返します。
+
+**制約**
+
+- `2^c|(m-1)`かつ`|a|+|b|-1<=2^c`なる`c`が存在する
+
+**計算量**
+
+- `O((|a|+|b|) log (|a|+|b|))`
+
+## convolution_ll
+
+`convolution`の`m`として十分大きな素数を用いることで対応できます。
+
+`1e15`を超える`NTT-Friendly`な素数の1例として、`1125900443713537 = 2^29×2097153+1`があります。
+
+# Verified
+
+- [C - 高速フーリエ変換](https://atcoder.jp/contests/atc001/tasks/fft_c)
+  - `m = 1012924417`
+    [1272ms](https://atcoder.jp/contests/atc001/submissions/17193829)
+  - `m = 1125900443713537`
+    [2448 ms](https://atcoder.jp/contests/atc001/submissions/17193739)
+
+# 参考
+
+- [本家 ACL のドキュメント Convolution](https://atcoder.github.io/ac-library/master/document_ja/convolution.html)

data/document_ja/dsu.md

@@ -0,0 +1,183 @@
+# DSU - Disjoint Set Union
+
+別名: Union Find (むしろ、こちらが有名)
+
+無向グラフに対して、
+
+- 辺の追加(2頂点の連結)
+- 2頂点が連結かの判定(別の頂点を通して行き来できることを含みます)
+
+をならし`O(α(n))`時間で高速に計算できます。
+
+
+また、内部的に、各連結成分ごとに代表元(代表となる頂点)を1つ持っています。
+
+辺の追加により連結成分がマージされるとき、元の連結成分の代表元のどちらかが新たな代表元になります。
+
+## 使い方
+
+```rb
+d = DSU.new(5)
+p d.groups      # => [[0], [1], [2], [3], [4]]
+p d.same?(2, 3) # => false
+p d.size(2)     # => 1
+
+d.merge(2, 3)
+p d.groups      # => [[0], [1], [2, 3], [4]]
+p d.same?(2, 3) # => true
+p d.size(2)     # => 2
+```
+
+## 特異メソッド
+
+### new(n = 0) -> DSU
+
+```rb
+d = DSU.new(n)
+```
+
+`n` 頂点, `0` 辺の無向グラフを生成します。
+
+頂点の番号は `0` 始まりで数え、`n - 1` までになります。
+
+**計算量**
+
+- `O(n)`
+
+
+**エイリアス**
+
+- `DSU`
+- `DisjointSetUnion`
+- `UnionFind`
+- `UnionFindTree`
+
+## インスタンスメソッド
+
+### merge(a, b) -> Integer
+
+```rb
+d.merge(a, b)
+```
+
+頂点 `a` と頂点 `b`を連結させます。
+`a`, `b` が既に連結だった場合はその代表元、非連結だった場合は連結させたあとの新たな代表元を返します。
+
+計算量 ならしO(α(n))
+
+**エイリアス**
+
+- `merge`
+- `unite`
+
+### same?(a, b) -> bool
+
+```rb
+d.same?(a, b)
+```
+
+頂点aと頂点bが連結であるとき `true` を返し、そうでないときは `false` を返します。
+
+**計算量**
+
+- ならし `O(α(n))`
+
+**エイリアス**
+
+- `same?`
+- `same`
+
+### leader(a) -> Integer
+
+```rb
+d.leader(a)
+```
+
+頂点 `a` の属する連結成分の代表(根)の頂点を返します。
+
+**計算量**
+
+- ならし `O(α(n))`
+
+**エイリアス**
+
+- `leader`
+- `root`
+- `find`
+
+### size(a) -> Integer
+
+```rb
+d.size(3)
+```
+
+頂点 `a` の属する連結成分の頂点数(サイズ)を返します。
+
+**計算量**
+
+- ならし `O(α(n))`
+
+### groups -> [[Integer]]
+
+```rb
+d.groups
+```
+
+グラフを連結成分に分けて、その情報を返します。
+
+返り値は2次元配列で、内側・外側ともに配列内の順番は未定義です。
+
+**計算量**
+
+- `O(n)`
+
+## Verified
+
+- [A \- Disjoint Set Union](https://atcoder.jp/contests/practice2/tasks/practice2_a)
+
+## 参考リンク
+
+- 当ライブラリ
+  - [当ライブラリの実装コード dsu.rb](https://github.com/universato/ac-library-rb/blob/master/lib/dsu.rb)
+  - [当ライブラリのテストコード dsu_test.rb](https://github.com/universato/ac-library-rb/blob/master/test/dsu_test.rb)
+- 本家ライブラリ
+  - [本家ACLのドキュメント dsu.md](https://github.com/atcoder/ac-library/blob/master/document_ja/dsu.md)
+  - [本家ACLのコード dsu.hpp](https://github.com/atcoder/ac-library/blob/master/atcoder/dsu.hpp)
+
+## その他
+
+### DSUよりUnionFindの名称の方が一般的では?
+
+UnionFindの方が一般的だと思いますが、本家ライブラリに合わせています。なお、`UnionFind`をクラスのエイリアスとしています。
+
+https://twitter.com/3SAT_IS_IN_P/status/1310122929284210689 (2020/9/27)
+
+Google Scholar によると
+- "union find": 8,550件
+- "union find data structure": 1,970件
+- "disjoint set union": 1,610 件
+- "disjoint set data structure": 1,030 件
+
+### メソッドのエイリアスについて
+
+本家ライブラリでは`same`ですが、真偽値を返すメソッドなのでRubyらしく`same?`をエイリアスとして持ちます。`same`はdeprecated(非推奨)ですので、`same?`を使って下さい。
+
+本家は`merge`ですが、`unite`もエイリアスに持ちます。`unite`は`UnionFind`に合っているし、蟻本などでも一般的にもよく使われています。
+
+本家は`leader`ですが、`UnionFind`に合っていて一般的であろう`find`や`root`もエイリアスとして持ちます。
+
+### mergeの返り値について
+
+本家ライブラリの実装に合わせて、`merge`メソッドは新たにマージしたか否かにかかわらず、代表元を返しています。
+
+ただ、新たに結合した場合には代表元を返して、そうでない場合は`nil`か`false`を返すような`merge?`(あるいは`merge!`)を入れようという案があります。Rubyに, true/false以外を返す` nonzero?`などもあるので、`merge?`という名称は良いと思います。
+
+### 実装の説明
+
+本家ライブラリに合わせて、内部のデータを`@parent_or_size`というインスタンス変数名で持っています。これは、根(代表元)となる要素の場合は、その連結成分の要素数(サイズ)に-1を乗じた値を返し、それ以外の要素の場合に、その要素の属する連結成分の代表(根)の番号を返します。
+
+なお、インスタンスが初期化され、まだどの頂点どうしも連結されていないときは、全ての頂点が自分自身を代表元とし、サイズ1の連結成分が頂点の数だけできます。そのため、初期化されたときは、内部の配列`@parent_or_size`は、要素が全て-1となります。
+
+### 変更履歴
+
+2020/10/25 [PR #64] メソッド`groups`(正確には内部で定義されている`groups_with_leader`)のバグを修正しました。