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

113, 523, 887, 1129, 1951, 2179, 5531, 5953, 8971, 10009, 10399, 10531, 10909, 13063, 13187, 13933, 13967, 14251, 14983, 16381, 16573, 17627, 18553, 18869, 27701, 27851, 29683, 32653, 34549, 37747, 40387, 41299, 42863, 45083, 45197, 46771, 46957, 47743, 47981, 49957

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

Extensions

Definition amended, offset changed and more terms from Georg Fischer, Apr 06 2022