OFFSET
1,2
MATHEMATICA
NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; p = 1; Do[ If[ Mod[p = NextPrim[p], n] == 10, Print[n]], {n, 1, 10^9}] (* Robert G. Wilson v, Feb 18 2004 *)
PROG
(Sage)
def A023152(max) :
terms = []
p = 2
for n in range(1, max+1) :
if (p - 10) % n == 0 : terms.append(n)
p = next_prime(p)
return terms
# Eric M. Schmidt, Feb 05 2013
(PARI) n=0; forprime(p=2, 1e9, if(Mod(p, n++)==10, print1(n", "))) \\ Charles R Greathouse IV, Apr 29 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Extended by Robert G. Wilson v, Feb 18 2004
a(9)-a(14) from Robert G. Wilson v, Feb 22 2006
First two terms inserted by Eric M. Schmidt, Feb 05 2013
a(15)-a(25) from Giovanni Resta, Feb 23 2020
STATUS
approved