|
| |
|
|
A244339
|
|
Expansion of (-2 * a(q) + 3*a(q^2) + 2*a(q^4)) / 3 in powers of q where a() is a cubic AGM theta function.
|
|
3
|
|
|
|
1, -4, 6, -4, 0, 0, 6, -8, 6, -4, 0, 0, 0, -8, 12, 0, 0, 0, 6, -8, 0, -8, 0, 0, 6, -4, 12, -4, 0, 0, 0, -8, 6, 0, 0, 0, 0, -8, 12, -8, 0, 0, 12, -8, 0, 0, 0, 0, 0, -12, 6, 0, 0, 0, 6, 0, 12, -8, 0, 0, 0, -8, 12, -8, 0, 0, 0, -8, 0, 0, 0, 0, 6, -8, 12, -4, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
0,2
|
|
|
COMMENTS
|
|
|
|
LINKS
|
|
|
|
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.
a(2*n) = A004016(n). a(2*n + 1) = -4 * A033762(n). a(3*n) = a(n). a(6*n + 1) = -4 * A097195(n). a(6*n + 5) = 0.
|
|
|
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
|
|
|
|
STATUS
|
approved
|
| |
|
|
