OFFSET
0,2
COMMENTS
Partial sums of A022629.
LINKS
FORMULA
G.f.: (1/(1 - x))*exp(Sum_{j>=1} Sum_{k>=1} (-1)^(j+1)*k^j*x^(j*k)/j).
MAPLE
b:= proc(n, i) option remember; `if`(i*(i+1)/2<n, 0,
`if`(n=0, 1, b(n, i-1)+`if`(i>n, 0, i*b(n-i, i-1))))
end:
a:= proc(n) option remember; `if`(n<0, 0, a(n-1)+b(n$2)) end:
seq(a(n), n=0..40); # Alois P. Heinz, Apr 13 2018
MATHEMATICA
nmax = 40; CoefficientList[Series[1/(1 - x) Product[(1 + k x^k), {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 40; CoefficientList[Series[1/(1 - x) Exp[Sum[Sum[(-1)^(j + 1) k^j x^(j k)/j, {k, 1, nmax}], {j, 1, nmax}]], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 13 2018
STATUS
approved