OFFSET
16,9
LINKS
Alois P. Heinz, Table of n, a(n) for n = 16..20000
Eric Weisstein's World of Mathematics, Ordered Factorization
FORMULA
MAPLE
b:= proc(n) option remember; series(x*(1+add(b(n/d),
d=numtheory[divisors](n) minus {1, n})), x, 5)
end:
a:= n-> coeff(b(n), x, 4):
seq(a(n), n=16..112); # Alois P. Heinz, Feb 22 2021
MATHEMATICA
b[n_] := b[n] = Series[x*(1 + Sum[b[n/d],
{d, Divisors[n] ~Complement~ {1, n}}]), {x, 0, 5}];
a[n_] := Coefficient[b[n], x, 4];
Table[a[n], {n, 16, 112}] (* Jean-François Alcover, Feb 28 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 22 2021
STATUS
approved