[논문 리뷰] On the Complexity of Real Root Isolation

Michael Sagraloff|arXiv (Cornell University)|2010. 11. 01.

Polynomial and algebraic computation참고 문헌 32인용 수 22

한 줄 요약

이 논문은 실수 계수를 가진 정수 다항식의 실근을 고립시키기 위한 새로운 결정론적 알고리즘을 제시한다. 이 알고리즘은 근사적인 데카르트 방법을 적용하며, 중간 다항식은 유한 정밀도로 계산된다. 다항식의 계수 크기나 유형과는 무관하게 실근 고립 복잡도는 근의 기하학적 성질에만 의존함을 보여, 정수 다항식에 대해 $ frac{O}(n^3 au^2)$의 비트 복잡도를 달성함으로써 기존 실용적 방법들보다 $n$ 배 향상된다.

ABSTRACT

We introduce a new approach to isolate the real roots of a square-free polynomial $F=\sum_{i=0}^n A_i x^i$ with real coefficients. It is assumed that each coefficient of $F$ can be approximated to any specified error bound. The presented method is exact, complete and deterministic. Due to its similarities to the Descartes method, we also consider it practical and easy to implement. Compared to previous approaches, our new method achieves a significantly better bit complexity. It is further shown that the hardness of isolating the real roots of $F$ is exclusively determined by the geometry of the roots and not by the complexity or the size of the coefficients. For the special case where $F$ has integer coefficients of maximal bitsize $τ$, our bound on the bit complexity writes as $ ilde{O}(n^3τ^2)$ which improves the best bounds known for existing practical algorithms by a factor of $n=deg F$. The crucial idea underlying the new approach is to run an approximate version of the Descartes method, where, in each subdivision step, we only consider approximations of the intermediate results to a certain precision. We give an upper bound on the maximal precision that is needed for isolating the roots of $F$. For integer polynomials, this bound is by a factor $n$ lower than that of the precision needed when using exact arithmetic explaining the improved bound on the bit complexity.

연구 동기 및 목표

• 실수 계수를 가진 정수 다항식의 실근을 고립시키기 위한 실용적이고 정확하며 결정론적인 알고리즘을 개발하는 것.
• 기존 실용적 방법들(예: 데카르트 또는 연분수 방법)보다 실근 고립의 비트 복잡도를 낮추는 것.
• 실근 고립의 복잡도가 계수 유형이나 크기와는 무관하게 근의 기하학적 성질—특히 $\Gamma_F$와 $\Sigma_F$—에 의해 결정됨을 보여주는 것.
• 정확한 산술 대신 근사 산술을 사용할 경우 worst-case 시나리오에서 증명 가능한 복잡도 향상이 이루어질 수 있음을 보여주는 것.
• 복잡한 근 고립을 위한 점근적 빠른 방법과 실용적 알고리즘 사이의 격차를 메우기 위해 실근에 대해 거의 최적의 복잡도를 달성하는 것.

제안 방법

• 각 재귀 단계에서 다항식 계수의 근사값을 사용하여 수정된 데카르트 방법을 적용한다.
• 각 구간 $I$에서 정밀도 한계 $\rho_I$를 유지하며, 근 고립에 필요한 정확도에 따라 동적으로 조정한다.
• 근이 고립되었는지 확인하기 위해 $\rho_I$-근사값 $g$가 정확히 한 개의 부호 변화를 가지는지 검증하는 인증 단계를 사용한다.
• 근사값 $\tilde{f}_I$를 필요한 정밀도까지만 계산함으로써 계산 비용을 줄이는 재귀적 분할 전략을 채택한다.
• 작은 구간에서는 변환된 다항식 $f_I(x) = f(a + w(I)x)$의 주계수들이 작아지므로, 근사치에서 차수를 감소시킬 수 있음을 활용한다.
• 계산적으로 효율적이고 실무적으로 데카르트의 정부호 법칙과 유사한 $\mathcal{T}'_{3/2}$-테스트 기반의 수정된 포함성 판별식을 사용한다.

실험 결과

연구 질문

• RQ1 데카르트 방법에서 근사 산술을 사용할 경우 실근 고립의 worst-case 비트 복잡도가 증명 가능하게 감소할 수 있는가?
• RQ2 실근 고립의 복잡도가 계수 크기보다 근의 기하학적 성질에 얼마나 의존하는가?
• RQ3 정확한 산술을 피하면서도 정확성과 완전성을 유지할 수 있는 복잡도 향상이 가능할까?
• RQ4 근 간격과 로그적 근 범위에 따라 근사값에 필요한 정밀도가 어떻게 변화하는가?
• RQ5 분할 과정 중에 근사 다항식의 차수를 줄일 수 있으며, 이는 성능 향상에 기여하는가?

주요 결과

• 정수 계수 다항식에 대해 이 알고리즘은 worst-case 비트 복잡도 $ frac{O}(n^3\tau^2)$를 달성하며, 기존 실용적 방법의 최고 알려진 복잡도보다 $n$ 배 향상된다.
• 실근 고립의 복잡도는 계수 유형이나 크기와는 무관하게 근의 기하학적 성질—즉, $\Gamma_F$와 $\Sigma_F$—에 의해 결정된다.
• 정확한 산술 대비 근사값에 필요한 최대 정밀도가 $n$ 배 감소하여, 복잡도 향상 결과가 설명된다.
• 근이 서로 가까운 다항식, 예를 들어 미그노트 다항식의 경우, 분할 과정에서 근사 다항식 $\tilde{f}_I$의 차수가 감소하여 계산 비용이 극적으로 감소한다.
• 고립된 구간 사이의 영역만 분할함으로써 인증 단계의 실행 시간을 거의 무시할 수 있을 정도로 줄일 수 있다.
• 알고리즘은 고전적인 데카르트 방법과의 구조적 유사성 덕분에 실용적이고 쉽게 구현할 수 있으며, 프로토타입 구현 결과 뛰어난 성능을 보였다.