[论文解读] Efficient Random-Walk Methods for Approximating Polytope Volume

Ioannis Z. Emiris, Vissarion Fisikopoulos|arXiv (Cornell University)|Dec 10, 2013

Computational Geometry and Mesh Generation参考文献 40被引用 6

一句话总结

本论文提出了一种高效实现Hit-and-run马尔可夫链蒙特卡洛方法的方案，用于近似高维凸多面体的体积，通过优化随机游走参数、迭代条件数调整以及利用射线追踪技术的快速边界预言机。该方法首次实现了对维度高达15的Birkhoff多面体体积的实用近似，其速度与精度显著优于现有精确与近似方法。

ABSTRACT

We experimentally study the fundamental problem of computing the volume of a convex polytope given as an intersection of linear inequalities. We implement and evaluate practical randomized algorithms for accurately approximating the polytope's volume in high dimensions (e.g. one hundred). To carry out this efficiently we experimentally correlate the effect of parameters, such as random walk length and number of sample points, on accuracy and runtime. Moreover, we exploit the problem's geometry by implementing an iterative rounding procedure, computing partial generations of random points and designing fast polytope boundary oracles. Our publicly available code is significantly faster than exact computation and more accurate than existing approximation methods. We provide volume approximations for the Birkhoff polytopes B_11,...,B_15, whereas exact methods have only computed that of B_10.

研究动机与目标

• 解决以H-多面体（半空间交集）表示的高维凸多面体体积计算问题，该问题为#P难问题，精确方法在维度超过15后即不可行。
• 通过优化随机游走长度和采样点数量等关键参数，改进现有随机近似算法，以提升准确率与效率。
• 开发并实现基于超平面查询的射线追踪技术的快速边界预言机，支持高维环境下的可扩展采样。
• 首次提供Birkhoff多面体B₁₁至B₁₅体积的实用近似，扩展了此前仅能精确计算B₁₀的范围。

提出的方法

• 采用Hit-and-run随机游走算法，通过选择随机方向并在多面体边界交线处均匀采样，生成近似均匀的样本。
• 通过射线追踪实现边界预言机，计算与面超平面的交点，降低随机游走步骤中成员查询的开销。
• 引入迭代条件数调整过程，改善多面体的条件数，从而提升马尔可夫链的混合时间。
• 通过经验相关性优化随机游走长度与采样点数量，平衡准确率与运行时间。
• 结合几何与算法技术，包括时空权衡及近似最近邻数据结构，加速预言机查询。
• 采用采样双重方法，利用极点查询与边界预言机之间的对偶性，并应用ε-core-set进行近似宽度估计。

实验结果

研究问题

• RQ1 能否在实践中通过随机算法实现高维H-多面体的高效且准确的体积近似，超越理论复杂度界限？
• RQ2 随机游走长度与采样点数量等关键参数如何影响高维体积估计中准确率与运行时间的权衡？
• RQ3 迭代条件数调整与快速边界预言机能否显著提升高维环境中Hit-and-run采样的收敛性与性能？
• RQ4 现有精确体积计算软件在多维扩展方面能延伸至多大维度？近似方法在实践中表现如何？
• RQ5 能否通过实用化的随机算法实现对先前难以处理的Birkhoff多面体（如B₁₁至B₁₅）的高置信度、低相对误差近似？

主要发现

• 所提出的实现方案在不到一小时内以平均误差≤2%计算出Birkhoff多面体B₁₀的体积，显著优于先前研究中耗时17年的结果。
• 首次对Birkhoff多面体B₁₁至B₁₅的体积进行了近似，提供区间估计，当ε=1时，B₁₀的误差为0.7%（ε=0.5时）。
• 该方法的Hit-and-run步数是[13]中Matlab实现的100倍，且在d>70时，cube-d的误差降低4至100倍。
• 代码在2秒内以小于1%的误差近似计算cross-10的体积，rh-15-60在3.44秒内完成，优于LRS、Normaliz和Qhull等精确工具（后者失败或超时）。
• 在d<80时，该实现比[13]快100倍以上，并在立方体与Birkhoff多面体上均实现了显著更高的准确率，尤其在高维中表现更优。
• 该方法首次实现了对B₁₀以外Birkhoff多面体体积的实用近似，计算结果优于目前已知的最佳渐近上界。