OFFSET
1,3
LINKS
R. J. Mathar, Table of n, a(n) for n = 1..399
MAPLE
read("transforms") :
n := 1;
for i from 1 do
p := A002113(i) ;
if digsum(p) = digsum(p^2) then
printf("%d %d\n", n, p) ;
n := n+1 ;
end if;
end do: # R. J. Mathar, Sep 09 2015
MATHEMATICA
Select[Range[0, 78*10^4], PalindromeQ[#]&&Total[IntegerDigits[#]] == Total[ IntegerDigits[ #^2]]&] (* Harvey P. Dale, Sep 26 2021 *)
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Patrick De Geest, Dec 15 2000
STATUS
approved