[论文解读] 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 向量及无界多面体的情形。
• 基于熵近似方法，为多指标运输多面体与列联表提供渐近公式。
• 建立一个基于最大熵分布与局部中心极限定理的严格理论框架，以保证估计的准确性。

提出的方法

• 在离散集 X ⊂ ℝⁿ（如 ℤⁿ₊ 或 {0,1}ⁿ）上构造一个概率测度 μ，使得 μ 在 P ∩ X 上为常数，且其期望位于 P 内。
• 将 μ 的构造表述为一个凸优化问题——在矩约束下进行熵最大化，可通过内点法求解。
• 利用最大熵分布通过中心极限定理近似 P ∩ X 上的计数测度，从而实现高斯型尾部估计。
• 应用傅里叶分析与复积分表示，将整数点数量表示为包含生成函数的振荡积分。
• 通过二次型的特征值估计与特征函数衰减性，控制积分近似中的误差。
• 将积分区域划分为内区、中区与外区，以控制误差项，并推导出具有显式误差率的渐近界。

实验结果

研究问题

• RQ1 是否可以使用最大熵分布来高精度地近似多面体中整数点的数量？
• RQ2 在离散集上进行熵最大化如何导致多面体中体积与点数的计算可处理的近似？
• RQ3 由最大熵分布导出的高斯近似在估计多面体上指数和时的准确性如何？
• RQ4 该方法能否为多指标运输多面体与列联表提供渐近公式？
• RQ5 当多面体无界或集合 X 为布尔立方体时，近似的误差界是什么？

主要发现

• 整数点数量 |P ∩ ℤⁿ| 近似为 e^{g(z)} / (2π)^d 乘以一个平行多面体上的积分，该积分通过熵与特征值相关的衰减估计进行有界控制。
• 该方法实现了形如 (1 + (2/5)απ²)^{-m_k} 的误差界，其中 m_k ≥ τ_k² / (ρ_k π²)，ρ_k 为二次型 ψ_k 的最大特征值。
• 对于积分的内区，当对所有 j 满足 |⟨a_j, t⟩| ≤ ϵ/(10σ(2ζ_j + 1)) 时，相位函数 f(t) 的误差被控制在 ϵ/30 以内。
• 通过设置 σ = 4d + 10 ln(1/ϵ)，确保在中心区域的高概率集中。
• 该方法在满足温和正则性条件下，为多维列联表的数量与运输多面体的体积提供了渐近公式。
• 该方法对无界多面体也具有鲁棒性，只要指数和或积分收敛即可。