This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A357080 #14 Sep 26 2022 20:35:20 %S A357080 0,1,80,162,243,476,486 %N A357080 Numbers k such that the sum of the digits of k multiplied by the sum of the digits of k^2 equals k. %C A357080 Suppose k has m digits, then the sum of the digits of k multiplied by the sum of the digits of k^2 is bounded by 9m times 9*(2m), which equals 162m^2. On the other hand, k is greater than 10^(m-1), which grows much faster than 162m^2. It follows that k can't have more than 4 digits. %e A357080 The sum of the digits of 80 is 8, the sum of the digits of 80^2 = 6400 is 10. The number 80 itself is 8*10. Thus, 80 is in this sequence. %t A357080 Select[Range[100000], # == Total[IntegerDigits[#]] Total[IntegerDigits[#^2]] &] %o A357080 (PARI) isok(k) = k == sumdigits(k)*sumdigits(k^2); \\ _Michel Marcus_, Sep 11 2022 %o A357080 (Python) %o A357080 def sd(n): return sum(map(int, str(n))) %o A357080 def ok(n): return sd(n) * sd(n*n) == n %o A357080 print([k for k in range(10**5) if ok(k)]) # _Michael S. Branicky_, Sep 11 2022 %Y A357080 Cf. A007953, A004159, A130181. %K A357080 nonn,base,fini,full %O A357080 1,3 %A A357080 _Tanya Khovanova_, Sep 10 2022