(PARI) {a(n)=local(A=1+x*y); for(i=1, n, A=1+sum(m=1, n, x^m*y*A^m/(1-y*A^(2*m)+x*O(x^n)+y*O(y^n)))); polcoeff(polcoeff(A, n, x), n, y)}
(PARI) {a(n)=local(A=1+x*y); for(i=1, n, A=1+sum(m=1, n, y^m*x*A^(2*m-1)/(1-x*A^(2*m-1)+x*O(x^n)+y*O(y^n)))); polcoeff(polcoeff(A, n, y), n, x)}