[논문 리뷰] The largest 5th pivot may be the root of a 61st degree polynomial

James Chen, Alan Edelman외 1명|arXiv (Cornell University)|2026. 02. 23.

Polynomial and algebraic computation인용 수 0

한 줄 요약

본 논문은 JuMP + Groebner basis + discriminant polynomial 방법을 도입하여 완전 피벗팅으로 가우스 소거의 정확한 최대 성장 인자를 결정하고, n=5의 최대가 61차 다항식의 해에 해당함을 보이며, n=6–8에 대한 결과를 제시한다.

ABSTRACT

This paper introduces a number of new techniques in the study of the famous question from numerical linear algebra: what is the largest possible growth factor when performing Gaussian elimination with complete pivoting? This question is highly complex, due to a complicated set of polynomial inequalities that need to be simultaneously satisfied. This paper introduces the JuMP + Groebner basis + discriminant polynomial approach as well as the use of interval arithmetic computations. Thus, we are introducing a marriage of numerical and exact mathematical computations. In 1988, Day and Peterson performed numerical optimization on $n=5$ with NPSOL and obtained a largest seen value of $4.1325...$. This same best value was reproduced by Gould with LANCELOT in 1991. We ran extensive comparable experiments with the modern software tool JuMP and also saw the same value $4.1325...$. While the combinatorial explosion of possibilities prevents us from knowing whether there may not be a larger maximum, we succeed in obtaining the exact mathematical value: the number $4.1325...$ is exactly the root of a 61st degree polynomial provided in this work, and is a maximum given the equality constraints seen by JuMP. In light of the numerics, we pose the conjecture that this lower bound is indeed the maximum. We also apply this technique to $n = 6$, $7$, and $8$. Furthermore, in 1969, an upper bound of $4\frac{17}{18}\approx 4.94$ was produced for the maximum possible growth for $n = 5$. We slightly lower this upper bound to $4.84$.

연구 동기 및 목표

• 완전 피벗팅이 적용된 가우스 소거에서 최대 성장 인자를 동기 부여하고 결정한다.
• 수치 최적화와 대수 기법을 결합한 정확하고 계산 가능한 프레임워크를 개발한다.
• n=5에서 실제 최대에 대한 엄밀한 상한/하한 분석과 추측을 제공한다.

제안 방법

• 동등 제약이 있는 상태에서 완전 피벗팅 하에서 성장 인자를 수치적으로 최적화하기 위해 JuMP를 사용한다.
• 수치 결과에서 타이트한 제약을 고정하고 남은 변수에 대해 대수적 소거를 사용하여 해를 구한다.
• 비선형 다항식에 의존하는 변수들을 제거하기 위해 Groebner 기저 계산을 적용한다.
• 성장 인자의 일변수 다항식으로 후보 최대값을 찾기 위해 판별식 다항식을 사용한다.
• n=5에서 최대 성장 인자에 해당하는 실수 해가 (4,5)에 위치하는 명시적 61차 다항식을 도출한다.
• 해당 방법을 n=6,7,8로 확장하여 그 크기에 대한 하한을 얻는다.

Figure 1. The 61st degree polynomial with only 12 Mathematica input cells

실험 결과

연구 질문

• RQ1 n=5일 때 완전 피벗팅이 적용된 가우스 소거에서 가장 큰 가능한 성장 인자는 무엇인가?
• RQ2 그 최대값이 정수 계수를 가진 고차 다항식의 근으로 정확히 특성화될 수 있는가?
• RQ3 수치 최적화와 정확한 대수적 방법을 어떻게 결합하여 성장 인자 최대값을 검증할 수 있는가?
• RQ4 이 접근법의 n>5(예: n=6,7,8)에 대한 시사점과 확장은 무엇인가?

주요 결과

• 완전 피벗팅 하에서 5x5 행렬의 최대 성장 인자는 [4,5]에 유일한 실근인 61차 다항식 P5(g)의 근이다.
• 61차 다항식은 JuMP + Groebner basis + discriminant 파이프라인과 인수분해를 통해 명시적으로 얻어지며(판별식은 (4,5)에서 실근이 하나인 g의 다항식을 제공한다).
• 수치적으로 관찰된 값 4.1325…은 정확히 이 근이며, 고정밀 보증이 있다(4.1325170786... 및 4.132517078632472854223346853277).
• Lancelot 및 JuMP 기반 수치 실험은 동일한 값을 재현하여 이를 저지한 하한으로 지지하며,저자들이 최대값이라고 추정하는 바(Conjecture 2.1)에 동의한다.
• 해당 방법은 n=6에서 정확한 하한(값 5), n=7에서 여섯차 다항식의 근, 그리고 n=8에서 값 8을 제공한다.
• 4.84의 n=5에 대한 향상된 상한은 구간 산술과 해석의 조합을 통해 얻어지며 Tornheim의 1968년 상한을 개선한다.

(a)