OFFSET
1,2
COMMENTS
All terms are zeroless (A052382).
LINKS
David Consiglio, Jr., Table of n, a(n) for n = 1..6253
EXAMPLE
2771 is a term of this sequence because 2*7*7*1 = 98 --> 9 + 8 = 17 --> 2771 / 17 = 163.
MAPLE
d:= n-> convert(n, base, 10):
q:= n-> (m-> m>0 and irem(n, add(j, j=d(m)))=0)(mul(i, i=d(n))):
select(q, [$1..500])[]; # Alois P. Heinz, May 29 2019
MATHEMATICA
Select[Range[500], DigitCount[#, 10, 0]==0&&Divisible[#, Total[ IntegerDigits[ Times@@IntegerDigits[#]]]]&] (* Harvey P. Dale, Jan 24 2021 *)
PROG
(Python)
def dprod(n):
....x = str(n)
....start = 1
....for q in x:
........start *= int(q)
....return start
def dsum(n):
....x = str(n)
....start = 0
....for q in x:
........start += int(q)
....return start
seq_1 = [n for n in range(1, 10000) if dprod(n) != 0 and n % (dsum(dprod(n))) == 0]
print(seq_1)
(PARI) spd(n) = my(d=digits(n)); sumdigits(vecprod(d)); \\ A128212
isok(n) = my(p=spd(n)); p && (n % p == 0); \\ Michel Marcus, May 29 2019
(Magma) [n:n in [1..450]| not 0 in Intseq(n) and IsIntegral(n/(&+Intseq((&*(Intseq(n))))))]; // Marius A. Burtea, May 31 2019
CROSSREFS
KEYWORD
nonn,base
AUTHOR
David Consiglio, Jr., May 29 2019
STATUS
approved