OFFSET
0,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..1035
FORMULA
a(n) ~ c * d^n / n^(2/3), where d = 10.3933629985595735031515117628403087010816839988881759248638104... and c = 0.42093748110527419326289922348630166534660617909266766696... - Vaclav Kotesovec, Feb 24 2024
EXAMPLE
G.f.: A(x) = 1 + 4*x + 32*x^2 + 240*x^3 + 2048*x^4 + 17920*x^5 + 163904*x^6 + 1526784*x^7 + 14473216*x^8 + 138743808*x^9 + 1342326528*x^10 + ...
RELATED SERIES.
The cube of 1/A(x) equals the g.f. A370018 which starts as
1/A(x)^3 = 1 - 12*x + 176*x^3 - 2752*x^6 + 43776*x^10 - 699392*x^15 + 11186176*x^21 + ... + (-4)^n * (2*4^n + 1)/3 * x^(n*(n+1)/2) + ...
and 1/A(x) equals the g.f. of A370019, which begins
1/A(x) = 1 - 4*x - 16*x^2 - 48*x^3 - 384*x^4 - 2816*x^5 - 24384*x^6 - 206336*x^7 - 1815552*x^8 - 16189440*x^9 + ... + A370019(n)*x^n + ...
PROG
(PARI) {a(n) = my(A);
A = sum(m=0, sqrtint(2*n+1), (-4)^m * (1 + 2*4^m)/3 * x^(m*(m+1)/2) +x*O(x^n))^(-1/3);
polcoeff(H=A, n)}
for(n=0, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 24 2024
STATUS
approved