OFFSET
1,1
COMMENTS
All the terms are congruent to 2 (mod 3).
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..2610
EXAMPLE
a(2) = 131 is prime: 131 + (1 + 3 + 1)^3 = 256 = 16^2.
a(3) = 863 is prime: 863 + (8 + 6 + 3)^3 = 5776 = 76^2.
MATHEMATICA
Select[Prime[Range[100000]], IntegerQ[Sqrt[# + Plus @@ (IntegerDigits[#])^3]] &]
PROG
(PARI) forprime(p=1, 10^6, if(issquare(sumdigits(p)^3 + p), print1(p, ", ")))
(Magma) [p: p in PrimesUpTo(2*10^6) | IsSquare(&+Intseq(p)^3 + p)] ;
CROSSREFS
KEYWORD
nonn,base
AUTHOR
K. D. Bajpai, Jun 26 2015
STATUS
approved