login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A000945
Euclid-Mullin sequence: a(1) = 2, a(n+1) is smallest prime factor of 1 + Product_{k=1..n} a(k).
(Formerly M0863 N0329)
101
2, 3, 7, 43, 13, 53, 5, 6221671, 38709183810571, 139, 2801, 11, 17, 5471, 52662739, 23003, 30693651606209, 37, 1741, 1313797957, 887, 71, 7127, 109, 23, 97, 159227, 643679794963466223081509857, 103, 1079990819, 9539, 3143065813, 29, 3847, 89, 19, 577, 223, 139703, 457, 9649, 61, 4357
OFFSET
1,1
COMMENTS
"Does the sequence ... contain every prime? ... [It] was considered by Guy and Nowakowski and later by Shanks, [Wagstaff 1993] computed the sequence through the 43rd term. The computational problem inherent in continuing the sequence further is the enormous size of the numbers that must be factored. Already the number a(1)* ... *a(43) + 1 has 180 digits." - Crandall and Pomerance
If this variant of Euclid-Mullin sequence is initiated either with 3, 7 or 43 instead of 2, then from a(5) onwards it is unchanged. See also A051614. - Labos Elemer, May 03 2004
Wilfrid Keller informed me that a(1)* ... *a(43) + 1 was factored as the product of two primes on Mar 09 2010 by the GNFS method. See the post in the Mersenne Forum for more details. The smaller 68-digit prime is a(44). Terms a(45)-a(47) were easy to find. Finding a(48) will require the factorization of a 256-digit number. See the b-file for the four new terms. - T. D. Noe, Oct 15 2010
On Sep 11 2012, Ryan Propper factored the 256-digit number by finding a 75-digit factor by using ECM. Finding a(52) will require the factorization of a 335-digit number. See the b-file for the terms a(48) to a(51). - V. Raman, Sep 17 2012
Needs longer b-file. - N. J. A. Sloane, Dec 18 2015
A056756 gives the position of the k-th prime in this sequence for each k. - Jianing Song, May 07 2021
Named after the Greek mathematician Euclid (flourished c. 300 B.C.) and the American engineer and mathematician Albert Alkins Mullin (1933-2017). - Amiram Eldar, Jun 11 2021
REFERENCES
Richard Crandall and Carl Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2001; see p. 6.
Richard Guy and Richard Nowakowski, Discovering primes with Euclid, Delta (Waukesha), Vol. 5, pp. 49-63, 1975.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Samuel S. Wagstaff, Jr., Computing Euclid's primes, Bull. Institute Combin. Applications, Vol. 8 (1993), pp. 23-32.
LINKS
Ryan Propper, Table of n, a(n) for n = 1..51 (first 47 terms from T. D. Noe)
Andrew R. Booker, On Mullin's second sequence of primes, Integers, Vol. 12, No. 6 (2012), pp. 1167-1177; arXiv preprint, arXiv:1107.3318 [math.NT], 2011-2013.
Andrew R. Booker, A variant of the Euclid-Mullin sequence containing every prime, arXiv preprint arXiv:1605.08929 [math.NT], 2016.
Andrew R. Booker and Sean A. Irvine, The Euclid-Mullin graph, Journal of Number Theory, Vol. 165 (2016), pp. 30-57; arXiv preprint, arXiv:1508.03039 [math.NT], 2015-2016.
Cristian Cobeli and Alexandru Zaharescu, Promenade around Pascal Triangle-Number Motives, Bull. Math. Soc. Sci. Math. Roumanie, Vol. 56(104), No. 1 (2013), pp. 73-98.
Keith Conrad, The infinitude of the primes, University of Connecticut, 2020.
C. D. Cox and A. J. van der Poorten, On a sequence of prime numbers, Journal of the Australian Mathematical Society, Vol. 8 (1968), pp. 571-574.
FactorDB, Status of EM51.
Richard Guy and Richard Nowakowski, Discovering primes with Euclid, Research Paper No. 260 (Nov 1974), The University of Calgary Department of Mathematics, Statistics and Computing Science.
Lucas Hoogendijk, Prime Generators, Bachelor Thesis, Utrecht University (Netherlands, 2020).
Robert R. Korfhage, On a sequence of prime numbers, Bull Amer. Math. Soc., Vol. 70 (1964), pp. 341, 342, 747. [Annotated scanned copy]
Evelyn Lamb, A Curious Sequence of Prime Numbers, Scientific American blog (2019).
Des MacHale, Infinitely many proofs that there are infinitely many primes, Math. Gazette, Vol. 97, No. 540 (2013), pp. 495-498.
Mersenne Forum, Factoring EM47.
Albert A. Mullin, Research Problem 8: Recursive function theory, Bull. Amer. Math. Soc., Vol. 69, No. 6 (1963), p. 737.
Thorkil Naur, Letter to N. J. A. Sloane, Aug 27 1991, together with copies of "Mullin's sequence of primes is not monotonic" (1984) and "New integer factorizations" (1983) [Annotated scanned copies]
Paul Pollack and Enrique Treviño, The Primes that Euclid Forgot, Amer. Math. Monthly, Vol. 121, No. 5 (2014), pp. 433-437. MR3193727; alternative link.
Daphne Stouthart, Euclid and the infinite number of missing primes, Bachelor Thesis, Utrecht Univ (Netherlands, 2024). See p. 1.
Samuel S. Wagstaff, Jr., Emails to N. J. A. Sloane, May 30 1991.
Samuel S. Wagstaff, Jr., Computing Euclid's primes, Bull. Institute Combin. Applications, Vol. 8 (1993), pp. 23-32. (Annotated scanned copy)
EXAMPLE
a(5) is equal to 13 because 2*3*7*43 + 1 = 1807 = 13 * 139.
MAPLE
a :=n-> if n = 1 then 2 else numtheory:-divisors(mul(a(i), i = 1 .. n-1)+1)[2] fi: seq(a(n), n=1..15);
# Robert FERREOL, Sep 25 2019
MATHEMATICA
f[1]=2; f[n_] := f[n] = FactorInteger[Product[f[i], {i, 1, n - 1}] + 1][[1, 1]]; Table[f[n], {n, 1, 46}]
PROG
(PARI) print1(k=2); for(n=2, 20, print1(", ", p=factor(k+1)[1, 1]); k*=p) \\ Charles R Greathouse IV, Jun 10 2011
(PARI) P=[]; until(, print(P=concat(P, factor(vecprod(P)+1)[1, 1]))) \\ Jeppe Stig Nielsen, Apr 01 2024
KEYWORD
nonn,nice,hard
STATUS
approved