OFFSET
1,1
LINKS
Daniel Mondot, Table of n, a(n) for n = 1..13275
MAPLE
For Maple code see A058906.
MATHEMATICA
base=3; Do[k=n; While[Apply[Plus, IntegerDigits[k, base]] n!=k&&k<250 n, k+=n]; If[k==250 n, Print[n]], {n, 1, 10^3}] (* Vincenzo Librandi, Sep 23 2017 *)
PROG
(Python)
from itertools import count, islice, combinations_with_replacement
from sympy.ntheory import digits
def A058899_gen(startvalue=1): # generator of terms >= startvalue
for n in count(max(startvalue, 1)):
for l in count(1):
if 2*l*n < 3**(l-1):
yield n
break
for d in combinations_with_replacement((0, 1, 2), l):
if (s:=sum(d)) > 0 and sorted(digits(s*n, 3)[1:]) == list(d):
break
else:
continue
break
CROSSREFS
KEYWORD
nonn,base
AUTHOR
N. J. A. Sloane, Jan 09 2001
STATUS
approved