A220220 Primes p of the form p = A161671(k) = A161671(k+1).

2, 7, 43, 311, 491, 827, 1367, 1693, 1733, 1741, 2089, 2239, 2927, 3343, 5231, 5743, 9319, 9521, 11177, 12611, 13249, 15511, 16661, 17989, 24083, 24611, 25679, 25841, 28723, 37861, 39199, 46663, 47279, 51659, 53281, 58031, 58309, 58549, 59723, 64091, 68041, 70051, 70913, 71261

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.

Example

The prime 2 is in the sequence because 2 = A161671(1) = A161671(2).

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

Cf. A141468, A161671, A214627.

Sequence in context: A235714 A146759 A303031 * A144059 A352293 A078676

Adjacent sequences: A220217 A220218 A220219 * A220221 A220222 A220223

Keyword

nonn,less

Author

Juri-Stepan Gerasimov, Dec 07 2012

Status

approved