[논문 리뷰] From Approximate Factorization to Root Isolation with Application to Cylindrical Algebraic Decomposition

Kurt Mehlhorn, Michael Sagraloff외 1명|arXiv (Cornell University)|2013. 01. 21.

Polynomial and algebraic computation참고 문헌 36인용 수 7

한 줄 요약

이 논문은 다중근을 가진 단변수 복소다항식에 대한 근 고립 알고리즘을 제안한다. 이때 근의 수가 알려져 있고 계수는 임의의 정밀도로 근사 가능하다는 가정을 한다. 근사 인수분해(특히 Pan의 알고리즘)를 활용하여, 반경이 2^{-\kappa} 이하인 원판 내에서 근을 고립하는 데 비트 복잡도 \tilde{O}(d^3 + d^2\tau + d\kappa)를 달성한다. 이는 실평면 대수곡선의 위상 계산 및 이변수 시스템 해법과 같은 원통형 대수적 분할 작업의 랜덤화된 복잡도 한계를 크게 향상시킨다.

ABSTRACT

We present an algorithm for isolating the roots of an arbitrary complex polynomial $p$ that also works for polynomials with multiple roots provided that the number $k$ of distinct roots is given as part of the input. It outputs $k$ pairwise disjoint disks each containing one of the distinct roots of $p$, and its multiplicity. The algorithm uses approximate factorization as a subroutine. In addition, we apply the new root isolation algorithm to a recent algorithm for computing the topology of a real planar algebraic curve specified as the zero set of a bivariate integer polynomial and for isolating the real solutions of a bivariate polynomial system. For input polynomials of degree $n$ and bitsize $τ$, we improve the currently best running time from $ O(n^{9}τ+n^{8}τ^{2})$ (deterministic) to $ O(n^{6}+n^{5}τ)$ (randomized) for topology computation and from $ O(n^{8}+n^{7}τ)$ (deterministic) to $ O(n^{6}+n^{5}τ)$ (randomized) for solving bivariate systems.

연구 동기 및 목표

• 다중근을 포함한 단변수 복소다항식의 근 고립 알고리즘을 개발하여, 서로 다른 근의 수가 알려져 있을 경우에 대응한다.
• 다항식 계수의 임의 정밀도 근사치를 제공하는 계수 오라클의 사용을 가능하게 한다.
• 계산대수기하학의 근본 문제들, 특히 원통형 대수적 분할을 위한 근 고립 및 관련 문제의 비트 복잡도를 감소시킨다.
• 실평면 대수곡선의 위상 계산 및 이변수 다항식 시스템 해법과 같은 최신 복잡도 한계를 향상시킨다.

제안 방법

• 알고리즘은 계수를 오라클이 제공하는 유리수 근사치를 통해 얻는 근사 인수분해를 서브루틴으로 사용한다.
• 서로 겹치지 않는 원판을 계산하여 각 원판이 정확히 하나의 서로 다른 근과 그 정확한 중복도를 포함하도록 근을 고립한다.
• 각 고립 원판의 반경은 \sigma_i / (64n) 이하로 제한되며, 여기서 \sigma_i 는 근 i 와 다른 근들 사이의 분리도이다.
• 근의 기하적 위치에 따라 정밀도 요구 조건을 적응적으로 조정하며, 핵심 섹션을 통해 분할 효과를 극대화한다.
• 결과 다항식의 정밀 분석을 적용하고, 근의 분리도 및 계수 크기에 따라 복잡도 한계를 설정한다.
• Pan의 알고리즘을 이용한 근사 인수분해와 결합함으로써, 근 고립 및 후속 응용 분야에서 근사 최적 비트 복잡도를 달성한다.

실험 결과

연구 질문

• RQ1 다중근을 포함한 다항식의 근 고립 알고리즘이, 서로 다른 근의 수가 알려져 있고 계수 오라클을 통해 임의 정밀도 근사치에 접근 가능한 경우에 효율적으로 작동할 수 있는가?
• RQ2 임의 정밀도 계수 접근 조건 하에서 다중근을 가진 단변수 다항식의 모든 근을 고립하는 데 최적의 비트 복잡도는 무엇인가?
• RQ3 근 고립은 곡선 위상 계산 및 이변수 시스템 해법과 같은 고수준 기하 계산에서 복잡도를 어떻게 감소시킬 수 있는가?
• RQ4 근의 분리도에 기반한 적응형 정밀도 요구 조건이 대수 알고리즘의 전체 복잡도 한계에 얼마나 기여할 수 있는가?
• RQ5 적응형 근 고립 방법을 사용하여 원통형 대수적 분할 작업의 복잡도를 크게 감소시킬 수 있는가?

주요 결과

• 정수 계수 다항식의 비트 크기가 \tau 인 경우, 반경이 2^{-\kappa} 이하인 고립 원판을 위한 근 고립에 \tilde{O}(d^3 + d^2\tau + d\kappa) 비트 연산이 소요된다.
• 차수 d 및 비트 크기 \tau 를 가진 다항식에 대해, 이 방법은 \tilde{O}(d^3 + d^2\tau) 의 복잡도로 근 고립을 달성하며, 이는 현재 알려진 최고의 복잡도 한계와 일치한다.
• 실평면 대수곡선의 위상 계산에 적용했을 경우, 복잡도는 결정론적 기존의 \tilde{O}(n^9\tau + n^8\tau^2) 에서 랜덤화된 \tilde{O}(n^6 + n^5\tau) 로 향상된다.
• 이변수 다항식 시스템 해법에 적용했을 경우, 복잡도는 결정론적 기존의 \tilde{O}(n^8 + n^7\tau) 에서 랜덤화된 \tilde{O}(n^6 + n^5\tau) 로 감소된다.
• 향상은 근의 분리도 및 기하적 구성에 따라 정밀도 요구 조건이 적응적으로 조정되어 강력한 분할 효과를 발휘하기 때문이다.
• 이 방법은 적응형 복잡도 분석이 대수기하 알고리즘에서 최악의 복잡도 한계를 훨씬 초월한 향상을 이끌 수 있음을 보여준다.