En la monografía de Edsger Dijkstra "Notes on Structured Programming", describe un sencillo programa imperativo para generar una matriz de la primera $n$ primos. Para cada primo $p_n$ encuentra el siguiente primo comprobando cada impar $j \gt p_n$ para ver si tiene un divisor primo $d$ en la gama $2 \lt d \le \sqrt j$ . El algoritmo extrae estos candidatos de la matriz de números primos que ha ido construyendo, que en cualquier iteración contiene los primos comprendidos entre $2$ y $p_n$ inclusivo, por lo que hay una suposición implícita es que $\sqrt j \le p_n$ (el algoritmo itera sobre la matriz de primos hasta $p_n$ siendo las únicas condiciones de terminación $d|j$ o $d > \sqrt j$ por lo que si $\sqrt j$ podría ser $\ge p_n$ corremos el riesgo de leer más allá del final de la matriz). Citando a Dijkstra:
Con toda humildad cito el comentario de Don Knuth sobre una versión anterior de este programa, en la que daba por sentado este hecho:
"¡Aquí es usted culpable de una grave omisión! Su programa hace uso de un profundo resultado de la teoría de números, a saber, que si $p_n$ denota el $n$ ésimo número primo siempre tienen $p_{n+1} < p_n^2$ ."
Así que tengo curiosidad por conocer la historia de este "resultado profundo" y la dificultad de su demostración.