(PARI)
E(p, m)={sum(k=0, m, sum(i=0, k, polcoef(p, i)*i!*(m-i)! )*x^k/(k!*(m-k)!))}
seq(n)={my(v=vector(n)); v[1]=1; for(n=2, n, my(p=sum(k=1, n-1, v[k]*v[n-k])); v[n]=E(intformal(p), binomial(n, 2))); vector(n, n, my(m=binomial(n, 2)); m!*polcoef(v[n], m))} \\
Andrew Howroyd, Jul 31 2024