OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..10000
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
Expansion of b(q) * (b(q) + 2*b(q^4)) / (3 * b(q^2)) in powers of q where b() is a cubic AGM theta function.
Expansion of psi(-q) * chi(-q)^3 * phi(q^3) * chi(q^3)^3 in powers of q where phi(), psi(), chi() are Ramanujan theta functions.
Expansion of eta(q)^4 * eta(q^4) * eta(q^6)^8 / (eta(q^2)^4 * eta(q^3)^4 * eta(q^12)^3) in powers of q.
Euler transform of period 12 sequence [ -4, 0, 0, -1, -4, -4, -4, -1, 0, 0, -4, -2, ...].
Moebius transform is period 12 sequence [ -4, 10, 0, -6, 4, 0, -4, 6, 0, -10, 4, 0, ...].
a(n) = -4 * b(n) where b(n) is multiplicative with b(2^e) = (1 - (-1)^e) * -3/4 if e>0, b(3^e) = 1, b(p^e) = e+1 if p == 1 (mod 6), b(p^e) = (1 + (-1)^e) / 2 if p == 5 (mod 6) with a(0) = 1.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 48^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A244375.
EXAMPLE
G.f. = 1 - 4*q + 6*q^2 - 4*q^3 + 6*q^6 - 8*q^7 + 6*q^8 - 4*q^9 - 8*q^13 + ...
MATHEMATICA
a[ n_] := If[ n < 1, Boole[n == 0], 2 Sum[ (-1)^(n/d) {2, -1, 0, 1, -2, 0}[[ Mod[ d, 6, 1] ]], {d, Divisors @ n}]];
a[ n_] := SeriesCoefficient[ QPochhammer[ q]^4 QPochhammer[ q^4] QPochhammer[ q^6]^8 / (QPochhammer[ q^2]^4 QPochhammer[ q^3]^4 QPochhammer[ q^12]^3), {q, 0, n}];
PROG
(PARI) {a(n) = if( n<1, n==0, 2 * sumdiv(n, d, (-1)^(n/d) * [0, 2, -1, 0, 1, -2][d%6 + 1]))};
(PARI) {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^4 * eta(x^4 + A) * eta(x^6 + A)^8 / (eta(x^2 + A)^4 * eta(x^3 + A)^4 * eta(x^12 + A)^3), n))};
(PARI) {a(n) = if( n<0, 0, polcoeff( 1 + 2 * sum(k=1, n, x^k / (1 + x^k) * [0, -2, 1, 0, -1, 2][k%6 + 1], x * O(x^n)), n))};
(PARI) {a(n) = if( n<0, 0, polcoeff( 1 + 2 * sum(k=1, n, x^k / (1 + x^k + x^(2*k)) * [3, -2, 1, -2][k%4 + 1], x * O(x^n)), n))};
(PARI) {a(n) = my(A); if( n<1, n==0, A = factor(n); -4 * prod( j=1, matsize(A)[1], if( p = A[j, 1], e = A[j, 2]; if( p==2, (1 - (-1)^e) * -3/4, if( p==3, 1, if( p%6 == 1, e+1, (1 + (-1)^e) / 2))))))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Jun 26 2014
STATUS
approved