OFFSET
1,2
COMMENTS
If k is a term, then 10 * k is a term. There are an infinite number of terms that are not divisible by 10. The numbers m = 24 * 10^(294 * k - 292) +1 are divisible by 7^3 = digsum(m)^3. Also, the numbers s = 491 * 10^(4624 * k - 4623) + 3, k >= 1, are divisible by 17^3 = digsum(s)^3. - Marius A. Burtea, Mar 18 2020
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
EXAMPLE
k=98415: sumdigits(98415)=27, q=98415=5*27*27*27.
MATHEMATICA
sud[x_] := Apply[Plus, IntegerDigits[x]] Do[s=sud[n]^3; If[IntegerQ[n/s], Print[n]], {n, 1, 10000}]
Select[Range[50000], Divisible[#, Total[IntegerDigits[#]]^3]&] (* Harvey P. Dale, Mar 22 2016 *)
PROG
(Magma) [k:k in [1..41000]| k mod &+Intseq(k)^3 eq 0]; // Marius A. Burtea, Mar 18 2020
(PARI) is(n)=n%sumdigits(n)^3==0 \\ Charles R Greathouse IV, Mar 19 2020
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Labos Elemer, Jun 14 2002
STATUS
approved