🔢 mathematics

The largest 5th pivot may be the root of a 61st degree polynomial

Dit artikel introduceert een hybride numeriek-exacte methode om aan te tonen dat de grootste 5e pivot bij Gauss-eliminatie met volledige pivoting exact de wortel is van een polynoom van de 61e graad, waarmee een nieuwe bovengrens van 4,84 wordt vastgesteld voor de groeifactor.

Oorspronkelijke auteurs: James Chen, Alan Edelman, John Urschel

Gepubliceerd 2026-03-17

📖 5 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: James Chen, Alan Edelman, John Urschel

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

De Grote Grootte-Explosie: Een Reis door de Wiskunde van het "Groeifactor"-Raadsel

Stel je voor dat je een gigantische, complexe puzzel probeert op te lossen. Je hebt een reeks getallen (een matrix) en je moet ze stap voor stap herschikken en vereenvoudigen om een antwoord te vinden. Dit proces heet Gaussische eliminatie, en het is de ruggengraat van bijna alle moderne computers die berekeningen doen, van het voorspellen van het weer tot het ontwerpen van vliegtuigen.

Maar er is een probleem: bij het oplossen van deze puzzels kunnen de getallen soms enorm groot worden. Stel je voor dat je begint met getallen tussen 0 en 1, maar halverwege de berekening springen ze plotseling naar 100 of 1000. Dit heet de groeifactor. Als deze factor te groot wordt, kan de computer in de war raken en fouten maken. Wiskundigen willen dus weten: Wat is het ergste scenario dat kan gebeuren? Hoe groot kunnen die getallen maximaal worden?

Dit papier van James Chen, Alan Edelman en John Urschel gaat over het vinden van dit maximale getal voor een specifieke, lastige versie van de puzzel (voor matrices van 5x5).

1. Het oude raadsel en de nieuwe aanpak

Sinds de jaren 60 weten wiskundigen de maximale groei voor kleine puzzels (3x3 en 4x4). Maar voor de 5x5 puzzel was het een mysterie. Computers hadden al vaak een antwoord gevonden van ongeveer 4,1325. Maar was dit het echte maximum? Of was er ergens een nog groter getal verborgen?

Vroeger probeerden wiskundigen dit met pure numerieke berekeningen (proberen en fouten). Maar dit papier introduceert een nieuwe, krachtige methode: een huwelijk tussen numerieke kracht en exacte wiskunde.

De Numerieke Kracht (JuMP): Denk hieraan als een super-snelle zoekrobot. Deze robot loopt door miljarden mogelijke puzzelconfiguraties en vindt de beste kandidaat. Hij zegt: "Kijk, hier krijg ik 4,1325!"
De Exacte Wiskunde (Groebner-basis & Discriminant): Nu komt de echte magie. De auteurs nemen de resultaten van de robot en zeggen: "Oké, laten we dit niet benaderen, maar bewijzen." Ze gebruiken een wiskundig gereedschap (Groebner-basis) om de complexe vergelijkingen te "ontwarren", alsof je een knoop in een touw oplost tot een rechte lijn.

2. De verrassende ontdekking: De 61e-graads polynoom

Wat ze vonden, was verbazingwekkend. Het getal 4,1325... is niet zomaar een willekeurig decimaal getal. Het is een exact wiskundig getal dat de oplossing is van een vergelijking met een 61e graad.

De Analogie:
Stel je voor dat je een sleutel zoekt die een deur opent. De meeste mensen denken dat de sleutel een simpele, rechte staaf is. Maar deze auteurs ontdekten dat de sleutel een ingewikkeld, 61-stukjes tellend mechanisme is.

Het getal 4,1325... is de enige echte sleutel die past in dat specifieke 61-delige slot.
Ze hebben de vergelijking (de polynoom) die dit slot beschrijft, volledig opgeschreven. Het is een reusachtige formule met 61 termen, maar het is de exacte waarheid.

Ze concluderen: "We vermoeden sterk dat dit het echte maximum is. Hoewel we niet 100% kunnen bewijzen dat er geen nog complexere, onontdekte configuratie bestaat (de 'combinatorische explosie' is te groot), hebben we wel bewezen dat binnen de logische grenzen die de computer ons gaf, dit het absolute maximum is."

3. Het plafond verlagen

Naast het vinden van de exacte ondergrens (4,1325...), hebben ze ook het plafond iets verlaagd.

Vroeger: Wiskundigen wisten zeker dat het getal nooit hoger zou gaan dan ongeveer 4,94.
Nu: Met een slimme combinatie van wiskundige analyse en "intervalrekenen" (een methode waarbij je niet één getal berekent, maar een veiligheidskader van getallen om zekerheid te hebben), hebben ze bewezen dat het maximum 4,84 is.

De Analogie:
Stel je voor dat je weet dat een bal niet hoger dan 5 meter kan springen. Door de grond te analyseren en de windrichting te meten, zeggen deze auteurs: "Eigenlijk kan die bal nooit hoger dan 4,84 meter springen." Het is een kleine verbetering, maar in de wiskunde is elke centimeter telt.

4. Wat betekent dit voor de toekomst?

De auteurs hebben hun methode ook toegepast op grotere puzzels (6x6, 7x7 en 8x8).

Voor 6x6 is het maximum 5.
Voor 7x7 is het weer een getal dat de oplossing is van een 6e-graads vergelijking.
Voor 8x8 is het 8.

Ze concluderen dat de reeks van maximale groeifactoren (1, 2, 2,25, 4, [61e-graads getal], 5, [6e-graads getal], 8) een heel vreemd en mooi wiskundig patroon vormt.

Samenvatting voor de leek

Dit papier is een feest van samenwerking tussen menselijke creativiteit en computerkracht.

De Robot vond de beste kandidaat.
De Wiskundigen namen die kandidaat en bouwden er een onbreekbaar bewijs omheen, waardoor ze ontdekten dat het antwoord een zeer complex, maar exact getal is (de wortel van een 61e-graads vergelijking).
Ze bewezen ook dat het antwoord niet zomaar "ergens rond de 4,9" kan liggen, maar strikt onder de 4,84 moet blijven.

Het is een eerbetoon aan Nick Higham, een legende in dit vakgebied, die altijd gefascineerd was door deze vraag. De auteurs zeggen: "We hadden je dit niet eerder kunnen vertellen, maar nu hebben we het eindelijk opgelost."

Kortom: Ze hebben het raadsel van de 5x5 matrix niet alleen opgelost, maar het antwoord in een prachtige, alhoewel complexe, wiskundige vorm gegoten.

1. Het Probleem

Het artikel richt zich op een klassiek en complex vraagstuk binnen de numerieke lineaire algebra: wat is de maximale groeifactor ( $g(A)$ ) die kan optreden bij het uitvoeren van Gauss-eliminatie met volledige pivoting (complete pivoting) voor een reële $n \times n$ matrix?

De groeifactor is gedefinieerd als:
$g(A) = \frac{\max_{i,j,k} |A^{(k)}_{ij}|}{\max_{i,j} |A_{ij}|}$
waarbij $A^{(k)}$ de submatrix is na $k-1$ stappen van eliminatie. Een kleine groeifactor impliceert numerieke stabiliteit.

Voor kleine dimensies $n$ zijn de waarden bekend:

$n=1$ : 1
$n=2$ : 2
$n=3$ : $2.25 $($ 2 \frac{1}{4}$)
$n=4$ : 4

Voor $n=5$ was de waarde decennia lang een open probleem. Sinds 1988 (Day en Peterson) en bevestigd in 1991 (Gould) en 2024 (met moderne software zoals JuMP), werd numeriek een waarde van ongeveer 4.1325 gevonden. Echter, of dit de exacte maximumwaarde is, en welke wiskundige vorm deze heeft, was onbekend. Bestaande bovengrenzen waren los (bijv. 4.94 uit 1969).

2. Methodologie

De auteurs introduceren een innovatieve hybride aanpak die numerieke optimalisatie combineert met exacte algebraïsche methoden. De kern van hun methode bestaat uit drie stappen:

Numerieke Optimalisatie (JuMP):
Ze gebruiken de Julia-pakketten JuMP (Julia for Mathematical Programming) om numerieke optimalisatie uit te voeren. Dit levert een kandidaat-matrix op met een hoge groeifactor en identificeert welke constraints (gelijkheden) "strak" (tight) zijn. In plaats van alle combinatorische mogelijkheden te doorzoeken, laten ze de numerieke solver de structuur van de gelijkheidsconstraints suggereren.
Gröbner-basis Berekeningen:
Zodra de structuur van de strakke constraints bekend is, worden deze omgezet in een stelsel van niet-lineaire polynoomvergelijkingen. De auteurs gebruiken Gröbner-basis technieken (geïmplementeerd in Mathematica) om variabelen uit dit stelsel te elimineren. Dit reduceert het probleem van het oplossen van een complex stelsel naar het vinden van de wortels van een univariaat polynoom.
Discriminant Polynomen:
Om het exacte maximum te vinden binnen de gevonden constraints, gebruiken ze de discriminant van het resulterende polynoom. De discriminant geeft de waarden van de parameter (de groeifactor $g$ ) waarvoor het polynoom meervoudige wortels heeft, wat correspondeert met de stationaire punten (extrema) van de functie.
Intervalrekening en Branch-and-Bound (voor bovengrenzen):
Om een nieuwe bovengrens te bewijzen, combineren de auteurs wiskundige analyse met intervalrekening. Ze construeren overdekkingen (covers) van de ruimte van $3 \times 3$ matrices met hoge groeifactor en tonen aan dat bepaalde combinaties van deze matrices onverenigbaar zijn binnen een $5 \times 5$ matrix met volledige pivoting.

3. Belangrijkste Bijdragen en Resultaten

A. De Exacte Waarde voor $n=5$

De meest opvallende bevinding is dat de numeriek gevonden waarde $4.1325...$ niet willekeurig is, maar een algebraïsch getal is.

De auteurs bewijzen dat deze waarde de unieke reële wortel is van een polynoom van graad 61 met gehele coëfficiënten (aangeduid als $P_5(g)$ ).
De exacte waarde is:
$g \approx 4.13251707863247285422...$
Ze concluderen dat deze waarde het maximum is binnen de set van matrices die voldoen aan de door JuMP geïdentificeerde gelijkheidsconstraints (Stelling 2.2).
Conjecture 2.1: De auteurs vermoeden dat deze waarde ook het absolute maximum is voor alle $5 \times 5$ matrices met volledige pivoting, hoewel een volledig bewijs voor de volledige combinatorische ruimte (vanwege de explosie van mogelijkheden) nog niet geleverd kan worden.

B. Resultaten voor $n=6, 7, 8$

De methode wordt ook toegepast op hogere dimensies:

$n=6$ : De maximale groeifactor is exact 5.
$n=7$ : De waarde is de wortel van een polynoom van graad 6 ( $g \approx 6.0569...$ ).
$n=8$ : De waarde is 8 (bereikt door Hadamard-matrices, maar ook door niet-Hadamard matrices).

C. Verbeterde Bovengrens voor $n=5$

De auteurs verbeteren de bestaande bovengrens voor $n=5$ .

De oude bovengrens (Tornheim, 1969) was $4 \frac{17}{18} \approx 4.94$ .
Met hun nieuwe techniek (combinatie van analyse en intervalrekening) verlagen ze deze naar 4.84 (Stelling 3.5).
Hoewel dit nog steeds boven de vermoedelijke ondergrens van 4.1325 ligt, is het een significante verkleining van de kloof.

4. Significatie en Toekomstperspectief

Wiskundige Structuur: Het artikel onthult dat de oplossing voor een fundamenteel numeriek probleem een verrassend complexe algebraïsche vorm heeft (een polynoom van graad 61). Dit illustreert de diepe verbinding tussen numerieke linear algebra en algebraïsche meetkunde.
Hybride Berekening: De paper demonstreert succesvol hoe numerieke optimalisatie (die snel maar benaderend is) kan worden gebruikt om de structuur van een probleem te vinden, waarna exacte symbolische methoden (Gröbner-basis) worden ingezet om een rigoureus wiskundig bewijs te leveren. Dit is een nieuwe standaard voor het oplossen van dergelijke problemen.
Eenheid van Getallen: De rij van maximale groeifactoren voor $n=1$ tot $8$ blijkt een vreemde wiskundige reeks te zijn: $1, 2, 2.25, 4, \text{(wortel graad 61)}, 5, \text{(wortel graad 6)}, 8$ .
Hulde: Het artikel is opgedragen aan de nagedachtenis van Nick Higham, een vooraanstaand expert op het gebied van numerieke linear algebra en stabiliteit, die dit probleem en de technieken in dit artikel fascinerend vond.

Samenvattend biedt dit artikel een doorbraak in het begrijpen van de groeifactor bij Gauss-eliminatie voor $n=5$ , transformeert een numeriek mysterie in een exact algebraïsch resultaat, en stelt een nieuwe methodologische weg voor voor het oplossen van complexe optimalisatieproblemen in de wiskunde.