OFFSET
0,1
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
Euler transform of period 36 sequence [ 0, 0, 0, 1, -1, 0, 0, -1, 1, 0, 0, 0, -1, 1, 0, 0, 0, -2, 0, 0, 0, 1, -1, 0, 0, 0, 1, -1, 0, 0, -1, 1, 0, 0, 0, -1, ...].
G.f.: Sum_{k in Z} (-1)^k * x^(18*k^2) * (x^(17*k + 4) + x^k).
a(n) = A133985(3*n). a(5*n + 1) = a(5*n + 3) = 0.
EXAMPLE
G.f. = 1 + x^4 - x^5 - x^17 - x^19 - x^39 + x^42 + x^70 + x^74 + x^110 + ...
G.f. = q + q^289 - q^361 - q^1225 - q^1369 - q^2809 + q^3025 + q^5041 + q^5329 + ...
MATHEMATICA
QP:= QPochhammer; a[n_]:= SeriesCoefficient[QP[q^36, q^36]*(QP[q^17, q^36]*QP[q^19, q^36] + q^4*QP[q, q^36]*QP[q^35, q^36]), {q, 0, n}]; Table[a[n], {n, 0, 50}] (* G. C. Greubel, Apr 04 2018 *)
PROG
(PARI) {a(n) = my(m); if( n<0, 0, if( issquare( 72*n + 1, &m), if( m%36 != 1, m=-m); if( m%36 == 1, m = (m-1) / 36, if( m%36 != 19, m=-m); if( m%36! = 19, return(0)); m = (m+17) / 36); (-1)^m))};
(PARI) {a(n) = my(A); if( n<0, 0, n*=3; A = x * O(x^n); polcoeff( eta(x + A) * eta(x^4 + A) * eta(x^6 + A)^5 / ( eta(x^2 + A) * eta(x^3 + A) * eta(x^12 + A) )^2, n))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Mar 11 2011
STATUS
approved