OFFSET
0,2
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..300
FORMULA
G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) 1 = Sum_{n>=0} (-x/2)^n * (A(x)^n + (-1)^n)^n.
(2) 1 = Sum_{n>=0} 2 * (-x)^n * A(x)^(n^2) / (2 - (-1)^n * x * A(x)^n)^(n+1).
a(n) = A325574(n)/2^n for n >= 0.
EXAMPLE
G.f.: A(x) = 1 + 2*x + 8*x^2 + 50*x^3 + 376*x^4 + 3124*x^5 + 27804*x^6 + 260496*x^7 + 2539616*x^8 + 25556330*x^9 + ...
such that
1 = 1 - (x/2)*(A(x) - 1) + (x/2)^2*(A(x)^2 + 1)^2 - (x/2)^3*(A(x)^3 - 1)^3 + (x/2)^4*(A(x)^4 + 1)^4 - (x/2)^5*(A(x)^5 - 1)^5 + (x/2)^6*(A(x)^6 + 1)^6 - (x/2)^7*(A(x)^7 - 1)^7 + ...
also,
1 = 2/(2 - x) - 2*x*A(x)/(2 + x*A(x))^2 + 2*x^2*A(x)^4/(2 - x*A(x)^2)^3 - 2*x^3*A(x)^9/(2 + x*A(x)^3)^4 + 2*x^4*A(x)^16/(2 - x*A(x)^4)^5 - 2*x^5*A(x)^25/(2 + x*A(x)^5)^6 + 2*x^6*A(x)^36/(2 - x*A(x)^6)^7 ... + 2*(-x)^n*A(x)^(n^2)/(2 - (-1)^n*x*A(x)^n)^(n+1) + ...
PROG
(PARI) {a(n) = my(A=[1]);
for(i=1, n, A=concat(A, 0); A[#A] = 2*polcoeff( sum(m=0, #A, (-x/2)^m * (Ser(A)^m + (-1)^m)^m ), #A)); A[n+1]}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 03 2023
STATUS
approved