OFFSET
1,1
COMMENTS
All the terms are congruent to 2 or 0 (mod 3).
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 131: 131 + (1 + 3 + 1)^3 = 256 = 16^2.
a(4) = 171: 171 + (1 + 7 + 1)^3 = 900 = 30^2.
MATHEMATICA
Select[Range[50000], IntegerQ[Sqrt[# +Plus@@(IntegerDigits[#])^3]]&]
PROG
(PARI) for(n = 1, 10^5, if(issquare(sumdigits(n)^3 + n), print1(n, ", ")))
(Magma) [ n: n in [1..10^4] | IsSquare(&+Intseq(n)^3 + n) ] ;
CROSSREFS
KEYWORD
nonn,base
AUTHOR
K. D. Bajpai, Jun 28 2015
STATUS
approved