(PARI) {a(n)=local(A=1+x, F=sum(k=0, n, binomial(2*k+1, k)/(2*k+1)*x^k)+x*O(x^n)); for(i=0, n, A=subst(F^2, x, serreverse(x/(A+x*O(x^n))))); polcoeff(A, n)}
(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1+A*serreverse(x/(A+x*O(x^n))))^2); polcoeff(A, n)}