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,1

Links

Table of n, a(n) for n=1..50.

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 *)

Prog

(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

Crossrefs

Cf. A003635, A052491, A058898-A058907.

Sequence in context: A133897 A086848 A364219 * A181459 A169907 A095493

Adjacent sequences: A058901 A058902 A058903 * A058905 A058906 A058907

Keyword

nonn,base

Author

N. J. A. Sloane, Jan 09 2001

Status

approved