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.

Previous Showing 11-18 of 18 results.

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

Views

Author

Michael Somos, May 29 2012

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Mathematica
    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 *)
  • PARI
    {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))};

Formula

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.

A008444 Theta series of A_4 lattice.

Original entry on oeis.org

1, 20, 30, 60, 60, 120, 40, 180, 150, 140, 130, 240, 180, 360, 120, 260, 220, 480, 210, 400, 360, 240, 360, 660, 200, 620, 240, 600, 540, 600, 240, 640, 630, 720, 320, 780, 420, 1080, 600, 480, 650, 840, 360, 1260, 720, 840, 440, 1380, 660, 860, 630, 640, 1080, 1560, 400
Offset: 0

Views

Author

Keywords

Comments

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

Examples

			G.f. = 1 + 20*x + 30*x^2 + 60*x^3 + 60*x^4 + 120*x^5 + 40*x^6 + 180*x^7 + ...
G.f. = 1 + 20*q^2 + 30*q^4 + 60*q^6 + 60*q^8 + 120*q^10 + 40*q^12 + 180*q^14 + 150*q^16 + 140*q^18 + 130*q^20 + 240*q^22 + 180*q^24 + 360*q^26 + 120*q^28 + 260*q^30 + 220*q^32 + 480*q^34 + 210*q^36 + 400*q^38 + 360*q^40 + 240*q^42 + 360*q^44 + 660*q^46 + 200*q^48 + 620*q^50 + ...
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 110.

Crossrefs

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

Programs

  • Magma
    L := Lattice("A", 4); A := ThetaSeries(L, 120); A;
    
  • Magma
    A := Basis( ModularForms( Gamma1(5), 2), 55) ; A[1] + 20*A[2] + 30*A[3]; /* Michael Somos, Nov 13 2014 */
  • Mathematica
    a[ n_] := With[ {u1 = QPochhammer[ x], u5 = QPochhammer[ x^5]}, SeriesCoefficient[ u1^5/u5 + 25 x u5^5/u1, {x, 0, n}]]; (* Michael Somos, Nov 13 2014 *)
    terms = 55; f[q_] = LatticeData["A4", "ThetaSeriesFunction"][-I Log[q]/Pi]; s = Series[f[q], {q, 0, 2 terms}]; CoefficientList[s, q^2][[1 ;; terms]] (* Jean-François Alcover, Jul 04 2017 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^5 / eta(x^5 + A) + 25 * x * eta(x^5 + A)^5 / eta(x + A), n))}; /* Michael Somos, Feb 06 2011 */
    

Formula

Expansion of f(-x)^5 / f(-x^5) + 25 * x * f(-x^5)^5 / f(-x) in powers of x where f() is a Ramanujan theta function. - Michael Somos, Feb 06 2011
Expansion of (1 / Pi) integral_{0 .. Pi/2} theta_3(z, q)^5 + theta_4(z, q)^5 dz in powers of q^2. - Michael Somos, Jan 01 2012
Coefficient of x^0 in the expansion f(x * q, q / x)^5 in powers of q^2 where f() is a Ramanujan theta function. - Michael Somos, Jan 01 2012
G.f. is a period 1 Fourier series which satisfies f(-1 / (5 t)) = 5^(3/2) (t/i)^2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A023916. - Michael Somos, Feb 06 2011
A023916(5*n) = a(n) for all n in Z.

A033717 Number of integer solutions to the equation x^2 + 2*y^2 + 4*z^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

Views

Author

Keywords

Comments

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

Examples

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

References

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

Crossrefs

Programs

  • Magma
    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 */
  • Mathematica
    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 *)
  • PARI
    {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 */
    
  • PARI
    {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 */
    

Formula

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

A246631 Number of integer solutions to x^2 + 2*y^2 + 2*z^2 = n.

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, 0, 6, 32, 32, 24, 48, 32, 16, 0
Offset: 0

Views

Author

Michael Somos, Aug 31 2014

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma0(8), 3/2), 80); A[1] + 2*A[2];
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^2]^2, {q, 0, n}];
  • PARI
    {a(n) = if( n<1, n==0, 2 * qfrep([ 1, 0, 0; 0, 2, 0; 0, 0, 2], n)[n])};
    
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^4 + A)^8 / (eta(x + A)^2 * eta(x^8 + A)^4), n))};
    

Formula

Theta series of quadratic form with Gram matrix [ 1, 0, 0; 0, 2, 0; 0, 0, 2 ].
Expansion of phi(q) * phi(q^2)^2 = phi(-q^4)^4 / phi(-q) in powers of q where phi() is a Ramanujan theta function.
Expansion of eta(q^2) * eta(q^4)^8 / (eta(q)^2 * eta(q^8)^4) in powers of q.
Euler transform of period 8 sequence [ 2, 1, 2, -7, 2, 1, 2, -3, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 4 (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A014455.
G.f.: theta_3(q) * theta_3(q^2)^2.
G.f.: Product{k>0} (1 - x^(2*k)) * (1 - x^(4*k))^8 / ((1 - x^k)^2 * (1 - x^(8*k))^4).
G.f.: Product{k>0} (1 + x^(2*k)) * (1 + x^k)^2 * (1 - x^(4*k))^3 / (1 + x^(4*k))^4.
a(n) = (-1)^floor((n+1) / 2) * A212885(n) = abs(A212885(n)).
a(n) = A033717(2*n). a(2*n) = A014455(n). a(2*n + 1) = 2 * A246811(n).
a(4*n) = A005875(n). a(4*n + 1) = 2 * A045834(n). a(4*n + 2) = 4 * A045828(n).
a(8*n) = A004015(n). a(8*n + 1) = 2 * A213022(n). a(8*n + 2) = 4 * A213625(n). a(8*n + 3) = 8 * A008443(n). a(8*n + 4) = 2 * A045826(n). a(8*n + 5) = 8 * A045831(n). a(8*n + 6) = 8 * A213624(n). a(8*n + 7) = 0.

A386315 Number of points in a face-centered cubic lattice that intersect a sphere of radius n centered on a point in the lattice.

Original entry on oeis.org

1, 12, 12, 36, 12, 84, 36, 108, 12, 108, 84, 132, 36, 180, 108, 252, 12, 204, 108, 228, 84, 324, 132, 300, 36, 444, 180, 324, 108, 372, 252, 396, 12, 396, 204, 756, 108, 468, 228, 540, 84, 492, 324, 516, 132, 756, 300, 588, 36, 780, 444, 612, 180, 660, 324
Offset: 0

Views

Author

Charles L. Hohn, Aug 15 2025

Keywords

Comments

For all n > 0, the points at the 4 90-degree rotations of [n, 0, 0] and the eight 90-degree rotations and vertical reflections of [n/2, n/2, n*sqrt(1/2)] form the 12 vertices of a cuboctahedron (shown in red in the Links example). Points that otherwise lie on an axis plane have 24-point symmetry (green), and all other points have 48-point symmetry (blue). Thus, a(n) for all n > 0 are odd multiples of 12.
The number of noncongruent points for each sphere radius n (points within or on vertices or edges of each symmetry region in the Links example) gives A387222, the number of those points that are primitive for radius n (darker colors) gives A387223 for odd sphere radii, and the total primitive count divided by 12 gives A278081 for odd sphere radii. Examples of nonprimitive points include [3, 0, 0] and [3/2, 3/2, 3*sqrt(1/2)] for a(3), which reduce to a(1) primitive points [1, 0, 0] and [1/2, 1/2, sqrt(1/2)] respectively.
Analog for the simple cubic lattice is A016725.

Examples

			a(3) = 36, which is the sum of 4 90-degree rotations of [3, 0, 0], 8 90-degree rotations and vertical reflections of [3/2, 3/2, 3*sqrt(1/2)] and [1, 0, 4*sqrt(1/2)], and 16 90-degree rotations and vertical and horizontal reflections of [5/2, 3/2, sqrt(1/2)].
		

Crossrefs

Programs

  • PARI
    a(n)={my(c=0); for(x=0, n, for(y=0, min(x, sqrtint(n^2-x^2)), for(o=0, 1, my(m=2*(n^2-(x+o/2)^2-(y+o/2)^2)); if(!issquare(m), next); my(z=sqrtint(m)); if(z>=0 && z%2==o, c+=2^(4-if(!z, 1)-if(x==y, 1)-if(!min(x, y) && !o, 1)-if(!vecmax([x, y, z, o]), 1)))))); c}
    
  • PARI
    a(n)={if(!n, return(1)); my(f=Vec(factor(n)), o=12, r=o); for(i=if(#f[1] && f[1][1]==2, 2, 1), #f[1], my(m=if(f[1][i]%8>=4, 2)); f[2][i]++; while(f[2][i]--, o=o*f[1][i]+r*m); r=o); o}

Formula

a(n) = A004015(n^2).
a(2*n) = a(n).
a(p*n) = p*a(n) where p is a prime and p mod 8 is in {1, 3}.
a(p*n) = p*a(n) + 2*a(n/p^c) where p is a prime, p mod 8 is in {5, 7}, and c is the count of prime factors p in n.

A214813 Maximal contact number of a subset of n balls from the face-centered cubic lattice.

Original entry on oeis.org

0, 1, 3, 6, 9, 12, 15, 18, 21
Offset: 1

Views

Author

N. J. A. Sloane, Jul 31 2012

Keywords

Comments

If S is an arrangement of non-overlapping balls of radius 1, the contact number of S is the number of pairs of balls that just touch each other.
a(13) >= 36 (take one ball and its 12 neighbors), so this is different from A008486.
If b(n) denotes the maximal contact number of any arrangement of n balls then it is conjectured that a(n) = b(n) for n <= 9. It is also known that b(10)>=25, b(11)>=29, b(12)>=33 and of course b(13) >= a(13) >= 36. [Bezdek 2012]
Note that Figure 1e of Bezdek's arxiv:1601.00145 shows at n=5 a sphere packing with 9 contacts on the hexagonal close package (!), not on the cubic close package (which equals the f.c.c.). [In Figure 1e there is one sphere that touches from above a set of 3 spheres in a middle layer right above the bottom sphere; so this needs the ABABA... layer structures of the h.c.p, and cannot be done with the ABCABC... layer structure of the f.c.c.] So Figure 1e is not demonstrating a(5)=9. The correct value for the f.c.c is apparently a(5)=8 (where two structures with 8 contacts exist.) - R. J. Mathar, Mar 13 2018

Crossrefs

A297331 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals, where column k is the expansion of (theta_3(q^(1/2))^k + theta_4(q^(1/2))^k)/2.

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 1, 4, 2, 0, 1, 12, 4, 0, 0, 1, 24, 6, 0, 0, 0, 1, 40, 24, 24, 4, 0, 0, 1, 60, 90, 96, 12, 8, 0, 0, 1, 84, 252, 240, 24, 24, 0, 0, 0, 1, 112, 574, 544, 200, 144, 8, 0, 2, 0, 1, 144, 1136, 1288, 1020, 560, 96, 48, 4, 0, 0, 1, 180, 2034, 3136, 3444, 1560, 400, 192, 6, 4, 0, 0
Offset: 0

Views

Author

Ilya Gutkovskiy, Dec 28 2017

Keywords

Examples

			Square array begins:
1,  1,  1,   1,    1,    1,  ...
0,  0,  4,  12,   24,   40,  ...
0,  2,  4,   6,   24,   90,  ...
0,  0,  0,  24,   96,  240,  ...
0,  0,  4,  12,   24,  200,  ...
0,  0,  8,  24,  144,  560,  ...
		

Crossrefs

Programs

  • Mathematica
    Table[Function[k, SeriesCoefficient[(EllipticTheta[3, 0, q^(1/2)]^k + EllipticTheta[4, 0, q^(1/2)]^k)/2, {q, 0, n}]][j - n], {j, 0, 11}, {n, 0, j}] // Flatten

Formula

G.f. of column k: (theta_3(q^(1/2))^k + theta_4(q^(1/2))^k)/2, where theta_() is the Jacobi theta function.

A319078 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, 0, 6, 32
Offset: 0

Views

Author

Michael Somos, Sep 09 2018

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma0(16), 3/2), 66); A[1] + 2*A[2] - 4*A[3] - 8*A[4];
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 4, 0, q] EllipticTheta[ 3, 0, q]^2, {q, 0, n}];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 4, 0, q^2]^2, {q, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^9 / (eta(x + A)^2 * eta(x^4 + A)^4), n))};
    

Formula

Expansion of eta(q^2)^9 / (eta(q)^2 * eta(q^4)^4) in powers of q.
Expansion of phi(q) * phi(-q^2)^2 = phi(-q^2)^4 / phi(-q) in powers of q.
Euler transform of period 4 sequence [2, -7, 2, -3, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 2^(11/2) (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A045834.
G.f. Product_{k>0} (1 - x^k)^3 * (1 + x^k)^5 / (1 + x^(2*k))^4.
a(n) = (-1)^n * A212885(n) = A083703(2*n) = A080965(2*n).
a(4*n) = a(n) * -A132429(n + 2) where A132429 is a period 4 sequence.
a(4*n) = A005875(n). a(4*n + 1) = 2 * A045834(n). a(4*n + 2) = -4 * A045828(n).
a(8*n) = A004015(n). a(8*n + 1) = 2 * A213022(n). a(8*n + 2) = -4 * A213625(n). a(8*n + 3) = -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.
Previous Showing 11-18 of 18 results.