[논문 리뷰] Telescopers for Rational and Algebraic Functions via Residues

Shaoshi Chen, Manuel Kauers외 1명|arXiv (Cornell University)|2012. 01. 10.

Polynomial and algebraic computation참고 문헌 17인용 수 9

한 줄 요약

이 논문은 잔류항 분석을 이용한 새로운 알고리즘을 제안하며, 유리함수와 대수함수에 대한 텔레스코퍼를 계산한다. 이 알고리즘은 m개 변수를 가진 유리함수에 대한 텔레스코퍼와 m−1개 변수를 가진 대수함수에 대한 텔레스코퍼 사이의 등가성을 확립한다. 이 방법은 미분 연산자를 통한 잔류항 제거를 활용하여, 세 변수를 가진 유리함수에 대해 더 빠른 알고리즘을 제공하고, 두 변수를 가진 대수함수에 대한 텔레스코퍼 순서에 대한 새로운 상계를 도출한다.

ABSTRACT

We show that the problem of constructing telescopers for functions of m variables is equivalent to the problem of constructing telescopers for algebraic functions of m -1 variables and present a new algorithm to construct telescopers for algebraic functions of two variables. These considerations are based on analyzing the residues of the input. According to experiments, the resulting algorithm for rational functions of three variables is faster than known algorithms, at least in some examples of combinatorial interest. The algorithm for algebraic functions implies a new bound on the order of the telescopers.

연구 동기 및 목표

• m개 변수를 가진 유리함수에 대한 텔레스코퍼 계산과 m−1개 변수를 가진 대수함수에 대한 텔레스코퍼 계산 사이의 등가성을 확립하기.
• 잔류항 기반 제거를 이용한 두 변수 대수함수에 대한 텔레스코퍼 계산을 위한 새로운 알고리즘 개발.
• 두 변수 대수함수에 대한 텔레스코퍼 순서에 대한 새로운 상한 도출.
• 세 변수 유리함수에 대한 창의적 텔레스코핑의 계산 효율성을 향상시키기 위해 대수함수 문제로의 환원을 통한 개선.

제안 방법

• 방법은 입력 함수의 잔류항을 유리함수 또는 대수함수 분야에서 분석하는 데 기반한다.
• 함수의 잔류항이 모두 0이면 그 함수가 도함수임을 보장하는 성질을 이용하며, 이는 잔류항 제거가 텔레스코퍼를 구성하는 핵심이 된다.
• 유리함수의 경우 Rothstein-Trager 결과식을 사용해 잔류항을 계산하고, 이 잔류항을 제거하기 위한 미분 연산자를 구성한다.
• 대수함수의 경우 문제를 리만 곡면으로 올리고, 함수 체 확장에서 잔류항을 분석하는 방식을 취한다.
• 잔류항 제거 이후 허미트 축소를 적용하여 증명 함수를 추출한다.
• 프로토타입 구현은 Koutschan의 HolonomicFunctions.m 패키지를 사용하며, 기존 방법과의 성능 비교를 수행한다.

실험 결과

연구 질문

• RQ1 m개 변수를 가진 유리함수에 대한 텔레스코퍼 계산 문제를 m−1개 변수를 가진 대수함수에 대한 텔레스코퍼 계산 문제로 환원할 수 있는가?
• RQ2 두 변수를 가진 대수함수의 잔류항은 어떤 구조를 가지며, 이를 미분 연산자를 통해 어떻게 제거할 수 있는가?
• RQ3 세 변수 유리함수에 대한 기존 창의적 텔레스코핑 알고리즘과 비교했을 때, 잔류항 기반 접근법의 성능은 어떠한가?
• RQ4 잔류항 기반 방법이 두 변수 대수함수에 대한 텔레스코퍼 순서에 대해 더 낮은 상한을 도출할 수 있는가?

주요 결과

• 세 변수 유리함수에 대한 알고리즘이 특정 조합론적 의미를 지닌 예제에서 기존 방법보다 뛰어난 성능을 보이며, 한 경우에서는 1.04초로 실행 시간을 단축시켰고, 가장 느린 대안 방법은 150시간 이상 소요되었다.
• 두 변수 대수함수에 대한 새로운 알고리즘은 텔레스코퍼 순서에 대한 새로운 상한을 도출한다: deg_Dt(L) ≤ [E:K(x)]·deg_x(B*) + 2·genus(E), 여기서 B*는 분모의 제곱-free 부분이다.
• 예제 2에서 새로운 방법(EA)은 1646.53초에 텔레스코퍼를 계산했고, 표준 방법(CK)은 40043.01초가 소요되어 뚜렷한 성능 향상을 보였다.
• 성능 향상은 Rothstein-Trager 결과식이 작은 성분들로 분해되는지 여부에 따라 달라지며, 만약 그렇다면 하나의 큰 유리함수 문제를 푸는 것보다 여러 개의 작은 대수함수 문제를 푸는 것이 더 효율적이다.
• 알고리즘의 1단계(잔류항 제거)가 대부분의 텔레스코퍼를 포괄할 경우, 2단계에서의 연관된 시스템의 복잡도가 감소하여 이 방법이 가장 효과적으로 작용한다.