A058905 - OEIS

%I #13 May 10 2023 12:18:15

%S 46,47,48,56,58,66,76,86,136,138,167,176,222,227,228,248,258,298,302,

%T 308,312,316,318,338,343,344,347,348,352,354,356,358,362,374,383,384,

%U 392,398,402,403,404,406,407,408,411,412,414,416,422,423

%N Inconsummate numbers in base 9: no number is this multiple of the sum of its digits (in base 9).

%p For Maple code see A058906.

%o (Python)

%o from itertools import count, islice, combinations_with_replacement

%o from sympy.ntheory import digits

%o def A058905_gen(startvalue=1): # generator of terms >= startvalue

%o     for n in count(max(startvalue,1)):

%o         for l in count(1):

%o             if l*n<<3 < 9**(l-1):

%o                 yield n

%o                 break

%o             for d in combinations_with_replacement(range(9),l):

%o                 if (s:=sum(d)) > 0 and sorted(digits(s*n,9)[1:]) == list(d):

%o                     break

%o             else:

%o                 continue

%o             break

%o A058905_list = list(islice(A058905_gen(),20)) # _Chai Wah Wu_, May 10 2023

%Y Cf. A003635, A052491, A058898-A058907.

%K nonn,base

%O 1,1

%A _N. J. A. Sloane_, Jan 09 2001