login
Expansion of q * f(-q, -q^4)^5 / f(-q)^3 in powers of q where f() is a Ramanujan theta function.
4

%I #21 Sep 08 2022 08:46:06

%S 1,-2,4,-3,1,2,-2,0,3,-2,2,-2,4,-6,4,1,-2,4,0,-3,2,-4,4,0,1,2,0,-4,0,

%T 2,2,-2,8,-6,-2,1,-2,0,6,0,2,-4,4,-6,3,2,-2,4,-3,-2,2,-2,4,0,2,0,0,0,

%U 0,-2,2,-4,4,-3,4,4,-2,-4,6,-6,2,0,4,-6,4,0,-4

%N Expansion of q * f(-q, -q^4)^5 / f(-q)^3 in powers of q where f() is a Ramanujan theta function.

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

%H G. C. Greubel, <a href="/A229802/b229802.txt">Table of n, a(n) for n = 1..1000</a>

%H Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>

%F Expansion of q * f(-q)^2 * (f(-q^5) / f(-q^2, -q^3))^5 = q * f(-q, -q^4)^2 * (f(-q^5) / f(-q^2, -q^3))^3 in powers of q where f() is a Ramanujan theta function. - _Michael Somos_, Jun 10 2014

%F Euler transform of period 5 sequence [ -2, 3, 3, -2, -2, ...].

%F Moebius transform is period 5 sequence [ 1, -3, 3, -1, 0, ...]. - _Michael Somos_, Jun 10 2014

%F G.f.: x * (Product_{k>0} (1 - x^k)^2) / (Product_{k>0} (1 - x^(5*k - 2)) * (1 - x^(5*k - 3)))^5.

%F a(5*n) = a(n).

%e G.f. = q - 2*q^2 + 4*q^3 - 3*q^4 + q^5 + 2*q^6 - 2*q^7 + 3*q^9 - 2*q^10 + ...

%t a[ n_] := SeriesCoefficient[ q QPochhammer[ q]^2 / (QPochhammer[ q^2, q^5] QPochhammer[ q^3, q^5])^5, {q, 0, n}]; (* _Michael Somos_, Jun 10 2014 *)

%t a[ n_] := If[ n < 1, 0, Sum[ Im[(I - 3) {1, I, -I, -1, 0}[[ Mod[ d, 5, 1] ]] ], {d, Divisors @ n}]]; (* _Michael Somos_, Jun 10 2014 *)

%o (PARI) {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( prod(k=1, n, (1 - x^k + A)^[ 2, 2, -3, -3, 2][k%5 + 1]), n))};

%o (PARI) {a(n) = if( n<1, 0, sumdiv(n, d, imag( (I - 3) * [ 0, 1, I, -I, -1][ d%5 + 1])))}; /* _Michael Somos_, Jun 10 2014 */

%o (Sage) ModularForms( Gamma1(5), 1, prec=70).1;

%o (Magma) Basis( ModularForms( Gamma1(5), 1), 70) [2] ; /* _Michael Somos_, Jun 10 2014 */

%Y Cf. A227216.

%K sign

%O 1,2

%A _Michael Somos_, Sep 29 2013