[論文レビュー] Recursive Method for the Solution of Systems of Linear Equations

Gennadi Malaschonok|arXiv (Cornell University)|Mar 29, 2017

Algebraic and Geometric Analysis参考文献 2被引用数 3

ひとこと要約

本稿では、行列乗算の計算量と同一の漸近的複雑度 $O(n^{\text{exponent of matrix multiplication}})$ を達成する、可換整域上での一次方程式系を解く再帰的アルゴリズムを提案する。この手法は、再帰的分割と行列式に基づく還元を用い、係数行列を対角形に変換することで、隅小行列が非ゼロかつゼロ因子でない場合に限って、最小限の算術演算でクラメルの公式を用いて効率的に解を得られる。

ABSTRACT

New solution method for the systems of linear equations in commutative integral domains is proposed. Its complexity is the same that the complexity of the matrix multiplication.

研究の動機と目的

• 可換整域上での一次方程式系を解く手法を開発し、行列乗算の複雑度と同一の漸近的複雑度を達成すること。
• 前進/後退代入やワンパス法といった既存手法を統一的に扱う再帰的フレームワークに一般化すること。
• 再帰的分解と行列式算術を活用することで、一次方程式系の解法の計算複雑度を $O(n^3)$ から行列乗算と同等の漸近的順序に低減すること。
• すべての可換環への適用を保証するため、隅小行列が非ゼロかつゼロ因子でないという条件を要求すること。

提案手法

• アルゴリズムは係数行列を再帰的に上部と下部の部分行列に分割し、各部分に再帰的還元を適用する。
• 変換中に構造的不変性を維持するため、小行列 $\delta^{k} = \det \mathbf{A}_{kk}^{k}$ と余因子行列 $\delta^{k}_{ij}$ を用いた再帰的行列式計算を用いる。
• 重要な恒等式 $\delta^{k} \cdot \hat{A}_{2}^{2} = \delta^{s} \cdot A_{2}^{2} - A_{1}^{2} \cdot G_{2}^{1}$ により、下部ブロックの変換が行われ、行列式構造が保存される。
• この手法により、行列は対角形 $\tilde{A} \rightarrow (\delta^{l}I_{l-k}, \hat{G})$ に還元され、ここで $\delta^{l}$ は先頭小行列、$\hat{G}$ は一般化された余因子を含む。
• 再帰的ステップを全行列が還元されるまで繰り返し、最終的な解は計算された小行列と余因子行列を用いたクラメルの公式から導出される。
• 分割が1行の場合には前進/後退代入が、最適な複雑度を達成するための二分分割の場合には、ワンパス法と同様の動作を示す。

実験結果

リサーチクエスチョン

• RQ1 再帰的アルゴリズムは、整域上での一次方程式系の解法において、行列乗算と同等の漸近的複雑度を達成できるか？
• RQ2 ワンパス法と前進/後退代入法を統一する単一の再帰的フレームワークで、一次方程式系の解法をどのように統合できるか？
• RQ3 一次方程式系の解法の複雑度を行列乗算と一致させるために、どのような再帰的分解戦略が必要か？
• RQ4 環と行列構造にどのような条件下で、このアルゴリズムが有効かつ効率的であるか？

主な発見

• アルゴリズムは、ストラッセンのアルゴリズムを用いることで $O(n^{\log_2 7})$ の計算複雑度を達成し、古典的な $O(n^3)$ の複雑度に比べて顕著に向上する。
• $n = 2^p$ かつ $m = n+1$ の場合、乗算/除算演算の数は $MD^{S}_{n,n+1} = \frac{7}{15}n^{\log_2 7} + n^2(\log_2 n - \frac{2}{3}) + n(2\log_2 n + \frac{1}{5})$ で与えられ、近似的に最適なスケーリングを示す。
• この手法はワンパス法と前進/後退代入法を一般化する：分割が1行の場合、これらの既知の技法に簡略化される。
• 解は計算された小行列 $\delta^n$ と余因子行列 $G$ を用いたクラメルの公式により得られ、最終的なシステムは $\delta^n I_n$ と $\hat{G}$ で表される。
• すべての可換整域で有効である。条件はすべての隅小行列 $\delta^k$ が非ゼロかつゼロ因子でないことである。
• 整数環 $\mathbb{Z}$ 上での数値例により正しさが確認された：$4 \times 4$ の方程式系に対して、行列は $\delta^4 = 27$ を持つ対角形に還元され、解ベクトル $x = (1, 2, -2, -1)^T$ が得られた。