OFFSET
1,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..400
FORMULA
G.f. A(x) = Sum_{n>=1} a(n)*x^n satisfies the following formulas.
(1.a) A(x)^5 = A( x*A(x)^4 * (1 + 5*A(x)) ).
(1.b) A(x)^25 = A( x*A(x)^24 * (1 + 5*A(x))*(1 + 5*A(x)^5) ).
(1.c) A(x)^125 = A( x*A(x)^124 * (1 + 5*A(x))*(1 + 5*A(x)^5)*(1 + 5*A(x)^25) ).
(1.d) A(x)^(5^n) = A( x*A(x)^(5^n-1) * Product_{k=0..n-1} (1 + 5*A(x)^(5^k)) ).
(2) A(x) = x * Product_{n>=0} (1 + 5*A(x)^(5^n)).
(3) A(x) = Series_Reversion( x / Product_{n>=0} (1 + 5*x^(5^n)) ).
EXAMPLE
G.f.: A(x) = x + 5*x^2 + 25*x^3 + 125*x^4 + 625*x^5 + 3130*x^6 + 15800*x^7 + 81625*x^8 + 443125*x^9 + 2609375*x^10 + 16984500*x^11 + 121023875*x^12 + ...
where A(x)^5 = A( x*A(x)^4 + 5*x*A(x)^5 ).
PROG
(PARI) {a(n) = my(A=[0, 1]); for(i=1, n, A=concat(A, 0);
F=Ser(A); A[#A] = polcoeff(subst(F, x, x*F^4 + 5*x*F^5) - F^5, #A+3) ); A[n+1]}
for(n=1, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 27 2024
STATUS
approved