cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A000946 Euclid-Mullin sequence: a(1) = 2, a(n+1) is the largest prime factor of 1 + Product_{k=1..n} a(k).

Original entry on oeis.org

2, 3, 7, 43, 139, 50207, 340999, 2365347734339, 4680225641471129, 1368845206580129, 889340324577880670089824574922371, 20766142440959799312827873190033784610984957267051218394040721
Offset: 1

Views

Author

Keywords

Comments

Cox and van der Poorten show that 5, 11, 13, 17, ... (A216227) are not members of this sequence. - Charles R Greathouse IV, Jul 02 2007
Booker's abstract claims: "We consider the second of Mullin's sequences of prime numbers related to Euclid's proof that there are infinitely many primes. We show in particular that it omits infinitely many primes, confirming a conjecture of Cox and van der Poorten."

References

  • R. K. Guy and R. Nowakowski, Discovering primes with Euclid, Delta (Waukesha), Vol. 5, pp. 49-63, 1975.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 5.
  • 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).

Crossrefs

Programs

  • Mathematica
    f[1] = 2; f[n_] := f[n] = FactorInteger[Product[f[i], {i, 1, n - 1}] + 1][[-1, 1]]; Table[f[n], {n, 1, 10}] (* Alonso del Arte, Jun 25 2011 based on the program given for A000945 *)
  • PARI
    gpf(n)=my(f=factor(n)[, 1]);f[#f];
    first(m)=my(v=vector(m));v[1]=2;for(i=2,m,v[i]=gpf(1+prod(j=1,i-1,v[j])));v; \\ Anders Hellström, Aug 14 2015

Extensions

Extended by Andrew R. Booker, Mar 13 2013