OFFSET
0,2
FORMULA
G.f. A(x) satisfies the following identities.
(1) 1 = Sum_{n>=0} ( (1+x)^n - 9*x*A(x) )^n * 4^n / 5^(n+1).
(2) 1 = Sum_{n>=0} (1+x)^(n^2) * 4^n / (5 + 36*x*A(x)*(1+x)^n)^(n+1).
EXAMPLE
G.f.: A(x) = 1 + 44*x + 7096*x^2 + 1926724*x^3 + 711117536*x^4 + 325957899584*x^5 + 176862173366416*x^6 + 110333447177205584*x^7 + ...
such that
1 = 1/5 + ((1+x) - 9*x*A(x))*4/5^2 + ((1+x)^2 - 9*x*A(x))^2*4^2/5^3 + ((1+x)^3 - 9*x*A(x))^3*4^3/5^4 + ((1+x)^4 - 9*x*A(x))^4*4^4/5^5 + ...
Also,
1 = 1/(5 + 36*x*A(x)) + (1+x)*4/(5 + 36*x*A(x)*(1+x))^2 + (1+x)^4*4^2/(5 + 36*x*A(x)*(1+x)^2)^3 + (1+x)^9*4^3/(5 + 36*x*A(x)*(1+x)^3)^4 + ...
PROG
(PARI) \p120
{A=vector(1); A[1]=1; for(i=1, 20, A = concat(A, 0);
A[#A] = round( Vec( sum(n=0, 1000, ( (1+x +x*O(x^#A))^n - 9*x*Ser(A) )^n * 4^n/5^(n+1)*1.)/36 ) )[#A+1]); A}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 10 2019
STATUS
approved