OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..2500
FORMULA
Expansion of q^(-5/3) * eta(q) * eta(q^6)^9 / (eta(q^2) * eta(q^3))^3 in powers of q.
Euler transform of period 6 sequence [-1, 2, 2, 2, -1, -4, ...].
a(n) = 1/6 * b(3*n + 5) where b() is multiplicative with b(2^e) = 2 - 2^e, b(3^e) = 0^e, b(p^e) = (p^(e+1) - 1) / (p - 1) if p>3.
EXAMPLE
G.f. = 1 - x + 2*x^2 + 3*x^4 - 2*x^5 + 4*x^6 + 5*x^8 - 5*x^9 + 8*x^10 + ...
G.f. = q^5 - q^8 + 2*q^11 + 3*q^17 - 2*q^20 + 4*q^23 + 5*q^29 - 5*q^32 + ...
MATHEMATICA
a[ n_] := SeriesCoefficient[ QPochhammer[ x] QPochhammer[ x^6]^9 / (QPochhammer[ x^2] QPochhammer[ x^3])^3, {x, 0, n}];
a[ n_] := SeriesCoefficient[ (QPochhammer[ -x^3] QPochhammer[ x^12])^3 / (QPochhammer[ -x] QPochhammer[ x^4]), {x, 0, n}];
a[ n_] := If[ n < 0, 0, Times @@ (Which[# == 2, 2 - 2^#2, # == 3, 1, True, (#^(#2 + 1) - 1) / (# - 1)] & @@@ FactorInteger[3 n + 5]) / 6]; (* Michael Somos, Jul 09 2018 *)
PROG
(PARI) {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A) * eta(x^6 + A)^9 / (eta(x^2 + A) * eta(x^3 + A))^3, n))};
(PARI) {a(n) = my(A, p, e); if( n<0, 0, n = 3*n + 5; A = factor(n); 1/6 * prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2, 2 - 2^e, p==3, 0, (p^(e+1) - 1) / (p - 1))))};
(Magma) Basis( ModularForms( Gamma0(18), 2), 210) [6]; /* Michael Somos, Jul 09 2018 */
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Nov 22 2013
STATUS
approved