OFFSET
1,3
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..300
FORMULA
G.f. satisfies:
(1) A(x - A(x^4 - x^10)) = x.
(2) A(x) = x + Sum_{n>=0} d^n/dx^n A(x^4-x^10)^(n+1) / (n+1)!.
(3) A(x) = x * exp( Sum_{n>=0} d^n/dx^n A(x^4-x^10)^(n+1)/x / (n+1)! ).
EXAMPLE
G.f.: A(x) = x + x^4 + 4*x^7 + 21*x^10 + 126*x^13 + 817*x^16 + 5574*x^19 + 39418*x^22 + 286286*x^25 + 2122491*x^28 + 15995696*x^31 + 122166551*x^34 +...
such that A(x) = x + A( A(x)^4 - A(x)^10 ).
RELATED SERIES.
A(x)^4 = x^4 + 4*x^7 + 22*x^10 + 136*x^13 + 901*x^16 + 6248*x^19 + 44758*x^22 + 328520*x^25 + 2457286*x^28 + 18659736*x^31 + 143455026*x^34 +...
A(x)^10 = x^10 + 10*x^13 + 85*x^16 + 690*x^19 + 5520*x^22 + 44002*x^25 + 351045*x^28 + 2808040*x^31 + 22537355*x^34 + 181530280*x^37 + 1467320874*x^40 +...
A(x^4 - x^10) = x^4 - x^10 + x^16 - 4*x^22 + 10*x^28 - 32*x^34 + 106*x^40 - 350*x^46 + 1211*x^52 - 4242*x^58 + 15083*x^64 - 54404*x^70 + 198114*x^76 +...
where Series_Reversion(A(x)) = x - A(x^4 - x^10).
PROG
(PARI) {a(n) = my(A=x); for(i=1, 3*n, A = x + subst(A, x, A^4 - A^10 +x*O(x^(3*n)))); polcoeff(A, 3*n-2)}
for(n=1, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 20 2016
STATUS
approved