OFFSET
1,2
COMMENTS
In other words, for n > 1, a(n) = smallest k not in a(1)..a(n-1) such that the partial sums of the number of divisors of a(j) divides the partial sum of a(j), 1 <= j < n.
Conjecture: 3 does not appear in the sequence.
LINKS
Michael De Vlieger, Log log scatterplot of a(n), n = 1..40000.
Michael De Vlieger, Annotated log log scatterplot of a(n), n = 1..5000, color coded to show primes in red, composite prime powers in gold, squarefree composites in green, products of composite prime powers in large light blue, and other numbers in dark blue.
EXAMPLE
a(2) = 2 since t(1) + d(2) | s(1) + 2 = 1+2 | 1+2 = 3 | 3,
a(3) = 7 since t(2) + d(7) | s(2) + 7 = 3+2 | 3+7 = 5 | 10,
a(4) = 8 since t(3) + d(8) | s(3) + 8 = 5+4 | 10+8 = 9 | 18,
a(5) = 12 since t(4) + d(12) | s(4) + 12 = 9+6 | 18+12 = 15 | 30,
a(6) = 27 since t(5) + d(27) | s(5) + 27 = 15+4 | 30+27 = 19 | 57, etc.
MATHEMATICA
nn = 2^10; c[_] = False; a[1] = s = t = 1; c[1] = True; u = 2; Do[k = u; While[Nand[! c[k], Divisible[(s + k), (t + Set[d, DivisorSigma[0, k]])]], k++]; Set[{a[n], c[k]}, {k, True}]; s += k; t += d; If[k == u, While[c[u], u++]], {n, 2, nn}]; Array[a, nn]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger and David James Sycamore, Dec 07 2022
STATUS
approved