OFFSET
1,1
EXAMPLE
a(2)=31 is prime: [3^3 + 1^3 = 27 + 1] = 28; [31^3 = 29791, 2+9+7+9+1] = 28.
a(4)=1321 is prime: [1^3 + 3^3 + 2^3 + 1^3 = 1 + 27 + 8 + 1] = 37; [31^3 = 2305199161, 2+3+0+5+1+9+9+1+6+1] = 37.
MATHEMATICA
Select[Prime[Range[30000]], Total[IntegerDigits[#]^3] == Plus @@ IntegerDigits[#^3] &]
PROG
(PARI) forprime(p=1, 30000, d=digits(p); if(sum(i=1, length(d), d[i]^3) == sumdigits(p^3), print1(p", ")));
CROSSREFS
KEYWORD
nonn,base
AUTHOR
K. D. Bajpai, Aug 17 2017
STATUS
approved