A284397 a(1) = 1, a(2) = 2; a(n) is the largest prime <= (a(a(n-1)) + a(n-a(n-2)-1)) for n > 2.

1, 2, 3, 3, 3, 5, 5, 5, 5, 5, 5, 7, 7, 7, 7, 7, 7, 7, 7, 11, 11, 7, 7, 11, 11, 11, 11, 11, 11, 11, 11, 13, 17, 17, 13, 13, 13, 17, 17, 13, 13, 17, 17, 17, 17, 17, 17, 17, 17, 19, 23, 19, 17, 23, 19, 19, 19, 23, 23, 19, 19, 23, 23, 19, 19, 23, 23, 23, 23, 23, 23, 23, 23, 23

Offset

1,2

Links

Altug Alkan, Table of n, a(n) for n = 1..10000

Altug Alkan, Alternative scatterplot of A284397

Altug Alkan, Scatterplot of 4*a(n)-n

Formula

a(1) = 1, a(2) = 2; a(n) = A007917(a(a(n-1)) + a(n-a(n-2)-1)) for n > 2.

Example

a(6) = 5 because a(a(5)) + a(6 - a(4) - 1) = a(3) + a(2) = 3 + 2 = 5 and A007917(5) = 5.

Mathematica

a[n_]:= If[n<3, n, Prime@ PrimePi[a[a[n - 1]] + a[n - a[n - 2] - 1]]]; Table[a[n], {n, 1, 50}] (* Indranil Ghosh, Mar 26 2017 *)

Prog

(PARI) a=vector(1000); a[1]=1; a[2]=2; for(n=3, #a, a[n] = precprime(a[a[n-1]]+a[n-a[n-2]-1])); a

Crossrefs

Cf. A005185, A007917, A055748, A284374, A284383.

Sequence in context: A372546 A136348 A353710 * A284374 A137847 A295820

Adjacent sequences: A284394 A284395 A284396 * A284398 A284399 A284400

Keyword

nonn

Author

Altug Alkan, Mar 26 2017

Status

approved