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.

A383746 Numbers k such that k divides the sum of the digits of k^(3k).

This page as a plain text file.
%I A383746 #19 May 13 2025 16:42:53
%S A383746 1,2,3,6,9,11,18,38,43,87,126,670,1098,2421,3588,4201,5114,5877,5922,
%T A383746 6048,11799,46119,46419,55098,55945,77439,91541,129624,153229,182402
%N A383746 Numbers k such that k divides the sum of the digits of k^(3k).
%e A383746 2 is a term since the sum of digits of 2^(3*2) is 64, which is divisible by 2.
%e A383746 3 is a term since the sum of digits of 3^(3*3) is 19683, which is divisible by 3.
%e A383746 1098 is a term since the sum of digits of 1098^(3*1098) is 45018, which is divisible by 1098.
%t A383746 Do[If[Mod[Plus @@ IntegerDigits[n^(3*n)], n] == 0, Print[n]], {n, 1, 10000}]
%o A383746 (Python)
%o A383746 from gmpy2 import digits, mpz
%o A383746 def ok(n): return n and sum(map(mpz, digits(n**(3*n))))%n == 0
%o A383746 print([k for k in range(1100) if ok(k)]) # _Michael S. Branicky_, May 08 2025
%Y A383746 Cf. A083282, A108859.
%K A383746 nonn,base,hard,more
%O A383746 1,2
%A A383746 _J.W.L. (Jan) Eerland_, May 08 2025
%E A383746 a(21)-a(23) and a(28)-a(30) from _Michael S. Branicky_, May 08 2025