A109564 prime(k) for those k where floor((2*(prime(k+1)-prime(k))*PrimePi(k) mod (8*k))/k) = m with m = 10.
1669, 2477, 2971, 3137, 3271, 4297, 4831, 5119, 5351, 5749, 6491, 6917, 7253, 7759, 7963, 8389, 8893, 10799, 11743, 12163, 17257, 18803, 19087, 20443, 21433, 22193, 22307, 22817, 24281, 27143, 28351, 29881, 30593, 32261, 38393, 39461, 45013, 45779, 48907, 49559
Offset: 1
Programs
-
Mathematica
a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 10, Prime[n], {}], {n, 1, 2000}]]
Extensions
Name edited, offset corrected and more terms added by Amiram Eldar, Jun 02 2025