[论文解读] On the Computation of Hilbert Bases and Extreme Rays of Cones

Raymond Hemmecke|ArXiv.org|Mar 12, 2002

Commutative Algebra and Its Applications参考文献 7被引用 5

一句话总结

本文提出了一种新颖的投影-提升算法，用于计算任意整数格子 $\Lambda \subseteq \mathbb{Z}^n$ 的幺半群 $\Lambda \cap \mathbb{R}^n_+$ 的极小生成集（Hilbert 基），并将其扩展以计算有理锥的极射线。该方法通过逐次投影到低维空间并回溯提升来逐步构建解，避免引入额外变量，并支持早期剪枝，从而得到的中间超集显著小于以往方法。其主要贡献是一种高效且无需变量的算法，在 Hilbert 基与极射线计算等任务中，展现出优于现有技术的计算效率与可扩展性。

ABSTRACT

In this paper we present a novel project-and-lift approach to compute the set of minimal generators of the semigroup $(Λ\cap\R^n_+,+)$ for lattices $Λ\subseteq\Z^n$. This problem class includes the computation of Hilbert bases of cones $\{z:Az=0,z\in\R^n_+\}$ for integer matrices $A$. A similar approach can be used to compute only the extreme rays of such cones. Finally, some combinatorial applications and computational experience are presented.

研究动机与目标

• 开发一种高效算法，用于计算任意格子 $\Lambda \subseteq \mathbb{Z}^n$ 的幺半群 $\Lambda \cap \mathbb{R}^n_+$ 的极小生成集（Hilbert 基）。
• 提供一种方法，用于计算由 $\{z : Az = 0, z \in \mathbb{R}^n_+\}$ 定义的 pointed 有理锥的极射线。
• 通过避免引入辅助变量，并在计算过程中实现对冗余向量的早期剔除，提升计算效率。
• 将该算法应用于组合问题，如计数魔方阵与 t-设计，并通过截断 Graver 基在整数规划中加以应用。

提出的方法

• 该算法采用投影-提升策略，递归地计算前 $j$ 个坐标投影的极小生成集 $H^+_j$。
• 在每一步 $j$，通过在从 $H^+_j$ 和格子生成元导出的精心选择的输入集 $F$ 上执行完成过程，计算 $H^+_{j+1}$。
• 对于 $j < s$，输入集 $F$ 包含下一个格子生成元及其负向量的投影，以应对非单射投影的情况。
• 引入一种 S-向量的选择策略，显著加速了收敛过程。
• 该方法避免引入额外变量，并允许在计算早期阶段剔除非极小或不可行的向量。
• 该算法以 C 语言实现为 MLP 软件，支持完整与截断 Hilbert 基的计算。

实验结果

研究问题

• RQ1 投影-提升算法是否能在不引入辅助变量的前提下，比现有方法更高效地计算 $\Lambda \cap \mathbb{R}^n_+$ 的 Hilbert 基？
• RQ2 能否使用类似的提升策略计算有理锥 $\{z : Az = 0, z \in \mathbb{R}^n_+\}$ 的极射线？
• RQ3 在计算极小生成集时，智能的 S-向量选择策略对完成过程性能有何影响？
• RQ4 在 Hilbert 基计算中，对冗余或不可行向量的早期剪枝在多大程度上能减小中间超集的规模？
• RQ5 该算法能否有效应用于组合问题，如计数魔方阵与 t-设计？

主要发现

• 投影-提升算法在不引入辅助变量的情况下，计算了任意格子 $\Lambda \subseteq \mathbb{Z}^n$ 的 $\Lambda \cap \mathbb{R}^n_+$ 的 Hilbert 基，显著降低了计算开销。
• 该算法成功计算了所有 6×6 魔方阵的 Hilbert 基，得到 522,347 个极小生成元，耗时约 10 天（在 1GHz PC 上，配备 4GB 内存）。
• 在完成过程中引入选择性 S-向量策略，显著加快了收敛速度。
• 该方法可高效计算截断 Hilbert 基，对整数规划与组合计数具有实用价值。
• 实现代码 MLP 已公开发布，并成功用于求解魔方阵计数等非平凡组合问题。
• 该算法还提供了一种高效计算对偶锥极射线的方法，其对应于原锥的面法向量。