OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1500
FORMULA
Expansion of eta(q^2)^5 * eta(q^3)^3 / (eta(q) * eta(q^6)^3) in powers of q.
Expansion of b(q^2) * (4*b(q^4) - b(q)) / 3 in powers of q where b() is a cubic AGM theta function.
Euler transform of period 6 sequence [ 1, -4, -2, -4, 1, -4, ...].
a(n) is multiplicative with a(2^e) = -3 if e>0, a(3^e) = 4 - 3^(e+1), a(p^e) = (p^(e+1) - 1) / (p - 1) if p>3.
Dirichlet g.f.: zeta(s) * zeta(s-1) * (1 - 2^(2-s)) * (1 - 2^(1-s)) * (1 - 3^(2-s)). - Amiram Eldar, Oct 28 2023
EXAMPLE
1 + q - 3*q^2 - 5*q^3 - 3*q^4 + 6*q^5 + 15*q^6 + 8*q^7 - 3*q^8 - ...
MATHEMATICA
eta[q_] := q^(1/24)*QPochhammer[q]; a:= CoefficientList[Series[eta[q^2]^5 *eta[q^3]^3/(eta[q]*eta[q^6]^3), {q, 0, 55}], q]; Table[a[[n]], {n, 1, 50}] (* G. C. Greubel, Jul 11 2018 *)
PROG
(PARI) {a(n) = if( n<0, 0, polcoeff( sum( k=1, sqrtint(n), 2 * (-x)^k^2, 1 + x * O(x^n))^4 - 9 * sum( k=1, sqrtint(n\3), 2 * (-x^3)^k^2, 1 + x * O(x^n))^4, n) / -8)}
(PARI) {a(n) = if( n<1, n==0, -(-1)^n * ( sumdiv( n, d, d * (1 - if( d%3==0, 3) - if( d%4==0, 1) + if(d%12==0, 3)))))}
(PARI) {a(n) = local(A, p, e); if( n<1, n==0, A = factor(n); prod( k=1, matsize(A)[1], if( p=A[k, 1], e=A[k, 2]; if( p==2, -3, if( p==3, 4 - 3^(e+1), (p^(e+1) - 1) / (p - 1))))))}
(PARI) {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^5 * eta(x^3 + A)^3 / eta(x + A) / eta(x^6 + A)^3, n))}
CROSSREFS
KEYWORD
sign,mult
AUTHOR
Michael Somos, Apr 21 2006
STATUS
approved