OFFSET
0,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..200
PROG
(PARI) \\ here U(n, k) is k-th column of A339650 as a vector.
R(n, k)={my(v=vector(n)); v[1]=k; for(n=2, n, v[n]=sum(j=1, (n-1)\2, v[j]*v[n-j]) + if(n%2, 0, binomial(v[n/2]+1, 2))); v}
U(n, k)={my(g=x*Ser(R(n, k))); Vec(1 + g + (subst(g + O(x*x^(n\3)), x, x^3) - g^3)/3) }
seq(n)={sum(k=0, n, U(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) )}
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Dec 14 2020
STATUS
approved