cp's OEIS Frontend

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.

A357080 Numbers k such that the sum of the digits of k multiplied by the sum of the digits of k^2 equals k.

This page as a plain text file.
%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