login
A078656
a(n) = prime(k) where k = n-th prime congruent to 1 mod 10.
1
31, 127, 179, 283, 353, 547, 739, 877, 1087, 1153, 1297, 1523, 1597, 1741, 1823, 2063, 2221, 2749, 2909, 3001, 3259, 3517, 3733, 3911, 4153, 4421, 4663, 4759, 4943, 5189, 5281, 5701, 5801, 6229, 6311, 6841, 7109
OFFSET
1,1
LINKS
MATHEMATICA
Prime[10Select[Range[110], PrimeQ[10# + 1] &] + 1]
Prime[#]&/@Select[Prime[Range[200]], NumberDigit[#, 0]==1&] (* Requires Mathematica version 12 or later *) (* Harvey P. Dale, May 31 2021 *)
PROG
(PARI) pip(n, m, r) = {sr=0; forprime(x=3, n, if(x%m == r, v=prime(x); sr+=1.0/v; print1(v" "); ) ); print(); print("m="m" r="r" sr="sr); }
CROSSREFS
Cf. A024912.
Sequence in context: A042886 A042888 A183836 * A095322 A127578 A333245
KEYWORD
easy,nonn
AUTHOR
Cino Hilliard, Dec 14 2002
EXTENSIONS
Edited by Robert G. Wilson v, Dec 17 2002
STATUS
approved