A109558 a(n) = If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 4, Prime[n]]

7, 13, 19, 47, 53, 61, 73, 83, 131, 359, 389, 401, 449, 479, 491, 691, 709, 787, 811, 829, 919, 1021, 1039, 1051, 1153, 1171, 1249, 1399, 1471, 1627, 1699, 1723, 1801, 1879, 2017, 2029, 2053, 2069, 2089, 2099, 2143, 2297, 2399, 2447, 2579, 2621

Offset

0,1

Links

Table of n, a(n) for n=0..45.

Mathematica

a = Flatten[Table[If[Floor[2*Mod[(Prime[n + 1] - Prime[n])*PrimePi[n]/n, 8]] == 4, Prime[n], {}], {n, 1, 400}]]

Crossrefs

Cf. A001359.

Sequence in context: A181938 A073648 A298737 * A108106 A231506 A165724

Adjacent sequences: A109555 A109556 A109557 * A109559 A109560 A109561

Keyword

nonn

Author

Roger L. Bagula, Jun 27 2005

Status

approved