%I #23 May 06 2022 13:13:51
%S 1,2,3,4,5,6,7,8,9,11,12,15,24,36,111,128,135,144,175,384,672,735,
%T 1111,1296,1575,11111,22176,42624,82944,111111,139968,688128,719712,
%U 1111111,1161216,1492992,2241792,2794176,4136832,4741632,6838272,11111111,12171264,13395375,13436928
%N Zuckerman numbers which when divided by product of their digits, give a quotient which is also a Zuckerman number.
%C Alternative Name: Zuckerman numbers k such that k/(product of digits of k) is also a Zuckerman number. - _Wesley Ivan Hurt_, Apr 26 2021
%C All positive repunits are terms (A002275).
%H Giovanni Resta, <a href="/A343681/b343681.txt">Table of n, a(n) for n = 1..160</a> (terms < 3*10^14)
%H Giovanni Resta, <a href="https://www.numbersaplenty.com/set/Zuckerman_number/">Zuckerman numbers</a>, Numbers Aplenty.
%e 24 is a Zuckerman number as 24/(2*4) = 3, 3/3 = 1 so 3 is also a Zuckerman number, and 24 is a term.
%e 1296 is a Zuckerman number as 1296/(1*2*9*6) = 12, 12/(1*2) = 4 so 12 is also a Zuckerman number and 1296 is a term.
%t zuckQ[n_] := (prod = Times @@ IntegerDigits[n]) > 0 && Divisible[n, prod]; Select[Range[10^6], zuckQ[#] && zuckQ[#/Times @@ IntegerDigits[#]] &] (* _Amiram Eldar_, Apr 26 2021 *)
%o (PARI) isz(n) = my(p=vecprod(digits(n))); p && !(n % p); \\ A007602
%o isok(n) = isz(n) && isz(n/vecprod(digits(n))); \\ _Michel Marcus_, Apr 26 2021
%Y Cf. A002275, A007602, A288069, A343680, A343682.
%Y Cf. A235507 (similar, with Niven numbers).
%K nonn,base
%O 1,2
%A _Bernard Schott_, Apr 26 2021
%E More terms from _David A. Corneth_, Apr 26 2021