OFFSET
0,6
LINKS
T. D. Noe, Table of n, a(n) for n=0..1004
R. E. Borcherds, A Siegel cusp form of weight 12 ...
FORMULA
Expansion of eta(q^2)^5 * eta(q^8)^12 / (eta(q) * eta(q^4))^2 in powers of q.
Euler transform of period 8 sequence [ 2, -3, 2, -1, 2, -3, 2, -13, ...]. - Michael Somos, Mar 06 2004
G.f.: x^4 * (Product_{k>0} (1 - x^(2*k))^5 * (1 - x^(8*k))^12 / ((1 - x^k) * (1 - x^(4*k)))^2). - Michael Somos, Mar 06 2004
EXAMPLE
G.f. = q^4 + 2*q^5 + 2*q^9 - 12*q^13 - 22*q^14 - 24*q^17 + 56*q^21 + ...
MATHEMATICA
max = 76; f[x_] := x^4*Product[ (1 - x^(2k))^5 (1 - x^(8k))^12/((1 - x^k) (1 - x^(4k)))^2, {k, 1, max}]; CoefficientList[ Series[ f[x], {x, 0, max}], x] (* Jean-François Alcover, Dec 12 2011, after Michael Somos *)
QP = QPochhammer; s = q^4*QP[q^2]^5*(QP[q^8]^12/(QP[q]*QP[q^4])^2) + O[q]^80; CoefficientList[s, q] (* Jean-François Alcover, Nov 24 2015 *)
PROG
(PARI) {a(n) = if( n<4, 0, n-=4; polcoeff( eta(x^8 + x * O(x^n))^12 * sum(k=1, sqrtint(n), 2*x^k^2, 1), n))}; /* Michael Somos, Mar 06 2004 */
(PARI) {a(n) = local(A); if( n<4, 0, n-=4; A = x * O(x^n); polcoeff(eta(x^2 + A)^5 * eta(x^8 + A)^12 / (eta(x + A)^2 * eta(x^4 + A)^2), n))}; /* Michael Somos, Mar 06 2004 */
CROSSREFS
KEYWORD
sign,nice
AUTHOR
STATUS
approved