[论文解读] A Polynomial-time Algorithm to Compute Generalized Hermite Normal Form of Matrices over Z[x]

Rui-Juan Jing, Chun-Ming Yuan+1 more|arXiv (Cornell University)|Jan 6, 2016

Matrix Theory and Algorithms参考文献 19被引用 8

一句话总结

本论文提出了首个在 ℤ[x] 上计算广义赫尔米特正规型（GHNF）的多项式时间算法，其等价于计算 ℤ[x]-模的约化格罗布纳基。该方法结合了优化的 F4 风格延拓策略、GHNF 元素的紧致度与高度界，以及高效的整数矩阵 HNF 计算，实现了最坏情况下的位复杂度 O(n^{30.753} d^{17.377} (h + log(n²d))^{4.001})。

ABSTRACT

In this paper, a polynomial-time algorithm is given to compute the generalized Hermite normal form for a matrix F over Z[x], or equivalently, the reduced Groebner basis of the Z[x]-module generated by the column vectors of F. The algorithm is also shown to be practically more efficient than existing algorithms. The algorithm is based on three key ingredients. First, an F4 style algorithm to compute the Groebner basis is adopted, where a novel prolongation is designed such that the coefficient matrices under consideration have polynomial sizes. Second, fast algorithms to compute Hermite normal forms of matrices over Z are used. Third, the complexity of the algorithm are guaranteed by a nice estimation for the degree and height bounds of the polynomials in the generalized Hermite normal form.

研究动机与目标

• 开发一种用于计算 ℤ[x] 上矩阵广义赫尔米特正规型（GHNF）的多项式时间算法，其等价于计算 ℤ[x]-模的约化格罗布纳基。
• 解决 ℤ[x] 上 GHNF 计算缺乏高效且复杂度可分析算法的问题，该环并非主理想整环，因此不支持标准赫尔米特正规型。
• 在理论上的多项式时间位复杂度与实际效率之间取得平衡，超越 Magma 和 Maple 等系统中现有的格罗布纳基实现。
• 建立 GHNF 中多项式及其变换矩阵 U 的紧致度与高度界，以支持严格的复杂度分析。
• 利用计算 ℤ 上赫尔米特正规型的快速算法作为整体算法中的计算原 primitive。

提出的方法

• 采用 F4 风格的格罗布纳基计算策略，提出一种新颖的延拓方法，仅对 xg₁,…,xg_{t−1} 进行延拓，而非所有 xgᵢ，从而控制系数矩阵的规模。
• 建立 GHNF 中多项式度与高度的紧致上界：度 ≤ nd，高度 ≤ 6n³d²(h+1+log(n²d))，其中 d 和 h 分别为输入多项式的最大度与高度。
• 证明满足 G = FU 的变换矩阵 U 的度在 n、d 和 h 上具有多项式有界性，这对复杂度分析至关重要。
• 使用现有的快速算法计算 ℤ 上的赫尔米特正规型，作为处理延拓向量系数矩阵的核心计算引擎。
• 在 Magma 和 Maple 中实现该算法，利用其内置的 HNF 命令，与标准格罗布纳基例程对比实际性能。
• 基于 ℤ-矩阵 HNF 计算的最佳已知位大小复杂度界进行复杂度分析，得出最坏情况下的复杂度为 O(n^{30.753} d^{17.377} (h + log(n²d))^{4.001})。

实验结果

研究问题

• RQ1 能否设计一种多项式时间算法，用于计算 ℤ[x] 上矩阵的广义赫尔米特正规型（GHNF），而 ℤ[x] 是一个非主理想整环？
• RQ2 如何将 F4 风格的格罗布纳基计算适配于 ℤ[x]-模，以在延拓过程中控制系数矩阵的增长？
• RQ3 在给定输入多项式度与高度界 d 和 h 的前提下，GHNF 中元素的最紧致可能的度与高度界是什么？
• RQ4 能否通过将高效的 ℤ-矩阵 HNF 算法作为子程序，实现一种实际可行且可证明为多项式时间的 ℤ[x]-模格罗布纳基计算算法？
• RQ5 与 Magma 和 Maple 等计算机代数系统中的现有格罗布纳基实现相比，该算法在实际应用中的表现如何？

主要发现

• 所提出的算法实现了最坏情况下的位大小计算复杂度 O(n^{30.753} d^{17.377} (h + log(n²d))^{4.001})，首次建立了 ℤ[x] 上 GHNF 计算的多项式时间算法。
• 该算法显著优于 Magma 中的格罗布纳基命令，尤其在高次多项式情况下：对于度为 150 的情况，Magma 的命令超时（>1000 秒），而新算法在约 105 秒内完成。
• 对于度为 200 的情况，新算法在 Magma 中耗时约 385 秒，而 Magma 的格罗布纳基命令在合理时间内未能终止。
• 在 Maple 中实现的该算法在大度数时（如 d=60）已快于 Magma，表明其在输入规模增大时具有优异的性能扩展性。
• 该算法的高效性归因于新颖的延拓策略，相比标准 F4 风格方法，显著减小了系数矩阵的规模。
• GHNF 元素的度与高度界比以往已知结果更紧致，其推导基于 Aschenbrenner 的工作，但针对 ℤ[x]-模格罗布纳基的特殊结构进行了适配。