[論文レビュー] Fast Gröbner Basis Computation for Boolean Polynomials

Franziska Hinkelmann, Elizabeth Arnold|arXiv (Cornell University)|Oct 13, 2010

Polynomial and algebraic computation参考文献 5被引用数 7

ひとこと要約

この論文では、多項式の単項式の2進表現を用いてビット単位の演算を実行することで、有理数上の記号計算を回避するC++アクセラレート型のMacaulay2パッケージであるBooleanGBを提示する。この手法により、生物学的モデルやスジドクなど、変数数が増加する応用分野において、標準的なMacaulay2実装と比較して最大50倍の高速化が達成され、計算時間を数分から5秒未塔に短縮する。

ABSTRACT

We introduce the Macaulay2 package BooleanGB, which computes a Gröbner basis for Boolean polynomials using a binary representation rather than symbolic. We compare the runtime of several Boolean models from systems in biology and give an application to Sudoku.

研究の動機と目的

• 係数と次数の増大に起因する、Boolean多項式の記号的Gröbner基底計算における性能ボトルネックを解消すること。
• 商環 $\mathbb{F}_2[x_1,\dots,x_n]/\langle x_i^2 + x_i \rangle$ 専用に最適化された、高速で効率的なBuchbergerのアルゴリズムの実装を開発すること。
• システム生物学やスジドクのような組合せパズルのような複雑な系におけるGröbner基底の実用的計算を可能にすること。
• Macaulay2のエージェントとしてネイティブなC++パッケージとして統合し、高性能な代数計算を実現すること。

提案手法

• Boolean多項式の単項式を2進文字列によって整数に表現し、高速なビット単位の演算を可能にする。
• 単項式の乗算と加算にビット単位のORを、可除性のチェックにビット単位のANDを用いる。
• 最小公倍数と除法をビット操作で実装する：LCMはOR、可除性がある場合は除法はXORで実行する。
• Buchbergerの第一および第二の基準をビット操作を用いて適用し、冗長なS多項式計算を削減する。
• 体多項式 $x_i^2 + x_i$ を暗黙的に符号化し、S多項式が多項式の多項式に残るという定理に依存する。
• アルゴリズムをMacaulay2のエンジンにC++ベースのパッケージとして統合し、`gbBoolean`関数を通じて公開する。

実験結果

リサーチクエスチョン

• RQ1 Macaulay2における記号的メソッドと比較して、Boolean多項式のビット表現がGröbner基底計算時間を著しく短縮できるか？
• RQ2 2進表現手法が、BooleanイデアルのGröbner基底計算における係数および次数の爆発を回避できるか？
• RQ3 実世界の生物学的モデルや組合せ問題において、BooleanGBパッケージのパフォーマンスは、標準的なMacaulay2のGröbner基底アルゴリズムと比較してどうなるか？
• RQ4 変数数が大幅に増加するスジドクのような大規模システムにも、この手法はスケーラブルか？
• RQ5 ビット単位の演算を用いることで、Boolean多項式イデアルのBuchbergerのアルゴリズムにおいて、正しさを保ちながら高いパフォーマンスを達成できるか？

主な発見

• BooleanGBパッケージは、標準的なlexicographic Gröbner基底計算と比較して、Macaulay2で最大50倍の高速化が達成され、実行時間比は0.017562079%にまで低下した。
• Shidoku問題において、BooleanGBは4.11442秒でGröbナ基底を計算したが、標準的なMacaulay2（lexicographic順序）では470.703秒を要した。
• GINsimリポジトリに収録されたすべてのテスト済み生物学的モデルにおいて、BooleanGBは現在のMacaulay2実装を上回り、ほとんどの場合で0.00002秒未塔の実行時間にまで短縮された。
• 単項式を64ビット整数として表現し、ビット単位の演算を実行することで、係数と次数の爆発を回避し、定数時間の算術演算を保証した。
• 変数数が16から64に増加したShidoku変換においても、lexicographic順序での完全なGröbner基底が正しく計算された。
• 多項式の多項式とフィールド多項式の間のS多項式が常に多項式のまま保たれるという理論的結果により、フィールド多項式の特別な符号化を省略でき、実装が簡素化された。