[論文レビュー] Efficient Computation of the Permanent of Block Factorizable Matrices

Kristan Temme, Paweł Wocjan|arXiv (Cornell University)|Aug 31, 2012

Matrix Theory and Algorithms参考文献 24被引用数 3

ひとこと要約

本稿では、L 個の 1×1 または 2×2 ブロックを持つブロック対角行列の積として表現可能な行列のパーマネントを計算する多項式時間アルゴリズムを提示する。パーマネントを量子回路遷移振幅に写像し、行列積状態（MPS）を活用することで、時間計算量 O(N·2^{3L²}) および空間計算量 O(N·2^{2L²}) を達成し、帯行列構造および帯行列構造の逆行列をもつ行列において、効率的な正確な計算が可能になる。

ABSTRACT

We present an efficient algorithm for computing the permanent for matrices of size N that can written as a product of L block diagonal matrices with blocks of size at most 2. For fixed L, the time and space resources scale linearly in N, with a prefactor that scales exponentially in L. This class of matrices contains banded matrices with banded inverse. We show that such a factorization into a product of block diagonal matrices gives rise to a circuit acting on a Hilbert space with a tensor product structure and that the permanent is equal to the transition amplitude of this circuit and a product basis state. In this correspondence, a block diagonal matrix gives rise to one layer of the circuit, where each block to a gate acting either on a single tensor component or on two adjacent tensor components. This observation allows us to adopt matrix product states, a computational method from condensed matter physics and quantum information theory used to simulate quantum systems, to evaluate the transition amplitude.

研究の動機と目的

• 一般ケースが #P 困難であることを踏まえ、その拡張としての構造的行列クラスに対する効率的かつ正確なパーマネント計算アルゴリズムの開発。
• 小規模なブロック対角成分への行列因数分解を活用して、パーマネントの計算困難性を克服すること。
• パーマネントと量子回路振幅の間の関係を確立し、量子多体物理学におけるテンソルネットワーク手法の適用を可能にすること。
• 帯行列構造および帯行列構造の逆行列をもつ行列が、このような因数分解を許容することを示し、効率的なパーマネント計算を可能にすること。

提案手法

• パーマネントは、各ブロック対角因数が 1 もしくは 2 ビットゲートをもつ回路レイヤーに対応する量子回路の遷移振幅として表現される。
• 行列積状態（MPS）を用いて、各回路レイヤー後の量子状態を表現・効率的に更新し、系の時間発展を追跡する。
• MPS トランスフォーメーションの順次更新が行われ、各レイヤーで結合次元が増加し、最大結合次元が O(2^{L²}) で抑えられる。
• 最終的なパーマネントは、変換された MPS を計算基底状態 |1⟩^⊗N とのスカラー積として計算され、MPS トランスフォーメーションの縮約によって評価される。
• 各因数内でのブロックの可換性を活用して、ゲートの並列適用を可能とし、回路の深さを低減し、効率を向上させる。
• 実行時間と空間の境界は、結合次元の増加とテンソル縮約のコストに基づき導出され、O(N·2^{3L²}) 時間および O(N·2^{2L²}) 空間が得られる。

実験結果

リサーチクエスチョン

• RQ1 特殊構造をもつ行列のパーマネントは、#P 困難であるにもかかわらず、多項式時間で正確に計算可能か？
• RQ2 小規模なブロック対角行列への因数分解が、量子インスパイアドのテンソルネットワーク手法を用いた効率的パーマネント計算を可能にするか？
• RQ3 行列積状態（MPS）を用いて、行列のパーマネントに対応する量子回路振幅を評価できるか？
• RQ4 サイズが 2 以下のブロックをもつブロック対角行列の積である行列のパーマネントを計算する計算量はどの程度か？
• RQ5 MPS 表現における結合次元の増加は、因数分解レイヤー数 L に対してどのようにスケーリングするか？

主な発見

• アルゴリズムは、ブロック因数分解可能な行列のパーマネントを時間 O(N·2^{3L²}) および空間 O(N·2^{2L²}) で計算可能であり、行列サイズ N に対して線形スケーリングを示す。
• 本手法は、行列積状態（MPS）を用いて、A = F₁F₂⋯Fₗ に対応する量子回路の遷移振幅を効率的に評価する。
• 結合次元 D_max は、最大で O(2^{L²}) まで増加し、これが全体の実行時間および空間計算量を決定する。
• このような因数分解を許容する任意の行列 A は、帯行列であり、帯幅は最大 L である。逆行列が存在する場合、その逆行列も同じ帯幅を持つ帯行列である。
• 本アルゴリズムは、帯行列構造および帯行列構造の逆行列をもつ行列に適用可能であり、三重帯行列や帯行列構造のトーペリッツ行列など、O(L) 個のブロック対角因数に分解可能な行列に適用可能である。
• 境界は最悪ケースの見積もりであるが、MPS 截断における多くの特異値が無視可能である場合には、実際の性能はさらに向上する可能性がある。