A263571 Expansion of f(x^2, x^2) * f(x, x^5) in powers of x where f(, ) is Ramanujan's general theta function.

1, 1, 2, 2, 0, 1, 0, 2, 3, 2, 2, 0, 0, 2, 0, 0, 3, 0, 4, 2, 0, 1, 0, 4, 2, 0, 2, 0, 0, 2, 0, 0, 2, 3, 2, 2, 0, 2, 0, 2, 3, 2, 2, 0, 0, 0, 0, 0, 4, 0, 2, 4, 0, 2, 0, 2, 1, 0, 6, 0, 0, 0, 0, 0, 2, 3, 2, 2, 0, 0, 0, 2, 4, 4, 2, 0, 0, 2, 0, 0, 2, 0, 0, 4, 0, 1, 0

Offset

0,3

Comments

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

Links

G. C. Greubel, Table of n, a(n) for n = 0..5000

Michael Somos, Introduction to Ramanujan theta functions, 2019.

Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.

Formula

Expansion of chi(x) * phi(x^2) * psi(-x^3) in powers of x where phi(), psi(), chi() are Ramanujan theta functions.

Expansion of q^(-1/3) * eta(q^3) * eta(q^4)^4 * eta(q^12) / (eta(q) * eta(q^6) * eta(q^8)^2) in powers of q.

a(n) = b(3*n + 1) where b() is multiplicative with b(2^e) = b(3^e) = (-1)^e, b(p^e) = e+1 if p == 1, 7 (mod 24), b(p^e) = (e+1) * (-1)^e if p == 5, 11 (mod 24), b(p^e) = (1 + (-1)^e) / 2 if p == 13, 17, 19, 23 (mod 24).

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

a(n) = A115660(3*n + 1) = A192013(3*n + 1) = A128581(6*n + 2).

a(2*n) = A261115(n). a(2*n + 1) = A263548(n). a(4*n + 1) = a(n). a(4*n + 3) = 2 * A128582(n).

a(8*n + 4) = a(8*n + 6) = 0. a(8*n) = A113780(n). a(8*n + 2) = 2 * A260089(n).

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/sqrt(6) = 1.282549... . - Amiram Eldar, Dec 28 2023

Example

G.f. = 1 + x + 2*x^2 + 2*x^3 + x^5 + 2*x^7 + 3*x^8 + 2*x^9 + 2*x^10 + ...
G.f. = q + q^4 + 2*q^7 + 2*q^10 + q^16 + 2*q^22 + 3*q^25 + 2*q^28 + ...

Mathematica

a[ n_] := If[ n < 0, 0, With[ {m = 3 n + 1}, Sum[ KroneckerSymbol[ 2, d] KroneckerSymbol[ -3, m/d], {d, Divisors[ m]}]]];
a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^2] EllipticTheta[ 3, 0, x^2] EllipticTheta[ 2, Pi/4, x^(3/2)] / (2^(1/2) x^(3/8)), {x, 0, n}];

Prog

(PARI) {a(n) = my(m); if( n<0, 0, m = 3*n + 1; sumdiv( m, d, kronecker( 2, d) * kronecker( -3, m/d)))};
(PARI) {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^3 + A) * eta(x^4 + A)^4 * eta(x^12 + A) / (eta(x + A) * eta(x^8 + A)^2), n))};

Crossrefs

Cf. A113780, A115660, A128581, A128582, A192013, A260089, A261115, A263548, A263577.

Cf. A000122, A000700, A010054, A121373.

Sequence in context: A264403 A038190 A279947 * A364060 A361292 A251690

Adjacent sequences: A263568 A263569 A263570 * A263572 A263573 A263574

Keyword

nonn,easy

Author

Michael Somos, Oct 21 2015

Status

approved