OFFSET
0,3
COMMENTS
Invert transform of A022629.
a(n) is the number of compositions of n where there are A022629(k) sorts of part k. - Joerg Arndt, Jan 24 2024
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
FORMULA
a(0) = 1; a(n) = Sum_{k=1..n} A022629(k)*a(n-k).
MATHEMATICA
nmax = 30; CoefficientList[Series[1/(2 - Product[(1 + k x^k), {k, 1, nmax}]), {x, 0, nmax}], x]
PROG
(Magma)
m:=80;
R<x>:=PowerSeriesRing(Integers(), m);
Coefficients(R!( 1/(2 - (&*[(1+j*x^j): j in [1..m+2]])) ));
(SageMath)
m=80;
def f(x): return 1/( 2 - product(1+j*x^j for j in range(1, m+3)) )
def A307063_list(prec):
P.<x> = PowerSeriesRing(QQ, prec)
return P( f(x) ).list()
A307063_list(m) # G. C. Greubel, Jan 24 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Mar 21 2019
STATUS
approved