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

37, 43, 67, 151, 157, 167, 173, 233, 251, 257, 263, 271, 331, 353, 367, 373, 383, 433, 443, 503, 541, 557, 563, 571, 587, 593, 601, 607, 647, 653, 683, 701, 719, 743, 761, 911, 929, 983, 1109, 1163, 1193, 1373, 1439, 1523, 1559, 1571, 1733, 1823, 1979

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

Extensions

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