[论文解读] Sparse Gröbner Bases: the Unmixed Case

Jean‐Charles Faugère, Pierre-Jean Spaenlehauer+1 more|arXiv (Cornell University)|Feb 28, 2014

Polynomial and algebraic computation参考文献 31被引用 17

一句话总结

本文引入了半群代数中的稀疏 Gröbner 基，以加速稀疏多项式系统的符号求解，尤其针对所有多项式共享同一 Newton 多面体的非混合情形。该文提出了 F5 和 FGLM 算法的稀疏变体，利用 Newton 多面体的组合性质，在实验中实现了超过 100 倍的加速，并给出了依赖于多面体体积、正则性及 Ehrhart 多项式的复杂度界。

ABSTRACT

Toric (or sparse) elimination theory is a framework developped during the last decades to exploit monomial structures in systems of Laurent polynomials. Roughly speaking, this amounts to computing in a \emph{semigroup algebra}, \emph{i.e.} an algebra generated by a subset of Laurent monomials. In order to solve symbolically sparse systems, we introduce \emph{sparse Gröbner bases}, an analog of classical Gröbner bases for semigroup algebras, and we propose sparse variants of the $F_5$ and FGLM algorithms to compute them. Our prototype "proof-of-concept" implementation shows large speed-ups (more than 100 for some examples) compared to optimized (classical) Gröbner bases software. Moreover, in the case where the generating subset of monomials corresponds to the points with integer coordinates in a normal lattice polytope $\mathcal P\subset\mathbb R^n$ and under regularity assumptions, we prove complexity bounds which depend on the combinatorial properties of $\mathcal P$. These bounds yield new estimates on the complexity of solving $0$-dim systems where all polynomials share the same Newton polytope (\emph{unmixed case}). For instance, we generalize the bound $\min(n_1,n_2)+1$ on the maximal degree in a Gröbner basis of a $0$-dim. bilinear system with blocks of variables of sizes $(n_1,n_2)$ to the multilinear case: $\sum n_i - \max(n_i)+1$. We also propose a variant of Fröberg's conjecture which allows us to estimate the complexity of solving overdetermined sparse systems.

研究动机与目标

• 通过半群代数中的单项式结构开发高效算法求解稀疏 Laurent 多项式系统。
• 解决经典 Gröbner 基方法在非密集单项式支持系统中应用时的计算瓶颈。
• 建立稀疏 Gröbner 基计算的复杂度界，其依赖于 Newton 多面体的内在组合不变量。
• 利用稀疏 Gröbner 基理论将双线性系统的已知次数界推广至多线性情形。
• 提出 Fröberg 猜想的变体，用于估计过度确定稀疏系统中的复杂度。

提出的方法

• 在由有限子集 $M \subset \mathbb{Z}^n$ 生成的仿射半群 $S$ 的代数 $k[S]$ 中定义稀疏 Gröbner 基。
• 提出 MatrixF5 算法的稀疏变体，用于在齐次代数 $k[S^{(h)}]$ 中计算稀疏 Gröbner 基，其中 $S^{(h)}$ 是由 $(\alpha,1)$（$\alpha \in M$）生成的齐次化半群。
• 开发 FGLM 算法的稀疏变体，通过单项式映射将稀疏 Gröbner 基转换为有理一元表示。
• 将“见证次数”——约化稀疏 Gröbner 基中最大稀疏次数的上界——作为关键复杂度指标。
• 利用 Newton 多面体 $\mathscr{P}$ 的组合不变量，包括其归一化体积 $\operatorname{vol}(\mathscr{P})$、Castelnuovo-Mumford 正则性 $\operatorname{reg}(k[S^{(h)}])$ 和 Ehrhart 多项式 $\operatorname{HP}_{\mathscr{P}}(\ell)$，推导复杂度估计。
• 在满足正则性假设的前提下，将算法应用于支持位于 $\mathscr{P} \cap \mathbb{Z}^n$ 的系统，确保正确性与高效性。

实验结果

研究问题

• RQ1 如何通过利用 Newton 多面体的结构，加速稀疏多项式系统的 Gröbner 基计算？
• RQ2 稀疏 Gröbner 基计算的复杂度界如何以正常格多面体 $\mathscr{P}$ 的组合不变量表示？
• RQ3 能否利用稀疏 Gröbner 基理论将双线性系统中经典的次数界 $\min(n_1,n_2)+1$ 推广至多线性情形？
• RQ4 稀疏 F5 和 FGLM 算法在实际中相比其经典版本有多大优势？
• RQ5 Fröberg 猜想的修改版本能否通过见证次数预测过度确定稀疏系统的求解复杂度？

主要发现

• 稀疏-MatrixF5 和稀疏-FGLM 算法在实际中实现了超过 100 倍的加速，部分示例中相比 FGb 中的经典 F5 提升超过 19,000 倍。
• 在 Newton 多面体为 $\mathscr{P}$ 的非混合情形下，见证次数由商环的 Hilbert 系列中第一个非正系数所界定，将双线性情形的界 $\min(n_1,n_2)+1$ 推广至多线性情形的 $\sum n_i - \max(n_i) + 1$。
• 稀疏 Gröbner 基计算的复杂度界依赖于归一化体积 $\operatorname{vol}(\mathscr{P})$、Castelnuovo-Mumford 正则性 $\operatorname{reg}(k[S^{(h)}]) = n+1 - \ell$ 和 Ehrhart 多项式 $\operatorname{HP}_{\mathscr{P}}(\ell)$。
• 所提出的 Fröberg 猜想变体预测：$\mathbb{C}[\mathscr{P}]$ 中的半正则序列产生 Hilbert 系列 $\operatorname{HS}_{\mathbb{C}[\mathscr{P}]/\langle f_1,\ldots,f_m \rangle}(t) = \left[ \operatorname{HS}_{\mathscr{P}}(t) \prod_{i=1}^m (1 - t^{d_i}) \right]_+$，且见证次数由第一个零系数的索引界定。
• 在单项式支持稀疏的少项式系统中，稀疏算法保持正确性，并实现高达 16,800 倍的加速，尽管由于非正常半群的存在，复杂度分析不适用。
• 通过稀疏 FGLM 获得的有理一元表示将解表示为一元多项式 $Q(T)$ 的根，其中单项式 $X^\alpha$ 映射为有理函数 $Q_\alpha(T)$，从而实现高效的解恢复。