OFFSET
0,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..200
PROG
(PARI)
InvEulerT(v)={dirdiv(Vec(log(1+x*Ser(v)), -#v), vector(#v, n, 1/n))}
b(k, r)={if(k%2, if(r%2, 0, my(j=r/2); k^j*(2*j)!/(j!*2^j)), sum(j=0, r\2, binomial(r, 2*j)*k^j*(2*j)!/(j!*2^j)))}
C(k, r)={sum(i=0, r, (-1)^i/i!/k^i)}
S(n, k)={sum(r=0, 2*n\k, if(k*r%2==0, x^(k*r/2)*b(k, r)*C(k, r)), O(x*x^n))}
seq(n)={concat([1], InvEulerT(Vec(-1 + prod(k=1, 2*n, S(n, k)))))}
CROSSREFS
KEYWORD
nonn,new
AUTHOR
Andrew Howroyd, Jan 28 2025
STATUS
approved