cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A008408 Theta series of Leech lattice.

Original entry on oeis.org

1, 0, 196560, 16773120, 398034000, 4629381120, 34417656000, 187489935360, 814879774800, 2975551488000, 9486551299680, 27052945920000, 70486236999360, 169931095326720, 384163586352000, 820166620815360, 1668890090322000
Offset: 0

Views

Author

Keywords

Examples

			G.f. = 1 + 196560*q^2 + 16773120*q^3 + 398034000*q^4 + 4629381120*q^5 + ...
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Third Edition, Springer-Verlag,1993, pp. 51, 134-135.
  • W. Ebeling, Lattices and Codes, Vieweg; 2nd ed., 2002, see p. 113.
  • N. J. A. Sloane, Seven Staggering Sequences, in Homage to a Pied Puzzler, E. Pegg Jr., A. H. Schoen and T. Rodgers (editors), A. K. Peters, Wellesley, MA, 2009, pp. 93-110.

Crossrefs

Programs

  • Magma
    // Theta series of the Leech lattice, from John Cannon, Dec 29 2006
    A008408Q := function(prec) M12 := ModularForms(Gamma0(1), 12); t1 := Basis(M12)[1]; T := PowerSeries(t1, prec); return Coefficients(T); end function; Q := A008408Q(1000); Q[678];
    
  • Magma
    Basis( ModularForms( Gamma0(1), 12), 30) [1] ; /* Michael Somos, Jun 09 2014 */
    
  • Maple
    with(numtheory); f := 1+240*add(sigma[ 3 ](m)*q^(2*m),m=1..50); t := q^2*mul((1-q^(2*m))^24,m=1..50); series(f^3-720*t,q,51);
  • Mathematica
    max = 17; f = 1 + 240*Sum[ DivisorSigma[3, m]*q^(2m), {m, 1, max}]; t = q^2*Product[(1 - q^(2m))^24, {m, 1, max}]; Partition[ CoefficientList[ Series[f^3 - 720t, {q, 0, 2 max}], q], 2][[All, 1]] (* Jean-François Alcover , Oct 14 2011, after Maple *)
    (* From version 6 on *) f[q_] = LatticeData["Leech", "ThetaSeriesFunction"][x] /. x -> -I*Log[q]/Pi; Series[f[q], {q, 0, 32}] // CoefficientList[#, q^2]& (* Jean-François Alcover, May 15 2013 *)
    a[ n_] := If[ n < 1, Boole[n == 0], SeriesCoefficient[(1 + 240 Sum[ q^k DivisorSigma[ 3, k], {k, n}])^3 - 720 q QPochhammer[ q]^24, {q, 0, n}]]; (* Michael Somos, Jun 09 2014 *)
  • PARI
    {a(n) = if( n<1, n==0, polcoeff( 1 + (sum(k=1, n, sigma(k,11)*x^k) - x*eta(x + O(x^n))^24) * 65520/691, n))}; /* Michael Somos, Oct 19 2006 */
    
  • PARI
    {a(n) = if( n<1, n==0, polcoeff( sum(k=1, n, 240*sigma(k,3)*x^k, 1 + x*O(x^n))^3 - 720*x*eta(x + O(x^n))^24, n))}; /* Michael Somos, Oct 19 2006 */
    
  • Python
    from sympy import divisor_sigma
    def A008408(n): return 65520*(divisor_sigma(n,11)-(n**4*divisor_sigma(n)-24*((m:=n+1>>1)**2*(0 if n&1 else (m*(35*m - 52*n) + 18*n**2)*divisor_sigma(m)**2)+sum((i*(i*(i*(70*i - 140*n) + 90*n**2) - 20*n**3) + n**4)*divisor_sigma(i)*divisor_sigma(n-i) for i in range(1,m)))))//691 if n else 1 # Chai Wah Wu, Nov 17 2022
  • Sage
    A = ModularForms( Gamma0(1), 12, prec=30) . basis() ; A[1] - 65520/691*A[0] # Michael Somos, Jun 09 2014
    

Formula

The simplest way to obtain this is to take the cube of the theta series for E_8 (A004009) and subtract 720 times the g.f. for the Ramanujan numbers (A000594).
This theta series is thus also the q-expansion of (7/12) E_4(z)^3 + (5/12) E_6(z)^2. Cf. A013973. - Daniel D. Briggs, Nov 25 2011
a(n) = 65520*(A013959(n) - A000594(n))/691, n >= 1. a(0) = 1. Expansion of the Theta series of the Leech lattice in powers of q^2. See the Conway and Sloane reference. - Wolfdieter Lang, Jan 16 2017

A305920 Inverse Euler transform of A008408.

Original entry on oeis.org

0, 196560, 16773120, -18919981080, -3292295086080, 2312547886368720, 640457437563740160, -302667453389051314200, -123005476312830648176640, 39529719620247267255853008, 23306082528463942764630528000, -4849033309391159571741461446680
Offset: 1

Views

Author

Seiichi Manyama, Jun 14 2018

Keywords

Comments

a(n) is a multiple of 24.

Examples

			1 + 196560*q^2 + 16773120*q^3 + 398034000*q^4 + ... = (1-q^2)^(-196560) * (1-q^3)^(-16773120) * (1-q^4)^18919981080 * ... .
		

Crossrefs

A108772 ( (Theta series of E_8)^(1/8) - (theta series of Leech lattice)^(1/24) ) / 30.

Original entry on oeis.org

0, 1, -369, -9408, 22853743, 4794451200, -3085237931328, -819499646151424, 384739275624398865, 158606464687370095617, -50143660099205286196800, -29592112124024539896414528, 6061193531453485412550560256, 5483327637634568394533944708352
Offset: 0

Views

Author

N. J. A. Sloane and Nadia Heninger, Jun 21 2005

Keywords

Crossrefs

Equals (A108091-A108093)/30.
Showing 1-3 of 3 results.