A350741 Records in A095258.

1, 3, 4, 6, 9, 27, 34, 68, 94, 235, 289, 578, 799, 1921, 9683, 16021, 27421, 54842, 69301, 138602, 434789, 1787371, 5179771, 5655149, 9653251, 10209853, 20419706, 43184409, 301039141, 611363527, 1274384647, 5084853899, 14906805553, 14946637163, 22591381313, 69291164983

Offset: 1

Author

Michael De Vlieger, Jan 23 2022

Keywords

• nonn

Links

• Michael De Vlieger, Table of n, a(n) for n = 1..161

Crossrefs

Cf. A095258.

Programs

Mathematica

c[_] = 0; j = c[1] = r = 1; s = 3; Prepend[Reap[Do[d = Divisors[s]; k = 1; While[c[d[[k]]] > 0, k++]; Set[k, d[[k]]]; Set[c[k], i]; If[k > r, r = k; Sow[r]]; j = k; s += k, {i, 2, 2100}] ][[-1, -1]], 1]

Python

from itertools import islice
from sympy import divisors
def A350741_gen(): # generator of terms
    bset, c, s = {1}, 1, 3
    yield 1
    while True:
        for d in divisors(s):
            if d not in bset:
                if d > c:
                    yield d
                    c = d
                bset.add(d)
                s += d
                break
A350741_list = list(islice(A350741_gen(),20)) # Chai Wah Wu, Jan 25 2022