[論文レビュー] A baby step-giant step roadmap algorithm for general algebraic sets

Saugata Basu, Marie-Françoise Roy+2 more|arXiv (Cornell University)|Jan 31, 2012

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

ひとこと要約

本稿では、実閉体上の一般代数的集合に対して、決定的ベビーステップ・ギアントステップルートマップアルゴリズムを提示し、$ d^{O(k\bar{d})} $ の複雑性を達成する。ここで $ d = \deg(Q) $ であり、$ k $ は次元である。無限小変形技術と次元を $ \sqrt{k} $ ずつ削減する再帰的スキームを組み合わせることで、一般性の仮定を回避し、従来の $ d^{O(k^2)} $ のルートマップ計算の境界を改善する。

ABSTRACT

Let $\mathrm{R}$ be a real closed field and $\mathrm{D} \subset \mathrm{R}$ an ordered domain. We give an algorithm that takes as input a polynomial $Q \in \mathrm{D}[X_1,\ldots,X_k]$, and computes a description of a roadmap of the set of zeros, $\mathrm{Zer}(Q,\mathrm{R}^k)$, of $Q$ in $\mathrm{R}^k$. The complexity of the algorithm, measured by the number of arithmetic operations in the ordered domain $\mathrm{D}$, is bounded by $d^{O(k \sqrt{k})}$, where $d = \mathrm{deg}(Q)\ge 2$. As a consequence, there exist algorithms for computing the number of semi-algebraically connected components of a real algebraic set, $\mathrm{Zer}(Q,\mathrm{R}^k)$, whose complexity is also bounded by $d^{O(k \sqrt{k})}$, where $d = \mathrm{deg}(Q)\ge 2$. The best previously known algorithm for constructing a roadmap of a real algebraic subset of $\mathrm{R}^k$ defined by a polynomial of degree $d$ has complexity $d^{O(k^2)}$.

研究の動機と目的

• 多項式 $ Q \in \mathrm{D}[X_1,\dots,X_k] $ によって定義される実代数的集合のルートマップを計算する決定的アルゴリズムを開発すること。ここで $ \mathrm{D} $ は順序付き整域である。
• 従来の確率的アルゴリズムが必要としていた一般座標の選択を避けること。これらの座標は、複雑性の境界内では決定的に計算することが困難である。
• ベビーステップ・ギアントステップ戦略を用いて、再帰の深さを $ k $ から $ \sqrt{k} $ に削減することで、ルートマップ構築の漸近的複雑性を向上させること。
• 複雑な代数的幾何学に依存せずに、無限小変形を用いて座標の良好な性質を保ちながら、正しさと効率性を維持すること。

提案手法

• 再帰的深さを $ \sqrt{k} $ に制限するため、各再帰レベルで次元を $ \sqrt{k} $ ずつ削減するベビーステップ・ギアントステップ再帰的スキームを用いる。
• 単一の無限小を用いて入力多様体に無限小変形を適用し、一般射影を必要とせずに元の座標が計算に適切なままであることを保証する。
• パラメータ化された有界代数的サンプリングと曲線の再パラメータ化に依存して、無限小変形下での曲線セグメントの極限を計算する。
• 変形された構造上で三角符号決定と制限付き消去を用いてルートマップ成分を計算し、半代数的性質を維持する。
• 適切にパラメータ化された曲線セグメントを用い、$ \lim_{\varepsilon} \gamma(\varepsilon) $ を用いて極限を計算することで、変形された設定から実際のルートマップを再構成する。
• 次数の上限と複雑性の見積もりは、先行アルゴリズム（例：[2, Algorithm 12.18]）から得られ、擬似還元および消去技術と組み合わせることで、全体の複雑性境界を保証する。

実験結果

リサーチクエスチョン

• RQ1 一般代数的集合の決定的ルートマップアルゴリズムを、従来の $ d^{O(k^2)} $ の境界を上回る複雑性で設計することは可能か？
• RQ2 従来のルートマップアルゴリズムにおける一般座標への依存を、複雑性境界を損なわずに排除することは可能か？
• RQ3 既存のアルゴリズムにおける再帰の深さ $ k $ を $ \sqrt{k} $ に削減することは、正しさと効率性を保ちながら可能か？
• RQ4 無限小変形を用いて、一般性の仮定を回避しながらも、同じ漸近的複雑性クラスを維持することは可能か？
• RQ5 $ k $ 変数の次数 $ d $ の多項式によって定義される実代数的集合のルートマップを計算する際の、最もタイトな複雑性境界は何か？

主な発見

• アルゴリズムは $ \mathrm{Zer}(Q,\mathrm{R}^k) $ のルートマップを、$ d = \deg(Q) \geq 2 $ に対して $ d^{O(k\sqrt{k})} $ の時間内に計算可能であり、従来の $ d^{O(k^2)} $ の境界に比べて顕著な改善を達成する。
• この方法は決定的であり、一般射影を必要とせず、単一の無限小変形によって座標の質を保つ。
• 複雑性は順序付き整域 $ \mathrm{D} $ 上での算術演算の回数として測定され、すべての実閉体に一様に成り立つ。
• ルートマップ構築により、代数的集合の各半代数的連結成分が、ルートマップと半代数的に連結な部分集合で交わることが保証される。
• このアルゴリズムにより、$ \mathrm{Zer}(Q,\mathrm{R}^k) $ の半代数的連結成分の数の計算が、$ d^{O(k\sqrt{k})} $ で効率的に可能である。
• 半代数的幾何学の技術と再帰的分解を組み合わせることで、従来の再帰的手法よりも優れた漸近的性能を達成する。