login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A130828
Primes p such that the sum of the digits of p^p is a prime.
0
5, 11, 19, 29, 37, 43, 89, 97, 113, 139, 269, 311, 337, 359, 367, 433
OFFSET
1,1
COMMENTS
Computed by Emeric Deutsch.
EXAMPLE
For 5^5 = 625, 6 + 2 + 5 = 13, which is a prime.
MAPLE
sd := proc (n) options operator, arrow: add(convert(n, base, 10)[j], j = 1 .. nops(convert(n, base, 10))) end proc:
a := proc (n) if isprime(sd(ithprime(n)^ithprime(n))) = true then ithprime(n) else end if end proc:
seq(a(n), n = 1 .. 90); # Emeric Deutsch, Jul 19 2007
CROSSREFS
Sequence in context: A038872 A141158 A239732 * A244241 A108151 A209188
KEYWORD
nonn,base,less
AUTHOR
J. M. Bergot, Jul 17 2007
STATUS
approved