[论文解读] Computation of the highest coefficients of weighted Ehrhart quasi-polynomials of rational polyhedra

Velleda Baldoni, Nicole Berline+3 more|arXiv (Cornell University)|Nov 7, 2010

Advanced Combinatorial Mathematics参考文献 20被引用 4

一句话总结

本文提出了一种基于局部近似的算法，高效计算有理多面体加权Ehrhart准多项式的最高次系数，其中格点按多项式权重计数。该方法在Barvinok方法的基础上进行扩展，适用于加权情形，利用生成函数与阶梯多项式表示，以闭式形式计算这些系数，在基础的Maple实现中性能与最先进的软件LattE macchiato相比仍具竞争力。

ABSTRACT

This article concerns the computational problem of counting the lattice points inside convex polytopes, when each point must be counted with a weight associated to it. We describe an efficient algorithm for computing the highest degree coefficients of the weighted Ehrhart quasi-polynomial for a rational simple polytope in varying dimension, when the weights of the lattice points are given by a polynomial function h. Our technique is based on a refinement of an algorithm of A. Barvinok [Computing the Ehrhart quasi-polynomial of a rational simplex, Math. Comp. 75 (2006), pp. 1449--1466] in the unweighted case (i.e., h = 1). In contrast to Barvinok's method, our method is local, obtains an approximation on the level of generating functions, handles the general weighted case, and provides the coefficients in closed form as step polynomials of the dilation. To demonstrate the practicality of our approach we report on computational experiments which show even our simple implementation can compete with state of the art software.

研究动机与目标

• 解决在组合学、优化与统计学中核心问题的有理多面体中带多项式权重的格点计数计算挑战。
• 将Barvinok的无权重Ehrhart系数计算方法扩展至一般加权情形，其中每个格点通过多项式函数 $ h(x) $ 获得一个值。
• 开发一种高效计算加权Ehrhart准多项式最高次系数的算法，避免计算所有系数的NP难问题。
• 提供一种方法，以阶梯多项式形式表示最高次系数作为缩放参数 $ n $ 的闭式表达，实现实际且可扩展的计算。
• 通过计算实验展示该方法的实用性，其性能与最先进的工具如LattE macchiato相比具有竞争力。

提出的方法

• 该方法采用Barvinok算法的改进版本，通过处理加权格点计数的生成函数，将其适配于处理多项式权重 $ h(x) $。
• 在单形仿射锥上采用局部近似策略，通过捕捉顶点附近行为的拼接函数，对加权和的生成函数进行近似。
• 该方法引入中间生成函数，将问题分解为可管理的部分，通过积分与基于体积的公式计算最高次项。
• 关键创新在于使用阶梯多项式表示周期性系数 $ E_m(n \bmod q) $，从而实现最高次系数的闭式符号计算。
• 该算法利用Barvinok分解的原始变体，将其分解为单模锥，从而在高维单形中提高效率。
• 该方法在Maple中实现，并通过随机格点单形作为测试用例，与LattE macchiato进行验证。

实验结果

研究问题

• RQ1 能否为带多项式权重的有理多面体高效计算加权Ehrhart准多项式的最高次系数？
• RQ2 如何将Barvinok的无权重算法扩展至一般多项式权重函数 $ h(x) $ 的情形，同时保持计算效率？
• RQ3 在单形锥上对生成函数进行局部近似，在多大程度上能准确表示加权Ehrhart准多项式中的最高次系数？
• RQ4 所得系数能否以缩放参数 $ n $ 的阶梯多项式形式表示为闭式，从而实现符号与算法操作？
• RQ5 新算法在计算完整Ehrhart多项式时，与现有最先进的软件如LattE macchiato相比性能如何？

主要发现

• 当 $ k_0 $ 固定且维度为输入的一部分时，该算法可在多项式时间内计算有理多面体加权Ehrhart准多项式的前 $ k_0 $ 个系数。
• 最高次系数 $ E_{d+M} $ 恰好等于积分 $ \int_{\mathfrak{p}} h(x) \, dx $，与已知的几何解释一致。
• 在计算实验中，新方法实现了显著的速度提升：在维度11时，仅用1681个CPU秒就计算出前三项系数，而使用LattE macchiato的原始方法计算完整Ehrhart多项式则超过30分钟。
• 尽管实现基于Maple（解释型语言），其性能在维度5及以上时仍优于LattE macchiato的原始方法，表明在编译语言中进一步优化具有巨大潜力。
• 当采用Barvinok分解的截断策略（例如，将行列式上限设为1000）时，性能优势更加显著，如LattE macchiato的Primal 1000变体所示。
• 结果证实，该新算法不仅理论合理，而且在实际应用中具有竞争力，尤其在完整Ehrhart计算变得不可行的高维单形中表现突出。