A023151 Numbers k such that prime(k) == 9 (mod k).

1, 2, 10, 11, 35, 37, 80, 100364, 251711, 251717, 251731, 251735, 251741, 251770, 4124456, 4124582, 27067096, 27067520, 69709706, 69709717, 69709723, 69709868, 69709931, 69709933, 465769825, 465769826, 465769831, 1208198548, 8179002130

Offset: 1

Author

David W. Wilson

Keywords

• nonn

Links

• Giovanni Resta, Table of n, a(n) for n = 1..46

Crossrefs

Cf. A092051, A023143, A023144, A023145, A023146, A023147, A023148, A023149, A023150, A023152.

Programs

Mathematica

NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; p = 1; Do[ If[ Mod[p = NextPrim[p], n] == 9, Print[n]], {n, 1, 10^9}] (* Robert G. Wilson v, Feb 18 2004 *)

Sage

def A023151(max) :
    terms = []
    p = 2
    for n in range(1, max+1) :
        if (p - 9) % n == 0 : terms.append(n)
        p = next_prime(p)
    return terms
# Eric M. Schmidt, Feb 05 2013

Extensions

More terms from Robert G. Wilson v, Feb 18 2004

a(25)-a(29) from Robert G. Wilson v, Feb 22 2006

First two terms inserted by Eric M. Schmidt, Feb 05 2013