[論文レビュー] Computing Least Fixed Points of Probabilistic Systems of Polynomials

Javier Esparza, Andreas Gaiser+1 more|arXiv (Cornell University)|Dec 21, 2009

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

ひとこと要約

本稿では、確率的多項式系（PSP）の最小不動点がすべて1のベクトルに等しいかどうかを決定する強多項式アルゴリズムを提示し、正確かつ非正確な算術のハイブリッドを用いて、最小不動点に対する信頼性があり線形収束する下界と上界を計算する、初めての実用的アルゴリズムを導入する。この手法は、数値的不安定性を効率的に扱いながらも正しさを保証し、分岐過程モデルにおける実験で、従来のLPベースの手法を上回る性能を発揮する。

ABSTRACT

We study systems of equations of the form X1 = f1(X1, ..., Xn), ..., Xn = fn(X1, ..., Xn), where each fi is a polynomial with nonnegative coefficients that add up to 1. The least nonnegative solution, say mu, of such equation systems is central to problems from various areas, like physics, biology, computational linguistics and probabilistic program verification. We give a simple and strongly polynomial algorithm to decide whether mu=(1, ..., 1) holds. Furthermore, we present an algorithm that computes reliable sequences of lower and upper bounds on mu, converging linearly to mu. Our algorithm has these features despite using inexact arithmetic for efficiency. We report on experiments that show the performance of our algorithms.

研究の動機と目的

• 分岐過程における絶滅確率や確率的プログラムの終了性を特定する上で極めて重要な、確率的多項式系（PSP）の最小不動点が(1,…,1)に等しいかどうかを決定する強多項式アルゴリズムの開発。
• 非正確な浮動小数点算術を用いても、真の値に線形に収束する信頼性のある下界と上界の列を計算する実用的アルゴリズムの設計。
• PSPに対するニュートン法の数値的不安定性を、正しさを保証する正確算術と効率性を求める浮動小数点算術を組み合わせることで克服。
• 爆発または非終了の閾値が一貫性の境界に位置する、ニュートロン分岐過程や類似モデルにおける臨界半径の効率的計算を可能にする。
• 実世界のPSPインスタンスを用いた実験を通じて、提案手法が従来のLPベースの手法に比べて実用的で優れていることを示す。

提案手法

• 一貫性の検査は、問題を一次方程式系に還元することで、正確な線形計画法の計算コストを回避する強多項式アルゴリズムを可能にする。
• 境界の計算では、ニュートン法に基づく新しい反復的手順を用いるが、正確算術を組み合わせることで正しさを保ちつつ、速度向上のために浮動小数点算術を活用する。
• 区間算術と正確な有理数算術を組み合わせることで、各反復で保証された下界と上界を維持する。
• 計算中に精度（MapleのDigits）を動的に調整し、正しさを保つために必要な場合にのみ増加させるため、効率性と信頼性のバランスを図る。
• 収縮写像の原理を適用することで、下界と上界が最小不動点に向かって線形収束することを保証する。
• 実験では、離散化されたニュートロン分岐過程と確率的文脈自由文法を用い、実際のPSPインスタンスにおける性能を評価する。

実験結果

リサーチクエスチョン

• RQ1 確率的多項式系（PSP）の一貫性（すなわち、その最小不動点が(1,…,1)に等しいかどうか）を強多項式時間で決定できるか？
• RQ2 非正確算術を用いても、信頼性があり線形に収束する下界と上界の列を計算する実用的アルゴリズムを設計できるか？
• RQ3 提案された一貫性アルゴリズムは、実世界のPSPにおいて、従来のLPベースの手法と比較して速度とスケーラビリティの面で優れているか？
• RQ4 固定点計算における確率的モデルの数値的効率性と正しさの両立を実現するため、ハイブリッドな正確/非正確算術をどの程度活用できるか？
• RQ5 ニュートロン分岐過程における臨界半径とは何か？また、提案されたアルゴリズムを用いて効率的に計算できるか？

主な発見

• 提案された一貫性アルゴリズムは、ニュートロン分岐過程インスタンスにおいて、LPベースの手法と比較して最大10倍の高速化を達成し、テストケースでは約10倍のスピードアップを示した。
• 変数数n=150のPSPにおいて、臨界半径は[2.981, 2.991]の区間に位置しており、ハリスの推定値2.9に非常に近い結果を得た。
• n=20, 50, 100の各ケースで、一貫性アルゴリズムは1秒未塔で完了したが、正確なLP手法はn=100のケースで最大258秒を要した。これにより、スケーラビリティの優位性が明確になった。
• 境界計算アルゴリズムは線形収束を達成し、有効なビット数が反復回数に比例して線形に増加した。すべてのテストケースで、最大2回の精度上昇（Digits=20に）が必要だった。
• すべてのテストインスタンスにおいて、最初の数回の反復内で上界が1未満に厳密に計算され、一貫性の欠如が証明された。これにより、正しさと効率性が確認された。
• 7変数の小さなシステムですら、浮動小数点ニュートン法で一般的に生じる誤った正の判定（μ < 1 なのに μ ≥ 1 と誤って報告）を、本手法は効果的に回避した。