A058904 - OEIS

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

(list; graph; refs; listen; history; text; internal format)

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