A105120 a(1) = 2; k(1) = 0; for n > 1: k(n) = smallest number j >= k(n-1) such that 2*a(n-1) + j is prime; a(n) = 2*a(n-1) + k(n).

2, 5, 11, 23, 47, 97, 197, 397, 797, 1597, 3203, 6421, 12889, 25841, 51749, 103567, 207227, 414553, 829211, 1658533, 3317177, 6634469, 13269059, 26538257, 53076679, 106153547, 212307299, 424614829, 849229907, 1698460067, 3396920419

Offset

1,1

Comments

Numbers k(n) are given in A105121.

a(n) appears to tend toward C*A055496(n), C~ 0.992521946129820000. - Bill McEachen, Feb 21 2022

Links

Robert G. Wilson v, Table of n, a(n) for n = 1..3000..

Example

a(10) = 1597; k(10) = 3; 2*1597 + j is not prime for 3 <= j < 9, but 2*1597 + 9 = 3203 is prime. Hence k(11) = 9 and a(11) = 3203.

Mathematica

a[1] = {2, 0}; a[n_] := a[n] = Block[{m = 2a[n - 1][[1]], k = a[n - 1][[2]]}, While[ !PrimeQ[m + k], k++ ]; {m + k, k}]; Table[ a[n][[1]], {n, 30}] (* Robert G. Wilson v, Apr 08 2005 *)

Prog

(PARI) print1(a=2, ", "); k=0; for(n=2, 31, j=k; while(!isprime(2*a+j), j++); k=j; print1(a=2*a+k, ", ")) \\ Klaus Brockhaus, Apr 08 2005

Crossrefs

k(n) is in A105121.

Sequence in context: A357292 A334276 A055496 * A084403 A261201 A055011

Adjacent sequences: A105117 A105118 A105119 * A105121 A105122 A105123

Keyword

nonn

Author

Yasutoshi Kohmoto, Apr 08 2005

Extensions

Edited, corrected and extended by Klaus Brockhaus and Robert G. Wilson v, Apr 08 2005

Status

approved