A321527 Expansion of x^3 * c(x^2) * c(x^4)^2 / (9 * c(x)) in powers of x where c() is a cubic AGM theta function.

0, 0, 0, 1, -1, 0, 2, 0, -3, 4, 0, 0, 1, 0, 0, 6, -7, 0, 8, 0, -6, 8, 0, 0, -1, 0, 0, 13, -8, 0, 12, 0, -15, 12, 0, 0, 7, 0, 0, 14, -18, 0, 16, 0, -12, 24, 0, 0, -5, 0, 0, 18, -14, 0, 26, 0, -24, 20, 0, 0, 6, 0, 0, 32, -31, 0, 24, 0, -18, 24, 0, 0, 5, 0, 0, 31

Offset

0,7

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

Number 124 of the 126 eta-quotients listed in Table 1 of Williams 2012.

Links

Antti Karttunen, Table of n, a(n) for n = 0..16384

Antti Karttunen, Data supplement: n, a(n) computed for n = 0..100000

Michael Somos, Introduction to Ramanujan theta functions

Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

K. S. Williams, Fourier series of a class of eta quotients, Int. J. Number Theory 8 (2012), no. 4, 993-1004.

Formula

Expansion of x^3 * (psi(x^3) * psi(x^6))^3 / (psi(x) * psi(x^2)) in powers of x where psi() is a Ramanujan theta function.

Expansion of x^3 * chi(-x) * f(-x^12)^6 / (chi(-x^3)^3 * f(-x^4)^2) in powers of x where chi(), f() are Ramanujan theta functions.

G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = (1/6) (t / i)^2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A321528.

a(n) = s(n/3) - s(n/4) - s(n/6) + s(n/12) where s(x) = sum of divisors of x for integer x else 0.

a(6*n + 1) = a(6*n + 5) = a(12*n + 2) = a(12*n + 10) = 0.

a(n) = A224226(n) if n>0. a(2*n) = -A229615(n). a(6*n + 3) = A008438(n). a(12*n + 6) = 2*A008438(n).

a(12*n + 3) = A112610(n). a(12*n + 4) = -A144614(n). a(12*n + 8) = -3*A033686(n). a(12*n + 9) = 4*A097723(n).

Example

G.f. = x^3 - x^4 + 2*x^6 - 3*x^8 + 4*x^9 + x^12 + 6*x^15 - 7*x^16 + ...

Mathematica

a[ n_] := SeriesCoefficient[ x^3 QPochhammer[ x, x^2] QPochhammer[ x^12]^6 / (QPochhammer[ x^3, x^6]^3 QPochhammer[ x^4]^2), {x, 0, n}];
a[ n_] := With[ {s = If[ FractionalPart @ # > 0, 0, DivisorSigma[1, #]] &}, If[ n < 1, 0, s[n/3] - s[n/4] - s[n/6] + s[n/12]]];
a[ n_] := If[ n < 1, 0, Sum[ d {0, 0, 4, -3, 0, 2, 0, -3, 4, 0, 0, 0}[[Mod[d, 12, 1]]] / 12, {d, Divisors[n]}]];

Prog

(PARI) {a(n) = if( n<1, 0, sumdiv( n, d, d * [0, 0, 0, 4, -3, 0, 2, 0, -3, 4, 0, 0][d%12 + 1] / 12))};
(PARI) {a(n) = my(s = x -> if( frac(x), 0, sigma(x))); if( n<1, 0, s(n/3) - s(n/4) - s(n/6) + s(n/12))};
(PARI) {a(n) = my(A); n-=3; if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A) * eta(x^6 + A)^3 * eta(x^12 + A)^6 / (eta(x^2 + A) * eta(x^3 + A)^3 * eta(x^4 + A)^2), n))};
(Magma) A := Basis( ModularForms( Gamma0(12), 2), 76); A[4] - A[5];

Crossrefs

Cf. A008438, A033686, A097723, A112610, A144614, A224226, A229615, A321528.

Sequence in context: A021835 A112476 A370724 * A224226 A153250 A102389

Adjacent sequences: A321524 A321525 A321526 * A321528 A321529 A321530

Keyword

sign,changed

Author

Michael Somos, Nov 12 2018

Status

approved