OFFSET
1,3
COMMENTS
Let sod(n) := digital sum of n (A007953); here we have sod(n^2) = sod(n^4).
Trivial cases:
(I) Powers of 10, as sod((10^k)^2) = sod((10^k)^4) = 1.
(II) If N is a term of sequence, then so is 10 * N.
REFERENCES
Hans Schubart, Einfuehrung in die klassische und moderne Zahlentheorie, Vieweg, Braunschweig, 1974.
EXAMPLE
sod(3^2) = sod(9) = 9 = sod(81) = sod(3^4), so 3 is a term.
sod(17^2) = sod(289) = 19 = sod(83521) = sod(17^4), so 17 is a term.
MATHEMATICA
Select[Range[0, 2000], Total[IntegerDigits[#^2]]==Total[IntegerDigits[#^4]]&] (* Harvey P. Dale, Jan 19 2011 *)
CROSSREFS
KEYWORD
base,nonn,changed
AUTHOR
Ulrich Krug (leuchtfeuer37(AT)gmx.de), Apr 25 2010
EXTENSIONS
Edited by D. S. McNeil, Nov 21 2010
a(43)-a(51) from Jason Yuen, Oct 13 2024
STATUS
approved