OFFSET
1,2
COMMENTS
See A354606 for further details.
PROG
(Python)
from sympy import divisor_count
from collections import Counter
from itertools import count, islice
def f(n): return divisor_count(n)
def agen():
n, an, fan, inventory = 1, 1, 1, Counter([1])
yield n
for n in count(2):
an = inventory[fan]
fan = f(an)
inventory.update([fan])
if an == 1: yield n
print(list(islice(agen(), 50))) # Michael S. Branicky, Jul 09 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Scott R. Shannon, Jul 09 2022
STATUS
approved