OFFSET
1,1
COMMENTS
There are also composites m = A161671(k) = A161671(k+1). An example is 65 = A161671(26) = A161671(27). - Michael De Vlieger, Mar 22 2022
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Scatterplot of A161671(n), n = 1..120, showing and labeling primes p in this sequence in red.
Michael De Vlieger, Scatterplot of A161671(n), n = 1..2^12, showing primes p in this sequence in red.
MATHEMATICA
f[n_] := FixedPoint[n + PrimePi@ # &, n + PrimePi@ n]; Reap[Do[(If[PrimeQ[#] && # == j, Sow[#]]; j = #) &[Prime[i] - f[i - 1] ], {i, 8500}] ][[-1, -1]] (* Michael De Vlieger, Mar 22 2022, after Robert G. Wilson v at A141468 *)
CROSSREFS
KEYWORD
nonn,less
AUTHOR
Juri-Stepan Gerasimov, Dec 07 2012
STATUS
approved