A109569 prime(k) for those k where floor((2*(prime(k+1)-prime(k))*PrimePi(k) mod (8*k))/k) = m with m = 15.
19609, 25471, 35617, 44293, 58831, 89689, 107377, 162143, 218287, 253159, 290249, 302329, 360091, 404597, 507217, 514967, 517639, 531383, 570253, 619397, 621143, 639757, 642281, 673669, 716173, 736279, 794249, 795349, 815729, 873787, 912649, 933073, 937253, 1013063
Offset: 1
Programs
-
Mathematica
a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 15, Prime[n], {}], {n, 1, 10000}]]
Extensions
Name edited, offset corrected and more terms added by Amiram Eldar, Jun 02 2025