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

2, 809, 821, 827, 857, 881, 1019, 1031, 1049, 1061, 1091, 1151, 1229, 1277, 1289, 1301, 1319, 1427, 1451, 1481, 1487, 1607, 1619, 1667, 1697, 1721, 1787, 1871, 1877, 1931, 1949, 1997, 2027, 2081, 2087, 2111, 2129, 2141, 2237, 2267, 2309, 2339, 2381

Offset: 1

Author

Roger L. Bagula, Jun 27 2005

Keywords

• nonn
• less

Crossrefs

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]] == 0, Prime[n], {}], {n, 1, 400}]]

Extensions

Offset changed and definition amended by Georg Fischer, Apr 06 2022