OFFSET
1,1
COMMENTS
a(6) > 10^5. - Michael S. Branicky, Sep 28 2024
MATHEMATICA
f[n_] := FromDigits[Take[IntegerDigits[n^n], -n]]; Do[If[PrimeQ[f[n]], Print[n]], {n, 8, 14949}] (* Ryan Propper, Apr 01 2007 *)
PROG
(Python)
from sympy import isprime
A126703_list = [n for n in range(1, 2000) if isprime(pow(n, n, 10**n))] # Chai Wah Wu, Mar 23 2018
CROSSREFS
KEYWORD
base,hard,more,nonn
AUTHOR
Shyam Sunder Gupta, Feb 11 2007
EXTENSIONS
a(5) from Ryan Propper, Apr 01 2007
STATUS
approved