A058904 Inconsummate numbers in base 8: no number is this multiple of the sum of its digits (in base 8).

42, 44, 51, 52, 60, 105, 109, 116, 124, 173, 177, 178, 181, 201, 205, 209, 210, 213, 214, 217, 233, 237, 241, 242, 245, 249, 250, 251, 254, 255, 269, 273, 277, 278, 282, 285, 287, 290, 298, 299, 300, 308, 336, 343, 348, 352, 397, 401, 402, 403

Offset: 1

Author

N. J. A. Sloane, Jan 09 2001

Keywords

• nonn
• base

Crossrefs

Cf. A003635, A052491, A058898-A058907.

Programs

Maple

For Maple code see A058906.

Mathematica

base=8; Do[k=n; While[Apply[Plus, IntegerDigits[k, base]] n!=k&&k<250n, k+=n]; If[k==250 n, Print[n]], {n, 1, 10^3}] (* Vincenzo Librandi, Sep 21 2017 *)

Python

from itertools import count, islice, combinations_with_replacement
def A058904_gen(startvalue=1): # generator of terms
    for n in count(max(startvalue,1)):
        for l in count(1):
            if 7*l*n < 1<<3*(l-1):
                yield n
                break
            for d in combinations_with_replacement(range(8),l):
                if (s:=sum(d)) > 0 and sorted(oct(s*n)[2:]) == list(map(str,d)):
                    break
            else:
                continue
            break
A058904_list = list(islice(A058904_gen(),20)) # Chai Wah Wu, May 09 2023