OFFSET
1,2
COMMENTS
a(12) > 2*10^10. - Giovanni Resta, May 23 2013
MATHEMATICA
Select[Range[10^6], (p = Prime[#]; PowerMod[#, #, p] == p - 1) &] (* Giovanni Resta, May 23 2013 *)
PROG
(Python)
from sympy import nextprime, prime
from itertools import count, islice
def agen(startn=1): # generator of terms
pn = prime(startn)
for n in count(startn):
if pow(n, n, pn) == pn - 1:
yield n
pn = nextprime(pn)
print(list(islice(agen(), 5))) # Michael S. Branicky, May 25 2023
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Alex Ratushnyak, May 21 2013
EXTENSIONS
a(6)-a(11) from Giovanni Resta, May 23 2013
a(12) from Michael S. Branicky, May 25 2023
STATUS
approved