[論文レビュー] Maximum entropy Gaussian approximation for the number of integer points and volumes of polytopes

Alexander Barvinok, J. A. Hartigan|ArXiv.org|Mar 30, 2009

Point processes and geometric inequalities参考文献 12被引用数 8

ひとこと要約

本稿では、多面体内の整数点の数および体積を推定するために、定常質量を多面体上に持つ離散集合（例：ℤⁿ₊ や {0,1}ⁿ）上の確率分布を構築することにより、最大エントロピー正規近似法を導入する。エントロピー最大化と中心極限定理の議論を活用することで、多インデックス輸送多面体および連続表の体積および数え上げに対して、効率的かつ計算可能な近似が得られる。

ABSTRACT

We describe a maximum entropy approach for computing volumes and counting integer points in polyhedra. To estimate the number of points from a particular set X in R^n in a polyhedron P in R^n, by solving a certain entropy maximization problem, we construct a probability distribution on the set X such that a) the probability mass function is constant on the intersection of P and X and b) the expectation of the distribution lies in P. This allows us to apply Central Limit Theorem type arguments to deduce computationally efficient approximations for the number of integer points, volumes, and the number of 0-1 vectors in the polytope. As an application, we obtain asymptotic formulas for volumes of multi-index transportation polytopes and for the number of multi-way contingency tables.

研究の動機と目的

• 多面体内の整数点の数および体積を推定する計算的に効率的な手法の開発。
• 高次元多面体における整数点の数え上げおよび体積計算の計算困難性の解決。
• 0-1ベクトルや非有界多面体を含む、多面体上での指数和および積分の推定への手法の拡張。
• エントロピーに基づく近似を用いた、多インデックス輸送多面体および連続表の漸近公式の確立。
• 推定精度のための、最大エントロピー分布および局所中心極限定理に基づく厳密な理論的枠組みの構築。

提案手法

• P ∩ X 上で定数で、期待値が P 内にあるような、離散集合 X ⊂ ℝⁿ（例：ℤⁿ₊ または {0,1}ⁿ）上に確率測度 μ を構築する。
• μ の構築を凸最適化問題として定式化する——モーメント制約下でのエントロピー最大化——これは内点法により解ける。
• 最大エントロピー分布を用いて、中心極限定理に基づき、P ∩ X 上の数え上げ測度をガウス型尾確率推定で近似する。
• 生成関数を含む振動的積分として整数点の数を表現するためのフーリエ解析および複素積分表現を適用する。
• 2次形式の固有値に基づく推定と特性関数の減衰を用いて、積分近似の誤差を評価する。
• 積分領域を内側、中間、外側の領域に分割し、誤差項を制御し、明示的な誤差率を伴う漸近的境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 最大エントロピー分布を用いることで、多面体内の整数点の数を高い精度で近似できるか？
• RQ2 離散集合上でのエントロピー最大化は、多面体内の体積および点数の計算可能な近似をどのように導くか？
• RQ3 最大エントロピー分布から導かれるガウス近似の、多面体上での指数和推定における精度はいかほどか？
• RQ4 この手法により、多インデックス輸送多面体および連続表のための漸近公式を得られるか？
• RQ5 多面体が非有界である場合、または集合 X がブール立方体である場合の近似誤差の境界は何か？

主な発見

• 整数点の数 |P ∩ ℤⁿ| は、e^{g(z)} / (2π)^d に平行六面体上での積分を乗じた形で近似され、積分はエントロピーおよび固有値に基づく減衰推定で有界化される。
• 誤差境界は、(1 + (2/5)απ²)^{-m_k} の形を取り、m_k ≥ τ_k² / (ρ_k π²) のとき成り立つ。ここで ρ_k は2次形式 ψ_k の最大固有値である。
• 積分の内側領域では、すべての j に対して |⟨a_j, t⟩| ≤ ϵ/(10σ(2ζ_j + 1)) が成り立つとき、位相関数 f(t) の誤差は ϵ/30 以下に抑えられる。
• 誤差は σ = 4d + 10 ln(1/ϵ) により制御され、中心領域における高確率集中を保証する。
• この手法により、ややいびつな正則性条件のもとで、多方向連続表の数および輸送多面体の体積の漸近公式が得られる。
• 指数和または積分が収束する限り、非有界多面体に対してもこのアプローチはロバストである。