A045828 -id:A045828 - cp's OEIS frontend

A005875 Theta series of simple cubic lattice; also number of ways of writing a nonnegative integer n as a sum of 3 squares (zero being allowed).

1, 6, 12, 8, 6, 24, 24, 0, 12, 30, 24, 24, 8, 24, 48, 0, 6, 48, 36, 24, 24, 48, 24, 0, 24, 30, 72, 32, 0, 72, 48, 0, 12, 48, 48, 48, 30, 24, 72, 0, 24, 96, 48, 24, 24, 72, 48, 0, 8, 54, 84, 48, 24, 72, 96, 0, 48, 48, 24, 72, 0, 72, 96, 0, 6, 96, 96, 24, 48, 96, 48, 0, 36, 48, 120

Offset: 0

N. J. A. Sloane

Number of ordered triples (i, j, k) of integers such that n = i^2 + j^2 + k^2.
The Madelung Coulomb energy for alternating unit charges in the simple cubic lattice is Sum_{n>=1} (-1)^n*a(n)/sqrt(n) = -A085469. - R. J. Mathar, Apr 29 2006
a(A004215(k))=0 for k=1,2,3,... but no other elements of {a(n)} are zero. - Graeme McRae, Jan 15 2007

			Order and signs are taken into account: a(1) = 6 from 1 = (+-1)^2 + 0^2 + 0^2, a(2) = 12 from 2 = (+-1)^2 + (+-1)^2 + 0^2; a(3) = 8 from 3 = (+-1)^2 + (+-1)^2 + (+-1)^2, etc.
G.f. =  1 + 6*q + 12*q^2 + 8*q^3 + 6*q^4 + 24*q^5 + 24*q^6 + 12*q^8 + 30*q^9 + 24*q^10 + ...

H. Cohen, Number Theory, Vol. 1: Tools and Diophantine Equations, Springer-Verlag, 2007, p. 317.
J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 107.
H. Davenport, The Higher Arithmetic. Cambridge Univ. Press, 7th ed., 1999, Chapter V.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 3, p. 109.
E. Grosswald, Representations of Integers as Sums of Squares. Springer-Verlag, NY, 1985, p. 54.
L. Kronecker, Crelle, Vol. LVII (1860), p. 248; Werke, Vol. IV, p. 188.
C. J. Moreno and S. S. Wagstaff, Jr., Sums of Squares of Integers, Chapman and Hall, 2006, p. 43.
T. Nagell, Introduction to Number Theory, Wiley, 1951, p. 194.
W. Sierpiński, 1925. Teorja Liczb. pp. 1-410 (p.61).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
H. J. S. Smith, Report on the Theory of Numbers, reprinted in Vol. 1 of his Collected Math. Papers, Chelsea, NY, 1979, see p. 338, Eq. (B').

T. D. Noe, Table of n, a(n) for n = 0..10000
George E. Andrews, The Bhargava-Adiga Summation and Partitions, 2016. See Lemma 2.1.
P. T. Bateman, On the representations of a number as the sum of three squares, Trans. Amer. Math. Soc. 71 (1951), 70-101.
S. Bhargava and C. Adiga, A basic bilateral series summation formula and its applications, Integral Transforms and Special Functions, 2 (1994), 165-184.
J. M. Borwein, K-K S. Choi, On Dirichlet series for sums of squares, Raman. J. 7 (2003) 95-127
S. K. K. Choi, A. V. Kumchev, and R. Osburn, On sums of three squares, arXiv:math/0502007 [math.NT], 2005.
M. Doring, J. Haidenbauer, U.-G. Meissner, and A. Rusetsky, Dynamical coupled-channel approaches on a momentum lattice, arXiv:1108.0676 [hep-lat], 2011.
J. A. Ewell, Recursive determination of the enumerator for sums of three squares, Internat. J. Math. and Math. Sci, 24 (2000), 529-532.
O. Fraser and B. Gordon, On representing a square as the sum of three squares, Amer. Math. Monthly, 76 (1969), 922-923.
M. D. Hirschhorn and J. A. Sellers, On Representations of a Number as a Sum of Three Squares, Discrete Mathematics 199 (1999), 85-101.
M. D. Hirschhorn and J. A. Sellers, On Representations Of A Number As A Sum Of Three Squares
A. Martinez Torres, L. R. Dai, C. Koren, D. Jido, and E. Oset, The KD, eta D_s interaction in finite volume and the D_{s^*0}(2317) resonance, arXiv:1109.0396 [hep-lat], 2011.
R. J. Mathar, Hierarchical Subdivision of the Simple Cubic Lattice, arXiv:1309.3705 [math.MG], 2013.
S. C. Milne, Infinite families of exact sums of squares formulas, Jacobi elliptic functions, continued fractions and Schur functions, Ramanujan J., 6 (2002), 7-149.
J. L. Mordell, The Representation Of Integers By Three Positive Squares, Michigan Math. J. 7(3): 289-290 (1960).
Eric T. Mortenson, A Kronecker-type identity and the representations of a number as a sum of three squares, arXiv:1702.01627 [math.NT], 2017.
G. Nebe and N. J. A. Sloane, Home page for this lattice
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Theta Series
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
Index entries for sequences related to sums of squares

Row d=3 of A122141 and of A319574, 3rd column of A286815.
Cf. A074590 (primitive solutions), A117609 (partial sums), A004215 (positions of zeros).
Analog for 4 squares: A000118.
x^2+y^2+k*z^2: A005875, A014455, A034933, A169783, A169784.
Cf. A000164, A004013, A004015, A008443, A045828, A045831, A045834, A213384.
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

# JacobiTheta3 is defined in A000122.
A005875List(len) = JacobiTheta3(len, 3)
A005875List(75) |> println # Peter Luschny, Mar 12 2018

Basis( ModularForms( Gamma1(4), 3/2), 75) [1]; /* Michael Somos, Jun 25 2014 */

(sum(x^(m^2),m=-10..10))^3; seq(coeff(%,x,n), n=0..50);
Alternative:
A005875list := proc(len) series(JacobiTheta3(0, x)^3, x, len+1);
seq(coeff(%, x, j), j=0..len-1) end: A005875list(75); # Peter Luschny, Oct 02 2018

SquaresR[3,Range[0,80]] (* Harvey P. Dale, Jul 21 2011 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q]^3, {q, 0, n}]; (* Michael Somos, Jun 25 2014 *)
a[ n_] := Length @ FindInstance[ n == x^2 + y^2 + z^2, {x, y, z}, Integers, 10^9]; (* Michael Somos, May 21 2015 *)
QP = QPochhammer; CoefficientList[(QP[q^2]^5/(QP[q]*QP[q^4])^2)^3 + O[q]^80, q] (* Jean-François Alcover, Nov 24 2015 *)

{a(n) = if( n<0, 0, polcoeff( sum( k=1, sqrtint(n), 2 * x^k^2, 1 + x * O(x^n))^3, n))};

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A)^5 / (eta(x + A) * eta(x^4 + A))^2)^3, n))}; /* Michael Somos, Jun 03 2012 */

{a(n) = my(G); if( n<0, 0, G = [ 1, 0, 0; 0, 1, 0; 0, 0, 1]; polcoeff( 1 + 2 * x * Ser( qfrep( G, n)), n))}; /* Michael Somos, May 21 2015 */

# uses Python code for A004018
from math import isqrt
def A005875(n): return A004018(n)+(sum(A004018(n-k**2) for k in range(1,isqrt(n)+1))<<1) # Chai Wah Wu, Jun 21 2024

Q = DiagonalQuadraticForm(ZZ, [1]*3)
Q.representation_number_list(75) # Peter Luschny, Jun 20 2014

A number n is representable as the sum of 3 squares iff n is not of the form 4^a (8k+7) (cf. A000378).
There is a classical formula (essentially due to Gauss):
For sums of 3 squares r_3(n): write (uniquely) -n=D(2^vf)^2, with D<0 fundamental discriminant, f odd, v>=-1. Then r_3(n) = 12L((D/.),0)(1-(D/2)) Sum_{d | f} mu(d)(D/d)sigma(f/d).
Here mu is the Moebius function, (D/2) and (D/d) are Kronecker-Legendre symbols, sigma is the sum of divisors function, L((D/.),0)=h(D)/(w(D)/2) is the value at 0 of the L function of the quadratic character (D/.), equal to the class number h(D) divided by 2 or 3 in the special cases D=-4 and -3. - Henri Cohen (Henri.Cohen(AT)math.u-bordeaux1.fr), May 12 2010
a(n) = 3*T(n) if n == 1,2,5,6 mod 8, = 2*T(n) if n == 3 mod 8, = 0 if n == 7 mod 8 and = a(n/4) if n == 0 mod 4, where T(n) = A117726(n). [Moreno-Wagstaff].
"If 12E(n) is the number of representations of n as a sum of three squares, then E(n) = 2F(n) - G(n) where G(n) = number of classes of determinant -n, F(n) = number of uneven classes." - Dickson, quoting Kronecker. [Cf. A117726.]
a(n) = Sum_{d^2|n} b(n/d^2), where b() = A074590() gives the number of primitive solutions.
Expansion of phi(q)^3 in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Oct 25 2006.
Euler transform of period 4 sequence [ 6, -9, 6, -3, ...]. - Michael Somos, Oct 25 2006
G.f.: (Sum_{k in Z} x^(k^2))^3.
a(8*n + 7) = 0. a(4*n) = a(n).
a(n) = A004015(2*n) = A014455(2*n) = A004013(4*n) = A169783(4*n). a(4*n + 1) = 6 * A045834(n). a(8*n + 3) = 8 * A008443(n). a(8*n + 5) = 24 * A045831(n). - Michael Somos, Jun 03 2012
a(4*n + 2) = 12 * A045828(n). - Michael Somos, Sep 03 2014
a(n) = (-1)^n * A213384(n). - Michael Somos, May 21 2015
a(n) = (6/n)*Sum_{k=1..n} A186690(k)*a(n-k), a(0) = 1. - Seiichi Manyama, May 27 2017
a(n) = A004018(n) + 2*Sum_{k=1..floor(sqrt(n))} A004018(n - k^2). - Daniel Suteu, Aug 27 2021
Convolution cube of A000122. Convolution of A004018 and A000122. - R. J. Mathar, Aug 03 2025

More terms from James Sellers, Aug 22 2000

A004015 Theta series of face-centered cubic (f.c.c.) lattice.

Original entry on oeis.org

1, 12, 6, 24, 12, 24, 8, 48, 6, 36, 24, 24, 24, 72, 0, 48, 12, 48, 30, 72, 24, 48, 24, 48, 8, 84, 24, 96, 48, 24, 0, 96, 6, 96, 48, 48, 36, 120, 24, 48, 24, 48, 48, 120, 24, 120, 0, 96, 24, 108, 30, 48, 72, 72, 32, 144, 0, 96, 72, 72, 48, 120, 0, 144, 12, 48, 48, 168, 48, 96

Offset: 0

N. J. A. Sloane

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

Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

			G.f. = 1 + 12*x + 6*x^2 + 24*x^3 + 12*x^4 + 24*x^5 + 8*x^6 + 48*x^7 + 6*x^8 + ...
G.f. = 1 + 12*q^2 + 6*q^4 + 24*q^6 + 12*q^8 + 24*q^10 + 8*q^12 + 48*q^14 + 6*q^16 + ...
From _Michael Somos_, Jan 05 2012: (Start)
a(2) = 6 since (1, -1, -1) is a solution to x^2 + y^2 + z^2 + x*y + x*z + y*z = 2 and the other 5 solutions are permutations and negations of this one.
a(2) = 6 since (1, 1, -1, -1) is a solution to x + y + z + w = 0 and x^2 + y^2 + z^2 + w^2 = 4 and the other 5 solutions are permutations of this one. (End)

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 113.
L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 263.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
N. J. A. Sloane and B. K. Teo, Theta series and magic numbers for close-packed spherical clusters, J. Chem. Phys. 83 (1985) 6520-6534.
L. V. Woodcock, Nature, Jan 09 1997, pp. 141-143.

N. J. A. Sloane, Table of n, a(n) for n = 0..9999
Charles L. Hohn, a(0) to a(15), animated
G. Nebe and N. J. A. Sloane, Home page for this lattice
N. J. A. Sloane, A portion of the f.c.c. lattice packing.
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
Eric Weisstein's World of Mathematics, Theta Series
Index entries for sequences related to f.c.c. lattice

Cf. A004013, A005875, A005901, A045828. A055039 gives the positions of the 0's in this sequence.

Cf. A000007, A000122, A004016, A008444, A008445, A008446, A008447, A008448, A008449 (Theta series of lattices A_0, A_1, A_2, A_4, ...)

L := Lattice("A",3); A := ThetaSeries(L, 140); A; /* Michael Somos, Nov 13 2014 */

A := Basis( ModularForms( Gamma1(8), 3/2), 70); A[1] + 12*A[2] + 6*A[3] + 24*A[4]; /* Michael Somos, Sep 08 2018 */

maxd := 201: temp0 := trunc(evalf(sqrt(maxd)))+2: a := 0: for i from -temp0 to temp0 do a := a+q^( (i+1/2)^2): od: th2 := series(a,q,maxd); a := 0: for i from -temp0 to temp0 do a := a+q^(i^2): od: th3 := series(a,q,maxd); th4 := series(subs(q=-q, th3),q,maxd); series((1/2)*(th3^3+th4^3),q,200);

a[n_] := SquaresR[3, 2n]; Table[a[n], {n, 0, 69}] (* Jean-François Alcover, Jul 12 2012 *)
a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q]^3 + EllipticTheta[ 4, 0, q]^3) / 2, {q, 0, 2 n}]; (* Michael Somos, May 24 2013 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q^2]^3 + 12 q QPochhammer[ q^4]^3 QPochhammer[ q^8]^2 / QPochhammer[ q^2]^2, {q, 0, n}]; (* Michael Somos, Nov 13 2014 *)
SquaresR[3,2*Range[0,70]] (* Harvey P. Dale, Jun 01 2015 *)

{a(n) = if( n<0, 0, n*=2; polcoeff( sum( k=1, sqrtint(n), 2 * x^k^2, 1 + x * O(x^n))^3, n))}; /* Michael Somos, Oct 25 2006 */

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^4 + A)^5 / eta(x^2 + A)^2 / eta(x^8 + A)^2)^3 + 12 * x * eta(x^4 + A)^3 * eta(x^8 + A)^2 / eta(x^2 + A)^2, n))}; /* Michael Somos, May 17 2008 */

{a(n) = if( n<1, n==0, 2 * qfrep( [2, 1, 1; 1, 2, 1; 1, 1, 2], n, 1)[n])}; /* Michael Somos, Jan 02 2012 */

from math import prod, isqrt
from sympy import factorint
def A004018(n): return prod(1 if p==2 else (e+1 if p&3==1 else (e+1)&1) for p, e in factorint(n).items())<<2 if n else 1
def A004015(n): return A004018(m:=n<<1)+(sum(A004018(m-k**2) for k in range(1,isqrt(m)+1))<<1) # Chai Wah Wu, Feb 24 2025

Expansion of phi(q^2)^3 + 12 * q * phi(q^2) * psi(q^4)^2 in powers of q where phi(), psi() are Ramanujan theta functions. - Michael Somos, Oct 25 2006
Expansion of (phi(q)^3 + phi(-q)^3) / 2 in powers of q^2 where phi() is a Ramanujan theta function. - Michael Somos, Oct 25 2006
Expansion of b(q) * phi(q^18) + c(q^3) * phi(q^2) in powers of q^3 where b(), c() are cubic AGM theta functions and phi() is a Ramanujan theta function. - Michael Somos, Oct 25 2006
Expansion of (theta_3(q)^3 + theta_4(q)^3) / 2 in powers of q^2.
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 2^(7/2) (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A004013.
a(n) = A005875(2*n).
G.f.: Sum_{i, j, k in Z} x^( i*i + j*j + k*k + i*j + i*k + j*k ). - Michael Somos, Jan 02 2012
From Michael Somos, Jan 05 2012: (Start)
Number of integer solutions to x^2 + y^2 + z^2 + x*y + x*z + y*z = n.
Number of integer solutions to x + y + z even and x^2 + y^2 + z^2 = 2 * n.
Number of integer solutions to x + y + z + w = 0 and x^2 + y^2 + z^2 + w^2 = 2 * n. (End)
a(2*n) = A005875(n). a(2*n+1) = 12 * A045828(n). - Michael Somos, Dec 28 2017

A259825 a(n) = 12*H(n) where H() is the Hurwitz class number.

Original entry on oeis.org

-1, 0, 0, 4, 6, 0, 0, 12, 12, 0, 0, 12, 16, 0, 0, 24, 18, 0, 0, 12, 24, 0, 0, 36, 24, 0, 0, 16, 24, 0, 0, 36, 36, 0, 0, 24, 30, 0, 0, 48, 24, 0, 0, 12, 48, 0, 0, 60, 40, 0, 0, 24, 24, 0, 0, 48, 48, 0, 0, 36, 48, 0, 0, 60, 42, 0, 0, 12, 48, 0, 0, 84, 36, 0, 0

Offset: 0

Michael Somos, Jul 05 2015

sign

Coefficients of q-expansion of Eisenstein series G_{3/2}(tau) multiplied by 12. - N. J. A. Sloane, Mar 16 2019

			G.f. = -1 + 4*x^3 + 6*x^4 + 12*x^7 + 12*x^8 + 12*x^11 + 16*x^12 + 24*x^15 + ...

Seiichi Manyama, Table of n, a(n) for n = 0..10000
Kathrin Bringmann and Jeremy Lovejoy, Overpartitions and class numbers of binary quadratic forms, arXiv:0712.0631 [math.NT], 2007. See page 5, equation (1.12).
Archer Clayton, Helen Dai, Tianyu Ni, Erick Ross, Hui Xue, and Jake Zummo, Non-repetition of second coefficients of Hecke polynomials, arXiv:2411.18419 [math.NT], 2024. See p. 21.
Wikipedia, Hurwitz class number.
Don Zagier, Modular Forms of One Variable, Notes based on a course given in Utrecht, 1991. See page 50.

Cf. A004024, A005869, A005876, A005877, A005878, A005884, A005886, A008443.
Cf. A045828, A045831, A045834, A058305, A058306, A130695, A181648, A185220.
Cf. A188569, A213022, A213617, A213624, A213625, A213627, A259827.
See also A306934-A306937.

terms = 100; gf[m_] := With[{r = Range[-m, m]}, -2 Sum[(-1)^k*x^(k^2 + k)/(1 + (-x)^k)^2, {k, r}]/EllipticTheta[3, 0, x] - 2 Sum[(-1)^k*x^(k^2 + 2 k)/(1 + x^(2 k))^2, {k, r}]/EllipticTheta[3, 0, -x]]; gf[terms // Sqrt // Ceiling] + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Apr 02 2017 *)
a[ n_] := If[ n<1, -Boole[n==0], With[{m = Floor[(-1 + Sqrt[1 + 4*n])/2]}, -2*SeriesCoefficient[ Sum[(-1)^k*x^(k^2 + k)/(1 + (-x)^k)^2, {k, -m-1,m}] / EllipticTheta[3, 0, x] + Sum[(-1)^k*x^(k^2 + 2*k)/(1 + x^(2*k))^2, {k, -m-2,m}]/ EllipticTheta[3, 0, -x], {x, 0, n}]]]; (* Michael Somos, Feb 04 2022 *)

PARI
```
{a(n) = 12 * qfbhclassno(n)};
```

{a(n) = my(D, f); 12 * if( n<1, (n==0)/-12, [D, f] = core(-n, 1); if( D%4>1 && !(f%2), D*=4; f/=2); if( D%4<2, qfbclassno(D) / max(1, D+6), 0) * sumdiv(f, d, moebius(d) * kronecker(D, d) * sigma(f/d)))};

a(n) = 12 * A058305(n) / A058306(n). a(4*n + 1) = a(4*n + 2) = 0. a(3*n + 4) = 6 * A259827(n).
a(4*n + 3) = 4 * A130695(n). a(8*n + 3) = A005886(n) = 2 * A005869(n) = 4 * A008443(n). a(12*n + 7) = 12 * A259655(n).
a(16*n + 4) = 6 * A045834(n) = 3 * A005876(n). a(16*n + 8) = 12 * A045828(n) = 6 * A005884(n) = 3 * A005877(n).
a(24*n + 3) = 4 * A213627(n). a(24*n + 7) = 12 * A185220(n). a(24*n + 11) = 12 * A213617(n). a(24*n + 19) = 12 * A181648(n). a(24*n + 23) = 12 * A188569(n+1).
a(32*n + 4) = 6 * A213022(n). a(32*n + 8) = 12 * A213625(n). a(32*n + 12) = 16 * A008443(n) = 8 * A005869(n) = 4 * A005886(n) = 2 * A005878(n). a(32*n + 20) = 24 * A045831(n) = 6 * A004024(n). a(32*n + 24) = 24 * A213624(n).
G.f.: -2 * (Sum_{k in Z} (-1)^k * x^(k*k + k) / (1 + (-x)^k)^2) / (Sum_{k in Z} x^k^2) - 2 * (Sum_{k in Z} (-1)^k * x^(k^2 + 2*k) / (1 + x^(2*k))^2) / (Sum_{k in Z} (-x)^k^2).
a(n) >= 0 if n > 0. - Michael Somos, Feb 04 2022

A014455 Theta series of quadratic form with Gram matrix [ 1, 0, 0; 0, 1, 0; 0, 0, 2 ]. Number of integer solutions to x^2 + y^2 + 2*z^2 = n.

Original entry on oeis.org

1, 4, 6, 8, 12, 8, 8, 16, 6, 12, 24, 8, 24, 24, 0, 16, 12, 16, 30, 24, 24, 16, 24, 16, 8, 28, 24, 32, 48, 8, 0, 32, 6, 32, 48, 16, 36, 40, 24, 16, 24, 16, 48, 40, 24, 40, 0, 32, 24, 36, 30, 16, 72, 24, 32, 48, 0, 32, 72, 24, 48, 40, 0, 48, 12, 16, 48, 56, 48, 32, 48, 16, 30, 64

Offset: 0

N. J. A. Sloane

This is the tetragonal P lattice (the classical holotype) of dimension 3.

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

			G.f. = 1 + 4*q + 6*q^2 + 8*q^3 + 12*q^4 + 8*q^5 + 8*q^6 + 16*q^7 + 6*q^8 + 12*q^9 + ...

John Cannon, Table of n, a(n) for n = 0..10000
G. Nebe and N. J. A. Sloane, Home page for this lattice
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

Cf. A004015, A005875, A005877, A008443, A033763, A045826, A045828, A156384, A213024, A246631.

A := Basis( ModularForms( Gamma1(8), 3/2), 40); A[1] + 4*A[2] + 6*A[3] + 8*A[4]; /* Michael Somos, Aug 31 2014 */

r[n_, z_] := Reduce[x^2 + y^2 + 2*z^2 == n, {x, y}, Integers]; a[n_] := Module[{rn0, rnz, k0, k}, rn0 = r[n, 0]; k0 = If[rn0 === False, 0, If[Head[rn0] === And, 1, Length[rn0]]]; For[k = 0; z = 1, z <= Ceiling[Sqrt[n/2]], z++, rnz = r[n, z]; If[rnz =!= False, k = If[Head[rnz] === And, k+1, k + Length[rnz]]]]; k0 + 2*k]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Oct 07 2013 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q]^2 EllipticTheta[ 3, 0, q^2], {q, 0, n}]; (* Michael Somos, Aug 31 2014 *)
QP = QPochhammer; s = QP[q^2]^8*(QP[q^4]/(QP[q]^4*QP[q^8]^2)) + O[q]^80; CoefficientList[s, q] (* Jean-François Alcover, Nov 25 2015, after Michael Somos *)

{a(n) = if( n<1, n==0, 2 * qfrep([ 1, 0, 0; 0, 1, 0; 0, 0, 2], n)[n])}; /* Michael Somos, Jul 05 2005 */

{a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^8 * eta(x^4 + A) / (eta(x + A)^4 * eta(x^8 + A)^2), n))}; /* Michael Somos, Jul 05 2005 */

Expansion of phi(q)^2 * phi(q^2) = psi(q)^4 / psi(q^4) in powers of q where phi(), psi() are Ramanujan theta functions. - Michael Somos, Apr 07 2012
Expansion of eta(q^2)^8 * eta(q^4) / (eta(q)^4 * eta(q^8)^2) in powers of q. - Michael Somos, Jul 05 2005
Euler transform of period 8 sequence [4, -4, 4, -5, 4, -4, 4, -3, ...]. - Michael Somos, Jul 07 2005
G.f.: theta_3(q)^2 * theta_3(q^2) = Product_{k>0} (1 - x^(2*k))^8 * (1 - x^(4*k)) / ((1 - x^k)^4 * (1 - x^(8*k))^2).
There is a classical formula (essentially due to Gauss): Write (uniquely) -2n=D(2^vf)^2, with D<0 fundamental discriminant, f odd, v>=-1. Then a(n)=12L((D/.),0)(1-(D/2))\sum_{d\mid f}\mu(d)(D/d)sigma(f/d) (the formula for A005875), except that the factor (1-(D/2)) has to be replaced by 1/3 if v=-1 and by 1 if v=0 (and kept if v>=1). Here mu() is the Moebius function, (D/2) and (D/d) are Kronecker-Legendre symbols, sigma() is the sum of divisors function, L((D/.),0)=h(D)/(w(D)/2) is the value at 0 of the L() function of the quadratic character (D/.), equal to the class number h(D) divided by 2 or 3 in the special cases D=-4 and -3. - Henri Cohen (Henri.Cohen(AT)math.u-bordeaux1.fr), May 12 2010
a(2*n) = a(8*n) = A005875(n). a(2*n + 1) = A005877(n) = 4 * A045828(n). a(4*n) = A004015(n). a(4*n + 2) = 2 * A045826(n). a(8*n + 4) = 12 * A045828(n). a(8*n + 7) = 16 * A033763(n). a(16*n + 6) = 8 * A008443(n). a(16*n + 14) = 0. - Michael Somos, Apr 07 2012
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 32^(1/2) (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A246631.

A004013 Theta series of body-centered cubic (b.c.c.) lattice.

Original entry on oeis.org

1, 0, 0, 8, 6, 0, 0, 0, 12, 0, 0, 24, 8, 0, 0, 0, 6, 0, 0, 24, 24, 0, 0, 0, 24, 0, 0, 32, 0, 0, 0, 0, 12, 0, 0, 48, 30, 0, 0, 0, 24, 0, 0, 24, 24, 0, 0, 0, 8, 0, 0, 48, 24, 0, 0, 0, 48, 0, 0, 72, 0, 0, 0, 0, 6, 0, 0, 24, 48, 0, 0, 0, 36, 0, 0, 56, 24, 0, 0, 0, 24, 0, 0, 72, 48, 0, 0, 0, 24, 0, 0

Offset: 0

N. J. A. Sloane

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

The bcc lattice is also called the coweight lattice A^*3 (or D^*_3), dual to the root lattice A_3 (or D_3, or the fcc lattice), or the permutohedral lattice A^*_3. - _Andrey Zabolotskiy, Mar 08 2020

			G.f. = 1 + 8*x^3 + 6*x^4 + 12*x^8 + 24*x^11 + 8*x^12 + 6*x^16 + 24*x^19 + 24*x^20 + ...
G.f. = 1 + 8*q^(3/2) + 6*q^2 + 12*q^4 + 24*q^(11/2) + 8*q^6 + 6*q^8 + 24*q^(19/2) + 24*q^10 + 24*q^12 + 32*q^(27/2) + ...

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 116.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

John Cannon, Table of n, a(n) for n = 0..10000
S. Ahlgren, The sixth, eighth, ninth and tenth powers of Ramanujan's theta function, Proc. Amer. Math. Soc., 128 (1999), 1333-1338; F_4(q).
R. J. Mathar, Hierarchical Subdivision of the Simple Cubic Lattice, arXiv preprint arXiv:1309.3705 [math.MG], 2013.
G. Nebe and N. J. A. Sloane, Home page for this lattice
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
Eric Weisstein's World of Mathematics, Theta Series
Index entries for sequences related to b.c.c. lattice

Cf. A004015, A005875, A008443, A045826, A045828, A045834, A213056.
Indices of nonzero terms are A004014.
Cf. A023916-A023936.
Cf. A004011, A008422, A008425, A008423, A008427, A008424, A008426.

Basis( ModularForms( Gamma0(8), 3/2), 90) [1]; /* Michael Somos, Sep 04 2014 */

M:=100; M1:=M*(M+1)/2; ph:=series(add(q^(k^2),k=-M..M),q,M1): ps:=series(add(q^(k*(k+1)/2),k=0..M),q,M1): t1:=series(subs(q=q^2, ph)^3, q,M1): t2:=series((2*sqrt(q))^3*subs(q=q^4, ps)^3,q,M1): t3:=seriestolist(series(subs(q=q^2,t1+t2),q,M1)): for n from 0 to nops(t3)-1 do lprint(n,t3[n+1]); od:

m = 13; m1 = m*((m + 1)/2); ph[q_] = Series[ Sum[ q^k^2, {k, -m, m}], {q, 0, m1}]; ps[q_] = Series[ Sum[ q^(k*((k + 1)/2)), {k, 0, m}], {q, 0, m1}]; t1[q_] = Normal[ Series[ ph[q^2]^3, {q, 0, m1}]]; t2[q_] = Normal[ Series[ (2*Sqrt[q])^3*ps[q^4]^3, {q, 0, m1}]]; CoefficientList[ Series[ t1[q^2] + t2[q^2], {q, 0, m1}], q] (* Jean-François Alcover, Dec 20 2011, translated from Maple *)
(* From version 6 on *) terms=91; f[q_] = LatticeData["BodyCenteredCubic", "ThetaSeriesFunction"][-I Log[q]/Pi]; CoefficientList[Simplify[f[q] + O[q]^terms, q>0], q][[1 ;; terms]] (* Jean-François Alcover, May 15 2013, updated Jul 08 2017 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, x^4]^3 + EllipticTheta[ 2, 0, x^4]^3, {x, 0, n}]; (* Michael Somos, May 24 2013 *)

{a(n) = if( n<0, 0, if( n%4==0, n/=4; polcoeff( sum(k=1, sqrtint(n), 2*x^k^2, 1 + x * O(x^n))^3, n), n%8==3, n\=8; 8*polcoeff( sum(k=0, (sqrtint(8*n+1) - 1)\2, x^((k^2 + k)/2), x * O(x^n))^3, n)))}; /* Michael Somos, Oct 25 2006 */

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^8 + A)^5 / eta(x^4 + A)^2 / eta(x^16 + A)^2)^3 + (2 * x * eta(x^16 + A)^2 / eta(x^8 + A))^3, n))}; /* Michael Somos, May 17 2008 */

subs(q=q^2, ph)^3+(2*sqrt(q))^3*subs(q=q^4, ps)^3, where ps = A010054 = Sum_{k=0..infinity} q^(k*(k+1)/2), ph = A000122 = Sum_{k=-infinity, infinity} q^(k^2).
Expansion of phi(q^4)^3 + 8 * q^3 * psi(q^8)^3 in powers of q where phi(), psi() are Ramanujan theta functions. - Michael Somos, Oct 25 2006
a(4*n + 1) = a(4*n + 2) = a(8*n + 7) = 0. a(4*n) = A005875(n).
Expansion of theta_3(q)^3 + theta_2(q)^3 in powers of q^(1/4).
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 2 (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A004015.
a(8*n) = A004015(n). a(8*n + 3) = 8 * A008443(n). a(8*n + 4) = 2 * A045826(n). - Michael Somos, Jul 19 2015
a(12*n + 4) = 6 * A213056(n). a(16*n + 4) = 6 * A045834(n). a(16*n + 8) = 12 * A045828(n).

A213384 Expansion of phi(-q)^3 in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -6, 12, -8, 6, -24, 24, 0, 12, -30, 24, -24, 8, -24, 48, 0, 6, -48, 36, -24, 24, -48, 24, 0, 24, -30, 72, -32, 0, -72, 48, 0, 12, -48, 48, -48, 30, -24, 72, 0, 24, -96, 48, -24, 24, -72, 48, 0, 8, -54, 84, -48, 24, -72, 96, 0, 48, -48, 24, -72, 0, -72, 96

Offset: 0

Michael Somos, Jun 10 2012

sign

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

			G.f. = 1 - 6*q + 12*q^2 - 8*q^3 + 6*q^4 - 24*q^5 + 24*q^6 + 12*q^8 - 30*q^9 + ...

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..2500 from G. C. Greubel)
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

Cf. A004015, A005875, A008443, A045826, A045834.

# JacobiTheta4 is defined in A002448.
A213384List(len) = JacobiTheta4(len, 3)
A213384List(63) |> println # Peter Luschny, Mar 12 2018

A := Basis( ModularForms( Gamma0(16), 3/2), 63); A[1] - 6*A[2] + 12*A[3] - 8*A[4]; /* Michael Somos, May 21 2015 */

a[ n_] := (-1)^n SquaresR[ 3, n]; (* Michael Somos, May 21 2015 *)
a[ n_] := (-1)^n Length @ FindInstance[ n == x^2 + y^2 + z^2, {x, y, z}, Integers, 10^9]; (* Michael Somos, May 21 2015 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 4, 0, q]^3, {q, 0, n}]; (* Michael Somos, May 21 2015 *)
a[ n_] := SeriesCoefficient[ (QPochhammer[ q]^2 / QPochhammer[ q^2])^3, {q, 0, n}]; (* Michael Somos, May 21 2015 *)

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

{a(n) = if( n<0, 0, polcoeff( sum( k=1, sqrtint(n), 2 * (-x)^k^2, 1 + x * O(x^n))^3, n))}; /* Michael Somos, May 21 2015 */

{a(n) = my(G); if( n<0, 0, G = [ 1, 0, 0; 0, 1, 0; 0, 0, 1]; (-1)^n * polcoeff( 1 + 2 * x * Ser( qfrep( G, n)), n))}; /* Michael Somos, May 21 2015 */

Expansion of (eta(q)^2 / eta(q^2))^3 in powers of q.
Euler transform of period 2 sequence [ -6, -3, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 2^(15/2) (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A008443.
G.f.: (Sum_{k in Z} (-1)^k * x^k^2)^3.
a(n) = (-1)^n * A005875(n). a(2*n) = A004015(n). a(2*n + 1) = -2 * A045826(n). a(4*n) = A005875(n). a(4*n + 1) = -6 * A045834(n). a(4*n + 2) = 12 * A045828(n). a(8*n + 3) = -8 * A008443(n). a(8*n + 7) = 0.

A212885 Expansion of phi(q) * phi(-q)^2 in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

1, -2, -4, 8, 6, -8, -8, 0, 12, -10, -8, 24, 8, -8, -16, 0, 6, -16, -12, 24, 24, -16, -8, 0, 24, -10, -24, 32, 0, -24, -16, 0, 12, -16, -16, 48, 30, -8, -24, 0, 24, -32, -16, 24, 24, -24, -16, 0, 8, -18, -28, 48, 24, -24, -32, 0, 48, -16, -8, 72, 0, -24, -32

Offset: 0

Michael Somos, May 29 2012

sign

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

			G.f. = 1 - 2*q - 4*q^2 + 8*q^3 + 6*q^4 - 8*q^5 - 8*q^6 + 12*q^8 - 10*q^9 + ...

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

Cf. A004015, A004024, A005875, A005877, A005878, A005887, A008443, A045828, A045834, A213624.

a[n_]:= SeriesCoefficient[EllipticTheta[3, 0, q]* EllipticTheta[3, 0, -q]^2, {q, 0, n}]; Table[a[n], {n, 0, 50}] (* G. C. Greubel, Nov 30 2017 *)

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

Expansion of phi(-x) * phi(-x^2)^2 = phi(-x^2)^4 / phi(x) in powers of x where phi() is a Ramanujan theta function.
Expansion of eta(q^2)^3 * eta(q)^2 / eta(q^4)^2 in powers of q.
Euler transform of period 4 sequence [-2, -5, -2, -3, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 32 (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A045828.
G.f.: Product_{k>0} (1 - x^(2*k))^3 * (1 - x^k)^2 / (1 - x^(4*k))^2.
a(4*n) = A005875(n). a(4*n + 1) = -2 * A045834(n). a(4*n + 2) = - A005877(n) = -4 * A045828(n).
a(8*n) = A004015(n). a(8*n + 3) = A005878(n) = 8 * A008443(n). a(8*n + 4)= A005887(n). a(8*n + 5) = -2 * A004024(n). a(8*n + 6) = -8 * A213624(n). a(8*n + 7) = 0.

A005877 Theta series of cubic lattice with respect to square.

Original entry on oeis.org

4, 8, 8, 16, 12, 8, 24, 16, 16, 24, 16, 16, 28, 32, 8, 32, 32, 16, 40, 16, 16, 40, 40, 32, 36, 16, 24, 48, 32, 24, 40, 48, 16, 56, 32, 16, 64, 40, 32, 32, 36, 40, 48, 48, 32, 48, 48, 16, 80, 40, 24, 80

Offset: 0

N. J. A. Sloane

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 107.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

G. C. Greubel, Table of n, a(n) for n = 0..1000
G. Nebe and N. J. A. Sloane, Home page for this lattice

Cf. A045828.

terms = 52; QP = QPochhammer; s = 4*(QP[q^2]^3*QP[q^4]^2)/QP[q]^2 + O[q]^terms; CoefficientList[s, q] (* Jean-François Alcover, Jul 05 2017 *)

a(n) = 4*A045828(n).

A033717 Number of integer solutions to the equation x^2 + 2y^2 + 4z^2 = n.

Original entry on oeis.org

1, 2, 2, 4, 4, 4, 8, 8, 6, 6, 8, 4, 8, 12, 0, 8, 12, 8, 10, 12, 8, 8, 24, 8, 8, 14, 8, 16, 16, 4, 0, 16, 6, 16, 16, 8, 12, 20, 24, 8, 24, 8, 16, 20, 8, 20, 0, 16, 24, 18, 10, 8, 24, 12, 32, 24, 0, 16, 24, 12, 16, 20, 0, 24, 12, 8, 16, 28, 16, 16, 48, 8, 30, 32, 8, 20, 24, 16, 0, 16, 24, 18

Offset: 0

N. J. A. Sloane

nonn

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

			G.f. = 1 + 2*q + 2*q^2 + 4*q^3 + 4*q^4 + 4*q^5 + 8*q^6 + 8*q^7 + 6*q^8 + ...

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p 102 eq 9.

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

Cf. A004015, A008443, A014455, A033763, A045826, A045828, A045831, A045834, A213022, A213622, A213624, A213625, A246811.

A := Basis( ModularForms( Gamma1(16), 3/2), 82); A[1] + 2*A[2] + 2*A[3] + 4*A[4] + 4*A[5] + 4*A[6] + 8*A[7] + 8*A[8] + 6*A[9] + 8*A[10] + 4*A[11]; /* Michael Somos, Sep 03 2014 */

a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^2] EllipticTheta[ 3, 0, q^4], {q, 0, n}]; (* Michael Somos, Sep 03 2014 *)

{a(n) = my(G); if( n<0, 0, G = [1, 0, 0; 0, 2, 0; 0, 0, 4]; polcoeff( 1 + 2 * x * Ser(qfrep( G, n)), n))}; /* Michael Somos, Sep 03 2014 */

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^3 * eta(x^4 + A) * eta(x^8 + A)^3 / (eta(x + A)^2 * eta(x^16 + A)^2), n))}; /* Michael Somos, Sep 03 2014 */

Expansion of phi(q) * phi(q^2) * phi(q^4) in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Sep 03 2014
Euler transform of period 16 sequence [2, -1, 2, -2, 2, -1, 2, -5, 2, -1, 2, -2, 2, -1, 2, -3, ...]. - Michael Somos, Sep 03 2014
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 8 (t/i)^(3/2) f(t) where q = exp(2 Pi i t). - Michael Somos, Sep 03 2014
a(2*n + 1) = 2 * A045828(n). a(4*n) = A014455(n). a(4*n + 1) = 2 * A213625(n). a(4*n + 2) = 2 * A246811(n). a(4*n + 3) = 4 * A213624(n). - Michael Somos, Sep 03 2014
a(8*n) = A005875(n). a(8*n + 1) = 2 * A213622(n). a(8*n + 2) = 2 * A045834(n). a(8*n + 7) = 8 * A033763(n). - Michael Somos, Sep 03 2014
a(16*n) = A004015(n). a(16*n + 2) = 2 * A213022(n). a(16*n + 6) = 8 *
A008443(n). a(16*n + 8) = 2 * A045826(n). a(16*n + 10) = 8 * A045831(n). a(16*n + 14) = 0. - Michael Somos, Sep 03 2014
G.f.: theta_3(q) * theta_3(q^2) * theta_3(q^4).

cp's OEIS Frontend

A005875 Theta series of simple cubic lattice; also number of ways of writing a nonnegative integer n as a sum of 3 squares (zero being allowed).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Julia

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A004015 Theta series of face-centered cubic (f.c.c.) lattice.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Formula

A259825 a(n) = 12*H(n) where H() is the Hurwitz class number.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A014455 Theta series of quadratic form with Gram matrix [ 1, 0, 0; 0, 1, 0; 0, 0, 2 ]. Number of integer solutions to x^2 + y^2 + 2*z^2 = n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A004013 Theta series of body-centered cubic (b.c.c.) lattice.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

A033717 Number of integer solutions to the equation x^2 + 2y^2 + 4z^2 = n.