[論文レビュー] Homotopy methods for multiplication modulo triangular sets

Alin Bostan, Muhammad F. I. Chowdhury+2 more|ArXiv.org|Jan 23, 2009

Polynomial and algebraic computation参考文献 22被引用数 8

ひとこと要約

本稿では、変形パラメータに関する形式的冪級数上での評価・補間技術を活用して、多項式の三角的集合に関する高速乗算のためのホモトピーに基づくアルゴリズムを提示する。特定の三角的集合の族、特に再帰的システムやコーシー加群に対して、$ O(\tilde{u}) $ の準線形時間計算量を達成する。ここで $ u $ は入力／出力サイズを表す。これは従来の手法に比べ顕著な改善をもたらす。

ABSTRACT

We study the cost of multiplication modulo triangular families of polynomials. Following previous work by Li, Moreno Maza and Schost, we propose an algorithm that relies on homotopy and fast evaluation-interpolation techniques. We obtain a quasi-linear time complexity for substantial families of examples, for which no such result was known before. Applications are given to notably addition of algebraic numbers in small characteristic.

研究の動機と目的

• 多変数および再帰的設定において、多項式乗算を三角的集合で行う際の高い計算コストに対処すること。
• 標準的手法が単項式の数に指数的ブロードキャストを引き起こす状況において、三角的集合に関する乗算の高速アルゴリズムを開発すること。
• 一変数の場合に限らないホモトピーに基づく評価・補間技術を、スパースまたは低次数構造を有する多変数三角的集合へと拡張すること。
• 小特異性における代数的数体系およびコーシー加群における乗算において、準線形時間計算量を提供すること。
• 変形に基づく補間処理中に生じる冪級数計算における精度問題を解消すること。

提案手法

• 変形パラメータ $ \eta $ を導入し、$ \mathbf{U} $ が基本環に根を持つように、新たな三角的集合 $ \mathbf{S} = \eta\mathbf{T} + (1-\eta)\mathbf{U} $ を構築する。
• ヘンゼルの補題を用いて、$ \mathbf{S} $ が $ \mathsf{R}[[\eta]] $ に根を持つことを保証し、形式的冪級数上での評価と補間を可能にする。
• $ \mathsf{R}[[\eta]] $ 上での高速評価および補間アルゴリズムを適用し、変形を介して $ AB \mod \langle \mathbf{T} \rangle $ の積を計算する。
• 直接的な多項式乗算と還元を回避するため、ホモトピー技術を活用し、単項式数の指数的ブロードキャストを回避する。
• 三角的集合 $ \mathbf{T} $ の単項式サポートに依存する $ \eta $-進精度を制御することで、冪級数計算における精度を制御する。
• 小特異性において、最終結果を回復するためにニュートンの恒等式とリフト技術を用い、準線形ビット計算量を保証する。

実験結果

リサーチクエスチョン

• RQ1 単項式サポートがスパースまたは次数が低い場合に、多変数三角的集合に関する乗算を準線形時間で行うことは可能か？
• RQ2 ホモトピーおよび変形技術を一変数から多変数三角的集合へと拡張することで、高速な評価・補間を実現できるか？
• RQ3 ホモトピー法において、正しくかつ効率的に計算を行うために、冪級数計算における精度制御にどのような要件が生じるか？
• RQ4 提案手法が、小特異性における代数的数体系およびコーシー加群において準線形複雑度を達成できるか？
• RQ5 ホモトピーに基づくアプローチにおける計算ボトルネックは何か。それらはどのように緩和できるか？

主な発見

• 形式 $ X_i^2 - 2X_{i-1} $ または $ X_i^2 - X_{i-1} $ の三角的集合に対して、$ \langle \mathbf{T} \rangle $ での乗算は $ O(\tilde{\delta}) $ 時間で実行可能であり、ここで $ \delta = 2^n $ は入力／出力サイズを表す。
• 次数 $ n $ の一変数多項式 $ F $ から導かれるコーシー加群に対して、$ \langle \mathbf{T} \rangle $ での乗算も $ O(\tilde{\delta}) $ であり、$ \delta = n! $ である。
• 直接乗算（コスト $ \delta^{\log_2 3} $）による指数的ブロードキャストを回避するため、ホモトピーと冪級数上での評価・補間技術を用いる。
• 小特異性において、$ p^\alpha $ を法とする冪和をリフトし、ニュートンの恒等式を用いて最終多項式を回復することで、準線形ビット計算量を達成する。
• $ \eta $-進展開における精度制御が、$ \mathbf{T} $ の単項式サポートに依存することが示され、効率的な計算が可能になる。
• 実験結果により、アルゴリズムの実用的効率が確認され、計測時間は準線形複雑度の推定と一致し、Magmaの結果に基づく実装よりも2桁の速度向上を達成した。