%I #34 Sep 09 2022 15:59:07
%S 1,720,179280,16954560,396974160,4632858720,34413301440,187477879680,
%T 814940600400,2975469665040,9486467837280,27053330840640,
%U 70485969919680,169930679355360,384163875688320,820167497170560,1668890801059920,3249626139960480,6096884624994960
%N Theta series of direct sum of 3 copies of E_8 lattice (the Niemeier lattice of type E_8^3).
%C Also the theta series for the Niemeier lattice of type E_8 D_16. - _Ben Mares_, Jul 17 2022
%D J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 123, 407.
%H Seiichi Manyama, <a href="/A008411/b008411.txt">Table of n, a(n) for n = 0..10000</a>
%F This series is the q-expansion of E_4(z)^3. Cf. A004009. - _Daniel D. Briggs_, Nov 25 2011
%F 691*a(n) - A029828(n) = 432000*A000594(n). - _Seiichi Manyama_, Jan 28 2017
%e G.f. = 1 + 720*q + 179280*q^2 + 16954560*q^3 + 396974160*q^4 + ...
%t a[ n_] := SeriesCoefficient[ With[ {t2 = EllipticTheta[ 2, 0, q]^4, t3 = EllipticTheta[ 3, 0, q]^4}, (t2^2 + 14 t2 t3 + t3^2)^3 ], {q, 0, n}]; (* _Michael Somos_, Jan 28 2017 *)
%t terms = 19; QP = QPochhammer; s = (QP[x]^24 + 256*x*QP[x^2]^24)^3 / (QP[x]*QP[x^2])^24 + O[x]^terms; CoefficientList[s, x] (* _Jean-François Alcover_, Jul 07 2017, adapted from PARI *)
%t terms = 19; E4[x_] = 1 + 240*Sum[k^3*x^k/(1 - x^k), {k, 1, terms}]; E4[x]^3 + O[x]^terms // CoefficientList[#, x]& (* _Jean-François Alcover_, Feb 28 2018 *)
%o (PARI) {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x + A)^24 + 256 * x * eta(x^2 + A)^24)^3 / (eta(x + A) * eta(x^2 + A))^24, n))}; /* _Michael Somos_, Jan 28 2017 */
%o (Magma) A := Basis( ModularForms( Gamma1(1), 12), 19); A[1] + 720*A[2]; /* _Michael Somos_, Jan 28 2017 */
%Y Cf. A000594, A004009, A029828, A280869.
%K nonn,easy
%O 0,2
%A _N. J. A. Sloane_