[論文レビュー] Algebraic osculation and factorization of sparse polynomials

Martin Weimann|ArXiv.org|Apr 1, 2009

Polynomial and algebraic computation参考文献 28被引用数 3

ひとこと要約

本稿は、多項式のニュートン多面体の幾何学的性質を活用して、代数的接触を用いたトロイカル多様体における代数的接触を活用し、留数基準を用いて既約因子を検出することで、スパースな二変数多項式の絶対的因数分解のための新規なアルゴリズムを提案する。主な貢献は、Galligo-Rupprecht法などの古典的手法と比較して再結合の回数を減らす、幾何学的認識を持つ決定的Henselの上昇法である。

ABSTRACT

We prove a theorem on algebraic osculation and we apply our result to the Computer Algebra problem of polynomial factorization. We consider X a smooth completion of the complex plane and D an effective divisor supported on the boundary of X. Our main result gives explicit conditions equivalent to that a given Cartier divisor on D extends to X. These osculation criterions are expressed with residues. We derive from this result a toric Hensel lifting which permits to compute the absolute factorization of a bivariate polynomial by taking in account the geometry of its Newton polytope. In particular, we reduce the number of possible recombinations when compared to the Galligo-Rupprecht algorithm.

研究の動機と目的

• 既存のLRアルゴリズムにおけるスパース多項式に対する一般座標変換の非効率性を是正すること。
• ニュートン多面体の幾何学的性質を活用して、絶対的因数分解における再結合ケースの数を削減すること。
• 確率的仮定や漸近的近似を避ける決定的アルゴリズムを構築すること。
• トロイカル解消を用いて、特異曲線に対しても留数基準の拡張を適用可能にする。
• 境界交差データとCartier除法子の拡張を因数分解に用いる理論的基盤を提供すること。

提案手法

• 曲線 $C$ を多項式 $f$ で定義するため、$\mathbb{C}^2$ の滑らかなトロイカルコンactification $X$ を用い、ニュートン多面体の情報を保持する。
• すべての $\Psi \in H^0(X, \Omega_X^2(D))$ に対して $\langle \gamma, \Psi \rangle = 0$ であることを条件として、境界 $D = \partial X$ 上のCartier除法子 $\gamma$ が $X$ に拡張可能かどうかを判定する。
• コhomological障害を計算するために、Grothendieck留数とDolbeault $\bar{\partial}$-解消を用いる。
• 因数分解問題を、$\operatorname{Card}(Q \cap \mathbb{Z}^2)$ 個の未知数をもつ、サイズ $2\operatorname{Vol}(Q) + \deg(\gamma' \cdot \partial X)$ のスパース線形方程式系に還元する。
• Minkowski和分解 $N_f = P + Q$ を用いて、ラインバンドルの分解 $\mathcal{O}_X(C) \simeq \mathcal{L}_P \otimes \mathcal{L}_Q$ を関連付ける。
• 高次隐函数計算（最大 $na^2$ 階まで）を含むトロイカルHenselの上昇法を適用し、因数の決定的検出を保証する。

実験結果

リサーチクエスチョン

• RQ1 留数に基づく代数的接触基準を用いて、スパースな二変数多項式の既約因子を検出できるか。
• RQ2 ニュートン多面体の幾何学的性質をどのように活用すれば、絶対的因数分解における再結合の数を削減できるか。
• RQ3 トロイカル表面の境界除法子 $D$ 上のCartier除法子が、表面全体に拡張可能となる正確な条件は何か。
• RQ4 本手法は、古典的Henselの上昇法やGalligo-Rupprechtアルゴリズムと比較して、計算量と精度の点でどのように異なるか。
• RQ5 境界上に非トロイカル特異点がある場合、それが因数分解における再結合の削減にどのように寄与できるか。

主な発見

• 境界除法子 $D = \partial X$ 上のCartier除法子 $\gamma$ が $X$ に拡張可能であるための必要十分条件は、すべての $\Psi \in H^0(X, \Omega_X^2(D))$ に対して $\langle \gamma, \Psi \rangle = 0$ であることであり、留数に基づく拡張基準が得られる。
• 本アルゴリズムは、一般座標ではなく多面体幾何学を用いることで、Galligo-Rupprecht法と比較して再結合の可能性を削減する。
• 二重次数-$(a,b)$ の多項式に対して、本手法は $b$ 階までと $a$ 階までのそれぞれ $a$ 個と $b$ 個の隐函数を計算するが、GRアルゴリズムは $a+b$ 個の関数を $a+b$ 階まで計算する。
• 本手法では、隐函数の必要精度が $na^2$ に達する場合があり、古典的Henselの上昇法の $2na = \deg(f)$ の境界を著しく上回る。
• フェースト多項式が平方自由である場合、本アルゴリズムは留数条件による決定的検出を備えたゼロ和型LRアルゴリズムに帰着する。
• 境界上に特異点を有する曲線に対しては、非トロイカル解消 $\widetilde{X}$ を用いることで定理1を適用可能となり、追加の特異除法子が再結合数をさらに削減する。