%I #25 Jan 03 2023 14:50:04
%S 1,2,0,0,4,4,0,0,5,4,0,0,4,8,0,0,8,6,0,0,8,4,0,0,5,12,0,0,12,8,0,0,8,
%T 8,0,0,4,12,0,0,16,8,0,0,12,8,0,0,9,14,0,0,12,16,0,0,8,4,0,0,12,16,0,
%U 0,16,16,0,0,16,8,0,0,8,20,0,0,16,8,0,0,17
%N Half of theta series of b.c.c. lattice with respect to long edge.
%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
%C The body-centered cubic (b.c.c. also known as D3*) lattice is the set of all triples [a, b, c] where the entries are all integers or all one half an odd integer. A long edge is centered at a triple with two integer entries and the remaining entry is one half an odd integer. - _Michael Somos_, May 31 2012
%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 From _Michael Somos_, May 31 2012: (Start)
%F Expansion of x * phi(x) * psi(x^4)^2 = x * psi(-x^2)^4 / phi(-x) in powers of x where phi(), psi() are Ramanujan theta functions.
%F Expansion of eta(q^2)^5 * eta(q^8)^4 / (eta(q)^2 * eta(q^4)^4) in powers of q.
%F Euler transform of period 8 sequence [ 2, -3, 2, 1, 2, -3, 2, -3, ...].
%F a(4*n) = a(4*n + 3) = 0. a(n) = A004025(n) / 2. a(4*n + 1) = A045834(n). a(4*n + 2) = 2 * A045828(n). (End)
%e q + 2*q^2 + 4*q^5 + 4*q^6 + 5*q^9 + 4*q^10 + 4*q^13 + 8*q^14 + 8*q^17 + ...
%t a[n_] := Module[{A = x*O[x]^n}, SeriesCoefficient[QPochhammer[x^2+A]^5 * (QPochhammer[x^8+A]^4 / (QPochhammer[x+A]^2*QPochhammer[x^4+A]^4)), {x, 0, n}]]; Table[a[n], {n, 0, 80}] (* _Jean-François Alcover_, Nov 05 2015, adapted from PARI *)
%o (PARI) {a(n) = local(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A)^5 * eta(x^8 + A)^4 / (eta(x + A)^2 * eta(x^4 + A)^4), n))} /* _Michael Somos_, May 31 2012 */
%Y Cf. A004025, A045828, A045834.
%K nonn
%O 1,2
%A _N. J. A. Sloane_