OFFSET
1,2
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Harshad Numbers
EXAMPLE
132 is in the sequence because 1 + 3 + 2 = 1*3*2 = 6 and 6 divides 132.
MATHEMATICA
Select[Range[11300000], Divisible[#1, (Plus @@ IntegerDigits[#1])] && (Plus @@ IntegerDigits[#1]) == (Times @@ IntegerDigits[#1]) &]
nQ[n_]:=With[{idn=IntegerDigits[n]}, Mod[n, Total[idn]]==0&&Total[idn]==Times@@idn]; Select[Range[112*10^5], nQ] (* Harvey P. Dale, Nov 02 2024 *)
PROG
(PARI) isok(n) = (d=digits(n)) && ((n % vecsum(d)) == 0) && (vecsum(d) == prod(k=1, #d, d[k])); \\ Michel Marcus, Jan 02 2017
CROSSREFS
KEYWORD
nonn,base,easy
AUTHOR
Ilya Gutkovskiy, Jan 01 2017
STATUS
approved