OFFSET
0,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..200
FORMULA
G.f.: A(x) = Sum_{n>=0} x^n*((1+x)^n + sqrt(A(x)))^n / (1 + x*sqrt(A(x))*(1+x)^n)^(n+1).
G.f.: A(x) = Sum_{n>=0} x^n*((1+x)^n - sqrt(A(x)))^n / (1 - x*sqrt(A(x))*(1+x)^n)^(n+1).
EXAMPLE
G.f.: A(x) = 1 + x + 2*x^2 + 3*x^3 + 8*x^4 + 17*x^5 + 61*x^6 + 223*x^7 + 1058*x^8 + 5355*x^9 + 29477*x^10 + 171068*x^11 + 1042319*x^12 + ...
Let B = sqrt(A(x)) then
A(x) = 1/(1 + x*B) + x*((1+x) + B)/(1 + x*B*(1+x))^2 + x^2*((1+x)^2 + B)^2/(1 + x*B*(1+x)^2)^3 + x^3*((1+x)^3 + B)^3/(1 + x*B*(1+x)^3)^4 + x^4*((1+x)^4 + B)^4/(1 + x*B*(1+x)^4)^5 + x^5*((1+x)^5 + B)^5/(1 + x*B*(1+x)^5)^6 + ...
and
A(x) = 1/(1 - x*B) + x*((1+x) - B)/(1 - x*B*(1+x))^2 + x^2*((1+x)^2 - B)^2/(1 - x*B*(1+x)^2)^3 + x^3*((1+x)^3 - B)^3/(1 - x*B*(1+x)^3)^4 + x^4*((1+x)^4 - B)^4/(1 - x*B*(1+x)^4)^5 + x^5*((1+x)^5 - B)^5/(1 - x*B*(1+x)^5)^6 + ...
PROG
(PARI) {a(n) = my(A=[1]); for(i=1, n,
A = Vec( sum(m=0, #A, x^m*((1+x)^m + Ser(A)^(1/2))^m/(1 + x*Ser(A)^(1/2)*(1+x)^m)^(m+1)) )); A[n+1]}
for(n=0, 35, print1(a(n), ", "))
(PARI) {a(n) = my(A=[1]); for(i=1, n,
A = Vec( sum(m=0, #A, x^m*((1+x)^m - Ser(A)^(1/2))^m/(1 - x*Ser(A)^(1/2)*(1+x)^m)^(m+1)) )); A[n+1]}
for(n=0, 35, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 23 2019
STATUS
approved