[논문 리뷰] Efficient Computation of the Permanent of Block Factorizable Matrices

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

Matrix Theory and Algorithms참고 문헌 24인용 수 3

한 줄 요약

이 논문은 블록 인수분해 가능 행렬—1×1 또는 2×2 블록을 가진 L개의 블록 대각 행렬의 곱으로 표현할 수 있는 행렬—의 영인수를 다항 시간 내에 계산하는 알고리즘을 제시한다. 영인수를 양자 회로 전이 진폭으로 매핑하고, 행렬 곱 상태(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-난이도를 초월하는 특정 구조를 가진 행렬의 영인수를 다항 시간 내에 정확하게 계산할 수 있는 효율적 정확 알고리즘을 개발하기 위해.
• 작은 블록을 가진 블록 대각 성분으로의 행렬 인수화를 통해 영인수의 계산 비가역성 문제를 해결하기 위해.
• 영인수와 양자 회로 진폭 사이의 연결 고리를 수립함으로써, 양자 다체 물리학의 텐서 네트워크 방법을 활용할 수 있도록 하기 위해.
• 띠형 구조와 띠형 역행렬을 가진 행렬이 그러한 인수분해를 허용함을 보여주어, 영인수의 효율적 계산이 가능함을 입증하기 위해.

제안 방법

• 영인수는 각 블록 대각 인수 요소가 단일 또는 이중 큐비트 게이트를 갖는 회로 레이어에 대응하는 양자 회로의 전이 진폭으로 표현된다.
• 행렬 곱 상태(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에 대해 선형 스케일링을 보인다.
• 이 방법은 영인수의 전이 진폭을 표현하는 양자 회로의 진화 A = F₁F₂⋯Fₗ에 대해 MPS를 효율적으로 적용한다.
• 최대 결합 차원 D_max는 O(2^{L²}) 이하로 증가하며, 이는 전체 실행 시간과 공간 복잡도를 결정한다.
• 이러한 인수분해를 허용하는 모든 행렬 A는 최대 L의 폭을 가진 띠형이어야 하며, 역행렬이 존재할 경우 동일한 폭을 가진 띠형이어야 한다.
• 이 알고리즘은 띠형 구조와 띠형 역행렬을 가진 행렬에 적용 가능하며, 삼중대각 행렬과 띠형 토플리츠 행렬을 포함하여, O(L)개의 블록 대각 인수로 분해 가능한 행렬에 대해 적용 가능하다.
• 이 복잡도 추정치는 최악의 경우를 기반으로 하며, 실질적으로 MPS 절단에서 많은 특이값이 무시 가능할 경우 성능 향상이 가능하다.