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.

Previous Showing 11-20 of 9094 results. Next

A004023 Indices of prime repunits: numbers k such that 11...111 (with k 1's) = (10^k - 1)/9 is prime.

Original entry on oeis.org

2, 19, 23, 317, 1031, 49081, 86453, 109297, 270343, 5794777, 8177207
Offset: 1

Views

Author

Keywords

Comments

People who search for repunit primes or repdigit primes may be looking for this entry.
The indices of primes with digital product (i.e., product of digits) equal to 1.
As of August 2014, only the first five repunits, through (10^1031-1)/9, have been proved prime. The next four repunits are known only to be probable primes and have not been proved to be prime. - Robert Baillie, Aug 17 2014
These indices p must also be prime. If p is not prime, say p = m*n, then 10^(m*n) - 1 = ((10^m)^n) - 1 => 10^m - 1 divides 10^(m*n) - 1. Since 9 divides 10^m - 1 or (10^m - 1)/9 = q, it follows q divides (10^p - 1)/9. This is a result of the identity, a^n - b^n = (a - b)(a^(n-1) + a^(n-2)*b + ... + b^(n-1)). - Cino Hilliard, Dec 23 2008
The numbers R_n = 11...111 = (10^n - 1)/9 with n in this sequence A004023, except for n = 2, are prime repunits in base ten, so they are prime Brazilian numbers belonging to A085104. [See Links: Les nombres brésiliens.] - Bernard Schott, Dec 24 2012
Search limit is 10800000, currently. - Serge Batalov, Jul 01 2021
On March 22 2022 the probable prime R49081 was proved to be a prime, and on May 15 2023 the probable prime R86453 was proved to be a prime. - Bassam Abdul-Baki, Dec 17 2024

Examples

			2 appears because the 2-digit repunit 11 is prime.
3 does not appear because 111 = 3 * 37 is not prime.
19 appears because the 19-digit repunit 1111111111111111111 is prime.
		

References

  • J. Brillhart et al., Factorizations of b^n +- 1. Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 2nd edition, 1985; and later supplements.
  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 19, pp 6, Ellipses, Paris 2008.
  • R. K. Guy, Unsolved Problems in Number Theory, Section A3.
  • Graham, Knuth and Patashnik, Concrete Mathematics, Addison-Wesley, 1994; see p 146 problem 22.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 60.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 235.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See entry 142857 at pp. 197-198.

Crossrefs

See A004022 for the actual primes.

Programs

  • Magma
    [p: p in PrimesUpTo(500) | IsPrime((10^p - 1) div 9)]; // Vincenzo Librandi, Nov 06 2014
    
  • Mathematica
    Select[Range[271000], PrimeQ[FromDigits[PadRight[{}, #, 1]]] &] (* Harvey P. Dale, Nov 05 2011 *)
    repUnsUpTo[k_] := ParallelMap[If[PrimeQ[#] && PrimeQ[(10^# - 1)/9], #, Nothing] &, Range[k]]; repUnsUpTo[5000] (* Mikk Heidemaa, Apr 24 2017 *)
  • PARI
    forprime(x=2,20000,if(ispseudoprime((10^x-1)/9),print1(x","))) \\ Cino Hilliard, Dec 23 2008
    
  • Python
    from sympy import isprime; {print(n, end = ', ') for n in range(1, 10**7) if isprime(n) and isprime(10**n//9)} # (Note that sympy.isprime is only a pseudo-primality test.) - Ya-Ping Lu, Dec 20 2021, edited by M. F. Hasler, Mar 28 2022

Extensions

a(6) = 49081 PRP found by Harvey Dubner - posting to Number Theory List (NMBRTHRY(AT)LISTSERV.NODAK.EDU) Sep 09, 1999; proved prime by Paul Underwood, Mar 21 2022.
a(7) = 86453 found using pfgw (a faster version of PrimeForm) on Oct 26 2000 by Lew Baxter (posting to Number Theory List), Oct 26, 2000; proved prime by Andreas Enge, May 16 2023.
a(8) = 109297 was apparently discovered independently by (in alphabetical order) Paul Bourdelais and Harvey Dubner around Mar 26-28 2007.
a(9) = 270343, was found Jul 11 2007 by Maksym Voznyy and Anton Budnyy, subsequently confirmed as a(9) (see Repunit Primes Project link) by Robert Price, Dec 14 2010
a(10) = 5794777 was found Apr 20 2021 by Ryan Propper and Serge Batalov
a(11) = 8177207 was found May 08 2021 by Ryan Propper and Serge Batalov

A029471 Numbers k that divide the (left) concatenation of all numbers <= k written in base 2 (most significant digit on left).

Original entry on oeis.org

1, 85, 145, 245, 1189, 356717, 19590671, 35741759, 791822369, 25313027035
Offset: 1

Views

Author

Keywords

Comments

No other terms below 3*10^10.

Crossrefs

Programs

  • Mathematica
    b = 2; c = {}; Select[Range[10^4], Divisible[FromDigits[c = Join[IntegerDigits[#, b], c], b], #] &] (* Robert Price, Mar 12 2020 *)
  • Python
    from itertools import count
    def a029471():
        total = 0
        power_of_two = 1
        index_of_two = 0
        length_of_string = 0
        for n in count(1):
            total += (n<Christian Perfect, Feb 07 2014
    
  • Python
    def concat_mod(base, k, mod): ...  # See A029479
    for k in range(1, 3*10**10):
      if concat_mod(2, k, k) == 0: print(k) # Jason Yuen, Mar 24 2024

Extensions

One more term from Larry Reeves (larryr(AT)acm.org), Dec 03 2001
Edited and updated by Larry Reeves (larryr(AT)acm.org), Apr 12 2002
a(7)-a(8) from Max Alekseyev, May 12 2011
a(9)-a(10) from Jason Yuen, Mar 24 2024

A057468 Numbers k such that 3^k - 2^k is prime.

Original entry on oeis.org

2, 3, 5, 17, 29, 31, 53, 59, 101, 277, 647, 1061, 2381, 2833, 3613, 3853, 3929, 5297, 7417, 90217, 122219, 173191, 256199, 336353, 485977, 591827, 1059503
Offset: 1

Views

Author

Robert G. Wilson v, Sep 09 2000

Keywords

Comments

Some of the larger entries may only correspond to probable primes.
The 1137- and 1352-digit values associated with the terms 2381 and 2833 have been certified prime with Primo. - Rick L. Shepherd, Nov 12 2002
Or, numbers k such that A001047(k) is prime. - Zak Seidov, Sep 17 2006
3^k - 2^k were proved prime for k = 3613, 3853, 3929, 5297, 7417 with Primo. - David Harrison, Jun 08 2011

Crossrefs

Cf. A058765, A000043 (Mersenne primes), A001047 (3^n-2^n).
Subset of A000040.

Programs

Extensions

a(20) = 90217 found by Mike Oakes, Feb 23 2001
Terms a(21) = 122219, a(22) = 173191, a(23) = 256199 were found by Mike Oakes in 2003-2005. Corresponding numbers of decimal digits are 58314, 82634, 122238.
a(24) = 336353 found by Mike Oakes, Oct 15 2007. It corresponds to a probable prime with 160482 decimal digits.
a(25) = 485977 found by Mike Oakes, Sep 06 2009; it corresponds to a probable prime with 231870 digits. - Mike Oakes, Sep 08 2009
a(26) = 591827 found by Mike Oakes, Aug 25 2009; it corresponds to a probable prime with 282374 digits.
a(27) = 1059503 found by Mike Oakes, Apr 12 2012; it corresponds to a probable prime with 505512 digits. - Mike Oakes, Apr 14 2012

A059801 Numbers k such that 4^k - 3^k is prime.

Original entry on oeis.org

2, 3, 7, 17, 59, 283, 311, 383, 499, 521, 541, 599, 1193, 1993, 2671, 7547, 24019, 46301, 48121, 68597, 91283, 131497, 148663, 184463, 341233
Offset: 1

Views

Author

Mike Oakes, Feb 23 2001

Keywords

Comments

Some of the larger entries may only correspond to probable primes.
The values corresponding to 1193 (719 digits) and 1993 (1200 digits) have been certified prime with Primo. - Rick L. Shepherd, Sep 10 2002
8 more terms found by Jean-Louis Charton during 2004 - 2006. Corresponding numbers of decimal digits are 14461, 27876, 28972, 41300, 54958, 79170, 89505, 111058, 205443. - Alexander Adamchuk, Dec 02 2006

Crossrefs

Programs

A001605 Indices of prime Fibonacci numbers.

Original entry on oeis.org

3, 4, 5, 7, 11, 13, 17, 23, 29, 43, 47, 83, 131, 137, 359, 431, 433, 449, 509, 569, 571, 2971, 4723, 5387, 9311, 9677, 14431, 25561, 30757, 35999, 37511, 50833, 81839, 104911, 130021, 148091, 201107, 397379, 433781, 590041, 593689, 604711, 931517, 1049897, 1285607, 1636007, 1803059, 1968721, 2904353, 3244369, 3340367
Offset: 1

Views

Author

Keywords

Comments

Some of the larger entries may only correspond to probable primes.
Since F(n) divides F(mn) (cf. A001578, A086597), all terms of this sequence are primes except for a(2) = 4 = 2 * 2 but F(2) = 1. - M. F. Hasler, Dec 12 2007
What is the next larger twin prime after F(4) = 3, F(5) = 5, F(7) = 13? The next candidates seem to be F(104911) or F(1968721) (greater of a pair), or F(397379), F(931517) (lesser of a pair). - M. F. Hasler, Jan 30 2013, edited Dec 24 2016, edited Sep 23 2017 by Bobby Jacobs
_Henri Lifchitz_ confirms that the data section gives the full list (49 terms) as far as we know it today of indices of prime Fibonacci numbers (including proven primes and PRPs). - N. J. A. Sloane, Jul 09 2016
Terms n such that n-2 is also a term are listed in A279795. - M. F. Hasler, Dec 24 2016
There are no Fibonacci numbers that are twin primes after F(7) = 13. Every Fibonacci prime greater than F(4) = 3 is of the form F(2*n+1). Since F(2*n+1)+2 and F(2*n+1)-2 are F(n+2)*L(n-1) and F(n-1)*L(n+2) in some order, and F(n+2) > 1, L(n-1) > 1, F(n-1) > 1, and L(n+2) > 1 for n > 3. - Bobby Jacobs, Sep 23 2017
These primes are occurring with about the same normalized frequency as Repunit primes (see Generalized Repunit Conjecture Ref). Assuming a base=1.618 (ratio of sequential terms), then the best fit coefficient is 0.60324 for the first 56 terms, which is already approaching Euler's constant 0.56145948. - Paul Bourdelais, Aug 23 2024

References

  • Clifford A. Pickover, Mazes for the Mind, St. Martin's Press, NY, 1992, p. 350.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 54.
  • Paulo Ribenboim, The Little Book of Big Primes, Springer-Verlag, NY, 1991, p. 178.
  • 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

Subsequence of A046022.
Column k=1 of A303215.

Programs

  • Mathematica
    Select[Range[10^4], PrimeQ[Fibonacci[#]] &] (* Harvey P. Dale, Nov 20 2012 *)
    (* Start ~ 1.8x faster than the above *)
    Select[Range[10^4], PrimeQ[#] && PrimeQ[Fibonacci[#]] &] (* Eric W. Weisstein, Nov 07 2017 *)
    Select[Prime[Range[PrimePi[10^4]]], PrimeQ[Fibonacci[#]] &] (* Eric W. Weisstein, Nov 07 2017 *)
    (* End *)
  • PARI
    v=[3,4]; forprime(p=5,1e5, if(ispseudoprime(fibonacci(p)), v=concat(v,p))); v \\ Charles R Greathouse IV, Feb 14 2011
    
  • PARI
    is_A001605(n)={n==4 || isprime(n) & ispseudoprime(fibonacci(n))}  \\ M. F. Hasler, Sep 29 2012

Formula

Prime(i) = a(n) for some n <=> A080345(i) <= 1. - M. F. Hasler, Dec 12 2007

Extensions

Additional comments from Robert G. Wilson v, Aug 18 2000
More terms from David Broadhurst, Nov 08 2001
Two more terms (148091 and 201107) from T. D. Noe, Feb 12 2003 and Mar 04 2003
397379 from T. D. Noe, Aug 18 2003
433781, 590041, 593689 from Henri Lifchitz submitted by Ray Chandler, Feb 11 2005
604711 from Henri Lifchitz communicated by Eric W. Weisstein, Nov 29 2005
931517, 1049897, 1285607 found by Henri Lifchitz circa Nov 01 2008 and submitted by Alexander Adamchuk, Nov 28 2008
1636007 from Henri Lifchitz March 2009, communicated by Eric W. Weisstein, Apr 24 2009
1803059 and 1968721 from Henri Lifchitz, November 2009, submitted by Alex Ratushnyak, Aug 08 2012
a(49)=2904353 from Henri Lifchitz, Jul 15 2014
a(50)=3244369 from Henri Lifchitz, Nov 04 2017
a(51)=3340367 from Henri Lifchitz, Apr 25 2018
a(52)-a(56) from Ryan Propper added by Paul Bourdelais, Aug 23 2024

A002981 Numbers k such that k! + 1 is prime.

Original entry on oeis.org

0, 1, 2, 3, 11, 27, 37, 41, 73, 77, 116, 154, 320, 340, 399, 427, 872, 1477, 6380, 26951, 110059, 150209, 288465, 308084, 422429
Offset: 1

Views

Author

Keywords

Comments

If n + 1 is prime then (by Wilson's theorem) n + 1 divides n! + 1. Thus for n > 2 if n + 1 is prime n is not in the sequence. - Farideh Firoozbakht, Aug 22 2003
For n > 2, n! + 1 is prime <==> nextprime((n+1)!) > (n+1)nextprime(n!) and we can conjecture that for n > 2 if n! + 1 is prime then (n+1)! + 1 is not prime. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 03 2004
The prime members are in A093804 (numbers n such that Sum_{d|n} d! is prime) since Sum_{d|n} d! = n! + 1 if n is prime. - Jonathan Sondow
150209 is also in the sequence, cf. the link to Caldwell's prime pages. - M. F. Hasler, Nov 04 2011

Examples

			3! + 1 = 7 is prime, so 3 is in the sequence.
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 116, p. 40, Ellipses, Paris 2008.
  • Harvey Dubner, Factorial and primorial primes, J. Rec. Math., 19 (No. 3, 1987), 197-203.
  • Richard K. Guy, Unsolved Problems in Number Theory, Section A2.
  • F. Le Lionnais, Les Nombres Remarquables, Paris, Hermann, 1983, p. 100.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 118.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 70.

Crossrefs

Cf. A002982 (n!-1 is prime), A064295. A088332 gives the primes.
Equals A090660 - 1.
Cf. A093804.

Programs

  • Magma
    [n: n in [0..800] | IsPrime(Factorial(n)+1)]; // Vincenzo Librandi, Oct 31 2018
    
  • Mathematica
    v = {0, 1, 2}; Do[If[ !PrimeQ[n + 1] && PrimeQ[n! + 1], v = Append[v, n]; Print[v]], {n, 3, 29651}]
    Select[Range[100], PrimeQ[#! + 1] &] (* Alonso del Arte, Jul 24 2014 *)
  • PARI
    for(n=0,500,if(ispseudoprime(n!+1),print1(n", "))) \\ Charles R Greathouse IV, Jun 16 2011
    
  • Python
    from sympy import factorial, isprime
    for n in range(0,800):
        if isprime(factorial(n)+1):
            print(n, end=', ') # Stefano Spezia, Jan 10 2019

Extensions

a(19) sent in by Jud McCranie, May 08 2000
a(20) from Ken Davis (kraden(AT)ozemail.com.au), May 24 2002
a(21) found by PrimeGrid around Jun 11 2011, submitted by Eric W. Weisstein, Jun 13 2011
a(22) from Rene Dohmen, Jun 09 2012
a(23) from Rene Dohmen, Jan 12 2022
a(24)-a(25) from Dmitry Kamenetsky, Jun 19 2024

A059802 Numbers k such that 5^k - 4^k is prime.

Original entry on oeis.org

3, 43, 59, 191, 223, 349, 563, 709, 743, 1663, 5471, 17707, 19609, 35449, 36697, 45259, 91493, 246497, 265007, 289937
Offset: 1

Views

Author

Mike Oakes, Feb 23 2001

Keywords

Comments

Some of the larger terms may only correspond to probable primes.
5^1663 - 4^1663, a 1163-digit number, has been certified prime with Primo. - Rick L. Shepherd, Nov 13 2002
4 more terms found by Predrag Minovic in 2004: 35449, 36697, 45259, 91493. Corresponding numbers of decimal digits are 24778, 25651, 31635, 63951. - Alexander Adamchuk, Dec 02 2006

Crossrefs

Programs

  • Mathematica
    Select[Range[1000], PrimeQ[5^# - 4^#] &] (* Alonso del Arte, Sep 09 2013 *)
  • PARI
    forprime(p=2,1e5,if(ispseudoprime(5^p-4^p),print1(p", "))) \\ Charles R Greathouse IV, Jun 10 2011

Extensions

New term 246497 found by Jean-Louis Charton in 2008 corresponding to a probable prime with 172295 digits - Jean-Louis Charton, Sep 02 2009
New term a(19) = 265007 found by Jean-Louis Charton, Feb 19 2013
a(20) = 289937 found by Jean-Louis Charton, Mar 15 2013

A062572 Numbers k such that 6^k - 5^k is prime.

Original entry on oeis.org

2, 5, 11, 13, 23, 61, 83, 421, 1039, 1511, 31237, 60413, 113177, 135647, 258413
Offset: 1

Views

Author

Mike Oakes, May 18 2001, May 19 2001

Keywords

Comments

The 809- and 1176-digit numbers associated with the terms 1039 and 1511 have been certified prime with Primo. - Rick L. Shepherd, Nov 15 2002

Examples

			2 is in the sequence because 6^2 - 5^2 = 36 - 25 = 11, which is prime.
3 is not in the sequence because 6^3 - 5^3 = 216 - 125 = 91 = 7 * 13, which is not prime.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[1000], PrimeQ[6^# - 5^#] &] (* Alonso del Arte, Sep 04 2013 *)
  • PARI
    forprime(p=2,1e4,if(ispseudoprime(6^n-5^n),print1(p", "))) \\ Charles R Greathouse IV, Jun 10 2011

Extensions

Edited by T. D. Noe, Oct 30 2008
Two more terms (31237 and 60413) found by Predrag Minovic in 2004 corresponding to probable primes with 24308 and 47011 digits. Jean-Louis Charton, Oct 06 2010
Two more terms (113177 and 135647) found by Jean-Louis Charton in 2009 corresponding to probable primes with 88069 and 105554 digits. Jean-Louis Charton, Oct 13 2010
a(15) from Jean-Louis Charton, Apr 08 2013

A000372 Dedekind numbers or Dedekind's problem: number of monotone Boolean functions of n variables, number of antichains of subsets of an n-set, number of elements in a free distributive lattice on n generators, number of Sperner families.

Original entry on oeis.org

2, 3, 6, 20, 168, 7581, 7828354, 2414682040998, 56130437228687557907788, 286386577668298411128469151667598498812366
Offset: 0

Views

Author

Keywords

Comments

A monotone Boolean function is an increasing function from P(S), the set of subsets of S, to {0,1}.
The count of antichains includes the empty antichain which contains no subsets and the antichain consisting of only the empty set.
a(n) is also equal to the number of upsets of an n-set S. A set U of subsets of S is an upset if whenever A is in U and B is a superset of A then B is in U. - W. Edwin Clark, Nov 06 2003
Also the number of simple games with n players in minimal winning form. - Fabián Riquelme, May 29 2011
The unlabeled case is A003182. - Gus Wiseman, Feb 20 2019
From Amiram Eldar, May 28 2021 and Michel Marcus, Apr 07 2023: (Start)
The terms were first calculated by:
a(0)-a(4) - Dedekind (1897)
a(5) - Church (1940)
a(6) - Ward (1946)
a(7) - Church (1965, verified by Berman and Kohler, 1976)
a(8) - Wiedemann (1991)
a(9) - Jäkel (2023)
a(9) - independently computed by Lennart Van Hirtum, Patrick De Causmaecker, Jens Goemaere, Tobias Kenter, Heinrich Riebler, Michael Lass, and Christian Plessl (2023)
(End)

Examples

			a(2)=6 from the antichains {}, {{}}, {{1}}, {{2}}, {{1,2}}, {{1},{2}}.
From _Gus Wiseman_, Feb 20 2019: (Start)
The a(0) = 2 through a(3) = 20 antichains:
  {}    {}     {}        {}
  {{}}  {{}}   {{}}      {{}}
        {{1}}  {{1}}     {{1}}
               {{2}}     {{2}}
               {{12}}    {{3}}
               {{1}{2}}  {{12}}
                         {{13}}
                         {{23}}
                         {{123}}
                         {{1}{2}}
                         {{1}{3}}
                         {{2}{3}}
                         {{1}{23}}
                         {{2}{13}}
                         {{3}{12}}
                         {{12}{13}}
                         {{12}{23}}
                         {{13}{23}}
                         {{1}{2}{3}}
                         {{12}{13}{23}}
(End)
		

References

  • Ian Anderson, Combinatorics of Finite Sets. Oxford Univ. Press, 1987, p. 38.
  • Jorge Luis Arocha, Antichains in ordered sets [in Spanish], Anales del Instituto de Matematicas de la Universidad Nacional Autonoma de Mexico, Vol. 27 (1987), pp. 1-21.
  • Joel Berman and Peter Koehler, Cardinalities of finite distributive lattices, Mitteilungen aus dem Mathematischen Seminar Giessen, Vol. 121 (1976), pp. 103-124.
  • Garrett Birkhoff, Lattice Theory, American Mathematical Society, Colloquium Publications, Vol. 25, 3rd ed., Providence, RI, 1967, p. 63.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 273.
  • Michael A. Harrison, Introduction to Switching and Automata Theory, McGraw Hill, NY, 1965, p. 188.
  • Donald E. Knuth, The Art of Computer Programming, Vol. 4A, Section 7.1.1, p. 79.
  • A. D. Korshunov, The number of monotone Boolean functions, Problemy Kibernet, No. 38, (1981), 5-108, 272. MR0640855 (83h:06013)
  • W. F. Lunnon, The IU function: the size of a free distributive lattice, in D. J. A. Welsh, editor, Combinatorial Mathematics and Its Applications. Academic Press, NY, 1971, pp. 173-181.
  • Saburo Muroga, Threshold Logic and Its Applications. Wiley, NY, 1971, pp. 38 and 214.
  • R. A. Obando, On the number of nondegenerate monotone boolean functions of n variables in an n-variable boolean algebra. In preparation.
  • 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).
  • Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001, p. 349.

Crossrefs

Programs

  • Mathematica
    nn=5;
    stableSets[u_,Q_]:=If[Length[u]===0,{{}},With[{w=First[u]},Join[stableSets[DeleteCases[u,w],Q],Prepend[#,w]&/@stableSets[DeleteCases[u,r_/;r===w||Q[r,w]||Q[w,r]],Q]]]];
    Table[Length[stableSets[Subsets[Range[n]],SubsetQ]],{n,0,nn}] (* Gus Wiseman, Feb 20 2019 *)
    Table[Total[Boole[Table[UnateQ[BooleanFunction[k, n]], {k, 0, 2^(2^n) - 1}]]], {n, 0, 4}] (* Eric W. Weisstein, Jun 27 2023 *)

Formula

The asymptotics can be found in the Korshunov paper. - Boris Bukh, Nov 07 2003
a(n) = Sum_{k=1..n} binomial(n,k)*A006126(k) + 2, i.e., this sequence is the inverse binomial transform of A006126, plus 2. E.g., a(3) = 3*1 + 3*2 + 1*9 + 2 = 20. - Rodrigo A. Obando (R.Obando(AT)computer.org), Jul 26 2004
From J. M. Aranda, Jun 12 2021: (Start)
a(n) = A132581(2^n) = A132581(2^n-2^m) + A132581(2^n-2^(n-m)) for n >= m >= 0.
a(n) = A132582(3*2^n -1) for n >= 0.
(End)

Extensions

a(8) from D. H. Wiedemann, personal communication, Nov 03 1990
Additional comments from Michael Somos, Jun 10 2002
a(9) from C. Jäkel added by Michel Marcus, Apr 04 2023

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

Original entry on oeis.org

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

Views

Author

Keywords

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
In Ribenboim 2004, a wrong value of a(8) is given, 6221271 instead of 6221671. - Stefano Spezia, Mar 27 2025

Examples

			a(5) is equal to 13 because 2*3*7*43 + 1 = 1807 = 13 * 139.
		

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.
  • 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).
  • Samuel S. Wagstaff, Jr., Computing Euclid's primes, Bull. Institute Combin. Applications, Vol. 8 (1993), pp. 23-32.

Crossrefs

Programs

  • 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}]
    nxt[{p_,a_}]:=With[{c=FactorInteger[p+1][[1,1]]},{p*c,c}]; Rest[NestList[nxt,{1,2},20][[;;,2]]] (* Harvey P. Dale, Feb 02 2025 *)
  • 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
Previous Showing 11-20 of 9094 results. Next