OFFSET
0,7
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..2000
FORMULA
a(n) = [x^n] 1 / (1 - Sum_{d|n, d < sqrt(n)} x^d).
MAPLE
a:= proc(n) option remember; uses numtheory; local b, l;
l, b:= select(x-> is(x<sqrt(n)), divisors(n)),
proc(m) option remember; `if`(m=0, 1,
add(`if`(j>m, 0, b(m-j)), j=l))
end; b(n)
end:
seq(a(n), n=0..45); # Alois P. Heinz, Sep 23 2022
MATHEMATICA
a[n_] := SeriesCoefficient[1/(1 - Sum[Boole[d < Sqrt[n]] x^d, {d, Divisors[n]}]), {x, 0, n}]; Table[a[n], {n, 0, 45}]
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Ilya Gutkovskiy, Sep 23 2022
STATUS
approved