OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..120
MATHEMATICA
a[n_] := Plus @@ DigitCount[Divisors[n], 2, 1]; am = -1; c = 0; seq={}; Do[a1 = a[n]; If[a1 > am, am = a1; c++; AppendTo[seq, n]], {n, 1, 10^4}]; seq (* Amiram Eldar, Dec 17 2019 *)
PROG
(Python) # uses imports and definitions in A093653
from itertools import count, islice
def f(n): return A093653(n)
def agen(r=0): yield from ((m, r:=fm)[0] for m in count(1) if (fm:=f(m))>r)
print(list(islice(agen(), 45))) # Michael S. Branicky, Feb 15 2023
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jason Earls, May 16 2004
STATUS
approved
