[论文解读] Software for Exact Integration of Polynomials over Polyhedra

Jesús A. De Loera, Brandon Dutra+4 more|arXiv (Cornell University)|Jul 31, 2011

Polynomial and algebraic computation参考文献 18被引用 8

一句话总结

本论文提出 LattE integrale，一种使用三角剖分和锥分解方法对凸多面体上的多项式进行精确积分的 C++ 软件实现。该方法在单纯形和一般多面体上实现了多项式时间的多项式函数积分计算，关键成果包括对投票理论多面体的精确体积计算——四候选人选举中概率为 12.27%，三候选人选举中确认为 12.33%，展示了在高维组合问题中的高效性与可扩展性。

ABSTRACT

We are interested in the fast computation of the exact value of integrals of polynomial functions over convex polyhedra. We present speed ups and extensions of the algorithms presented in previous work. We present the new software implementation and provide benchmark computations. The computation of integrals of polynomials over polyhedral regions has many applications; here we demonstrate our algorithmic tools solving a challenge from combinatorial voting theory.

研究动机与目标

• 开发一种高性能、精确的软件工具，用于计算具有有理系数的凸多面体区域上多项式的积分。
• 通过提升速度并同时支持三角剖分和锥分解方法，扩展先前的精确积分算法。
• 解决组合投票理论中的一个难题，即在多数制和排序选择制下，不同选举结果出现的概率问题。
• 展示精确积分在高维多面体（如四人和五人候选人选举产生的多面体）中的实际可扩展性。
• 提供一种计算框架，利用 Ehrhart 理论和体积计算，高效提取拟多项式的首项系数。

提出的方法

• 该软件使用两种主要分解策略：基于 Brion–Lasserre–Lawrence–Varchenko 理论的多面体三角剖分和带符号的锥分解。
• 通过基于指数积分恒等式的闭式公式，计算在单形锥上对线性形式幂的积分。
• 该方法利用了具有有理系数的单纯形上多项式积分可通过生成函数和赋值理论在多项式时间内计算的事实。
• 对于一般多面体，软件将定义域分解为单形或锥，分别计算各部分的积分，并使用包含-排除原理或带符号分解规则合并结果。
• 该实现使用高效的有理数算术数据结构，并支持以标准多面体格式输入，包含线性不等式和等式。
• 关键优化在于通过直接计算多面体在 n=1 时的体积，仅计算 Ehrhart 拟多项式的首项系数，从而避免完整拟多项式枚举。

实验结果

研究问题

• RQ1 在高维空间中，是否能高效地对非满维多面体上的多项式进行精确积分？
• RQ2 在精确积分中，三角剖分与锥分解方法在速度和数值稳定性方面如何比较？
• RQ3 在大规模选民群体中，对于四名和五名候选人，多数制和排序选择制下产生不同胜者的精确概率是多少？
• RQ4 是否可以通过体积计算直接获得 Ehrhart 拟多项式的首项系数，而无需完整枚举？
• RQ5 精确积分工具在处理具有 24 个或更多变量的组合投票理论复杂问题时，其可扩展性达到何种程度？

主要发现

• 软件 LattE integrale 在 n=1 时成功计算出四候选人投票多面体的精确体积为 2988379676768359 / 7552997065814637134660504411827077120000。
• 通过锥分解，该软件在不到一分钟内计算出四名候选人下多数制与排序选择制下不同选举结果的概率为 12.27%。
• 该方法通过直接计算对应多面体的体积并应用体积比，确认了先前计算出的三名候选人选举中 12.33% 的概率。
• 四名候选人情况下，Ehrhart 拟多项式的首项系数通过在 n=1 时多面体的体积直接提取，避免了完整拟多项式计算。
• 对于五名候选人，该多面体有 120 个变量；LRS 枚举 1250 万个顶点后被中止，表明该方法已超出当前计算能力的极限。
• 该软件表明，通过体积计算实现的精确积分比完整拟多项式枚举更高效，尤其在高维投票多面体中优势显著。