OFFSET
1,2
COMMENTS
Starting with a(5), positions of 13 in A004648. - corrected by Eric M. Schmidt, Feb 05 2013
MATHEMATICA
NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; p = 1; Do[ If[ Mod[p = NextPrim[p], n] == 13, Print[n]], {n, 10^9}] (* Robert G. Wilson v, Feb 22 2006 *)
PROG
(Sage)
def A116658(max) :
terms = []
p = 2
for n in range(1, max+1) :
if (p - 13) % n == 0 : terms.append(n)
p = next_prime(p)
return terms
# Eric M. Schmidt, Feb 05 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Feb 21 2006
EXTENSIONS
a(24)-a(34) from Robert G. Wilson v, Feb 22 2006
First four terms inserted by Eric M. Schmidt, Feb 05 2013
STATUS
approved