|
|
A058904
|
|
Inconsummate numbers in base 8: no number is this multiple of the sum of its digits (in base 8).
|
|
1
|
|
|
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
|
|
|
MAPLE
|
|
|
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
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|