OFFSET
-1,3
LINKS
G. C. Greubel, Table of n, a(n) for n = -1..1000
Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], 2015-2016.
FORMULA
Expansion of (eta(q^3) * eta(q^15))^2 / (eta(q) * eta(q^5) * eta(q^9) * eta(q^45)) in powers of q.
Euler transform of period 45 sequence [ 1, 1, -1, 1, 2, -1, 1, 1, 0, 2, 1, -1, 1, 1, -2, 1, 1, 0, 1, 2, -1, 1, 1, -1, 2, 1, 0, 1, 1, -2, 1, 1, -1, 1, 2, 0, 1, 1, -1, 2, 1, -1, 1, 1, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u*v * (u*v + 3) - (u+v) * (u^2 + u*v + v^2).
G.f. is a period 1 Fourier series which satisfies f(-1 / (45 t)) = f(t) where q = exp(2 Pi i t).
a(n) = A058684(n) unless n=0.
a(n) ~ exp(4*Pi*sqrt(n/5)/3) / (5^(1/4) * sqrt(6) * n^(3/4)). - Vaclav Kotesovec, Oct 14 2015
EXAMPLE
1/q + 1 + 2*q + q^2 + 3*q^3 + 4*q^4 + 5*q^5 + 6*q^6 + 7*q^7 + 11*q^8 + ...
MATHEMATICA
nmax=60; CoefficientList[Series[Product[(1-x^(3*k))^2 * (1-x^(15*k))^2 / ((1-x^k) * (1-x^(5*k)) * (1-x^(9*k)) * (1-x^(45*k))), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 14 2015 *)
eta[q_] := q^(1/24)*QPochhammer[q]; a:= CoefficientList[Series[q*(eta[q^3]*eta[q^15])^2/(eta[q]*eta[q^5]*eta[q^9]*eta[q^45]), {q, 0, 60}], q]; Table[a[[n]], {n, 1, 50}] (* G. C. Greubel, Jun 20 2018 *)
PROG
(PARI) {a(n) = local(A); if( n<-1, 0, n++; A = x * O(x^n); polcoeff( (eta(x^3 + A) * eta(x^15 + A))^2 / (eta(x + A) * eta(x^5 + A) * eta(x^9 + A) * eta(x^45 + A)), n))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael Somos, May 24 2013
STATUS
approved