[論文レビュー] 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 の次数と高さに対するきめ細かな上限を確立し、厳密な複雑度解析を可能にすること。
• ℤ 上のヘリトゥルム正規形式（HNF）を計算する高速アルゴリズムを、全体のアルゴリズム内での計算プリミティブとして活用すること。

提案手法

• グレブナー基底計算のための F4 スタイル戦略を採用し、すべての xgᵢ ではなく xg₁,…,xg_{t−1} のみを延長することで、係数行列のサイズを制御する新しい延長法を導入する。
• GHNF に含まれる多項式の次数および高さに対するきめ細かな上限を確立する：次数 ≤ nd および高さ ≤ 6n³d²(h+1+log(n²d)) で、ここで d および h は入力多項式の最大次数および最大高さを表す。
• G = FU を満たす変換行列 U の次数が n, d, h に関して多項式的に有界であることを証明する。これは複雑度解析において極めて重要である。
• 既存の ℤ 行列 HNF 計算の高速アルゴリズムを、延長されたベクトルの係数行列を処理するためのコア計算エンジンとして使用する。
• Magma および Maple に実装し、それらの組み込み HNF コマンドを用いて、標準的なグレブナー基底ルーチンとの実用的性能を評価する。
• ℤ 行列 HNF 計算の最良のビットサイズ複雑度境界に基づく複雑度解析を適用し、最悪ケースの複雑度が O(n^{30.753} d^{17.377} (h + log(n²d))^{4.001}) であることを導出する。

実験結果

リサーチクエスチョン

• RQ1 主イデアル整域でない ℤ[x] 上の行列の一般化ヘリトゥルム正規形式（GHNF）を計算する多項式時間アルゴリズムを設計できるか？
• RQ2 F4 スタイルのグレブナー基底計算を ℤ[x]-加群に適応させることで、延長過程における係数行列のサイズの増大をどのように制御できるか？
• RQ3 入力多項式の上限 d および h が与えられたとき、ℤ[x] 上の行列の 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 要素の次数および高さの上限は、従来のものよりもきめ細かく、アッシュェンブレナーの研究手法を応用したが、ℤ[x]-加群グレブナー基底の特殊構造に適合させたものである。