login
Zuckerman numbers which divided by the product of their digits give integers which are also divisible by the product of their digits, and so on, until result is 1.
0

%I #48 Mar 31 2023 09:17:53

%S 1,2,3,4,5,6,7,8,9,12,15,24,36,128,135,144,175,384,672,735,1296,1575,

%T 82944,139968,1492992,27869184

%N Zuckerman numbers which divided by the product of their digits give integers which are also divisible by the product of their digits, and so on, until result is 1.

%C Repunits >= 11 (A002275) are not in the sequence because, as they are fixed points of this map, they don't fit the definition.

%C Question: is this sequence finite as the similar sequence with Niven numbers (A114440) that has 15095 terms?

%C No other terms up to 2*10^9. - _Michel Marcus_, Apr 27 2021

%C From _David A. Corneth_, Apr 27 2021: (Start)

%C Terms are 7-smooth. Any prime factor > 7 will not be divided away by dividing by product of digits.

%C Any number k > a(26)*10^163 with product of digits vp > 0 has k/vp > a(26) so it suffices to check all candidates <= a(26)*10^163. Doing so gives no more terms so this sequence is finite and full. (End)

%C The number of steps needed to reach 1, has a maximum of 3, which occurs for n = 21, 23..26. - _A.H.M. Smeets_, Apr 29 2021

%H Giovanni Resta, <a href="https://www.numbersaplenty.com/set/Zuckerman_number/">Zuckerman numbers</a>, Numbers Aplenty.

%e The integer 1296 is divisible by the product of its digits as 1296/(1*2*9*6) = 12, then 12/(1*2) = 6 and 6/6 = 1; hence, 1296 is a term of this sequence.

%t f[n_] := If[(prod = Times @@ IntegerDigits[n]) > 0 && Divisible[n, prod], n/prod, 0]; Select[Range[10^5], FixedPointList[f, #][[-1]] == 1 &] (* _Amiram Eldar_, Apr 27 2021 *)

%o (PARI) isz(n) = my(p=vecprod(digits(n))); p && !(n % p); \\ A007602

%o isok(n) = if (n==1, return(1)); my(m=n); until(m==1, if (isz(m), my(nm = m/vecprod(digits(m))); if (nm==m, return (0), m = nm), return(0))); return(1); \\ _Michel Marcus_, Apr 27 2021

%o (Python)

%o def proddigit(n):

%o p = 1

%o while n > 0:

%o n, p = n//10, p*(n%10)

%o return p

%o n, a = 1, 1

%o while n > 0:

%o aa, pa = a, proddigit(a)

%o while pa > 1 and aa%pa == 0 and aa > 1:

%o aa = aa//pa

%o pa = proddigit(aa)

%o if aa == 1:

%o print(n,a)

%o n = n+1

%o a = a+1 # _A.H.M. Smeets_, Apr 29 2021

%Y Cf. A007602, A288069.

%Y Cf. A114440 (similar for Harshad numbers).

%Y Subsequence of A002473 and of A343681.

%K nonn,base,fini,full

%O 1,2

%A _Bernard Schott_, Apr 27 2021

%E a(26) from _Michel Marcus_, Apr 27 2021