A109567 prime(k) for those k where floor((2*(prime(k+1)-prime(k))*PrimePi(k) mod (8*k))/k) = m with m = 13.

15683, 28229, 45893, 48679, 59281, 60539, 74959, 86869, 100853, 106033, 107509, 114493, 117443, 118973, 121789, 153001, 153191, 160253, 178697, 181303, 189067, 190921, 191353, 195599, 203057, 204173, 214993, 220063, 221101, 227303, 231719, 233777, 237409, 243709

Offset: 1

Author

Roger L. Bagula, Jun 27 2005

Keywords

• nonn
• less

Crossrefs

Cf. A000040, A000720, A001359.

Cf. A109556-A109569 for m = 2..15.

Programs

Mathematica

a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 13, Prime[n], {}], {n, 1, 10000}]]

Extensions

Name edited, offset corrected and more terms added by Amiram Eldar, Jun 02 2025