OFFSET
1,2
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..71
MATHEMATICA
p = 1; Do[ If[ Mod[p = NextPrime[p], n] == 7, Print[n]], {n, 1, 10^9}] (* Robert G. Wilson v, Feb 18 2004 *)
PROG
(Sage)
def A023149(max) :
terms = []
p = 2
for n in range(1, max+1) :
if (p - 7) % 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
a(27)-a(31) from Robert G. Wilson v, Feb 22 2006
First four terms inserted by Eric M. Schmidt, Feb 05 2013
STATUS
approved