A089208 Duplicate of A083282.
1, 64, 19683, 16777216, 30517578125, 101559956668416, 558545864083284007
Offset: 1
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.
Triangle begins: 1; 1, 4; 1, 8, 27; 1, 16, 81, 256; 1, 32, 243, 1024, 3125; 1, 64, 729, 4096, 15625, 46656; ...
a089072 = flip (^) a089072_row n = map (a089072 n) [1..n] a089072_tabl = map a089072_row [1..] -- Reinhard Zumkeller, Mar 18 2013
[k^n: k in [1..n], n in [1..12]]; // G. C. Greubel, Nov 01 2022
Column[Table[k^n, {n, 8}, {k, n}], Center] (* Alonso del Arte, Nov 14 2011 *)
flatten([[k^n for k in range(1,n+1)] for n in range(1,12)]) # G. C. Greubel, Nov 01 2022
2 is a term since the sum of digits of 2^(3*2) is 64, which is divisible by 2. 3 is a term since the sum of digits of 3^(3*3) is 19683, which is divisible by 3. 1098 is a term since the sum of digits of 1098^(3*1098) is 45018, which is divisible by 1098.
Do[If[Mod[Plus @@ IntegerDigits[n^(3*n)], n] == 0, Print[n]], {n, 1, 10000}]
from gmpy2 import digits, mpz def ok(n): return n and sum(map(mpz, digits(n**(3*n))))%n == 0 print([k for k in range(1100) if ok(k)]) # Michael S. Branicky, May 08 2025
Comments