OFFSET
1,1
COMMENTS
This allows larger terms of A051254 (which triple in digits each entry) to be given. Like A051254, currently requires Riemann Hypothesis to show sequence continues.
Currently a(11)=66768 generates only a probable prime number. - Arkadiusz Wesolowski, May 28 2011
Likewise a(12) and a(13) generate only a probable prime numbers, as well as being conditional on a(11) and a(12) being proved primes. Minimality of a(12)-a(13) is exhaustively tested. - Serge Batalov, Aug 06 2013
a(14) = 8436308 is found by Ryan Propper and Serge Batalov, Apr 29 2024, but a few remaining gaps below this value were being double-checked. The double-check is now complete (see GitHub link). - Ryan Propper and Serge Batalov, May 24 2024.
REFERENCES
T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 8.
LINKS
Chris K. Caldwell, Mills' Theorem - a generalization.
Chris K. Caldwell, The List of Largest Known Primes, The 11th Mills' prime
Chris K. Caldwell and Yuanyou Cheng, Determining Mills' Constant and a Note on Honaker's Problem, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.1.
Sergey Batalov, Mills' b(14) search report
Henri & Renaud Lifchitz, PRP Records and PRP Records, search:Mills
W. H. Mills, A prime-representing function, Bull. Amer. Math. Soc., Vol. 53 (1947), p. 604.
Eric Weisstein's World of Mathematics, Mills' Constant
Eric Weisstein's World of Mathematics, Mills' Prime
E. M. Wright, A class of representing functions, J. London Math. Soc., Vol. 29 (1954) pp. 63-71.
FORMULA
b(1) = 2; b(n+1) = nextprime(b(n)^3); a(n) = b(n+1)-b(n)^3;
EXAMPLE
The Mills' primes (given in A051254) are 2, 2^3+3 = 11, (2^3+3)^3+30 = 11^3+30 = 1361, ((2^3+3)^3+30)^3+6 = 1361^3+6 = 2521008887, etc. The terms added at each step yield this sequence. They are the least positive integers which added to the cube of the preceding prime yield again a prime, cf. formula. - M. F. Hasler, Jul 22 2013
MATHEMATICA
B[1] = 2; B[n_] := B[n] = NextPrime[B[n - 1]^3]; Table[B[n + 1] - B[n]^3, {n, 7}] (* Robert Price, Jun 09 2019 *)
PROG
(PARI) p=2; until(, np=nextprime(p^3); print1(np-p^3, ", "); p=np) \\ Jeppe Stig Nielsen, Apr 22 2020
CROSSREFS
KEYWORD
more,nonn,hard
AUTHOR
Chris K. Caldwell, Jun 22 2005
EXTENSIONS
a(9)-a(11) from Caldwell and Cheng, Aug 29 2005
Corrected by T. D. Noe, Sep 24 2008
a(12) (which generates a PRP) from Serge Batalov, Jul 19 2013
a(13) (which generates a PRP) from Serge Batalov, Aug 06 2013
a(14) (which generates a PRP) from Ryan Propper and Serge Batalov, May 24 2024
STATUS
approved