OFFSET
0,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..300
FORMULA
G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies the following formulas.
(1) A(x)^3 = Sum_{n>=0} 11^(n*(n-1)/2) * (1 + 11^(2*n+1))/12 * x^(n*(n+1)/2).
(2) A(x)^3 = Product_{n>=1} (1 + 11^(n-1)*x^n) * (1 - 11^n*x^n) * (1 + 11^(n+1)*x^n), by the Jacobi triple product identity.
(3) A(x) = Product_{n>=1} F( 11^(n-1)*x^n ), where F(x) = (1 + 111*x - 1221*x^2 - 1331*x^3)^(1/3).
a(n) ~ (-1)^(n+1) * c * 11^(2*n) / n^(4/3), where c = 0.2588865455859866840901787578907966... - Vaclav Kotesovec, Feb 27 2024
EXAMPLE
G.f.: A(x) = 1 + 37*x - 1369*x^2 + 133632*x^3 - 9888768*x^4 + 845367083*x^5 - 78838949376*x^6 + 7721334144755*x^7 - 776624602305024*x^8 + ...
where the cube of g.f. A(x) yields the series
A(x)^3 = 1 + 111*x + 147631*x^3 + 2161452161*x^6 + 348104014265601*x^10 + 616687495357008127151*x^15 + 12017494675541950940487123311*x^21 + ... + 11^(n*(n-1)/2) * (1 + 11^(2*n+1))/12 * x^(n*(n+1)/2) + ...
The cube of g.f. A(x) also equals the infinite product
A(x)^3 = (1 + x)*(1 - 11*x)*(1 + 11^2*x) * (1 + 11*x^2)*(1 - 11^2*x^2)*(1 + 11^3*x^2) * (1 + 11^2*x^3)*(1 - 11^3*x^3)*(1 + 11^4*x^3) * (1 + 11^3*x^4)*(1 - 11^4*x^4)*(1 + 11^5*x^4) * ...
PROG
(PARI) {a(n) = my(A); A = prod(m=1, n+1, (1 + 11^(m-1)*x^m) * (1 - 11^m*x^m) * (1 + 11^(m+1)*x^m) +x*O(x^n))^(1/3); polcoeff(H=A, n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Feb 25 2024
STATUS
approved