[論文レビュー] A bound on the minimum of a real positive polynomial over the standard simplex

Saugata Basu, Richard Leroy+1 more|ArXiv.org|Feb 19, 2009

Advanced Optimization Algorithms Research参考文献 5被引用数 12

ひとこと要約

本稿では、実代数幾何学的手法を用いて、標準単体上での正の多変数多項式の最小値に対する明示的かつ非退化な下界を提示する。結果式とシルベスター行列を用いた単変数結果の一般化により、結果式多項式にコーシーの不等式を適用することで、非退化性の仮定を必要とせず、次元 $k$、次数 $d$、係数のビットサイズ $ au$ のみに依存する計算可能な下界を確立する。

ABSTRACT

We consider the problem of bounding away from 0 the minimum value m taken by a polynomial P of Z[X_1,...,X_k] over the standard simplex, assuming that m>0. Recent algorithmic developments in real algebraic geometry enable us to obtain a positive lower bound on m in terms of the dimension k, the degree d and the bitsize of the coefficients of P. The bound is explicit, and obtained without any extra assumption on P, in contrast with previous results reported in the literature.

研究の動機と目的

• 標準単体 $\Delta \subset \mathbb{R}^k$ 上での正の多項式 $P \in \mathbb{Z}[X_1,\dots,X_k]$ の最小値 $m > 0$ に対する明示的下界を提供すること。
• 従来の手法が要請していた非退化性や特別な仮定を排除すること。
• 次元 $k$、次数 $d$、係数のビットサイズ $ au$ のみに依存する下界 $m_{k,d,\tau}$ を導出すること。
• アルゴリズム的実代数幾何学的手法を用いて、明示的かつ計算可能で普遍的な下界を提供することで、先行研究を改善すること。

提案手法

• 最小値が達成される臨界点を特定するため、単変数の場合の一般化として、結果式 $R(Z) = \operatorname{Res}_T(P(T) - Z, P'(T))$ を分析する。
• シルベスター行列の構造とハダマールの行列式上限を用いて、結果式多項式 $R(Z)$ の係数を制御する。
• 結果式多項式にコーシーの不等式を適用し、$m$（正の最小根）の下界を導出する。
• 単変数アプローチを単体の各面に沿って有理関数の単変数表現（RUR）を用いて多変数多項式へ拡張する。
• RUR多項式の次数および係数のビットサイズ（$d_u, \tau_u, d_{P,u}, \tau_{P,u}$）を制御し、結果式のサイズを制御する。
• 再帰的上限 $D, \rho, \rho'$ を用いて、最終的な下界 $m_{k,d,\tau} = \left[2^{\rho'+1}\sqrt{dD+1}\right]^{-D}\left[2^{\rho}\sqrt{D+1}\right]^{-dD}$ を構築する。

実験結果

リサーチクエスチョン

• RQ1 非退化性の仮定を必要とせず、標準単体上での正の多項式の最小値に対する普遍的かつ明示的な下界を導出可能か？
• RQ2 標準単体上での正の多項式の最小値は、次数 $d$、次元 $k$、係数のビットサイズ $ au$ とどのように関係するか？
• RQ3 アルゴリズム的実代数幾何学的手法を用いて得られる、$m = \min_{\Delta} P$ の最もタイトな明示的下界は何か？
• RQ4 キャンニーやデ・ルアとサントスの先行研究と比較して、提案された下界は指数関数的成長やパラメータ依存性においてどのように定量的に優れているか？

主な発見

• 本稿では、$D = 2d(2d-1)^{k-1}$、$\rho$、$\rho'$ が $k,d,\tau$ を用いた再帰的式で定義される明示的下界 $m_{k,d,\tau} = \left[2^{\rho'+1}\sqrt{dD+1}\right]^{-D}\left[2^{\rho}\sqrt{D+1}\right]^{-dD}$ を完全に確立する。
• 非退化性や横断的条件を要件としないが、キャンニーのギャップ定理やロジャシエヴィッチに基づく下界とは異なり、単体上でのすべての正の多項式に適用可能である。
• 非退化性の仮定を排除し、$k$、$d$、$ au$ に明示的な依存性を提供することで、先行研究を改善し、定量的推定 $\frac{1}{m_{k,d,\tau}} \leq \left(2^{\tau}\right)^{2^{k+3}d^{k+1}k}2^{2^{k+7}d^{k+2}k^{2}}$ を得る。
• 未知の普遍定数を含むデ・ルアとサントスの一般推定値よりも、指数部を明示的にすることでよりタイトな下界を達成する。
• 系 $\partial P_\sigma / \partial Y_i = 0$ を単一の式 $Q = \sum_i (\partial P_\sigma / \partial Y_i)^2 = 0$ に置き換えることで、退化ケースをカバー可能となる。
• 単変数の場合には $m > \frac{3^{d/2}}{2^{(2d-1)\tau}(d+1)^{2d-1/2}}$ という下界が得られ、これは[BCR]の最近の結果よりもわずかに優れている。