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 21-30 of 35 results. Next

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Nov 02 2005

Keywords

Comments

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

Examples

			G.f. = x + x^2 + x^4 + 2*x^7 + x^8 + 2*x^13 + 2*x^14 + x^16 + 2*x^19 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, If[ Mod[n, 3] == 0, 0, DivisorSum[ n, KroneckerSymbol[ -12, #] &]]]; (* Michael Somos, Jul 30 2015 *)
    a[ n_] := SeriesCoefficient[ x QPochhammer[ x^9]^3 / QPochhammer[ x^3] + x^2 QPochhammer[ x^18]^3 / QPochhammer[ x^6], {x, 0, n}]; (* Michael Somos, Jul 30 2015 *)
  • PARI
    {a(n) = if( n<1, 0, if( n%3, sumdiv(n,d, kronecker(-12, d))))};
    
  • PARI
    {a(n) = if( n<1, 0, direuler(p=2, n, if( p==3, 1, 1 / ((1 - X) * (1 - kronecker(-12, p)*X))))[n])}
    
  • 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==2, 1, p==3, 0, p%6==1, e+1, !(e%2))))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^2 + A)^2 * eta(x^9 + A) * eta(x^18 + A) / (eta(x + A) * eta(x^6 + A)), n))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x^9 + A)^3 / eta(x^3 + A) + x * eta(x^18 + A)^3 / eta(x^6 + A), n))};

Formula

Euler transform of period 18 sequence [ 1, -1, 1, -1, 1, 0, 1, -1, 0, -1, 1, 0, 1, -1, 1, -1, 1, -2, ...].
Moebius transform is period 18 sequence [ 1, 0, -1, 0, -1, 0, 1, 0, 0, 0, -1, 0, 1, 0, 1, 0, -1, 0, ...].
a(n) is multiplicative with a(2^e) = 1, a(3^e) = 0^e, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1 + (-1)^e)/2 if p == 5 (mod 6).
a(3*n) = 0, a(2*n) = a(n).
G.f.: Sum_{k>0} x^(6*k - 5) / (1 - x^(6*k - 5)) - x^(6*k - 1) / (1 - x^(6*k - 1)) - x^(18*k - 15) / (1 - x^(18*k - 15)) + x^(18*k - 6) / (1 - x^(18*k - 6)).
G.f.: Sum_{k>0} x^k * (1 - x^(2*k)) * (1 - x^(4*k)) * (1-x^(10*k)) / (1 - x^(18*k)).
Expansion of (c(q) + c(q^2))/3 in powers of q^(1/3) where c(q) is a cubic AGM theta function.
a(3*n + 1) = A033687(n). a(6*n + 1) = A097195(n). - Michael Somos, Jul 30 2015
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(3*sqrt(3)) = 0.604599... (A073010). - Amiram Eldar, Oct 15 2022

A138806 Expansion of (theta_3(q) * theta_3(q^27) + theta_2(q) * theta_2(q^27) - 1) / 2 in powers of q.

Original entry on oeis.org

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

Views

Author

Michael Somos, Mar 30 2008

Keywords

Comments

Half the number of integer solutions to x^2 + x*y + 7*y^2 = n. - Jianing Song, Nov 20 2019

Examples

			q + q^4 + 2*q^7 + 3*q^9 + 2*q^13 + q^16 + 2*q^19 + q^25 + 3*q^27 + ...
		

Crossrefs

Cf. A138805 (number of integer solutions to x^2 + x*y + 7*y^2 = n).
Similar sequences: A096936, A113406, A110399.

Programs

  • Mathematica
    f[p_, e_] := If[Mod[p, 6] == 1, e + 1, (1 + (-1)^e)/2]; f[2, e_] := 1 - Mod[e, 2]; f[3, e_] := 3; f[3, 1] = 0; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 07 2023 *)
  • PARI
    {a(n) = if( n<1, 0, if( n%3 == 2, 0, if( n%3==1, sumdiv(n, d, kronecker(-3, d)), if( n%9==0, 3 * sumdiv(n/9, d, kronecker(-3, d))))))}
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv(n, d, kronecker(-3, d)) - if( n%3==0, sumdiv(n/3, d, [0, 1, -1, -3, 1, -1, 3, 1, -1][d%9+1])))}
    
  • PARI
    {a(n) = if( n<1, 0, qfrep([2, 1; 1, 14], n, 1)[n])}

Formula

a(n) is multiplicative and a(3^e) = 3 if e>1, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1 + (-1)^e) / 2 if p == 5 (mod 6).
a(3*n + 2) = a(4*n + 2) = 0.
G.f.: (Sum_{i,j} x^(i*i + i*j + 7*j*j) - 1) / 2.
A138805(n) = 2 * a(n) unless n=0. A033687(n) = a(3*n + 1). A097195(n) = a(6*n + 1). A123884(n) = a(12*n + 1). 2 * A121361(n) = a(12*n + 7).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(3*sqrt(3)) = 0.604599... (A073010). - Amiram Eldar, Nov 16 2023

A123863 Expansion of (c(q^3) - c(q^6) - 2*c(q^12)) / 3 in powers of q where c(q) is a cubic AGM theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Oct 14 2006

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^4 + 2*q^7 - q^8 + 2*q^13 - 2*q^14 - q^16 + 2*q^19 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, Times @@ (Which[ # < 4, {1, -1, 0}[[Mod[#, 3, 1]]], Mod[#, 6] == 1, #2 + 1, True, (1 + (-1)^#2) / 2] & @@@ FactorInteger @ n)]; (* Michael Somos, Aug 03 2015 *)
    a[ n_] := SeriesCoefficient[ x EllipticTheta[ 2, 0, x^(9/2)] EllipticTheta[ 2, Pi/4, x^(1/2)] EllipticTheta[ 4, 0, x^18] / (2^(3/2) x^(5/4) QPochhammer[ x^6]), {x, 0, n}]; (* Michael Somos, Aug 03 2015 *)
  • 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^18 + A)^4 / (eta(x^2 + A) * eta(x^6 + A) * eta(x^9 + A) * eta(x^36 + A)), n))};
    
  • 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==2, -1, p==3, 0, p%6==1, e+1, !(e%2))))};

Formula

Expansion of (a(x) - a(x^2) - a(x^3) - 2*a(x^4) + a(x^6) + 2*a(x^12)) / 6 in powers of x where a() is a cubic AGM theta function. - Michael Somos, Aug 03 2015
Expansion of psi(-x) * psi(-x^9) * phi(x^9) / f(-x^6) in powers of x where phi(), psi(), f() are Ramanujan theta functions. - Michael Somos, Aug 03 2015
Expansion of eta(q) * eta(q^4) * eta(q^18)^4 / (eta(q^2) * eta(q^6) * eta(q^9) * eta(q^36)) in powers of q.
Euler transform of period 36 sequence [ -1, 0, -1, -1, -1, 1, -1, -1, 0, 0, -1, 0, -1, 0, -1, -1, -1, -2, -1, -1, -1, 0, -1, 0, -1, 0, 0, -1, -1, 1, -1, -1, -1, 0, -1, -2, ...].
a(n) is multiplicative with a(2^e) = -1 if e>0, a(3^e) = 0^e, a(p^e) = e+1 if p == 1 (mod 6), a(p^e) = (1 + (-1)^e)/2 if p == 5 (mod 6).
a(3*n) = a(6*n + 5) = 0.
a(2*n) = -A113448(n). a(6*n + 2) = -A033687(n).
a(3*n + 1) = A227696(n). a(6*n + 1) = A097195(n). a(12*n + 1) = A123884(n). a(12*n + 7) = 2 * A121361(n). - Michael Somos, Aug 03 2015

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Jan 22 2008

Keywords

Comments

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

Examples

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

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, DivisorSum[ n, (Mod[#, 2] - 4 Boole[Mod[#, 8] == 4]) KroneckerSymbol[ -3, n/#] &]]; (* Michael Somos, Oct 12 2015 *)
    a[ n_] := If[ n < 1, 0, Times @@ (Which[# == 1 || # == 3, 1, # == 2, If[#2 < 2, -1, -3 (-1)^#2], Mod[#, 6] == 1, #2 + 1, True, 1 - Mod[#2, 2]] & @@@ FactorInteger@n)]; (* Michael Somos, Oct 12 2015 *)
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( eta(x + A) * eta(x^3+A) * eta(x^4 + A)^4 * eta(x^24 + A)^2 / ( eta(x^2 + A) * eta(x^8+A) * eta(x^12+A) )^2, n))};
    
  • 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==2, if( e<2, -1, -3 * (-1)^e), p==3, 1, p%6>1, !(e%2), e+1)))};
    
  • PARI
    {a(n) = if( n<1, 0, sumdiv(n, d, ((d%2) -4 * (d%8==4)) * kronecker(-3, n/d)))};

Formula

Expansion of eta(q) * eta(q^3) * eta(q^4)^4 * eta(q^24)^2 / (eta(q^2) * eta(q^8) * eta(q^12))^2 in powers of q.
Euler transform of period 24 sequence [ -1, 1, -2, -3, -1, 0, -1, -1, -2, 1, -1, -2, -1, 1, -2, -1, -1, 0, -1, -3, -2, 1, -1, -2, ...].
a(n) is multiplicative with a(2) = -1, a(2^e) = -3 * (-1)^e if e>1, 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. is a period 1 Fourier series which satisfies f(-1 / (24 t)) = 12^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A123484.
G.f.: x * Product_{k>0} (1 - x^k)^2 * (1 + x^(2*k))^2 * (1 + x^k + x^(2*k)) * (1 - x^(4*k) + x^(8*k))^2.
Moebius transform is period 24 sequence [ 1, -2, 0, -2, -1, 0, 1, 6, 2, -1, 0, 1, -2, 0, -6, -1, 0, 1, 2, 0, 2, -1, 0, ...].
a(2*n) = A244375(n). a(2*n + 1) = A033762(n). a(3*n) = a(n). a(3*n + 1) = A122861(n).
a(4*n) = -3 * A093829(n). a(4*n + 1) = A112604(n). a(4*n + 2) = -A033762(n). a(4*n + 3) = A112605(n).
a(6*n + 1) = A097195(n). a(6*n + 5) = 0.
Expansion of q * f(-q, -q) * f(q^2, q^10) / f(-q, -q^5)^2 in powers of q where f(, ) is Ramanujan's general theta function. - Michael Somos, Oct 12 2015
Sum_{k=1..n} abs(a(k)) ~ (Pi*sqrt(3)/4) * n. - Amiram Eldar, Jan 28 2024

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

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

Original entry on oeis.org

1, -4, 6, -4, 0, 0, 6, -8, 6, -4, 0, 0, 0, -8, 12, 0, 0, 0, 6, -8, 0, -8, 0, 0, 6, -4, 12, -4, 0, 0, 0, -8, 6, 0, 0, 0, 0, -8, 12, -8, 0, 0, 12, -8, 0, 0, 0, 0, 0, -12, 6, 0, 0, 0, 6, 0, 12, -8, 0, 0, 0, -8, 12, -8, 0, 0, 0, -8, 0, 0, 0, 0, 6, -8, 12, -4, 0
Offset: 0

Views

Author

Michael Somos, Jun 26 2014

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

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, Boole[n == 0], 2 Sum[ (-1)^(n/d) {2, -1, 0, 1, -2, 0}[[ Mod[ d, 6, 1] ]], {d, Divisors @ n}]];
    a[ n_] := SeriesCoefficient[ QPochhammer[ q]^4 QPochhammer[ q^4] QPochhammer[ q^6]^8 / (QPochhammer[ q^2]^4 QPochhammer[ q^3]^4 QPochhammer[ q^12]^3), {q, 0, n}];
  • PARI
    {a(n) = if( n<1, n==0, 2 * sumdiv(n, d, (-1)^(n/d) * [0, 2, -1, 0, 1, -2][d%6 + 1]))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^4 * eta(x^4 + A) * eta(x^6 + A)^8 / (eta(x^2 + A)^4 * eta(x^3 + A)^4 * eta(x^12 + A)^3), n))};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 + 2 * sum(k=1, n, x^k / (1 + x^k) * [0, -2, 1, 0, -1, 2][k%6 + 1], x * O(x^n)), n))};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 + 2 * sum(k=1, n, x^k / (1 + x^k + x^(2*k)) * [3, -2, 1, -2][k%4 + 1], x * O(x^n)), n))};
    
  • PARI
    {a(n) = my(A);  if( n<1, n==0, A = factor(n); -4 * prod( j=1, matsize(A)[1], if( p = A[j,1], e = A[j,2]; if( p==2, (1 - (-1)^e) * -3/4, if( p==3, 1, if( p%6 == 1, e+1, (1 + (-1)^e) / 2))))))};

Formula

Expansion of b(q) * (b(q) + 2*b(q^4)) / (3 * b(q^2)) in powers of q where b() is a cubic AGM theta function.
Expansion of psi(-q) * chi(-q)^3 * phi(q^3) * chi(q^3)^3 in powers of q where phi(), psi(), chi() are Ramanujan theta functions.
Expansion of eta(q)^4 * eta(q^4) * eta(q^6)^8 / (eta(q^2)^4 * eta(q^3)^4 * eta(q^12)^3) in powers of q.
Euler transform of period 12 sequence [ -4, 0, 0, -1, -4, -4, -4, -1, 0, 0, -4, -2, ...].
Moebius transform is period 12 sequence [ -4, 10, 0, -6, 4, 0, -4, 6, 0, -10, 4, 0, ...].
a(n) = -4 * b(n) where b(n) is multiplicative with b(2^e) = (1 - (-1)^e) * -3/4 if e>0, b(3^e) = 1, b(p^e) = e+1 if p == 1 (mod 6), b(p^e) = (1 + (-1)^e) / 2 if p == 5 (mod 6) with a(0) = 1.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 48^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A244375.
a(2*n) = A004016(n). a(2*n + 1) = -4 * A033762(n). a(3*n) = a(n). a(6*n + 1) = -4 * A097195(n). a(6*n + 5) = 0.
Sum_{k=1..n} abs(a(k)) ~ (2*Pi/sqrt(3)) * n. - Amiram Eldar, Jun 08 2025

A246650 Expansion of phi(x) * chi(-x) * psi(x^3) in powers of x where phi(), psi(), chi() are Ramanujan theta functions.

Original entry on oeis.org

1, 1, -2, 0, 2, -3, -2, 0, 1, 2, -2, 0, 2, 0, -2, 0, 3, 2, 0, 0, 2, -3, -2, 0, 2, 2, -2, 0, 0, 0, -4, 0, 2, 1, -2, 0, 2, -6, 0, 0, 1, 2, -2, 0, 4, 0, -2, 0, 0, 2, -2, 0, 2, 0, -2, 0, 3, 2, -2, 0, 2, 0, 0, 0, 2, 3, -2, 0, 0, -6, -2, 0, 4, 0, -2, 0, 2, 0, 0, 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 + x - 2*x^2 + 2*x^4 - 3*x^5 - 2*x^6 + x^8 + 2*x^9 - 2*x^10 + ...
G.f. = q + q^4 - 2*q^7 + 2*q^13 - 3*q^16 - 2*q^19 + q^25 + 2*q^28 - ...
		

Crossrefs

Programs

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

Formula

Expansion of q^(-1/3) * eta(q^2)^4 * eta(q^6)^2 / (eta(q) * eta(q^3) * eta(q^4)^2) in powers of q.
a(2*n) = A129451(n). a(4*n) = A123884(n). a(4*n + 1) = A122861(n). a(4*n + 2) = -2 * A121361(n). a(4*n + 3) = 0.
a(8*n) = A131961(n). a(8*n + 1) = A097195(n). a(8*n + 2) = -2 * A131962(n). a(8*n + 4) = 2 * A131963(n). a(8*n + 6) = -2 * A131964(n).
a(16*n + 1) = A123884(n). a(16*n + 5) = -3 * A033687(n). a(16*n + 9) = 2 * A121361(n). a(16*n + 13) = 0.

A253625 Expansion of psi(q^2) * f(-q, q^2)^2 / f(-q, -q^5) in powers of q where psi(), f() are Ramanujan theta functions.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jan 06 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. = 1 - 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), 81); A[1] - A[2] + 3*A[3] - A[4] + 3*A[5];
  • Mathematica
    a[ n_] := If[ n < 1, Boole[n == 0], (-1)^ n Sum[(-1)^ Quotient[ d, 3] {1, 1, 0}[[ Mod[d, 3, 1] ]] {1, 2}[[ Mod[n/d, 2, 1] ]], {d, Divisors @ n}]];
    a[ n_] := SeriesCoefficient[ QPochhammer[ q] QPochhammer[ q^3] (QPochhammer[ -q^3, q^6] QPochhammer[ -q^2, q^2])^4, {q, 0, n}];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q^3]^2 EllipticTheta[ 2, 0, q]^2 / (EllipticTheta[ 2, 0, q^(1/2)] EllipticTheta[ 2, 0, q^(3/2)]), {q, 0, n}];
  • PARI
    {a(n) = if( n<1, n==0, (-1)^n * sumdiv(n, d, (-1)^(d\3) * (d%3>0) * (2-(n\d)%2)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A) * eta(x^4 + A)^4 * eta(x^6 + A)^8 / (eta(x^2 + A)^4 * eta(x^3 + A)^3 * eta(x^12 + A)^4), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, n==0, A = factor(n); - prod( k=1, matsize(A)[1], if( p=A[k, 1], e=A[k, 2]; if( p==2, -3, if( p==3, 1, if( p%6 == 1, e+1, 1-e%2))))))};
    

Formula

Expansion of psi(q^2)^2 * phi(q^3)^2 / (psi(q) * psi(q^3)) = f(-q) * f(-q^3) * (chi(q^3) / chi(-q^2))^4 in powers of q where phi(), psi(), chi(), f() are Ramanujan theta functions.
Expansion of (-a(q) - 3*a(q^2) + 4*a(q^4)) / 6 = b(q^4) * (b(q) + 2*b(q^4)) / (3*b(q^2)) in powers of q where a(), b() are cubic AGM theta functions.
Expansion of eta(q) * eta(q^4)^4 * eta(q^6)^8 / (eta(q^2)^4 * eta(q^3)^3 * eta(q^12)^4) in powers of q.
Euler transform of period 12 sequence [ -1, 3, 2, -1, -1, -2, -1, -1, 2, 3, -1, -2, ...].
Moebius transform is period 12 sequence [ -1, 4, 0, 0, 1, 0, -1, 0, 0, -4, 1, 0, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 48^(-1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A253623.
a(n) = -b(n) where b() is multiplicative with b(2^e) = -3 if e>0, b(3^e) = 1, b(p^e) = e+1 if p == 1 (mod 6), b(p^e) = (1 + (-1)^e)/2 if p == 5 (mod 6).
G.f.: 1 + Sum_{k>0} (3 - (k mod 2)*4) * (x^k + x^(3*k)) / (1 + x^(2*k) + x^(4*k)).
G.f.: Product_{k>0} (1 - x^k) * (1 - x^(3*k)) * (1 + x^(3*k))^4 / (1 - x^(2*k) + x^(4*k))^4.
a(n) = (-1)^n * A253626(n). a(2*n) = A107760(n). a(2*n + 1) = - A033762(n). a(3*n) = a(n). a(3*n + 1) = - A122861(n). a(4*n + 1) = - A112604(n). a(4*n + 2) = 3 * A033762(n). a(4*n + 3) = - A112605(n).
a(6*n + 1) = - A097195(n). a(6*n + 2) = 3 * A033687(n). a(6*n + 5) = 0. a(12*n + 1) = - A123884(n). a(12*n + 7) = -2 * A121361(n). a(12*n + 10) = 0.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=0..m} a(k) = Pi/(2*sqrt(3)) = 0.906899... (A093766). - Amiram Eldar, Jun 08 2025

Extensions

Typo in formula fixed by Colin Barker, Jan 08 2015

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Sep 23 2013

Keywords

Comments

Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Rogers and Zudilin (2011) page 6: "This identity can be verified by eliminating b(q) with b(q^{1/3}) - b(q) = 3c(q^3) - c(q)."
The zeros of the g.f. A(q) where q = exp(2 Pi i t) are of the form t = (m/2 + sqrt(-3)/18) / n where m is an odd integer and n is in A004611. For example, (1/2 + sqrt(-3)/18) / 1, (1/2 + sqrt(-3)/18) / 7, (5/2 + sqrt(-3)/18) / 13.

Examples

			G.f. = q - 3*q^3 + q^4 + 2*q^7 - 3*q^12 + 2*q^13 + q^16 + 2*q^19 - 6*q^21 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(27), 1), 85); A[2] - 3*A[4] + A[5] + 2*A[8] - 3*A[13] + 2*A[14] + A[15]; /* Michael Somos, Jun 16 2015 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ (QPochhammer[ q^3]^4 - QPochhammer[ q^9] QPochhammer[ q]^3) / (3 QPochhammer[ q^3] QPochhammer[ q^9]), {q, 0, n}]; (* Michael Somos, Jun 16 2015 *)
    a[ n_] := SeriesCoefficient[ q (QPochhammer[ q^9]^4 - 3 q^2 QPochhammer[ q^3] QPochhammer[ q^27]^3) / (QPochhammer[ q^3] QPochhammer[ q^9]), {q, 0, n}]; (* Michael Somos, Jun 16 2015 *)
    f[p_, e_] := If[Mod[p, 3] == 1, e+1, (1 + (-1)^e) / 2]; f[3, 1] = -3; f[3, e_] := 0; a[n_] := Times @@ f @@@ FactorInteger[n]; a[0] = 0; a[1] = 1; Array[a, 100, 0] (* Amiram Eldar, Sep 04 2023 *)
  • 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==3, -3 * (e==1), p%3==1, e+1, !(e%2))))};
    
  • PARI
    {a(n) = my(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x^9 + A)^4 - 3 * x^2 * eta(x^3 + A) * eta(x^27 + A)^3) / (eta(x^3 + A) * eta(x^9 + A)), n))};
    

Formula

Expansion of c(q^3) / 3 - c(q^9) in powers of q where c() is a cubic AGM theta function.
Expansion of (a(q) - 4*a(q^3) + 3*a(q^9)) / 6 in powers of q where a() is a cubic AGM theta function.
Expansion of (eta(q^3)^4 - eta(q)^3 * eta(q^9)) / (3 * eta(q^3) * eta(q^9)) in powers of q.
a(n) is multiplicative with a(3) = -3, a(3^e) = 0 if e>1, a(p^e) = e+1 if p == 1 (mod 3), a(p^e) = (1 + (-1)^e) / 2 if p == 2 (mod 3).
G.f. is a period 1 Fourier series which satisfies f(-1 / (27 t)) = 27^(1/2) (t/i) f(t) where q = exp(2 Pi i t).
a(3*n + 2) = a(4*n + 2) = a(9*n) = a(9*n + 6) = 0. a(3*n + 1) = A033687(n). a(9*n + 3) = -3 * A033687(n).
From Michael Somos, Jun 16 2015: (Start)
a(4*n) = a(n). a(6*n + 1) = A097195(n). a(12*n + 1) = A123884(n). a(12*n + 7) = 2 * A121361(n).
a(n) = Sum_{d|n} A259024(n/d) * [ 0, 1, 0, -2, 0, 1][mod(d, 6) + 1]. (End)

A253626 Expansion of psi(q^2) * f(q, q^2)^2 / f(q, q^5) in powers of q where psi(), f() are Ramanujan theta functions.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jan 06 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. = 1 + q + 3*q^2 + q^3 + 3*q^4 + 3*q^6 + 2*q^7 + 3*q^8 + q^9 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(12), 1), 86); A[1] + A[2] + 3*A[3] + A[4] + 3*A[5];
  • Mathematica
    a[ n_] := If[ n < 1, Boole[n == 0], Sum[ (-1)^ Quotient[ d, 3] {1, 1, 0}[[ Mod[d, 3, 1] ]] {1, 2}[[ Mod[n/d, 2, 1] ]], {d, Divisors @ n}]];
    a[ n_] := SeriesCoefficient[ QPochhammer[ -q] QPochhammer[ -q^3] (QPochhammer[ q^3, q^6] QPochhammer[ -q^2, q^2])^4, {q, 0, n}];
  • PARI
    {a(n) = if( n<1, n==0, sumdiv(n, d, kronecker(-12, d) + if(d%2, 0, 2 * kronecker(-12, d/2))))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^3 + A)^3 * eta(x^4 + A)^3 / ( eta(x + A) * eta(x^2 + A) * eta(x^6 + A) * eta(x^12 + A) ), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<1, n==0, A = factor(n); prod( k=1, matsize(A)[1], if( p=A[k, 1], e=A[k, 2]; if( p==2, 3, if( p==3, 1, if( p%6 == 1, e+1, 1-e%2))))))};
    

Formula

Expansion of psi(q^2)^2 * phi(-q^3)^2 / (psi(-q) * psi(-q^3)) = f(q) * f(q^3) * (chi(-q^3) / chi(-q^2))^4 in powers of q where phi(), psi(), chi(), f() are Ramanujan theta functions.
Expansion of (a(q) + 3*a(q^2) + 2*a(q^4)) / 6 = b(q^4) * (-b(q) + 4*b(q^4)) / (3*b(q^2)) in powers of q where a(), b() are cubic AGM theta functions.
Expansion of eta(q^3)^3 * eta(q^4)^3 / ( eta(q) * eta(q^2) * eta(q^6) * eta(q^12) ) in powers of q.
Euler transform of period 12 sequence [ 1, 2, -2, -1, 1, 0, 1, -1, -2, 2, 1, -2, ...].
Moebius transform is period 12 sequence [ 1, 2, 0, 0, -1, 0, 1, 0, 0, -2, -1, 0, ...].
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).
a(n) is multiplicative with a(0) = 1, a(2^e) = 3 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.: 1 + Sum_{k>0} (3 - (k mod 2)*2) * (q^k + q^(3*k)) / (1 + q^(2*k) + q^(4*k)).
G.f.: Product_{k>0} (1 - q^(3*k))^3 * (1 - q^(4*k))^3 / ( (1 - q^k) * (1 - q^(2*k)) * (1 - q^(6*k)) * (1 - q^(12*k)) ).
a(n) = (-1)^n * A253625(n). a(2*n) = A107760(n). a(2*n + 1) = A033762(n). a(3*n) = a(n). a(3*n + 1) = A122861(n). a(4*n + 1) = A112604(n). a(4*n + 2) = 3 * A033762(n). a(4*n + 3) = A112605(n).
a(6*n + 1) = A097195(n). a(6*n + 2) = 3 * A033687(n). a(6*n + 5) = 0. a(12*n + 1) = A123884(n). a(12*n + 7) = 2 * A121361(n). a(12*n + 10) = 0.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/sqrt(3) = 1.813799... (A093602). - Amiram Eldar, Jan 21 2024
Previous Showing 21-30 of 35 results. Next