A128616 Expansion of q * psi(q^3) * psi(q^5) in powers of q where psi() is a Ramanujan theta function.

1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 2, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 2, 0, 0, 2, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 2, 0, 0, 1, 0, 0, 0, 0, 0, 1, 2, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 2, 0, 0, 2, 0, 1, 0, 0, 0, 2, 0, 0

Offset: 1

Author

Michael Somos, Mar 13 2007

Keywords

• nonn

Comments

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

Examples

			G.f. = x + x^4 + x^6 + x^9 + x^10 + x^15 + x^16 + 2*x^19 + x^24 + x^25 + 2*x^31 + ...
		

References

• B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 377, Entry 9(iv).

Links

• G. C. Greubel, Table of n, a(n) for n = 1..10000
• Michael Somos, Introduction to Ramanujan theta functions
• Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

Crossrefs

Cf. A082451, A121362.

Programs

Mathematica

a[ n_] := If[ n < 1, 0, DivisorSum[ n, KroneckerSymbol[ -60, #] + KroneckerSymbol[ 20, #] KroneckerSymbol[ -3, n/#] &] / 2]; (* Michael Somos, Nov 12 2015 *)
a[ n_] := SeriesCoefficient[ q(QPochhammer[ q^6] QPochhammer[ q^10])^2 / (QPochhammer[ q^3] QPochhammer[ q^5]), {q, 0, n}]; (* Michael Somos, Nov 12 2015 *)

PARI

{a(n) = if( n<1, 0, sumdiv(n, d, kronecker(-60, d) + kronecker(20, d) * kronecker(-3, n/d) )/2)};

PARI

{a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x^6 + A) * eta(x^10 + A))^2 / (eta(x^3 + A) * eta(x^5 + A)), n))};

Formula

Expansion of (eta(q^6) * eta(q^10))^2 / (eta(q^3) * eta(q^5)) in powers of q.

Euler transform of period 30 sequence [ 0, 0, 1, 0, 1, -1, 0, 0, 1, -1, 0, -1, 0, 0, 2, 0, 0, -1, 0, -1, 1, 0, 0, -1, 1, 0, 1, 0, 0, -2, ...].

For n>0, n in A028957 equivalent to a(n) nonzero. If a(n) nonzero, a(n) = A082451(n) and a(n) = A121362(n).

a(n) = (A082451(n) + A121362(n))/2.

G.f.: x * Product_{k>0} (1 - x^(3*k)) * (1 - x^(5*k)) * (1 + x^(6*k))^2 * (1 + x^(10*k))^2.