[논문 리뷰] Rigid continuation paths I. Quasilinear average complexity for solving polynomial systems

Pierre Lairez|arXiv (Cornell University)|2017. 11. 09.

Polynomial and algebraic computation참고 문헌 46인용 수 10

한 줄 요약

이 논문은 랜덤 가우시안 다항식 시스템을 quasi-linear 평균 복잡도로 해결하기 위해 강성 연속 경로를 사용하는 새로운 수치적 연속 방법을 제안한다. 선형 경로가 아닌 방정식의 강성 운동을 따라 해를 추적함으로써, 평균적으로 $ O(n^4D^2) $개의 연속 단계를 달성하며, 이는 이전의 $ O((\text{입력 크기})^{3/2 + o(1)}) $ bound를 $ O((\text{입력 크기})^{1 + o(1)}) $로 향상시킨다. 이는 스말레의 17번 문제를 최적의 효율성으로 해결한다.

ABSTRACT

How many operations do we need on the average to compute an approximate root of a random Gaussian polynomial system? Beyond Smale's 17th problem that asked whether a polynomial bound is possible, we prove a quasi-optimal bound $ ext{(input size)}^{1+o(1)}$. This improves upon the previously known $ ext{(input size)}^{\frac32 +o(1)}$ bound. The new algorithm relies on numerical continuation along \emph{rigid continuation paths}. The central idea is to consider rigid motions of the equations rather than line segments in the linear space of all polynomial systems. This leads to a better average condition number and allows for bigger steps. We show that on the average, we can compute one approximate root of a random Gaussian polynomial system of~$n$ equations of degree at most $D$ in $n+1$ homogeneous variables with $O(n^5 D^2)$ continuation steps. This is a decisive improvement over previous bounds that prove no better than $\sqrt{2}^{\min(n, D)}$ continuation steps on the average.

연구 동기 및 목표

• 랜덤 다항식 시스템을 풀기 위한 준선형 평균 복잡도를 달성함으로써 스말레의 17번 문제를 해결하기 위해.
• 선형 호모토피 경로의 한계를 극복하기 위해, 조건 수와 단계 크기를 향상시키기 위해 방정식의 강성 운동을 도입하기 위해.
• 근사 해를 계산하기 위해 필요한 연속 단계 수를 $ O((\text{입력 크기})^{3/2 + o(1)}) $에서 $ O((\text{입력 크기})^{1 + o(1)}) $로 줄이기 위해.
• 랜덤 밀도 다항식 시스템의 한 근을 찾는 데 있어 최적의 평균 복잡도를 가지는 결정적 알고리즘을 제공하기 위해.

제안 방법

• 선형 경로가 아닌 다항식 시스템의 강성 운동(회전 및 이동)을 고려함으로써 강성 연속 경로를 도입하기 위해.
• 시스템이 강성 변환을 통해 진화하는 새로운 경로 추적 전략을 정의함으로써, 더 나은 평균 조건 수 $ \mu(F_t, \zeta_t) $를 도출하기 위해.
• 연속 단계의 $ \mu $-추정을 사용: $ K(F,G,\zeta) \leq C \int_0^1 \frac{\mu(F_t, \zeta_t)}{\| \dot{F}_t \| + \| \dot{\zeta}_t \|} dt $, 그러나 개선된 경로 기하학을 적용하기 위해.
• 가우시안 선형 사상과 커널 샘플링을 사용하여 $ H $에 속하는 랜덤 가우시안 시스템 $ G $와 알려진 영점 $ \zeta \in \mathbb{P}(\mathbb{C}^{n+1}) $를 샘플링하기 위해.
• 프로젝티브 뉴턴 반복을 적용하여 강성 경로를 따라 영점을 추적함으로써, 개선된 조건 수 덕분에 더 큰 단계 크기로 수렴 보장하기 위해.
• 자연적 확률 측도에 따라 시스템 공간의 단위 구 위에서 적분함으로써 평균 복잡도를 분석하기 위해.

실험 결과

연구 질문

• RQ1 랜덤 다항식 시스템을 풀기 위한 평균 연속 단계 수를 입력 크기에 대해 준선형으로 줄일 수 있는가?
• RQ2 선형 보간 대신 시스템의 강성 운동을 사용하면 평균 조건 수가 향상되고 더 큰 단계 크기가 가능해지는가?
• RQ3 강성 경로를 따라 $ \mu $-추정을 효과적으로 적용하여 단계 수에 대한 근접 최적의 경계를 확보할 수 있는가?
• RQ4 랜덤 시스템과 알려진 영점을 다항식 시간 내에 샘플링할 수 있는가? 이는 연속의 효율적 초기화를 가능하게 한다.

주요 결과

• 랜덤 가우시안 다항식 시스템의 한 근을 근사적으로 계산하기 위해 필요한 평균 연속 단계 수는 $ O(n^4D^2) $이며, 이는 입력 크기 $ N $에 대해 준선형이다.
• 이것은 이전에 알려진 최선의 평균 복잡도 경계인 $ O(N^{3/2 + o(1)}) $를 초월하여 $ O(N^{1 + o(1)}) $ 복잡도로 향상시켰다.
• 강성 연속 경로의 사용은 선형 경로에 비해 평균 조건 수를 감소시켜 더 크고 안정적인 단계 크기를 가능하게 한다.
• 알고리즘은 최적의 복잡도 경계 $ O(N^{1 + o(1)}) $를 달성하며, 스말레의 17번 문제에 결정적이고 효율적인 방법으로 답을 제시한다.
• 이 방법은 랜덤 시스템과 알려진 영점을 효율적으로 샘플링할 수 있는 새로운 샘플링 기법에 의존하며, 이는 연속 프레임워크와 호환된다.
• 이론적 분석은 기대 단계 수가 $ O(n^4D^2) $로 스케일링되며, 해의 수와 무관하게 최적의 평균 복잡도임을 확인한다.