トーリックイデアルのグレブナ基底を求めるアルゴリズムの研究.pptVIP

トーリックイデアルの　　グレブナ基底を求める　　　　アルゴリズム – F4およびF5 – 東京大学　情報理工学系研究科 中山 裕貴 December 15, 2003 発表の背景 グレブナ基底の計算方法 従来はBuchbergerの算法およびその亜種（改良は主に、不要なペアを除く手法について） F4,F5とは Faugèreによる、グレブナ基底計算の新しい算法 F4[Faugère ‘99] ??? 多項式同士の計算を、まとめて行列の形で行う F5[Faugère ’02] ??? 新しい多項式を追加するとき、以前の計算結果を部分的に保存しておく この２つの算法により、今まで解けなかったサイズの問題が解けたと言われている 発表の内容 今回は、対象をトーリックイデアルの場合に限定、アルゴリズムを専用に特化させる F4アルゴリズムを改良し、高速化?メモリ節約 このとき、簡約操作を非常に効率よく行うことができる F5アルゴリズムを新たに実装 上記２つのアルゴリズムについて、２つの異なる項順序で実行時間を計測 実行時間の解析のため、計算途中におけるsugarや次数、およびペア数を比較する アルゴリズムはasir上で実装?評価 トーリックイデアルとは 多項式イデアル　　　　　　　　　　　を考える 　　　　ならば、 　　および多項式 　について、 トーリックイデアル　　は行列　　で定まり を満たす二項式イデアル (係数は1,-1) 単項間の順序関係 任意の単項m1,m2に対し、順序関係を定める 整列順序であり、最小元は 　　　　ならば、 多変数の場合、順序は一意ではない 辞書式順序 全次数逆辞書式順序 グレブナ基底とは イデアルI に含まれる多項式に対し、その先頭項で生成されるイデアル　　 を考える 一般に、イデアルI の生成元　　　　　　について しかし、生成元　　　　　がグレブナ基底のときまたそのときに限って グレブナ基底を求めるBuchbergerのアルゴリズム 入力: イデアルの生成元 出力: グレブナ基底 do { 　　　　 ; Gの任意のペア(i ,j )に対し、Spol(gi ,gj )を計算; 　Spol(gi ,gj )をGの要素で簡約したものをg’ とする; g’ が0でなければ、Gに加える; } while(　　　　) Buchbergerアルゴリズムの主な改良法 F4アルゴリズムの概要 複数のS多項式を計算しておく それらの簡約に必要な元の候補を選ぶ 簡約時に加減算だけで済むよう、単項式倍しておく Symbolic Preprocessing F4アルゴリズムの概要 S多項式?簡約に使う元を行列で表す 行が多項式、列が単項の種類を表す 簡約 = 行列の対角化 (多項式の定数倍?加減算) トーリックイデアルに対するF4アルゴリズム トーリックイデアル F5アルゴリズムの概要 新たに多項式を生成したとき、生成に使った多項式の情報を保持しておく Buchberger, 改良したF4, F5のベンチマーク – 実験の環境 入力: 行列　　　　　　　　　　のトーリックイデアルの生成元 項順序 全次数逆辞書式順序 辞書式順序 計算代数システムasir上での実装 ３つのアルゴリズムの比較 組み込みのBuchbergerアルゴリズム (modular計算はしない) 今回改良したF4アルゴリズム F5アルゴリズム　　 CPU : UltraSPARC-II 360MHz, メモリ : 2GB 　 ベンチマーク (n=10～40,全次数逆辞書式順序) ベンチマーク (n=10～40,辞書式順序) F4,F5の計算時間増大の原因 辞書式順序でF4,F5の計算時間が増える理由 生成される多項式の次数(sugar)が増大 A20の場合、基底の最大次数は全次数逆辞書式の場合　??? 3辞書式の場合　??? 20 F4の場合、簡約行列のサイズが冗長になる? F5の場合、不要なペア?Ruleの数が増大? A20を例に実験を行った結果???次ページ F4による A20のグレブナ基底計算の効率 全次数逆辞書式順序　??? F5 によるA20のグレブナ基底計算の効率 全次数逆辞書式順序 ???　 実験結果 組み込み関数とF4,F5の比較 全次数逆辞書式順序では、F4,F5の方が高速 辞書式順序では逆に時間がかかる F4の場合、わずかな多項式を得るために必要以上の計算を行っている F5の場合、次数が大きくなったときに、ペアの選択ストラテジーがうまく働かない F4について、簡約時に生成される行列をより小さくする F5について、sugarを使う方式に書き換える