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-10 of 17 results. Next

A033762 Product t2(q^d); d | 3, where t2 = theta2(q) / (2 * q^(1/4)).

Original entry on oeis.org

1, 1, 0, 2, 1, 0, 2, 0, 0, 2, 2, 0, 1, 1, 0, 2, 0, 0, 2, 2, 0, 2, 0, 0, 3, 0, 0, 0, 2, 0, 2, 2, 0, 2, 0, 0, 2, 1, 0, 2, 1, 0, 0, 0, 0, 4, 2, 0, 2, 0, 0, 2, 0, 0, 2, 2, 0, 0, 2, 0, 1, 0, 0, 2, 2, 0, 4, 0, 0, 2, 0, 0, 0, 3, 0, 2, 0, 0, 2, 0, 0, 2, 0, 0, 3, 2, 0
Offset: 0

Views

Author

Keywords

Comments

Number of solutions of 8*n + 4 = x^2 + 3*y^2 in positive odd integers. - Michael Somos, Sep 18 2004
Half the number of integer solutions of 4*n + 2 = x^2 + y^2 + z^2 where 0 = x + y + z and x and y are odd. - Michael Somos, Jul 03 2011
Given g.f. A(x), then q^(1/2) * 2 * A(q) is denoted phi_1(z) where q = exp(Pi i z) in Conway and Sloane.
Half of theta series of planar hexagonal lattice (A2) with respect to an edge.
Bisection of A002324. Number of ways of writing n as a sum of a triangular plus three times a triangular number [Hirschhorn]. - R. J. Mathar, Mar 23 2011
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).

Examples

			G.f. = 1 + x + 2*x^3 + x^4 + 2*x^6 + 2*x^9 + 2*x^10 + x^12 + x^13 + 2*x^15 + ...
G.f. = q + q^3 + 2*q^7 + q^9 + 2*q^13 + 2*q^19 + 2*q^21 + q^25 + q^27 + 2*q^31 + ...
a(6) = 2 since 8*6 + 4 = 52 = 5^2 + 3*3^2 = 7^2 + 3*1^2.
		

References

  • Burce C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, 1991, see p. 223 Entry 3(i).
  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, 1999, p. 103. See Eq. (13).
  • Nathan J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 78, Eq. (32.27).

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(12), 1), 202); A[2] + A[4]; /* Michael Somos, Jul 25 2014 */
  • Mathematica
    a[ n_] := If[ n < 0, 0, DivisorSum[ 2 n + 1, Mod[(3 - #)/2, 3, -1] &]]; (* Michael Somos, Jul 03 2011 *)
    QP = QPochhammer; s = (QP[q^2]*QP[q^6])^2/(QP[q]*QP[q^3]) + O[q]^100; CoefficientList[s, q] (* Jean-François Alcover, Nov 27 2015, adapted from PARI *)
    a[ n_] := If[ n < 1, Boole[n == 0], Times @@ (Which[# < 2, 0^#2, Mod[#, 6] == 5, 1 - Mod[#2, 2], True, #2 + 1] & @@@ FactorInteger@(2 n + 1))]; (* Michael Somos, Mar 06 2016 *)
    %t A033762 a[ n_] := SeriesCoefficient[ (1/4) x^(-1/2) EllipticTheta[ 2, 0, x^(1/2)] EllipticTheta[ 2, 0, x^(3/2)], {x, 0, n}]; (* Michael Somos, Mar 06 2016 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A) * eta(x^6 + A))^2 / (eta(x + A) * eta(x^3 + A)), n))}; /* Michael Somos, Sep 18 2004 */
    
  • PARI
    {a(n) = if( n<0, 0, n = 2*n + 1; sumdiv( n, d, kronecker( -12, d) * (n / d % 2)))}; /* Michael Somos, Nov 04 2005 */
    
  • PARI
    {a(n) = if( n<0, 0, n = 8*n + 4; sum( j=1, sqrtint( n\3), (j%2) * issquare(n - 3*j^2)))} /* Michael Somos, Nov 04 2005 */
    
  • PARI
    {a(n) = if( n<0, 0, sumdiv(2*n + 1, d, kronecker(-3, d)))}; /* Michael Somos, Mar 06 2016 */
    

Formula

Expansion of q^(-1/2) * (eta(q^2) * eta(q^6))^2 / (eta(q) * eta(q^3)) in powers of q. - Michael Somos, Apr 18 2004
Expansion of q^(-1) * (a(q) - a(q^4)) / 6 in powers of q^2 where a() is a cubic AGM theta function. - Michael Somos, Oct 24 2006
Expansion of psi(x) * psi(x^3) in powers of x where psi() is a Ramanujan theta function. - Michael Somos, Jul 03 2011
Euler transform of period 6 sequence [ 1, -1, 2, -1, 1, -2, ...]. - Michael Somos, Apr 18 2004
From Michael Somos, Sep 18 2004: (Start)
Given g.f. A(x), then B(x) = (x * A(x^2))^2 satisfies 0 = f(B(x), B(x^2), B(x^4)) where f(u, v, w) = v^3 + 4*u*v*w + 16*v*w^2 - 8*w*v^2 - w*u^2.
a(n) = b(2*n + 1) where b() is multiplicative with b(2^e) = 0^e, b(3^e) = 1, b(p^e) = (1 + (-1)^e) / 2 if p==5 (mod 6) otherwise b(p^e) = e+1. (Clarification: the g.f. A(x) is not the primary function of interest, but rather B(x) = x * A(x^2), which is an eta-quotient and is the generating function of a multiplicative sequence.)
G.f.: (Sum_{j>0} x^((j^2 - j) / 2)) * (Sum_{k>0} x^(3(k^2 - k) / 2)) = Product_{k>0} (1 + x^k) * (1 - x^(2*k)) * (1 + x^(3*k)) * (1 - x^(6*k)).
G.f.: Sum_{k>=0} a(k) * x^(2*k + 1) = Sum_{k>0} x^k * (1 - x^k) * (1 - x^(4*k)) * (1 - x^(5*k)) / (1 - x^(12*k)). (End)
G.f.: s(4)^2*s(12)^2/(s(2)*s(6)), where s(k) := subs(q=q^k, eta(q)), where eta(q) is Dedekind's function, cf. A010815. [Fine]
G.f.: Sum_{k>=0} a(k) * x^(2*k + 1) = Sum_{k>0} x^k / (1 + x^k + x^(2*k)) - x^(4*k) / (1 + x^(4*k) + x^(8*k)). - Michael Somos, Nov 04 2005
a(n) = A002324(2*n + 1) = A035178(2*n + 1) = A091393(2*n + 1) = A093829(2*n + 1) = A096936(2*n + 1) = A112298(2*n + 1) = A113447(2*n + 1) = A113661(2*n + 1) = A113974(2*n + 1) = A115979(2*n + 1) = A122860(2*n + 1) = A123331(2*n + 1) = A123484(2*n + 1) = A136748(2*n + 1) = A137608(2*n + 1). A005881(n) = 2*a(n).
6 * a(n) = A004016(6*n + 3). - Michael Somos, Mar 06 2016
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(2*sqrt(3)) = 0.906899... (A093766). - Amiram Eldar, Nov 23 2023

Extensions

Corrected by Charles R Greathouse IV, Sep 02 2009

A123884 Expansion of phi(x) * phi(-x^6) / chi(-x^2) in powers of x where phi(), chi() are Ramanujan theta functions.

Original entry on oeis.org

1, 2, 1, 2, 3, 2, 2, 0, 2, 2, 1, 4, 0, 2, 3, 2, 2, 0, 4, 2, 2, 0, 0, 2, 1, 4, 2, 2, 2, 2, 3, 2, 0, 2, 2, 2, 2, 0, 2, 4, 4, 0, 0, 0, 1, 2, 4, 0, 2, 4, 2, 2, 1, 6, 0, 2, 2, 0, 0, 2, 4, 2, 0, 2, 2, 0, 4, 0, 4, 2, 1, 2, 0, 2, 4, 0, 0, 2, 2, 4, 3, 4, 0, 2, 2, 2, 2
Offset: 0

Views

Author

Michael Somos, Oct 17 2006

Keywords

Comments

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

Examples

			G.f. = 1 + 2*x + x^2 + 2*x^3 + 3*x^4 + 2*x^5 + 2*x^6 + 2*x^8 + 2*x^9 + x^10 + ...
G.f. = q + 2*q^13 + q^25 + 2*q^37 + 3*q^49 + 2*q^61 + 2*q^73 + 2*q^97 + 2*q^109 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x]^2 EllipticTheta[ 4, 0, x^6] / EllipticTheta[ 4, 0, x^2], {x, 0, n}]; (* Michael Somos, Oct 01 2015 *)
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, x] EllipticTheta[ 4, 0, x^6] QPochhammer[ -x^2, x^2], {x, 0, n}]; (* Michael Somos, Oct 01 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^4 * eta(x^6 + A)^2 / (eta(x + A)^2 * eta(x^4 + A) * eta(x^12 + A)), n))};

Formula

Expansion of q^(-1/12) * eta(q^2)^4 * eta(q^6)^2 / (eta(q)^2 * eta(q^4) * eta(q^12)) in powers of q.
Euler transform of period 12 sequence [ 2, -2, 2, -1, 2, -4, 2, -1, 2, -2, 2, -2, ...].
a(n) = A093829(12*n + 1).
a(n) = (-1)^n * A248886(n). a(2*n) = A131961(n). a(2*n + 1) = 2 * A131963(n). - Michael Somos, Oct 01 2015

A121361 Expansion of f(x^1, x^5) * psi(x^2) in powers of x where psi(), f() are Ramanujan theta functions.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 2, 2, 1, 1, 0, 1, 0, 1, 2, 0, 1, 1, 0, 2, 0, 2, 1, 0, 1, 1, 1, 1, 2, 1, 0, 1, 2, 1, 0, 0, 1, 1, 1, 1, 0, 0, 2, 1, 2, 0, 1, 1, 1, 2, 1, 1, 0, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 3, 0, 0, 1, 0, 1, 0, 0, 2, 1, 1, 1, 1, 1, 2, 0, 1, 0, 2, 2, 1, 3, 0, 0, 0, 1, 0, 0
Offset: 0

Views

Author

Michael Somos, Jul 16 2006

Keywords

Comments

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

Examples

			G.f. = 1 + x + x^2 + x^3 + x^5 + x^6 + 2*x^7 + x^8 + x^10 + x^11 + ...
G.f. = q^7 + q^19 + q^31 + q^43 + q^67 + q^79 + 2*q^91 + q^103 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x^1, x^6] QPochhammer[ -x^5, x^6] QPochhammer[ x^6] EllipticTheta[ 2, 0, x] / (2 x^(1/4)), {x, 0, n}]; (* Michael Somos, Sep 02 2014 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^3 + A) * eta(x^4 + A) * eta(x^12 + A) / (eta(x + A) * eta(x^6 + A)), n))};

Formula

Expansion of q^(-7/12) * eta(q^2) * eta(q^3) * eta(q^4) * eta(q^12) /
(eta(q) * eta(q^6)) in powers of q.
Euler transform of period 12 sequence [ 1, 0, 0, -1, 1, 0, 1, -1, 0, 0, 1, -2, ...].
2*a(n) = A093829(12*n + 7).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(2*sqrt(3)) = 0.906899... (A093766). - Amiram Eldar, Jan 20 2025

A035178 a(n) = Sum_{d|n} Kronecker(-12, d) (= A134667(d)).

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 2, 1, 1, 0, 0, 1, 2, 2, 0, 1, 0, 1, 2, 0, 2, 0, 0, 1, 1, 2, 1, 2, 0, 0, 2, 1, 0, 0, 0, 1, 2, 2, 2, 0, 0, 2, 2, 0, 0, 0, 0, 1, 3, 1, 0, 2, 0, 1, 0, 2, 2, 0, 0, 0, 2, 2, 2, 1, 0, 0, 2, 0, 0, 0, 0, 1, 2, 2, 1, 2, 0, 2, 2, 0, 1, 0, 0, 2, 0, 2, 0, 0, 0, 0, 4, 0, 2, 0, 0, 1, 2, 3, 0, 1, 0, 0, 2, 2, 0
Offset: 1

Views

Author

Keywords

Comments

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

Examples

			G.f. = q + q^2 + q^3 + q^4 + q^6 + 2*q^7 + q^8 + q^9 + q^12 + 2*q^13 + 2*q^14 + ...
		

References

  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, p. 346.

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(6), 1), 88); B := (A[1] - 1) / 3 + A[2]; B; /* Michael Somos, Aug 04 2015 */
  • Mathematica
    a[ n_] := If[ n < 1, 0, Sum[ KroneckerSymbol[ -12, d], { d, Divisors[ n]}]]; (* Michael Somos, Jun 24 2011 *)
    a[ n_] := If[ n < 1, 0, Times @@ (Which[ # < 5, 1, Mod[#, 6] == 5, 1 - Mod[#2, 2], True, #2 + 1 ] & @@@ FactorInteger@n)]; (* Michael Somos, Aug 04 2015 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, 0, q^(1/2)]^3 / EllipticTheta[ 2, 0, q^(3/2)] - 4) / 12, {q, 0, n}]; (* Michael Somos, Aug 04 2015 *)
    a[n_] := DivisorSum[n, KroneckerSymbol[-12, #]&]; Array[a, 105] (* Jean-François Alcover, Dec 01 2015 *)
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, kronecker( -12, d)))}; /* Michael Somos, Apr 18 2004 */
    
  • PARI
    {a(n) = if( n<1, 0, direuler( p=2, n, 1 / ((1 - X) * (1 - kronecker( -12, p) * X))) [n])}; /* Michael Somos, Jun 24 2011 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^3 + A) * eta(x^2 + A)^6 / (eta(x^6 + A)^2 * eta(x + A)^3) - 1) / 3, n))}; /* Michael Somos, Aug 11 2009 */
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k, ]; if( p<5, 1, p%6==5, 1-e%2, 1+e)))}; /* Michael Somos, Aug 04 2015 */
    

Formula

Moebius transform is period 6 sequence [ 1, 0, 0, 0, -1, 0, ...]. - Michael Somos, Feb 14 2006
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = (u1 - u2) * (u1 - u2 - u3 + u6) - (u2 -u6) * (1 + 3*u6). - Michael Somos, May 29 2005
Dirichlet g.f.: zeta(s) * L(chi,s) where chi(n) = Kronecker( -12, n). Sum_{n>0} a(n) / n^s = Product_{p prime} 1 / ((1 - p^-s) * (1 - Kronecker( -12, p) * p^-s)). - Michael Somos, Jun 24 2011
a(n) is multiplicative with a(p^e) = 1 if p=2 or p=3, a(p^e) = 1+e if p == 1 (mod 6), a(p^e) = (1 + (-1)^e)/2 if p == 5 (mod 6).
G.f.: Sum_{k>0} (x^k + x^(3*k)) / (1 + x^(2*k) + x^(4*k)) = Sum_{k>=0} x^(6*k + 1) / (1 - x^(6*k + 1)) - x^(6*k + 5) / (1 - x^(6*k + 5)). - Michael Somos, Feb 14 2006
a(n) = |A093829(n)| = -(-1)^n * A137608(n) = a(2*n) = a(3*n). a(6*n + 1) = A097195(n). a(6*n + 5) = 0.
From Michael Somos, Aug 11 2009: (Start)
3 * a(n) = A107760(n) unless n=0. a(2*n + 1) = A033762(n). a(3*n + 1) = A033687(n). a(4*n + 1) = A112604(n). a(4*n + 3) = A112605(n).
a(8*n + 1) = A112606(n). a(8*n + 3) = A112608(n). a(8*n + 5) = 2 * A112607(n). a(8*n + 7) = 2 * A112608(n). a(12*n + 1) A123884(n). a(12*n + 7) = 2 * A121361(n).
a(24*n + 1) = A131961(n). a(24*n + 7) = 2 * A131962(n). a(24*n + 13) = 2 * A131963(n). a(24*n + 19) = 2 * A131964(n). (End)
Expansion of (psi(q)^3 / psi(q^3) - 1) / 3 in powers of q where psi() is a Ramanujan theta function. - Michael Somos, Aug 04 2015
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(2*sqrt(3)) = 0.906899... (A093766). - Amiram Eldar, Nov 16 2023

Extensions

Definition edited by Michael Somos, Aug 11 2009

A113447 Expansion of i * theta_2(i * q^3)^3 / (4 * theta_2(i * q)) in powers of q^2.

Original entry on oeis.org

1, 1, 1, -1, 0, 1, 2, 1, 1, 0, 0, -1, 2, 2, 0, -1, 0, 1, 2, 0, 2, 0, 0, 1, 1, 2, 1, -2, 0, 0, 2, 1, 0, 0, 0, -1, 2, 2, 2, 0, 0, 2, 2, 0, 0, 0, 0, -1, 3, 1, 0, -2, 0, 1, 0, 2, 2, 0, 0, 0, 2, 2, 2, -1, 0, 0, 2, 0, 0, 0, 0, 1, 2, 2, 1, -2, 0, 2, 2, 0, 1, 0, 0, -2, 0, 2, 0, 0, 0, 0, 4, 0, 2, 0, 0, 1, 2, 3, 0, -1, 0, 0, 2, 2, 0
Offset: 1

Views

Author

Michael Somos, Nov 02 2005

Keywords

Comments

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

Examples

			G.f. = q + q^2 + q^3 - q^4 + q^6 + 2*q^7 + q^8 + q^9 - q^12 + 2*q^13 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(24), 1), 106); A[2] + A[3] + A[4] - A[5] + A[7] + 2*A[8] + A[9] + A[10]; /* Michael Somos, May 07 2015 */
  • Mathematica
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, {1, 0, 0, -2, -1, 0, 1, 2, 0, 0, -1, 0}[[Mod[#, 12, 1]]] &]]; (* Michael Somos, Jan 31 2015 *)
  • PARI
    {a(n) = if( n<1, 0, -(-1)^max( 1, valuation( n, 2)) * sumdiv(n, d, kronecker( -12, d)))};
    
  • PARI
    {a(n) = if( n<1, 0, direuler( p=2, n, if( p==2, 1 + X / (1 + X), 1 / ((1 - X) * (1 - kronecker( -12, p) * X))))[n])};
    
  • PARI
    {a(n) = local(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^3 + A)^3 * eta(x^12 + A)^3 / (eta(x + A) * eta(x^4 + A) * eta(x^6 + A)^3), n))};
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv(n, d, [ 0, 1, 0, 0, -2, -1, 0, 1, 2, 0, 0,-1][d%12 + 1]))}; /* Michael Somos, May 07 2015 */
    

Formula

Expansion of (eta(q^2) * eta(q^3)^3 * eta(q^12)^3) / (eta(q) * eta(q^4) * eta(q^6)^3) in powers of q.
Euler transform of period 12 sequence [1, 0, -2, 1, 1, 0, 1, 1, -2, 0, 1, -2, ...].
Moebius transform is period 12 sequence [1, 0, 0, -2, -1, 0, 1, 2, 0, 0, -1, 0, ...].
a(n) is multiplicative and a(2^e) = -(-1)^e if e>0, a(3^e) = 1, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1+(-1)^e)/2 if p == 5 (mod 6).
G.f.: Sum_{k>0} x^(6*k - 5) / (1 - x^(6*k - 5)) - x^(6*k - 1) / (1 - x^(6*k - 1)) - 2 * x^(12*k - 8) / (1 - x^(12*k - 8)) + 2 * x^(12*k - 4) / (1 - x^(12*k-4)).
G.f.: Sum_{k>0} x^k * (1 - x^(3*k))^2 / (1 + x^(4*k) + x^(8*k)).
G.f.: x * Product_{k>0} (1 - x^k) / (1 - x^(4*k - 2)) * ((1 - x^(12*k - 6)) / (1 - x^(3*k)))^3.
Expansion of theta_2(i * q^3)^3 / (4 * theta_2(i * q)) in powers of q^2.
Expansion of q * psi(-q^3)^3 / psi(-q) in powers of q where psi() is a Ramanujan theta function.
Expansion of (c(q) * c(q^4)) / (3 * c(q^2)) in powers of q where c() is a cubic AGM theta function.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = (4/3)^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g(t) is the g.f. for A132973.
a(n) = -(-1)^n * A093829(n). - Michael Somos, Jan 31 2015
Convolution inverse of A133637.
a(3*n) = a(n). a(6*n + 5) = a(12*n + 10) = 0. |a(n)| = A035178(n).
a(2*n) = A093829(n). a(2*n + 1) = A033762(n).
a(4*n + 1) = A112604(n). a(4*n + 3) = A112605(n).
a(6*n + 1) = A097195(n). a(6*n + 2) = A033687(n).
a(8*n + 1) = A112606(n). a(8*n + 3) = A112608(n). a(8*n + 5) = 2 * A112607(n). a(8*n + 6) = A112605(n). a(8*n + 7) = 2 * A112609(n).
a(12*n + 1) = A123884(n). a(12*n + 7) = 2 * A121361(n).
a(24*n + 1) = A131961(n). a(24*n + 7) = 2 * A131962(n). a(24*n + 13) = 2 * A121963(n). a(24*n + 19) = 2 * A131964(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(6*sqrt(3)) = 0.604599... (A073010). - Amiram Eldar, Nov 23 2023

A123633 Expansion of (c(q^2)/c(q))^3 in powers of q where c() is a cubic AGM theta function.

Original entry on oeis.org

1, -3, 3, 5, -18, 15, 24, -75, 57, 86, -252, 183, 262, -744, 522, 725, -1998, 1365, 1852, -4986, 3336, 4436, -11736, 7719, 10103, -26322, 17067, 22040, -56682, 36306, 46336, -117867, 74700, 94378, -237744, 149277, 186926, -466836, 290706, 361126, -895014, 553224
Offset: 1

Views

Author

Michael Somos, Oct 03 2006, Jan 21 2009

Keywords

Comments

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).
In the arXiv:2305.13951 paper on page 21 is this: "The q-expansion of y coincides with the sequence A123633 in the OEIS". - Michael Somos, May 26 2023

Examples

			G.f. = q - 3*q^2 + 3*q^3 + 5*q^4 - 18*q^5 + 15*q^6 + 24*q^7 - 75*q^8 + 57*q^9 + ...
		

Crossrefs

Programs

  • Magma
    M := Basis(ModularForms(Gamma1(6), 1), 43); M1 := M[1]; M2 := M[2]; A := M2/(M1 + 2*M2); A; /* Michael Somos, May 26 2023 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ q / (QPochhammer[ q^3, q^6]^3 / QPochhammer[ q, q^2])^3, {q, 0, n}]; (* Michael Somos, Feb 19 2015 *)
    a[ n_] := SeriesCoefficient[ q (Product[ 1 - q^k, {k, 1, n, 2}] / Product[ 1 - q^k, {k, 3, n, 6}]^3)^3, {q, 0, n}]; (* Michael Somos, Feb 19 2015 *)
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x + A) / eta(x^2 + A))^3 * (eta(x^6 + A) / eta(x^3 + A))^9, n))};
    

Formula

Expansion of q / (chi(-q^3)^3 / chi(-q))^3 in powers of q where chi() is a Ramanujan theta function.
Euler transform of period 6 sequence [ -3, 0, 6, 0, -3, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v)= u^2 - v - u*v * (6 + 8*v).
G.f.: x * (Product_{k>0} (1 - x^(2*k - 1)) / (1 - x^(6*k - 3))^3 )^3.
G.f. is a period 1 Fourier series which satisfies f(-1 / (6 t)) = (1 / 8) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A128642.
A128636(n) = a(n) unless n = 0. Convolution inverse of A105559.
Convolution cube of A092848.
Convolution with A123330 is A093829. - Michael Somos, May 26 2023

A228447 Expansion of q * (psi(q^3) * psi(q^6)) / (psi(q) * phi(q)) in powers of q where phi(), psi() are Ramanujan theta functions.

Original entry on oeis.org

1, -3, 7, -15, 30, -57, 104, -183, 313, -522, 852, -1365, 2150, -3336, 5106, -7719, 11538, -17067, 25004, -36306, 52280, -74700, 105960, -149277, 208951, -290706, 402127, -553224, 757158, -1031166, 1397744, -1886151, 2534316, -3391254, 4520112, -6002007
Offset: 1

Views

Author

Michael Somos, Oct 26 2013

Keywords

Comments

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

Examples

			G.f. = q - 3*q^2 + 7*q^3 - 15*q^4 + 30*q^5 - 57*q^6 + 104*q^7 - 183*q^8 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ (1/4) EllipticTheta[ 2, 0, q^(3/2)]^3 / (EllipticTheta[ 2, 0, q^(1/2)] EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^3]), {q, 0, n}]
  • PARI
    {a(n) = local(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x + A)^3 * eta(x^4 + A)^2 * eta(x^6 + A) * eta(x^12 + A)^2 / (eta(x^2 + A)^7 * eta(x^3 + A)), n))}

Formula

Expansion of q * (psi(q^3)^3 / psi(q)) / (phi(q) * phi(q^3)) in powers of q where phi(), psi() are Ramanujan theta functions.
Expansion of eta(q)^3 * eta(x^4)^2 * eta(x^6) * eta(x^12)^2 / (eta(x^2)^7 * eta(x^3)) in powers of q.
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^2 * (1 - 3*v) - v * (1 - 4*v) * (1 - 3*u)^2.
a(n) = -(-1)^n * A187100(n). a(2*n) = -3 * A128638(n).
Convolution inverse is A187145. Convolution with A033716 is A093829.

A229615 Expansion of q^2 * psi(q^3)^6 / psi(q)^2 in powers of q where psi() is a Ramanujan theta function.

Original entry on oeis.org

1, -2, 3, 0, -1, 0, 7, -8, 6, 0, 1, 0, 8, -12, 15, 0, -7, 0, 18, -16, 12, 0, 5, 0, 14, -26, 24, 0, -6, 0, 31, -24, 18, 0, -5, 0, 20, -28, 42, 0, -8, 0, 36, -48, 24, 0, 13, 0, 31, -36, 42, 0, -25, 0, 56, -40, 30, 0, 6, 0, 32, -64, 63, 0, -12, 0, 54, -48, 48, 0
Offset: 2

Views

Author

Michael Somos, Sep 26 2013

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Magma
    Basis( ModularForms( Gamma0(6), 2), 70)[3] /* Michael Somos, Mar 05 2023 */
  • Mathematica
    a[ n_] := If[n < 1, 0, Sum[ {0, 1, -2, 1, 0, 0}[[ Mod[d, 6, 1]]] n/d, {d, Divisors[n]}]];
    a[ n_] := If[n < 1, 0, Sum[ {0, 1/2, -2/3, 1/2, 0, 0}[[ Mod[d, 6, 1]]] d, {d, Divisors[n]}]];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, 0, q^(3/2)]^6 / EllipticTheta[ 2, 0, q^(1/2)]^2 / 16, {q, 0, n}];
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, n/d * [0, 0, 1, -2, 1, 0][d%6 + 1]))};
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, d * [0, 0, 1/2, -2/3, 1/2, 0][d%6 + 1]))};
    
  • PARI
    {a(n) = my(A); if( n<2, 0, n -= 2; A = x * O(x^n); polcoeff( (eta(x + A) * eta(x^6 + A)^6 / (eta(x^2 + A)^2 * eta(x^3 + A)^3))^2, n))};
    
  • Sage
    ModularForms( Gamma0(6), 2, prec=70).2;
    

Formula

Expansion of (a(q) - a(q^2))^2 / 36 = c(q^2)^4 / (9 * c(q)^2) in powers of q where a(), c() are cubic AGM theta functions.
Expansion of ((eta(q) * eta(q^6)^6) / (eta(q^2)^2 * eta(q^3)^3))^2 in powers of q.
Euler transform of period 6 sequence [ -2, 2, 4, 2, -2, -4, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (6 t)) = (1/12) (t / i)^2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A229616.
G.f.: sum_{k>0} x^(6*k-4) / (1 - x^(6*k-4))^2 - 2 * x^(6*k-3) / (1 - x^(6*k-3))^2 + x^(6*k-2) / (1 - x^(6*k-2))^2.
G.f.: sum_{k>0} (3*k-2) * x^(6*k-4) / (1 - x^(6*k-4)) - (4*k-2) * x^(6*k-3) / (1 - x^(6*k-3)) + (3*k-1) * x^(6*k-2) / (1 - x^(6*k-2)).
a(6*n + 1) = a(6*n + 5) = 0. a(6*n + 2) = A144614(n). a(6*n + 3) = -2 * A008438(n). a(6*n + 4) = 3 * A033686(n).
Convolution square of A093829.

A112298 Expansion of (a(q) - 3*a(q^2) + 2*a(q^4)) / 6 in powers of q where a() is a cubic AGM theta function.

Original entry on oeis.org

1, -3, 1, 3, 0, -3, 2, -3, 1, 0, 0, 3, 2, -6, 0, 3, 0, -3, 2, 0, 2, 0, 0, -3, 1, -6, 1, 6, 0, 0, 2, -3, 0, 0, 0, 3, 2, -6, 2, 0, 0, -6, 2, 0, 0, 0, 0, 3, 3, -3, 0, 6, 0, -3, 0, -6, 2, 0, 0, 0, 2, -6, 2, 3, 0, 0, 2, 0, 0, 0, 0, -3, 2, -6, 1, 6, 0, -6, 2, 0, 1, 0, 0, 6, 0, -6, 0, 0, 0, 0, 4, 0, 2, 0, 0, -3, 2, -9, 0, 3, 0, 0, 2, -6, 0
Offset: 1

Views

Author

Michael Somos, Sep 02 2005

Keywords

Comments

Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

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

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(12), 1), 106); A[2] - 3*A[3] + A[4] + 3*A[5]; /* Michael Somos, Jan 17 2015 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ q QPochhammer[ q, q^2]^3 QPochhammer[ -q^6, q^6]^3 EllipticTheta[ 4, 0, q^2] EllipticTheta[ 2, 0, q^(3/2)] / (2 q^(3/8)), {q, 0, n}]; (* Michael Somos, Jan 17 2015 *)
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, JacobiSymbol[ -3, n/#] {1, -2, 1, 0}[[Mod[#, 4, 1]]] &]]; (* Michael Somos, Jan 17 2015 *)
  • PARI
    {a(n) = if( n<1, 0, sumdiv(n, d, kronecker(-3, n/d)*[0, 1, -2, 1][d%4 + 1]))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x + A) * eta(x^12 + A))^3/ (eta(x^2 + A) * eta(x^3 + A) * eta(x^4 + A) * eta(x^6 + A)), n))};
    

Formula

From Michael Somos, Jan 17 2015: (Start)
Expansion of b(q) * (b(q^4) - b(q)) / (3*b(q^2)) in powers of q where b() is a cubic AGM theta function.
Expansion of q * chi(-q)^3 * phi(-q^2) * psi(q^3) / chi(-q^6)^3 in powers of q where phi(), psi(), chi() are Ramanujan theta functions.
Expansion of q * phi(-q)^2 * psi(q^6)^2 / (psi(-q) * psi(-q^3)) in powers of q where phi(), psi() are Ramanujan theta functions.
Expansion of q * f(q) * f(-q, -q^5)^4 / f(q^3)^3 in powers of q where f() is a Ramanujan theta function. (End)
Expansion of (eta(q) * eta(q^12))^3 / (eta(q^2) * eta(q^3) * eta(q^4) * eta(q^6)) in powers of q.
Euler transform of period 12 sequence [ -3, -2, -2, -1, -3, 0, -3, -1, -2, -2, -3, -2, ...].
Moebius transform is period 12 sequence [ 1, -4, 0, 6, -1, 0, 1, -6, 0, 4, -1, 0, ...].
Multiplicative with a(2^e) = 3(-1)^e if e>0, a(3^e)=1, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1 + (-1)^e)/2 if p == 2 (mod 6).
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 12^(1/2) (t/i) f(t) where q = exp(2 Pi i t).
G.f.: Sum_{k>0} Kronecker(-3, k) * x^k * (1 - x^k)^2 / (1 - x^(4*k)).
a(n) = -(-1)^n * A244375(n). a(6*n + 5) = 0, a(3*n) = a(n).
a(2*n) = -3 * A093829(n). a(2*n + 1) = A033762(n). a(3*n + 1) = A129576(n). a(4*n + 1) = A112604(n). a(4*n + 3) = A112605(n). a(6*n + 1) = A097195(n). a(6*n + 2) = -3 * A033687(n).
Sum_{k=1..n} abs(a(k)) ~ (Pi/sqrt(3)) * n. - Amiram Eldar, Jan 23 2024

A112848 Expansion of eta(q)*eta(q^2)*eta(q^18)^2/(eta(q^6)*eta(q^9)) in powers of q.

Original entry on oeis.org

1, -1, -2, 1, 0, 2, 2, -1, -2, 0, 0, -2, 2, -2, 0, 1, 0, 2, 2, 0, -4, 0, 0, 2, 1, -2, -2, 2, 0, 0, 2, -1, 0, 0, 0, -2, 2, -2, -4, 0, 0, 4, 2, 0, 0, 0, 0, -2, 3, -1, 0, 2, 0, 2, 0, -2, -4, 0, 0, 0, 2, -2, -4, 1, 0, 0, 2, 0, 0, 0, 0, 2, 2, -2, -2, 2, 0, 4, 2, 0, -2, 0, 0, -4, 0, -2, 0, 0, 0, 0, 4, 0, -4, 0, 0, 2, 2, -3, 0, 1, 0, 0, 2, -2, 0
Offset: 1

Views

Author

Michael Somos, Sep 22 2005

Keywords

Crossrefs

Cf. A033687, A033762, A092829, A093829, A097195, A248897, A255648 (absolute values).

Programs

  • Mathematica
    QP = QPochhammer; s = QP[q]*QP[q^2]*(QP[q^18]^2/(QP[q^6]*QP[q^9])) + O[q]^100; CoefficientList[s, q] (* Jean-François Alcover, Nov 25 2015 *)
    f[p_, e_] := If[Mod[p, 6] == 1, e+1, (1+(-1)^e)/2]; f[2, e_] := (-1)^e; f[3, e_]:= -2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jan 28 2024 *)
  • PARI
    {a(n)=if(n<1, 0, if(n%3==0, n/=3; -2,1)* sumdiv(n,d,kronecker(-12,d) -if(d%2==0, 2*kronecker(-3,d/2))))}
    
  • PARI
    {a(n)=local(A); if (n<1, 0, n--; A=x*O(x^n); polcoeff( eta(x+A)*eta(x^2+A)*eta(x^18+A)^2/ eta(x^6+A)/eta(x^9+A), n))}

Formula

Euler transform of period 18 sequence [ -1, -2, -1, -2, -1, -1, -1, -2, 0, -2, -1, -1, -1, -2, -1, -2, -1, -2, ...].
Moebius transform is period 18 sequence [1, -2, -3, 2, -1, 6, 1, -2, 0, 2, -1, -6, 1, -2, 3, 2, -1, 0, ...].
Multiplicative with a(2^e) = (-1)^e, a(3^e) = -2 if e>0, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1+(-1)^e)/2 if p == 5 (mod 6).
G.f.: Sum_{k>0} Kronecker(-3, k) x^k(1-x^(2k))^2/(1-x^(6k)) = x Product_{k>0} (1-x^k)(1-x^(2k))(1+x^(9k))(1+x^(6k)+x^(12k)).
a(3n) = -2*A092829(n). a(3n+1) = A093829(3n+1) = A033687(n). a(3n+2) = A093829(3n+2). a(6n)/2 = A093829(n). a(6n+1) = A097195(n). a(6n+3) = -2*A033762(n). a(6n+5) = 0.
Sum_{k=1..n} abs(a(k)) ~ c * n, where c = 2*Pi/(3*sqrt(3)) = 1.209199... (A248897). - Amiram Eldar, Jan 23 2024
Showing 1-10 of 17 results. Next