OFFSET
1,1
COMMENTS
This is to the exponent 4 as A182404 is to the exponent 2.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
21 is in the sequence because sum of digits 2+1= 3 is prime, and sum of the 4th powers of the digits 2^4+1^4=17 is a prime.
MATHEMATICA
Select[Range[350], AllTrue[{Total[IntegerDigits[#]], Total[ IntegerDigits[ #]^4]}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Apr 01 2019 *)
PROG
(PARI) dspow(n, b, k)=my(s); while(n, s+=(n%b)^k; n\=b); s
select(n->isprime(sumdigits(n))&&isprime(dspow(n, 10, 4)), vector(10^3, i, i)) \\ Charles R Greathouse IV, May 11 2012
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Jonathan Vos Post, May 10 2012
STATUS
approved