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-6 of 6 results.

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

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

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, Apr 17 2004

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
    Basis( ModularForms( Gamma1(6), 1), 90) [2]; /* Michael Somos, Jul 02 2014 */
  • Mathematica
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, {1, -2, 0, 2, -1, 0} [[ Mod[#, 6, 1]]] &]];
    QP = QPochhammer; s = (QP[q]*QP[q^6]^6)/(QP[q^2]^2*QP[q^3]^3) + O[q]^105; CoefficientList[s, q] (* Jean-François Alcover, Nov 30 2015, adapted from PARI *)
  • PARI
    {a(n) = if( n<1, 0, polcoeff( sum( k=0, n, x^k * (1 - x^k)^2 / (1 + x^(2*k) + x^(4*k)), x * O(x^n)), n))};
    
  • PARI
    {a(n) = local(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x + A) * eta(x^6 + A)^6 / (eta(x^2 + A)^2 * eta(x^3 + A)^3), n))};
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, kronecker( -12, d) - if( d%2==0, 2 * kronecker( -3, d/2) ) ))}; /* Michael Somos, May 29 2005 */
    
  • Sage
    ModularForms( Gamma1(6), 1, prec=90).1; # Michael Somos, Sep 27 2013
    

Formula

Expansion of (a(q) - a(q^2)) / 6 = c(q^2)^2 / (3 * c(q)) in powers of q where a(), c() are cubic AGM functions. - Michael Somos, Sep 06 2007
Expansion of (eta(q) * eta(q^6)^6) / (eta(q^2)^2 * eta(q^3)^3) in powers of q.
Euler transform of period 6 sequence [ -1, 1, 2, 1, -1, -2, ...].
Moebius transform is period 6 sequence [ 1, -2, 0, 2, -1, 0, ...] = A112300. - Michael Somos, Jul 16 2006
Multiplicative with a(p^e) = (-1)^e if p=2; a(p^e) = 1 if 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. is a period 1 Fourier series which satisfies f(-1 / (6 t)) = 12^(-1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is g.f. for A122859.
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = w * (u + v)^2 - v * (v + w) * (v + 4*w).
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u2 * (u2 - u3 - 4*u6) - (u3 + u6) * (u1 - 3*u3 - 3*u6).
G.f.: Sum_{k>0} (x^k - 2 * x^(2*k) + 2 * x^(4*k) - x^(5*k)) / (1 - x^(6*k)) = x * Product_{k>0} ((1 - x^k) * (1 - x^(6*k))^6) / ((1 - x^(2*k))^2 * (1 - x^(3*k))^3).
a(n) = -(-1)^n * A113447(n). - Michael Somos, Jan 31 2015
a(2*n) = -a(n). a(3*n) = a(n). a(6*n + 5) = 0.
A035178(n) = |a(n)|. A033762(n) = a(2*n + 1). A033687(n) = a(3*n + 1).
a(4*n + 1) = A112604(n). a(4*n + 3) = A112605(n). a(6*n + 1) = A097195(n). a(8*n + 1) = A112606(n). a(8*n + 3) = A112608(n). a(8*n + 5) = 2 * A112607(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(6*sqrt(3)) = 0.302299894039... . - Amiram Eldar, Nov 21 2023

A132973 Expansion of psi(-q)^3 / psi(-q^3) in powers of q where psi() is a Ramanujan theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Sep 07 2007

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. = 1 - 3*q + 3*q^2 - 3*q^3 + 3*q^4 + 3*q^6 - 6*q^7 + 3*q^8 - 3*q^9 + 3*q^12 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, Pi/4, q^(1/2)]^3 / EllipticTheta[ 2, Pi/4, q^(3/2)]/2, {q, 0, n}]; (* Michael Somos, May 26 2013 *)
  • PARI
    {a(n) = if( n<1, n==0, 3 * (-1)^n * sumdiv(n, d, kronecker(-12, d)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^3 * eta(x^4 + A)^3 * eta(x^6 + A) / (eta(x^2 + A)^3 * eta(x^3 + A) * eta(x^12 + A )), n))};

Formula

Expansion of b(q^2)^2 / b(-q) = b(q) * b(q^4) / b(q^2) in powers of q where b() is a cubic AGM theta function.
Expansion of (a(q^2) + 2 * a(q^4) - a(q)) / 2 = (c(q)^2 - 5 * c(q) * c(q^4) + 4 * c(q^4)^2) / (3 * c(q^2)) in powers of q where a(), c() are cubic AGM theta functions. - Michael Somos, May 26 2013
Expansion of eta(q)^3 * eta(q^4)^3 * eta(q^6) / (eta(q^2)^3 * eta(q^3) * eta(q^12)) in powers of q.
Euler transform of period 12 sequence [ -3, 0, -2, -3, -3, 0, -3, -3, -2, 0, -3, -2, ...].
Moebius transform is period 12 sequence [ -3, 6, 0, 0, 3, 0, -3, 0, 0, -6, 3, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 108^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g(t) is the g.f. for A113447.
G.f.: Product_{k>0} (1 - x^k)^3 * (1 + x^(2*k))^3 / ((1 - x^(3*k)) * (1 + x^(6*k))).
G.f.: 1 + 3 * Sum_{k>0} (-1)^k * (x^k + x^(3*k)) / (1 + x^k + x^(2*k)).
G.f.: 1 + 3 * ( Sum_{k>0} x^(6*k-5) / ( 1 + x^(6*k-5) ) - x^(6*k-1) / ( 1 + x^(6*k-1) )).
a(n) = (-1)^n * A107760(n). Convolution inverse of A132974.
a(2*n) = A107760(n). a(2*n + 1) = -3 * A033762(n). a(3*n) = A132973(n). a(3*n + 1) = -3 * A227696(n). - Michael Somos, Oct 31 2015
a(6*n + 1) = -3 * A097195(n). a(6*n + 2) = 3 * A033687(n). a(6*n + 5) = 0. - Michael Somos, Oct 31 2015

A137608 Expansion of (1 - psi(-q)^3 / psi(-q^3)) / 3 in powers of q where psi() is a Ramanujan theta function.

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, Jan 29 2008

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

  • Mathematica
    a[ n_] := If[ n < 1, 0, -(-1)^n DivisorSum[n, KroneckerSymbol[ -12, #] &]]; (* Michael Somos, May 06 2015 *)
    a[ n_] := SeriesCoefficient[ (4 + EllipticTheta[ 2, Pi/4, q^(1/2)]^3 / EllipticTheta[ 2, Pi/4, q^(3/2)]) / 6, {q, 0, n}]; (* Michael Somos, May 06 2015 *)
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, {1, -2, 0, 0, -1, 0, 1, 0, 0, 2, -1, 0}[[Mod[#, 12, 1]]] &]]; (* Michael Somos, May 07 2015 *)
  • PARI
    {a(n) = if( n<1, 0, -(-1)^n * sumdiv(n, d, kronecker(-12, d)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (1 - eta(x + A)^3 * eta(x^4 + A)^3 * eta(x^6 + A) / (eta(x^2 + A)^3 * eta(x^3 + A) * eta(x^12 + A))) / 3, n))}; /* Michael Somos, May 06 2015 */

Formula

Expansion of (1 - b(q^2)^2 / b(-q) ) / 3 in powers of q where b() is a cubic AGM function.
Moebius transform is period 12 sequence [ 1, -2, 0, 0, -1, 0, 1, 0, 0, 2, -1, 0, ...].
a(n) is multiplicative with a(2^e) = -1 unless 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} (-1)^k * (x^k + x^(3*k)) / (1 + x^k + x^(2*k)).
G.f.: ( Sum_{k>0} x^(6*k-5) / ( 1 + x^(6*k-5) ) - x^(6*k-1) / ( 1 + x^(6*k-1) )).
a(n) = -(-1)^n * A035178(n). -3 * a(n) = A132973(n) unless n = 0.
a(2*n) = -A035178(n). a(2*n + 1) = A033762(n). a(3*n) = a(n). a(3*n + 1) = A227696(n).
a(4*n + 1) + A112604(n). a(4*n + 3) = A112605(n). a(6*n + 1) = A097195(n). a(6*n + 5) = 0.
a(8*n + 1) = A112606(n). a(8*n + 3) = A112608(n). a(8*n + 5) = 2 * A112607(n-1). 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 * A131963(n). a(24*n + 19) = 2 * A131964(n).

A226132 Expansion of - c(-q) * c(q^2) / 9 in powers of q where c() is a cubic AGM theta function.

Original entry on oeis.org

1, -1, 3, -1, 6, -3, 8, -1, 9, -6, 12, -3, 14, -8, 18, -1, 18, -9, 20, -6, 24, -12, 24, -3, 31, -14, 27, -8, 30, -18, 32, -1, 36, -18, 48, -9, 38, -20, 42, -6, 42, -24, 44, -12, 54, -24, 48, -3, 57, -31, 54, -14, 54, -27, 72, -8, 60, -30, 60, -18, 62, -32, 72
Offset: 1

Views

Author

Michael Somos, May 27 2013

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Number 91 of the 126 eta-quotients listed in Table 1 of Williams 2012.

Examples

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

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, Sum[ If[ OddQ[d] && ! Divisible[ n/d, 3], -d (-1)^(n/d), 0], {d, Divisors[ n]}]];
    a[ n_] := If[ n < 2, Boole[n == 1], Times @@ (Which[ # == 2, -1, # == 3, #^#2, True, (#^(#2 + 1) - 1) / (# - 1)] & @@@ FactorInteger[n])];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q^3]^3 / EllipticTheta[ 3, 0, q] EllipticTheta[ 2, 0, I q^(3/2)]^3 / EllipticTheta[ 2, 0, I q^(1/2)] / (4 (-1)^(1/4)), {q, 0, n}];
  • PARI
    {a(n) = if( n<1, 0, sumdiv( n, d, if(d%2 && (n/d)%3, -d * (-1)^(n/d))))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, 0, A = factor(n); prod( k= 1, matsize(A)[1], if( p=A[k, 1], e=A[k, 2]; if( p==2, -1, if( p==3, p^e, (p^(e+1) - 1) / (p - 1))))))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x + A) * eta(x^4 + A) * eta(x^6 + A)^12 / (eta(x^2 + A)^4 * eta(x^3 + A)^3 * eta(x^12 + A)^3), n))};

Formula

Expansion of (a(q) - a(q^2)) * (a(q^2) + 2 * a(q^4)) / 18 = c(q^2)^4 / (9 * c(q) * c(q^4)) = (b(-q) * b(q^2) - a(-q) * a(q^2)) / 9 in powers of q where a(), b(), c(q) are cubic AGM theta functions.
Expansion of q * (phi(q^3)^3 / phi(q)) * (ps(-q^3)^3 / psi(-q)) in powers of q where phi(), psi() are Ramanujan theta functions.
Expansion of eta(q) * eta(q^4) * eta(q^6)^12 / (eta(q^2)^4 * eta(q^3)^3 * eta(q^12)^3) in powers of q.
Euler transform of period 12 sequence [ -1, 3, 2, 2, -1, -6, -1, 2, 2, 3, -1, -4, ...].
Multiplicative with a(2^e) = -1 if e>0, a(3^e) = 3^e, a(p^e) = (p^(e+1) - 1) / (p - 1) if p>3.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 4/3 (t/i)^2 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A226139.
G.f.: Sum_{k>0 not 3|k} x^k / (1 - (-x)^k)^2 = Sum_{k>0 not 2|k} k * x^k * (1 - x^k) / (1 + x^(3*k)).
G.f.: x * Product_{k>0} (1 - x^k) * (1 - x^(4*k)) * (1 - x^(3*k))^6 * (1 + x^(3*k))^9 / ((1 - x^(2*k))^4 * (1 + x^(6*k))^3).
a(2*n) = - A121443(n). a(2*n + 1) = A185717(n).
a(n) = -(-1)^n * A121443(n). Convolution of A113447 and A113973.
From Amiram Eldar, Sep 12 2023: (Start)
Dirichlet g.f.: (1 - 1/2^(s-1))^2 * (1 - 1/3^s) * zeta(s-1) * zeta(s).
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/54 = 0.18277045... . (End)

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Aug 04 2015

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. = x + x^2 - 2*x^3 - x^4 - 2*x^6 + 2*x^7 + x^8 - 2*x^9 + 2*x^12 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(36), 1), 80); A[2] + A[3] - 2*A[4] - A[5] - 2*A[7] + 2*A[8] + A[9] - 2*A[10] + 2*A[13] + 2*A[14] + 2*A[15] - A[17] - 2*A[19] - 4*A[20];
  • Mathematica
    a[ n_] := If[ n < 1, 0, Sum[ {1, 1, 0, -1, -1, 0}[[Mod[ d, 6, 1]]] {1, 0, -2, 0, 1, 0}[[Mod[ n/d, 6, 1]]], {d, Divisors @ n}]]
    a[ n_] := If[ n < 1, 0, Times @@ (Which[ # == 1, 1, # == 2, -(-1)^#2, # == 3, -2, Mod[#, 6] == 5, 1 - Mod[#2, 2], True, #2 + 1] & @@@ FactorInteger @ n)];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, Pi/4, q^(1/2)] EllipticTheta[ 2, Pi/4, q^(9/2)] EllipticTheta[ 3, 0, q] / (2 q^(1/4) QPochhammer[ q^6]), {q, 0, n}];
  • PARI
    {a(n) = if( n<1, 0, sumdiv(n, d, [0, 1, 1, 0, -1, -1][d%6 + 1] * [0, 1, 0, -2, 0, 1][n\d%6 + 1]))};
    
  • PARI
    {a(n) = if( n<1, 0, A = factor(n); prod(k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2, -(-1)^e, p==3, -2, p%6==5, 1-e%2, e+1)))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A)^4 * eta(x^9 + A) * eta(x^36 + A) / (eta(x + A) * eta(x^4 + A) * eta(x^6 + A) * eta(x^18 + A)), n))};
    

Formula

Expansion of (a(q) + a(q^2) - 3*a(q^3) - 2*a(q^4) - 3*a(q^6) + 6*a(q^12)) / 6 in powers of q where a() is a cubic AGM theta function.
Expansion of q * phi(q) * psi(-q) * psi(-q^9) / f(-q^6) in powers of q where phi(), psi(), f() are Ramanujan theta functions.
Expansion of eta(q^2)^4 * eta(q^9) * eta(q^36) / (eta(q) * eta(q^4) * eta(q^6) * eta(q^18)) in powers of q.
Euler transform of period 36 sequence [ 1, -3, 1, -2, 1, -2, 1, -2, 0, -3, 1, -1, 1, -3, 1, -2, 1, -2, 1, -2, 1, -3, 1, -1, 1, -3, 0, -2, 1, -2, 1, -2, 1, -3, 1, -2, ...].
Moebius transform is period 36 sequence [ 1, 0, -3, -2, -1, 0, 1, 2, 0, 0, -1, 6, 1, 0, 3, -2, -1, 0, 1, 2, -3, 0, -1, -6, 1, 0, 0, -2, -1, 0, 1, 2, 3, 0, -1, 0, ...].
a(n) is multiplicative with a(2^e) = -(-1)^e if e>0, 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. is a period 1 Fourier series which satisfies f(-1 / (36 t)) = 108^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A123863.
a(2*n) = A112848(n). a(2*n + 1) = A123530(n). a(3*n) = -2 * A113447(n). a(3*n + 1) = A227696(n).
a(4*n) = - A112848(n). a(4*n + 1) = A253243(n). a(4*n + 2) = A123530(n). a(4*n + 3) = -2 * A246838(n).
a(6*n) = -2 * A093829(n). a(6*n + 1) = A097195(n). a(6*n + 2) = A033687(n). a(6*n + 3) = -2 * A033762(n). a(6*n + 5) = 0.
a(8*n + 1) = A260941(n). a(8*n + 2) = A253243(n). a(8*n + 3) = -2 * A260943(n). a(8*n + 4) = - A123530(n). a(8*n + 5) = 2 * A260942(n). a(8*n + 6) = -2 * A246838(n). a(8*n + 7) = 2 * A260944(n).
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-6 of 6 results.