A226059 Expansion of eta(q) * eta(q^9) * eta(q^21)^2 / (eta(q^3)^2 * eta(q^7) * eta(q^63)) in powers of q.

1, -1, -1, 2, -2, -1, 5, -3, -4, 8, -5, -6, 16, -8, -11, 23, -15, -16, 39, -21, -26, 58, -35, -39, 92, -51, -58, 132, -77, -85, 194, -108, -125, 276, -156, -174, 393, -218, -245, 542, -304, -336, 755, -417, -467, 1026, -573, -627, 1401, -770, -853, 1870

Offset

-1,4

Links

G. C. Greubel, Table of n, a(n) for n = -1..1500

Formula

Euler transform of period 63 sequence [ -1, -1, 1, -1, -1, 1, 0, -1, 0, -1, -1, 1, -1, 0, 1, -1, -1, 0, -1, -1, 0, -1, -1, 1, -1, -1, 0, 0, -1, 1, -1, -1, 1, -1, 0, 0, -1, -1, 1, -1, -1, 0, -1, -1, 0, -1, -1, 1, 0, -1, 1, -1, -1, 0, -1, 0, 1, -1, -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 / (63 t)) = 1 / f(t) where q = exp(2 Pi i t).

Example

1/q - 1 - q + 2*q^2 - 2*q^3 - q^4 + 5*q^5 - 3*q^6 - 4*q^7 + 8*q^8 - 5*q^9 + ...

Mathematica

eta[q_] := q^(1/24)*QPochhammer[q]; b := eta[q]*eta[q^9]*eta[q^21]^2/ (eta[q^3]^2*eta[q^7]*eta[q^63]); a:= CoefficientList[Series[q*b , {q, 0, 100}], q]; Table[a[[n]], {n, 1, 80}] (* G. C. Greubel, Jul 03 2018 *)

Prog

(PARI) {a(n) = local(A); if( n<-1, 0, n++; A = x * O(x^n); polcoeff( eta(x + A) * eta(x^9 + A) * eta(x^21 + A)^2 / (eta(x^3 + A)^2 * eta(x^7 + A) * eta(x^63 + A)), n))}

Crossrefs

Sequence in context: A184050 A367094 A324798 * A353738 A127742 A110438

Adjacent sequences: A226056 A226057 A226058 * A226060 A226061 A226062

Keyword

sign

Author

Michael Somos, May 24 2013

Status

approved