A005875 -id:A005875 - cp's OEIS frontend

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.

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.

A224442 Numbers that are the sum of three squares (square 0 allowed) in exactly two ways.

Original entry on oeis.org

9, 17, 18, 25, 26, 27, 29, 33, 34, 36, 38, 45, 49, 51, 53, 57, 59, 61, 62, 68, 69, 72, 73, 75, 77, 82, 83, 85, 94, 97, 100, 102, 104, 105, 106, 107, 108, 109, 116, 118, 123, 130, 132, 136, 138, 139, 141, 144, 147, 152, 154, 155, 157, 158, 165, 177, 180, 187

Offset: 1

Wolfdieter Lang, Apr 08 2013

nonn

These are the numbers for which A000164(a(n)) = 2.
a(n) is the n-th largest number which has a representation as sum of three integer squares (square 0 allowed), in exactly two ways, if neither the order of terms nor the signs of the numbers to be squared are taken into account. The multiplicity for a(n) with order and signs taken into account is A005875(a(n)).
This sequence is a proper subsequence of A000378.

			a(1) = 9 = 0^2 + 0^2 + 3^2 = 1^2 + 2^2 + 2^2, and 9 is the smallest number m with A000164(m) = 2 for m >= 0. The multiplicity with order and signs taken into account is 2*3 + 8*3 = 30 = A005875(9).
The two representations [a,b,c] for a(n), n = 1, ..., 10, are
n=1,   9 = [0, 0, 3], [1, 2, 2],
n=2,  17 = [0, 1, 4], [2, 2, 3],
n=3,  18 = [0, 3, 3], [1, 1, 4],
n=4,  25 = [0, 0, 5], [0, 3, 4],
n=5,  26 = [0, 1, 5], [1, 3, 4],
n=6,  27 = [1, 1, 5], [3, 3, 3],
n=7,  29 = [0, 2, 5], [2, 3, 4],
n=8,  33 = [1, 4, 4], [2, 2, 5],
n=9,  34 = [0, 3, 5], [3, 3, 4],
n=10, 36 = [0, 0, 6], [2, 4, 4].

Alois P. Heinz, Table of n, a(n) for n = 1..1000

Cf. A000164, A005875, A000378, A094942 (one way), A224443 (three ways).

b:= proc(n, i, t) option remember; `if`(n=0, 1, `if`(t*i^22, 3, min(3, b(n, i-1, t)+
      `if`(i^2>n, 0, b(n-i^2, i, t-1))))))
    end:
a:= proc(n) option remember; local k;
      for k from 1 +`if`(n=1, 0, a(n-1))
      while b(k, isqrt(k), 3)<>2 do od; k
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, Apr 09 2013

Select[ Range[0, 200], Length[ PowersRepresentations[#, 3, 2]] == 2 &] (* Jean-François Alcover, Apr 09 2013 *)

This sequence gives the increasingly ordered elements of the set {m integer | m = a^2 + b^2 + c^2, a, b and c integers with 0 <= a <= b <= c, and m has exactly two such representation}.

The sequence gives the increasingly ordered members of the set {m integer | A000164(m) = 2, m >= 0}.

A236933 Number of integer solutions to a^2 + b^2 + c^2 + 12*d^2 = n.

Original entry on oeis.org

1, 6, 12, 8, 6, 24, 24, 0, 12, 30, 24, 24, 10, 36, 72, 16, 18, 96, 84, 24, 48, 108, 72, 48, 40, 78, 168, 32, 12, 168, 120, 48, 60, 144, 96, 48, 78, 84, 216, 64, 24, 240, 144, 24, 48, 168, 144, 96, 70, 114, 252, 64, 84, 312, 240, 48, 120, 252, 168, 120, 32

Offset: 0

N. J. A. Sloane, Feb 16 2014

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..10000

Different from A005875.

For number of solutions to a^2+b^2+c^2+k*d^2=n for k=1,2,3,4,5,6,7,8,12 see A000118, A236928, A236926, A236923, A236930, A236931, A236932, A236927, A236933.

Maple
```
See A236924.
```

G.f.: theta_3(q)^3*theta_3(q^12), where theta_3() is the Jacobi theta function. - Ilya Gutkovskiy, Aug 03 2018

A016725 Number of integer solutions to x^2+y^2+z^2 = n^2, allowing zeros and distinguishing signs and order.

Original entry on oeis.org

1, 6, 6, 30, 6, 30, 30, 54, 6, 102, 30, 78, 30, 78, 54, 150, 6, 102, 102, 126, 30, 270, 78, 150, 30, 150, 78, 318, 54, 174, 150, 198, 6, 390, 102, 270, 102, 222, 126, 390, 30, 246, 270, 270, 78, 510, 150, 294, 30, 390, 150, 510, 78, 318, 318, 390, 54, 630, 174, 366

Offset: 0

csvcjld(AT)nomvst.lsumc.edu

Hurwitz found a formula for a(n). See the paper by Olds.

			1 + 6*x + 6*x^2 + 30*x^3 + 6*x^4 + 30*x^5 + 30*x^6 + 54*x^7 + 6*x^8 + ...

T. D. Noe, Table of n, a(n) for n = 0..10000
Michael Gilleland, Some Self-Similar Integer Sequences
Werner Hürlimann, Exact and Asymptotic Evaluation of the Number of Distinct Primitive Cuboids, Journal of Integer Sequences, Vol. 18 (2015), Article 15.2.5.
Jean Lagrange, Décomposition d'un entier en somme de carrés et fonction multiplicative, Séminaire Delange-Pisot-Poitou. Théorie des nombres, 14 no. 1 (1972-1973), Exp. No. 1, 5 p.
C. D. Olds, On the representations, N_3(n^2), Bull. Amer. Math. Soc. 47 (1941), 499-503.
Eric Weisstein's World of Mathematics, Sum of Squares Function

Cf. A005875.

Column k=3 of A302996.

for n from 0 to 60 do s:=0: for x from -n to n do for y from -n to n do for z from -n to n do if (x^2+y^2+z^2) = n^2 then s:=s+1 fi od od od: printf("%d, ",s) od: # C. Ronaldo (aga_new_ac(AT)hotmail.com), Dec 13 2004

Mathematica
```
SquaresR[3, Range[0,100]^2]
```

{a(n) = if( n<1, n==0, polcoeff( sum( k=1, n, 2 * x^k^2, 1 + x * O(x^n^2))^3, n^2))} /* Michael Somos, Nov 18 2011 */

{a(n) = local(A, p, e); if( n<1, n==0, A = factor(n); 6 * prod( k=1, matsize(A)[1], if( p = A[k, 1], e = A[k, 2]; if( p==2, 1, p^e + if( p%4 == 1, 0, 2 * (p^e - 1) / (p - 1))))))} /* Michael Somos, Nov 18 2011 */

a(n) = 6 * b(n) if n>0 where b(n) is multiplicative with b(2^e) = 1, b(p^e) = p^e if p == 1 (mod 4), b(p^e) = p^e + 2 * (p^e - 1) / (p - 1) if p == 3 (mod 4). - Michael Somos, Nov 18 2011
a(n) = A005875(n^2).
a(n) = [x^(n^2)] theta_3(x)^3, where theta_3() is the Jacobi theta function. - Ilya Gutkovskiy, Apr 20 2018

Revised description from C. Ronaldo (aga_new_ac(AT)hotmail.com), Dec 13 2004

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).

A213022 Expansion of phi(x)^2 * psi(x) in powers of x where phi(), psi() are Ramanujan theta functions.

Original entry on oeis.org

1, 5, 8, 5, 8, 16, 9, 8, 16, 8, 17, 24, 8, 16, 16, 13, 24, 16, 16, 24, 32, 13, 8, 32, 8, 24, 40, 16, 25, 24, 24, 24, 32, 16, 16, 40, 17, 32, 32, 16, 40, 48, 16, 16, 32, 21, 48, 32, 16, 24, 40, 32, 24, 56, 24, 45, 40, 16, 32, 24, 32, 40, 48, 16, 32, 64, 25, 24

Offset: 0

Michael Somos, Jun 03 2012

nonn

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

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(0) = 1 since the norm squared of point [0, 0, 0] with respect to [0, 0, 1/4] is 1/16 = 1/16 + 1/2*0.
a(1) = 5 since the norm squared of points [-1/2, -1/2, -1/2], [-1/2, 1/2, -1/2], [0, 0, -1], [1/2, -1/2, -1/2], [1/2, 1/2, -1/2] with respect to [0, 0, 1/4] is 9/16 = 1/16 + 1/2*1.
1 + 5*x + 8*x^2 + 5*x^3 + 8*x^4 + 16*x^5 + 9*x^6 + 8*x^7 + 16*x^8 + 8*x^9 + ...
q + 5*q^9 + 8*q^17 + 5*q^25 + 8*q^33 + 16*q^41 + 9*q^49 + 8*q^57 + 16*q^65 + ...

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. A005875.

CoefficientList[QPochhammer[q^2]^12/(QPochhammer[q]^5*QPochhammer[q^4]^4) + O[q]^70, q] (* Jean-François Alcover, Nov 05 2015 *)

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

Expansion of q^(-1/8) * eta(q^2)^12 / (eta(q)^5 * eta(q^4)^4) in powers of q.
Expansion of q^(-1/16) times theta series of b.c.c. lattice with respect to point [0, 0, 1/4] in powers of q^(1/2).
Euler transform of period 4 sequence [ 5, -7, 5, -3, ...].
6 * a(n) = A005875(8*n + 1).

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.

A094739 Numbers m such that 4^k*m, for integer k >= 0, are numbers having a unique partition into three squares.

Original entry on oeis.org

1, 2, 3, 5, 6, 10, 11, 13, 14, 19, 21, 22, 30, 35, 37, 42, 43, 46, 58, 67, 70, 78, 91, 93, 115, 133, 142, 163, 190, 235, 253, 403, 427

Offset: 1

T. D. Noe, May 24 2004

Lehmer's paper has an erroneous version of this sequence. He omits 163 and includes 162 (which has 4 partitions) and 182 (which has 3 partitions). Lehmer conjectures that there are no more terms. Note that squares are allowed to be zero.
From Wolfdieter Lang, Aug 27 2020: (Start)
Another name is: Integers not divisible by 4 that are uniquely represented as x^2 + y^2 + z^2 with integers 0 <= x <= y <= z.
This sequence of 33 numbers is complete. See Arno, Theorem 8, p. 332, where 19 is missing, as observed by Kaplansky, Remark 2.1. (a) - (c), p. 87.
All positive integers represented uniquely as sum of three squares of nonnegative numbers, ignoring order and signs, are given by 4^k*a(n), for integer k >= 0 and n = 1 .. 33. See Arno, also p. 322, with some known results, and Kaplansky's Remark 2.1.(c). (End)

			The unique partitions of m*4^k into three squares are,
for m = 1:
1 = 1^2 + 0^2 + 0^2;
4 = 2^2 + 0^2 + 0^2;
16 = 4^2 + 0^2 + 0^2;
...
for m = 163:
163 = 9^2 + 9^2 + 1^2;
163*4 = 18^2 + 18^2 + 2^2;
163*16 = 36^2 + 36^2 + 4^2;
...

Steven Arno, The imaginary quadratic fields of class number 4, Acta Arithmetica 60.4 (1992) 321 - 334.
Irving Kaplansky, Integers Uniquely Represented by Certain Ternary Forms, in "The Mathematics of Paul Erdős I", Ronald. L. Graham and Jaroslav Nešetřil (Eds.), Springer, 1997, pp. 86 - 94.
D. H. Lehmer, On the Partition of Numbers into Squares, The American Mathematical Monthly, Vol. 55, No.8, October 1948, pp. 476-481.

Cf. A005875 (number of ways of writing n as the sum of three squares), A094740 (n having a unique partition into three positive squares).

lim=100; nLst=Table[0, {lim^2}]; Do[n=a^2+b^2+c^2; If[n>0 && n0&]

Keyword full added by Wolfdieter Lang, Aug 27 2020

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.

A290733 Number of compact partitions of n where each partition is counted with a certain weight.

Original entry on oeis.org

0, -1, 2, -1, 0, -3, 3, 2, 0, -3, 1, -2, -1, 0, 5, 3, -2, -4, 1, -2, 1, -3, -1, 4, 2, 1, 6, -3, -3, -6, 1, 2, -2, -1, 2, -4, 3, 4, 4, 3, 2, -8, -1, -2, -1, -4, 0, 4, -2, -1, 4, -3, 3, 0, 7, 1, 3, 2, -6, -6, -5, -4, 4, 2, -2

Offset: 0

N. J. A. Sloane, Aug 10 2017

sign

See Andrews (2016) for the definition of the particular weight used here.

4*a(n) + 2*A290734(n) = (-1)^n*A005875(n) for n > 0.

George E. Andrews, The Bhargava-Adiga Summation and Partitions, 2016. See Lemma 2.1.

Cf. A005875, A290734.

M:=101;
B:=proc(a,q,n) local j,t1; global M;
t1:=1;
for j from 0 to M do
t1:=t1*(1-a*q^j)/(1-a*q^(n+j));
od;
t1; end;
# c_0
t2:=add((-1)^m*q^m*B(-q,q,m-1)/(1+q^m), m=1..M):
series(t2,q,M);
seriestolist(%);

See Maple program for g.f.

cp's OEIS Frontend

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

A224442 Numbers that are the sum of three squares (square 0 allowed) in exactly two ways.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A236933 Number of integer solutions to a^2 + b^2 + c^2 + 12*d^2 = n.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Formula

A016725 Number of integer solutions to x^2+y^2+z^2 = n^2, allowing zeros and distinguishing signs and order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Formula

A213022 Expansion of phi(x)^2 * psi(x) in powers of x where phi(), psi() are Ramanujan theta functions.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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