(PARI)
P(n, k) = k*(2*n+2*k-4)!*(2*n+k-1)!/((n+k-1)!*(n+k-2)!*n!*(n+k)!);
F(K, N=23) = {
my(x='x + O('x^(K+1)), t='t + O('t^(N+1)),
r='t*Ser(vector(N, n, sqr(binomial(2*n, n)/(n+1))), 't),
p=x^3*Ser(apply(k->Ser(vector(N, n, P(n-1, k)), 't), [3..K])),
s=serreverse(t*(1+r)), f=subst(subst(p, 't, s), 'x, 'x*s/'t));
Vec(polcoeff(f, K));
};