OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
257 is in the sequence because (1) it is prime and (2) the sum of its digits (2 + 5 + 7 = 14) plus the sum of the digits of 257^2 (66049 -> 6 + 6 + 0 + 4 + 9 = 25) plus the sum of the digits of 257^3 (16974593 -> 1 + 6 + 9 + 7 + 4 + 5 + 9 + 3 = 44) equals 14 + 25 + 44 = 83, which is prime.
MAPLE
filter:= n -> isprime(n) and isprime(convert(convert(n, base, 10), `+`)+convert(convert(n^2, base, 10), `+`)+convert(convert(n^3, base, 10), `+`)):
select(filter, [seq(i, i=5..10000, 6)]); # Robert Israel, Jul 28 2020
MATHEMATICA
Select[Prime[Range[300]], PrimeQ[Total[IntegerDigits[#]] + Total[IntegerDigits[#^2]] + Total[IntegerDigits[#^3]]]&]
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Harvey P. Dale, Feb 16 2011
STATUS
approved