OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
PROG
(PARI) SumD(x)= { local(s=0); while (x>9, s+=x%10; x\=10); return(s + x) } { n=0; for (m=1, 10^9, if (isprime(SumD(m^m)), write("b066233.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, Feb 07 2010
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Jason Earls, Dec 19 2001
STATUS
approved