OFFSET
1,3
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..200
FORMULA
O.g.f. A(x) satisfies: [x^n] exp( 4*n * A(x) ) = (n + 3) * [x^(n-1)] exp( 4*n * A(x) ) for n>=1.
a(n) ~ c * n! * n^7, where c = 0.00000132855349... - Vaclav Kotesovec, Mar 20 2018
EXAMPLE
O.g.f.: A(x) = x + x^2 + 6*x^3 + 50*x^4 + 520*x^5 + 6312*x^6 + 86080*x^7 + 1288704*x^8 + 20862720*x^9 + 361454720*x^10 + ...
where
A(x) = x * (1 - 3*x*A'(x)) / (1 - 4*x*A'(x)).
RELATED SERIES.
exp(A(x)) = 1 + x + 3*x^2/2! + 43*x^3/3! + 1369*x^4/4! + 69561*x^5/5! + 4991371*x^6/6! + 471516403*x^7/7! + 56029153713*x^8/8! + 8112993527089*x^9/9! + ... + A300988(n)*x^n/n! + ...
A'(x) = 1 + 2*x + 18*x^2 + 200*x^3 + 2600*x^4 + 37872*x^5 + 602560*x^6 + 10309632*x^7 + 187764480*x^8 + 3614547200*x^9 + ...
PROG
(PARI) {a(n) = my(A=x); for(i=1, n, A = x*(1-3*x*A')/(1-4*x*A' +x*O(x^n))); polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) /* [x^n] exp( 4*n * A(x) ) = (n + 3) * [x^(n-1)] exp( 4*n * A(x) ) */
{a(n) = my(A=[1]); for(i=1, n+1, A=concat(A, 0); V=Vec(Ser(A)^(4*(#A-1))); A[#A] = ((#A+2)*V[#A-1] - V[#A])/(4*(#A-1)) ); polcoeff( log(Ser(A)), n)}
for(n=1, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 17 2018
STATUS
approved