OFFSET
1,3
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..5000
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
Expansion of (1 - eta(q)^2 * eta(q^2) / eta(q^4)) / 2 in powers of q.
G.f.: (1 - Product_{k>0} (1 - x^k)^2 / (1 + x^(2*k)) )/2.
EXAMPLE
G.f. = q + q^2 - 2*q^3 - q^4 + 2*q^6 - q^8 + 3*q^9 - 2*q^11 - 2*q^12 - q^16 + ...
MATHEMATICA
a[ n_] := If[ n < 1, 0, (-1)^Quotient[ n - 1, 2] DivisorSum[n, KroneckerSymbol[-2, #] &]]; (* Michael Somos, Jun 09 2015 *)
a[ n_] := SeriesCoefficient[ (1 - EllipticTheta[ 4, 0, q] EllipticTheta[ 4, 0, q^2]) / 2 , {q, 0, n}]; (* Michael Somos, Jun 09 2015 *)
a[ n_] := SeriesCoefficient[ (1 - QPochhammer[ q]^2 QPochhammer[ q^2] / QPochhammer[ q^4]) / 2 , {q, 0, n}]; (* Michael Somos, Jun 09 2015 *)
PROG
(PARI) {a(n) = if( n<1, 0, (-1)^((n-1)\2) * sumdiv(n, d, kronecker( -2, d)))};
(PARI) {a(n) = local(A); if( n<1, 0, A = x * O(x^n); polcoeff( (1 - eta(x + A)^2 * eta(x^2 + A) / eta(x^4 + A)) / 2, n))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Mar 30 2007
STATUS
approved