[论文解读] Barvinok's Rational Functions: Algorithms and Applications to Optimization, Statistics, and Algebra

Yoshida, Ruriko|ArXiv.org|Jun 15, 2004

Commutative Algebra and Its Applications参考文献 53被引用 4

一句话总结

本文提出基于Barvinok有理函数的算法，以高效计算有理凸多面体中的格点数量，并求解混合整数规划问题。通过将生成函数编码为短有理函数，该方法实现了在固定维度下对环理想、格伯纳基、Ehrhart系列及整数规划最优解的多项式时间计算，实验结果通过LattE软件套件得到验证。

ABSTRACT

The main theme of this dissertation is the study of the lattice points in a rational convex polyhedron and their encoding in terms of Barvinok's short rational functions. The first part of this thesis looks into theoretical applications of these rational functions to Optimization, Statistics, and Computational Algebra. The main theorem on Chapter 2 concerns the computation of the \emph{toric ideal} $I_A$ of an integral $n imes d$ matrix $A$. We encode the binomials belonging to the toric ideal $I_A$ associated with $A$ using Barvinok's rational functions. If we fix $d$ and $n$, this representation allows us to compute a universal Gröbner basis and the reduced Gröbner basis of the ideal $I_A$, with respect to any term order, in polynomial time. We derive a polynomial time algorithm for normal form computations which replaces in this new encoding the usual reductions of the division algorithm. Chapter 3 presents three ways to use Barvinok's rational functions to solve Integer Programs. The second part of the thesis is experimental and consists mainly of the software package { t LattE}, the first implementation of Barvinok's algorithm. We report on experiments with families of well-known rational polytopes: multiway contingency tables, knapsack type problems, and rational polygons. We also developed a new algorithm, {\em the homogenized Barvinok's algorithm} to compute the generating function for a rational polytope. We showed that it runs in polynomial time in fixed dimension. With the homogenized Barvinok's algorithm, we obtained new combinatorial formulas: the generating function for the number of $5 imes 5$ magic squares and the generating function for the number of $3 imes 3 imes 3 imes 3$ magic cubes as rational functions.

研究动机与目标

• 通过Barvinok有理函数编码，开发高效算法以计算有理凸多面体中的格点数量。
• 利用有理函数计算环理想及其格伯纳基，实现固定维度下的多项式时间复杂度。
• 将这些方法应用于混合整数规划、统计学（如列联表）以及代数几何问题。
• 在LattE软件系统中实现并评估这些算法，以检验其实际性能与可扩展性。
• 为Ehrhart系列计算、多面体上的优化以及整数规划的精确解提供计算框架。

提出的方法

• 利用Barvinok将多面体中格点的生成函数表示为短有理函数，通过其分解为单模锥的方法。
• 使用锥的带符号分解，表达编码Ehrhart系列与计数函数的有理函数。
• 实现同伦化的Barvinok算法与BBS（基于基的二分搜索）算法，以计算格点数量并优化多面体。
• 应用挖掘算法，通过迭代改进可行解，利用有理函数表示求解混合整数规划问题。
• 利用计算机代数系统与LattE中的C++实现，计算Ehrhart多项式、生成函数与最优值。
• 通过有理函数的泰勒级数展开，提取在缩放下对应格点数量的系数。

实验结果

研究问题

• RQ1 对于固定维度，Barvinok的有理函数能否在多项式时间内计算环理想的一个通用格伯纳基？
• RQ2 利用有理函数编码与锥分解，有理多面体中格点计数的计算效率如何？
• RQ3 BBS与挖掘算法在求解混合整数规划问题时，相较于传统方法的性能表现如何？
• RQ4 即使完整化简代价高昂，能否通过有理函数表示精确计算多面体的Ehrhart系列？
• RQ5 有理函数编码在多大程度上可实现对背包问题与运输多面体的精确优化？

主要发现

• 对于固定的 $ n $ 和 $ d $，与 $ n \times d $ 整数矩阵 $ A $ 相关的环理想 $ I_A $ 可通过Barvinok的有理函数编码，实现其通用与约化格伯纳基的多项式时间计算。
• LattE成功计算出 $ 4 \times 4 $ 幻方多面体的Ehrhart系列为 $ \frac{t^8 + 4t^7 + 18t^6 + 36t^5 + 50t^4 + 36t^3 + 18t^2 + 4t + 1}{(-1+t)^4(-1+t^2)^4} $，其系数与已知计数结果在 $ t^{15} $ 以内一致。
• 24胞体恰好包含33个格点，其中仅1个位于其内部，该结果通过标准与同伦化Barvinok算法均得到验证。
• 通过挖掘算法求得背包问题"cuww1"的最优解为 $ (7334, 0, 0, 0, 0) $，最优值为1,562,142，且经BBS算法验证。
• BBS算法共执行26次迭代，使用125,562个单模锥，计算最优值，展现出高计算成本下的可扩展性。
• 内存高效的泰勒展开模式可在2秒内完成Ehrhart系列最多15项的计算，结果与先前研究一致。