Un contraejemplo es $n_1 := 6561101970383$ con $$ \log_{10} \left( (n_1/e)^{n_1} \sqrt{2\pi n_1} \right) = 81244041273652.999999999999995102483 - \phantom; , $$ pero $$ \log_{10} (n_1!) = 81244041273653.000000000000000618508 + \phantom;. $$ Si he calculado bien, $n_1$ es el primer contraejemplo, y el único hasta $10^{13}$ . El cálculo debe alcanzar $10^{15}$ en algún momento de la próxima semana, con una probabilidad de alrededor de $1 - \exp(-\frac16) \sim 15\%$ de encontrar un $n_2$ .
El cálculo (en gp/pari ) necesitó aquí unas 40 horas de CPU, que se redujeron a 4 horas al ejecutarse en paralelo en 10 de los 12 cabezales de alhambra.math.harvard.edu . Esto fue no calculando $\log_{10} (n!)$ con suficiente precisión para cada $n \leq 10^{13}$ lo que habría llevado cientos de veces más tiempo. El problema de encontrar valores casi integrales de $\log_{10} (n!)$ es un caso especial del "dilema del fabricante de mesas" (Wikipedia atribuye esta feliz acuñación a William Kahan); en este caso, la técnica de aproximación lineal sugerida por Lefèvre al final de la página 15 de su diapositivas lleva tiempo $\tilde O(N^{2/3})$ para encontrar todos los ejemplos con $n < N$ . Eso es lo que se está ejecutando en alhambra ahora.
Por el camino aparecieron algunos términos más de la secuencia A177901: $252544447$ , $1430841730$ , $5042264463$ , $31774693500$ , $40752166709$ , $46787073630$ , $129532358256$ , $421559495894$ , $2418277169072$ , $6105111564681$ , y entonces $n_1 = 6561101970383$ que podría incluso convertirse en el último plazo hasta el $10^{15}$ porque $\log_{10} (n_1!)$ está tan cerca de un número entero (aproximadamente $9$ veces más cerca de lo necesario para nuestro propósito). [ EDITAR Es la última legislatura $<10^{14}$ pero no $10^{15}$ ver más abajo]. El término $252544447$ fue reportado en math.se #8323 por Byron Schmuland [ EDITAR y unos meses antes por David Cantrell en sci.math], aunque todavía no se ha publicado en OEIS. Los siguientes parecen ser nuevos, y los publicaré en OEIS en breve.
Kamenetsky tenía razón al sugerir que la aproximación debería fallar a veces: en base 10, esperamos que $n$ sea un contraejemplo con una probabilidad aproximada de $1/cn$ con $c = 12 \log 10$ por lo que, en promedio, cada rango $[N, 10^{12}N]$ debería tener alrededor de uno. Así pues, no es de extrañar que el primero (pasado $n=1!$ ) resulta tener $13$ dígitos. Esta heurística es también la fuente de la estimación $1-\exp(-\frac16)$ para la probabilidad de otro contraejemplo en $ [10^{13}, 10^{15}]$ .
ACTUALIZACIÓN El cálculo ya ha pasado $10^{14}$ sin encontrar un nuevo contraejemplo. Sin embargo, sí encontró un nuevo término para la secuencia OEIS un poco más allá de $10^{14}$ : $n=125291661119688$ con $\log_{10}(n!)$ cercano pero justo por debajo del entero más próximo $1711938609606982$ (donde un contraejemplo debe estar un poco por encima), y tampoco tan cerca como $1/(12n)$ - la diferencia es de $1/(8.4n)$ .
Ya que estoy: Debería haber mencionado que el cálculo gp/pari también encontró (en un minuto o dos) todos los términos en $[10^4,10^8]$ de la OEIS, lo que da cierta credibilidad a los nuevos resultados; y agradezco a Gerry Myerson que me llamara la atención sobre esta cuestión con su edición de hace unas dos semanas.
