[논문 리뷰] Partial fraction decompositions and an algorithm for computing the vector partition function

Todor Milev|arXiv (Cornell University)|2009. 10. 26.

Advanced Combinatorial Mathematics참고 문헌 12인용 수 3

한 줄 요약

이 논문은 포인트드 다면체 콘에 대해 벡터 분할 함수를 쿼asi다항식으로 계산하기 위한 두 가지 알고리즘을 제시한다. 이는 라우렌트 다항식의 링에서 부분 분수 분해를 사용한다. 논문은 예외적인 리 대수계에 대해 코스타ント 분할 함수의 주기 상한을 명시적으로 제공하며, 각각 E6: 6, E7: 12, E8: 60, F4: 12, G2: 6이다. 더불어 더 빠른 알고리즘은 공개된 C++ 프로그램으로 구현되어 있으며, LGPLv3.0 라이선스 하에 vectorpartition.sourceforge.net 에서 이용 가능하다.

ABSTRACT

This paper gives an exposition of well known results on vector partition functions. The exposition is based on works of M. Brion, A. Szenes and M. Vergne and is geared toward explicit computer realizations. In particular, the paper presents two algorithms for computing the vector partition function with respect to a finite set of vectors $I$ as a quasipolynomial over a finite set of pointed polyhedral cones. We use the developed techniques to relate a result of P. Tumarkin and A. Felikson (and present an independent proof in the particular case of finite-dimensional root systems) to give bounds for the periods of the Kostant partition functions of $E_6$, $E_7$, $E_8$, $F_4$, $G_2$ (the periods are divisors of respectively 6, 12, 60, 12, 6). The first of the described algorithms has been realized and is publicly available under the Library General Public License v3.0 at {http://vectorpartition.sourceforge.net/}. We include (non-unique) partial fraction decompositions for the generating functions of the Kostant partition function for $A_2$, $A_3$, $A_4$, $B_2$, $B_3$, $C_2$, $C_3$, $G_2$ in the appendix.

연구 동기 및 목표

• 다각형 콘 위에서 벡터 분할 함수를 쿼اسي다항식으로 계산하기 위한 명시적이고 계산 가능성이 있는 알고리즘을 개발하는 것.
• 입력 벡터의 질량 n-1 부분집합이 생성하는 초평면에 평행한 벽을 가진 쿼اسي다항식의 성질에 대한 구조적 증명을 제공하는 것.
• 부분 분수 분해 기법을 사용하여 예외적인 리 대수에 대한 코스타ント 분할 함수의 주기 상한을 유도하는 것.
• 벡터 분할 함수와 그 침대 분할을 계산하기 위한 고성능 C++ 프로그램을 구현하고 공개하는 것.
• 고전적 루트 체계(A2–A4, B2–B3, C2–C3, G2)에서 코스타ント 분할 함수의 생성함수에 대한 비유일한 부분 분수 분해를 제시하는 것.

제안 방법

• 복소수 단위근을 피하기 위해 계산 효율성을 높이기 위해, 분모가 (1−x^α)^l 형태인 ℚ[x₁,…,xₙ, x₁⁻¹,…,xₙ⁻¹]의 링에서 부분 분수 분해를 사용한다.
• 초기 부분 분수 분해를 생성하기 위해 연장된 쎄네스-버그네 공식(보조정리 2.1(c))을 적용한다.
• 단계적 단순화를 위해 반복적으로 적용하는 감소 알고리즘(섹션 6.1.1)을 사용하며, 항등식 1/(1−x^α) = (1+x^α)/(1−x^{2α})를 이용해 짧은 루트 분모를 제거한다.
• 두 단계 알고리즘을 활용한다: 첫 번째로 유리 함수를 단순화하기 위한 재귀적 감소 과정을 수행하고, 두 번째로 중복된 짧은 루트를 제거하기 위한 변환을 수행한다.
• 다항식, 유리 산술, 해시 테이블, 다면체 콘 분할을 위한 템플릿 클래스를 사용한 맞춤형 C++ 구현을 사용한다.
• 정확도 검증을 위해 체크섬을 사용한다: 각 변환 전후로 유리수 값을 대입하여 프로그래밍 오류를 탐지한다.

실험 결과

연구 질문

• RQ1 벡터 분할 함수 P_I(γ)는 유한한 수의 닫힌, 포인트드 다면체 콘 위에서 어떤 구조를 가지는가?
• RQ2 복소수 단위근을 피하면서, 벡터 분할 함수의 생성함수에 대해 부분 분수 분해를 체계적으로 계산할 수 있는 방법은 무엇인가?
• RQ3 예외적인 루트 체계 E6, E7, E8, F4, G2에 대해 코스타ント 분할 함수의 정확한 주기는 무엇인가?
• RQ4 벡터 분할 함수와 그 침대 분할을 계산하기 위한 효율적이고 공개 가능한 알고리즘이 구현될 수 있는가?
• RQ5 제안된 두 알고리즘의 성능는 어떻게 비교되며, 고전적 루트 체계에 대해 어떤 수정이 효율성을 향상시키는가?

주요 결과

• 벡터 분할 함수 P_I(γ)는 I의 질량 n-1 부분집합이 생성하는 초평면에 평행한 벽을 가진 유한한 수의 닫힘, 포인트드 다면체 콘 위에서 쿼اسي다항식이다.
• 섹션 2와 3의 알고리즘(부분 분수 분해 기반)은 섹션 4의 기본 알고리즘보다 실질적으로 훨씬 빠르다.
• E6, E7, E8, F4, G2에 대한 코스타ント 분할 함수의 주기는 각각 6, 12, 60, 12, 6 이하이며, 이 값들의 약수이다.
• 섹션 6.1.1의 수정된 알고리즘은 고전적 루트 체계에 특화되어 있어, 섹션 2.3의 일반 알고리즘보다 최적화된 감소 규칙 덕분에 성능이 뛰어나다.
• 벡터 분할 함수 프로그램의 구현은 http://vectorpartition.sourceforge.net/ 에서 LGPLv3.0 라이선스 하에 공개되어 있다.
• A2, A3, A4, B2, B3, C2, C3, G2에 대해 코스타ント 분할 함수의 생성함수에 대한 비유일한 부분 분수 분해가 명시적으로 계산되어 부록에 제시되어 있다.