A144717 a(n) = smallest positive integer > a(n-1) such that 2*a(1)*a(2)*...*a(n) + 1 is prime.

1, 2, 3, 5, 7, 8, 9, 11, 12, 14, 17, 20, 24, 30, 34, 44, 72, 85, 86, 92, 115, 122, 125, 132, 142, 150, 161, 162, 181, 186, 198, 224, 248, 252, 282, 283, 290, 307, 319, 321, 344, 350, 376, 445, 476, 567, 623, 676, 682, 704, 741, 749, 786, 803, 806, 893, 1014, 1046

Offset

1,2

Links

Jon E. Schoenfield, Table of n, a(n) for n = 1..505 (lists all terms < 10^5)

Example

a(1)=1 because a(0) is not defined and 2*1 + 1 = 3 is prime;
a(2)=2 because 2*1*2 + 1 = 5 is prime;
a(3)=3 because 2*1*2*3 + 1 = 13 is prime;
a(4) is not 4 because 2*1*2*3*4 + 1 = 49 is not prime, but a(4)=5 works because 2*1*2*3*5 + 1 = 61 is prime.

Mathematica

k = 2; a = {}; Do[If[PrimeQ[k n + 1], k = k n; AppendTo[a, n]], {n, 1, 3000}]; a (* Artur Jasinski *)
nxt[{p_, a_}]:=Module[{k=a+1}, While[!PrimeQ[p*k+1], k++]; {p*k, k}]; NestList[ nxt, {2, 1}, 60][[All, 2]] (* Harvey P. Dale, Aug 18 2021 *)

Prog

(Python)
from sympy import isprime
from itertools import count, islice
def agen(): # generator of terms
    an, p = 1, 2
    while True:
        yield an
        an = next(k for k in count(an+1) if isprime(p*k+1))
        p *= an
print(list(islice(agen(), 58))) # Michael S. Branicky, Jan 13 2023

Crossrefs

Cf. A046966, A046972, A144718, A144722, A144723, A144724, A144725, A144726, A144727, A144728, A144729, A144730, A144731.

Sequence in context: A350799 A246438 A158791 * A202728 A244016 A186042

Adjacent sequences: A144714 A144715 A144716 * A144718 A144719 A144720

Keyword

nonn,nice

Author

Artur Jasinski, Sep 19 2008

Extensions

Edited by N. J. A. Sloane, Sep 21 2017 following suggestions from Richard C. Schroeppel

Status

approved