[논문 리뷰] Fast (Multi-)Evaluation of Linearly Recurrent Sequences: Improvements and Applications

Martin Ziegler|ArXiv.org|2005. 11. 08.

Numerical Methods and Algorithms참고 문헌 26인용 수 4

한 줄 요약

이 논문은 다항식 다중점 평가와 빠른 행렬 산술을 활용하여 선형 재귀 수열의 개별적 및 병렬 평가를 최적화한 알고리즘을 제시한다. 수직 다항식 평가에 대해 균일한 $Ó(\sqrt{n} \cdot \log n)$ 복잡도를 달성하고, 다중 평가에 대해서는 $Ó(\sqrt{n} \cdot \log n + \ell \cdot \log(n/\ell))$의 복잡도를 확보하며, 초함수 및 고차 다항식 계수 계산에의 적용 가능성을 제시한다.

ABSTRACT

For a linearly recurrent vector sequence P[n+1] = A(n) * P[n], consider the problem of calculating either the n-th term P[n] or L<=n arbitrary terms P[n_1],...,P[n_L], both for the case of constant coefficients A(n)=A and for a matrix A(n) with entries polynomial in n. We improve and extend known algorithms for this problem and present new applications for it. Specifically it turns out that for instance * any family (p_n) of classical orthogonal polynomials admits evaluation at given x within O(n^{1/2} log n) operations INDEPENDENT of the family (p_n) under consideration. * For any L indices n_1,...,n_L <= n, the values p_{n_i}(x) can be calculated simultaneously using O(n^{1/2} log n + L log(n/L)) arithmetic operations; again this running time bound holds uniformly. * Every hypergeometric (or, more generally, holonomic) function admits approximate evaluation up to absolute error e>0 within O((log(1/e)^{1/2} loglog(1/e)) -- as opposed to O(log(1/e)) -- arithmetic steps. * Given m and a polynomial p of degree d over a field of characteristic zero, the coefficient of p^m to term X^n can be computed within O(d^2 M(n^{1/2})) steps where M(n) denotes the cost of multiplying two degree-n polynomials. * The same time bound holds for the joint calculation of any L<=n^{1/2} desired coefficients of p^m to terms X^{n_i}, n_1,...,n_L <= n.

연구 동기 및 목표

• 상수 또는 다항계수를 가진 선형 재귀 수열 평가를 위한 기존 알고리즘을 향상 및 일반화하기 위해.
• 임의의 인덱스에서 수열 항을 근사 최적의 시간 복잡도로 효율적으로 다중 평가할 수 있도록 하기 위해.
• 기호 계산의 기본 문제들, 특히 수직 다항식 평가 및 고차 다항식의 계수 추출에 개선된 알고리즘을 적용하기 위해.
• 홀로노믹 및 초함수를 포함한 다양한 수열 유형에 대해 균일한 복잡도 한계를 설정하기 위해.
• 새로운 알고리즘이 고전적 방법을 상수 인자 오버헤드만으로 일반화하면서도 점근적 최적성은 유지함을 보여주기 위해.

제안 방법

• 선형 재귀 수열 평가의 비용을 줄이기 위해 빠른 다항식 다중점 평가 및 빠른 푸리에 변환을 활용한다.
• 상수 계수 재귀의 고차 인덱스 항을 $Ó(k^3 \cdot \log n)$ 연산 내에서 행렬 거듭제곱 및 반복 제곱 기법을 적용하여 계산한다.
• 동반 행렬의 구조에 기반한 분할 정복 전략을 활용해 동시에 여러 개의 비연속 항을 계산한다.
• 유리 함수 산술과 크래머의 법칙을 적용하여 홀로노믹 수열의 합과 곱에서 계수의 차수 한계를 분석한다.
• 수직 다항식 평가 문제를 다항계수를 가진 선형 재귀로 환원하여 균일한 복잡도 한계를 확보한다.
• 다항식 행렬의 역행렬의 항이 유리 함수 차수에 의해 유계임을 이용해 합과 곱의 수열에 대한 재귀 관계를 효율적으로 계산할 수 있음을 활용한다.

실험 결과

연구 질문

• RQ1 고전적 방법보다 더 빠르게 $\mathcal{O}(n)$ 연산을 초과하지 않고 선형 재귀 수열의 $n$번째 항을 고급 대수 기법을 통해 평가할 수 있는가?
• RQ2 선형 재귀 수열의 $\ell$개 비연속 항 평가에 대해 최적의 복잡도는 무엇이며, 다양한 수열 유형 간에 균일하게 달성 가능한가?
• RQ3 균일한 복잡도 한계를 확보하면서 재귀 기반 방법을 활용해 고전적 수직 다항식의 평가를 어떻게 가속화할 수 있는가?
• RQ4 예를 들어 $p^m$에서의 특정 계수 계산은 출력에 민감한 복잡도로 수행될 수 있는가?
• RQ5 다항계수를 가진 홀로노믹 수열의 합과 곱을 지배하는 재귀 관계의 차수 및 복잡도 한계는 무엇인가?

주요 결과

• 모든 고전적 수직 다항식 가족의 $n$번째 항은 주어진 $x$에서 $\mathcal{O}(\sqrt{n} \cdot \log n)$의 산술 연산으로 평가될 수 있으며, 이는 특정 다항식 가족에 따라 달라지지 않는다.
• 모든 수직 다항식 가족에 대해 $n_i \leq n$인 임의의 $\ell$개 항 $p_{n_i}(x)$의 평가를 동시에 $\mathcal{O}(\sqrt{n} \cdot \log n + \ell \cdot \log(n/\ell))$의 산술 연산으로 수행할 수 있으며, 이는 모든 수열 유형에 대해 균일하다.
• 절대 오차 $\epsilon > 0$ 이내로 초함수(또는 홀로노믹) 함수를 근사 평가하기 위해서는 $\mathcal{O}(\sqrt{\log(1/\epsilon)} \cdot \log\log(1/\epsilon))$의 산술 단계가 필요하며, 이는 고전적 $\mathcal{O}(\log(1/\epsilon))$의 bound를 향상시킨다.
• $p$가 차수 $d$인 다항식일 때, $p^m$의 $X^n$ 계수는 $\mathcal{O}(d^2 \cdot M(\sqrt{n}))$의 단계 내에서 계산될 수 있으며, 여기서 $M(n)$은 다항식 곱셈 비용을 나타낸다.
• $\ell \leq \sqrt{n}$인 경우, $p^m$의 지정된 차수 $n_i \leq n$에서의 계수를 공동으로 계산하는 데에도 $\mathcal{O}(d^2 \cdot M(\sqrt{n}))$의 단계가 소요되며, 이는 단일 계수의 bound와 일치한다.
• 다항계수의 차수 $d$를 가진 두 홀로노믹 수열의 합과 곱은 각각 깊이 최대 $k + \ell$ 및 $k \cdot \ell$인 재귀 관계를 가지며, 유리 계수의 차수는 각각 $\mathcal{O}(d(k + \ell)^2)$ 및 $\mathcal{O}(d k^2 \ell^2)$이다.