OFFSET
1,2
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..46
EXAMPLE
204475053103 = prime(8179002124) and 204475053103 = 25*8179002124 + 3.
MATHEMATICA
NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; p = 1; Do[ If[ Mod[p = NextPrim[p], n] == 3, Print[n]], {n, 1, 10^9}] (* Robert G. Wilson v, Feb 18 2004 *)
Select[Range[100000], Mod[Prime[#] - 3, #] == 0 &] (* T. D. Noe, Feb 05 2013 *)
PROG
(Sage)
def A023145(max) :
terms = []
p = 2
for n in range(1, max+1) :
if (p - 3) % n == 0 : terms.append(n)
p = next_prime(p)
return terms
# Eric M. Schmidt, Feb 05 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Robert G. Wilson v, Feb 18 2004
2 more terms from Giovanni Resta, Feb 22 2006
a(29) from Robert G. Wilson v, Feb 22 2006
First two terms inserted by Eric M. Schmidt, Feb 05 2013
Terms a(30) and beyond from Giovanni Resta, Feb 23 2020
STATUS
approved