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 11 results. Next

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Mar 11 2007

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^3 + x^4 - 2*x^5 + 2*x^7 - 2*x^10 + 3*x^12 + ...
G.f. = q - q^3 - 2*q^5 + 2*q^7 + q^9 - 2*q^11 + 2*q^15 - 2*q^21 + 3*q^25 + ...
		

Crossrefs

Programs

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

Formula

a(n) = b(2*n + 1) where b() is multiplicative with b(2^e) = 0^e, b(3^e) = (-1)^e, b(p^e) = e+1 if p == 1, 7 (mod 24), b(p^e) = (e+1)* (-1)^e if p == 5, 11 (mod 24), b(p^e) = (1+(-1)^e)/2 if p == 13, 17, 19, 23 (mod 24).
Euler transform of period 24 sequence [ -1, -2, 0, 0, -1, -1, -1, -1, 0, -2, -1, -2, -1, -2, 0, -1, -1, -1, -1, 0, 0, -2, -1, -2, ...].
a(12*n + 6) = a(12*n + 8) = a(12*n + 9) = a(12*n + 11)= 0.
G.f.: Product_{k>0} (1 - x^(8*k)) * (1 - x^(12*k))^2 / ((1 + x^k) * (1 + x^(2*k))^2 * (1 - x^(3*k)) * (1 + x^(12*k))).
G.f.: Sum_{k>=0} a(k) * x^(2*k + 1) = Sum_{k>0} (x^k - x^(3*k)) / (1 + x^(4*k)) * Kronecker(-12, k) = Sum_{k>0} (x^k + x^(3*k)) / (1 + x^(2*k) + x^(4*k)) * Kronecker(2, k).
a(n) = A128581(2*n + 1) = A115660(2*n + 1). a(3*n + 2) = -2 * A128582(n). a(12*n) = A113780(n).
a(2*n) = A259668(n). a(3*n + 1) = - A128580(n). - Michael Somos, Jul 12 2015

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Mar 11 2007

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 - 2*q^5 - q^6 + 2*q^7 + q^8 + q^9 - 2*q^10 - ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 1, 0, -(-1)^n DivisorSum[ n, KroneckerSymbol[ #, 8] KroneckerSymbol[ n/#, 3] &]]; (* Michael Somos, Nov 15 2015 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q^2] EllipticTheta[ 4, 0, q^3] - EllipticTheta[ 4, 0, q] EllipticTheta[ 3, 0, q^6])/2, {q, 0, n}]; (* Michael Somos, Nov 15 2015 *)
    a[ n_] := SeriesCoefficient[ q QPochhammer[ -q] QPochhammer[ q^24] QPochhammer[ q^4, q^8] QPochhammer[ q^3, -q^3], {q, 0, n}]; (* Michael Somos, Nov 15 2015 *)
  • PARI
    {a(n) = -(-1)^n * sumdiv(n, d, kronecker(d, 8) * kronecker(n/d, 3))}
    
  • PARI
    {a(n) = my(A =  x * O(x^n)); polcoeff( eta(x^2 + A)^3 * eta(x^3 + A) * eta(x^12 + A) * eta(x^24 + A) / (eta(x + A) * eta(x^6 + A)^2 * eta(x^8 + A)), n)}
    
  • PARI
    A128581(n)={prod(i=1,matsize(n=factor(n))[1], if(12>n[i,1]%24, if(bittest(12,n[i,1]),(-1)^(n[i,1]+n[i,2]-1), if(bittest(n[i,2],0)&&n[i,1]%6>1,-n[i,2]-1,n[i,2]+1)), !bittest(n[i,2],0)))} \\ a(p^e) as given in formula. - M. F. Hasler, May 07 2018

Formula

Expansion of eta(q^2)^3 * eta(q^3) * eta(q^12) * eta(q^24) / (eta(q) * eta(q^6)^2 * eta(q^8)) in powers of q.
Euler transform of period 24 sequence [ 1, -2, 0, -2, 1, -1, 1, -1, 0, -2, 1, -2, 1, -2, 0, -1, 1, -1, 1, -2, 0, -2, 1, -2, ...].
Multiplicative with a(2^e) = -(-1)^e if e>0, a(3^e) = (-1)^e, a(p^e) = e+1 if p == 1, 7 (mod 24), a(p^e) = (e+1)(-1)^e if p == 5, 11 (mod 24), a(p^e) = (1 + (-1)^e)/2 if p == 13, 17, 19, 23 (mod 24).
a(n) = -(-1)^n * A115660(n). a(2*n) = A115660(n). a(2*n + 1) = A128580(n).
abs(a(n)) = A000377(n) = A192013(n). - M. F. Hasler, May 07 2018

A259668 Expansion of psi(-x)^2 * psi(x^3)^2 / (phi(-x^4) * psi(-x^6)) in power of x where phi(), psi() are Ramanujan theta functions.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jul 02 2015

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^5 + 3*x^6 - 2*x^7 + 2*x^8 - 2*x^11 + 3*x^12 + ...
G.f. = q - 2*q^5 + q^9 - 2*q^21 + 3*q^25 - 2*q^29 + 2*q^33 - 2*q^45 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ x^8] QPochhammer[x^12, x^24] QPochhammer[ x^6] (QPochhammer[ x, x^6] QPochhammer[ x^5, x^6])^2, {x, 0, n}];
    a[ n_] := SeriesCoefficient[ (EllipticTheta[2, 0, x^(3/2)] EllipticTheta[ 2, Pi/4, x^(1/2)])^2 / (2^(5/2) x^(1/4) EllipticTheta[ 4, 0, x^4] EllipticTheta[ 2, Pi/4, x^3]), {x, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^6 + A)^3 * eta(x^8 + A) * eta(x^12 + A) / (eta(x^2 + A)^2 * eta(x^3 + A)^2 * eta(x^24 + A)), n))};

Formula

Expansion of f(-x^8) * f(-x, -x^5)^2 / psi(-x^6) in powers of x where psi(), f() are Ramanujan theta functions.
Euler transform of period 24 sequence [ -2, 0, 0, 0, -2, -1, -2, -1, 0, 0, -2, -2, -2, 0, 0, -1, -2, -1, -2, 0, 0, 0, -2, -2, ...].
a(n) = A128580(2*n) = A134177(2*n) = A115660(4*n) = A128581(4*n).
a(6*n + 1) = -2 * A113780(n).

A263548 Expansion of f(x, x) * f(x^2, x^10) in powers of x where f(, ) is Ramanujan's general theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Oct 20 2015

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 + 2*x^4 + 2*x^6 + 2*x^9 + x^10 + 4*x^11 + ...
G.f. = q^2 + 2*q^5 + q^8 + 2*q^11 + 2*q^14 + 2*q^20 + 2*q^29 + q^32 + ...
		

Crossrefs

Programs

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

Formula

Expansion of phi(x) * chi(x^2) * psi(-x^6) in powers of x where phi(), psi(), chi() are Ramanujan theta functions.
Expansion of q^(-2/3) * eta(q^2)^4 * eta(q^6) * eta(q^24) / (eta(q)^2 * eta(q^8) * eta(q^12)) in powers of q.
Euler transform of period 24 sequence [ 2, -2, 2, -2, 2, -3, 2, -1, 2, -2, 2, -2, 2, -2, 2, -1, 2, -3, 2, -2, 2, -2, 2, -2, ...].
a(n) = A000377(3*n + 2) = A192013(3*n + 2) = A115660(6*n + 4).

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Jul 22 2015

Keywords

Comments

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

Examples

			G.f. = 1 + x + 3*x^3 + 2*x^4 + 3*x^6 + 2*x^9 + x^10 + 2*x^12 + 4*x^13 + ...
G.f. = q + q^9 + 3*q^25 + 2*q^33 + 3*q^49 + 2*q^73 + q^81 + 2*q^97 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, DivisorSum[ 8 n + 1, KroneckerSymbol[ -6, #] &]];
    a[ n_] := If[ n < 0, 0, Times @@ (Which[ # <= 3, Mod[#, 2], Mod[#, 24] > 12, 1 - Mod[#2, 2], True, (#2  + 1) KroneckerSymbol[ 3, #]^#2] & @@@ FactorInteger @ (8 n + 1))];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, x^3] EllipticTheta[ 2, 0, x^(1/2)] / (2 x^(1/8)), {x, 0, n}];
  • PARI
    {a(n) = if( n<0, 0, sumdiv( 8*n + 1, d, kronecker( -6, d)))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<0, 0, factor(8*n + 1); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2, 0, p==3, 1, p%24>12, !(e%2), (e+1) * kronecker(3, p)^e)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^2 * eta(x^6 + A)^5 / (eta(x + A) * eta(x^3 + A)^2 * eta(x^12 + A)^2), n))};

Formula

Expansion of q^(-1/8) * eta(q^2)^2 * eta(q^6)^5 / (eta(q) * eta(q^3)^2 * eta(q^12)^2) in powers of q.
Euler transform of period 12 sequence [ 1, -1, 3, -1, 1, -4, 1, -1, 3, -1, 1, -2, ...].
a(n) = A259668(2*n) = A128580(4*n) = A129402(4*n) = A134177(4*n) = A190615(4*n) = A115660(8*n + 1) = A128581(8*n + 1) = A192013(8*n + 1).

A261115 Expansion of f(x, x) * f(x^4, x^8) in powers of x where f(,) is Ramanujan's general theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Aug 08 2015

Keywords

Comments

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

Examples

			G.f. = 1 + 2*x + 3*x^4 + 2*x^5 + 3*x^8 + 4*x^9 + 2*x^12 + 2*x^13 + 2*x^16 + ...
G.f. = q + 2*q^7 + 3*q^25 + 2*q^31 + 3*q^49 + 4*q^55 + 2*q^73 + 2*q^79 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, x] EllipticTheta[ 4, 0, x^12] / QPochhammer[ x^4, x^8], {x, 0, n}];
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^5 * eta(x^8 + A) * eta(x^12 + A)^2 / (eta(x + A)^2 * eta(x^4 + A)^3 * eta(x^24 + A)), n))};

Formula

Expansion of q^(-1/6) * eta(q^2)^5 * eta(q^8) * eta(q^12)^2 / (eta(q)^2 * eta(q^4)^3 * eta(q^24)) in powers of q.
Euler transform of period 24 sequence [ 2, -3, 2, 0, 2, -3, 2, -1, 2, -3, 2, -2, 2, -3, 2, -1, 2, -3, 2, 0, 2, -3, 2, -2, ...].
a(n) = (-1)^n * A260110(n) = A128580(3*n) = A129402(3*n) = A115660(6*n + 1) = A128581(6*n + 1) = A192013(6*n + 1).
a(4*n) = A113780(n). a(4*n + 1) = 2 * A260089(n). a(4*n + 2) = a(4*n + 3) = 0.

A263571 Expansion of f(x^2, x^2) * f(x, x^5) in powers of x where f(, ) is Ramanujan's general theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Oct 21 2015

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^3 + x^5 + 2*x^7 + 3*x^8 + 2*x^9 + 2*x^10 + ...
G.f. = q + q^4 + 2*q^7 + 2*q^10 + q^16 + 2*q^22 + 3*q^25 + 2*q^28 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, With[ {m = 3 n + 1}, Sum[ KroneckerSymbol[ 2, d] KroneckerSymbol[ -3, m/d], {d, Divisors[ m]}]]];
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^2] EllipticTheta[ 3, 0, x^2] EllipticTheta[ 2, Pi/4, x^(3/2)] / (2^(1/2) x^(3/8)), {x, 0, n}];
  • PARI
    {a(n) = my(m); if( n<0, 0, m = 3*n + 1; sumdiv( m, d, kronecker( 2, d) * kronecker( -3, m/d)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^3 + A) * eta(x^4 + A)^4 * eta(x^12 + A) / (eta(x + A) * eta(x^8 + A)^2), n))};

Formula

Expansion of chi(x) * phi(x^2) * psi(-x^3) in powers of x where phi(), psi(), chi() are Ramanujan theta functions.
Expansion of q^(-1/3) * eta(q^3) * eta(q^4)^4 * eta(q^12) / (eta(q) * eta(q^6) * eta(q^8)^2) in powers of q.
a(n) = b(3*n + 1) where b() is multiplicative with b(2^e) = b(3^e) = (-1)^e, b(p^e) = e+1 if p == 1, 7 (mod 24), b(p^e) = (e+1) * (-1)^e if p == 5, 11 (mod 24), b(p^e) = (1 + (-1)^e) / 2 if p == 13, 17, 19, 23 (mod 24).
G.f. is a period 1 Fourier series which satisfies f(-1 / (72 t)) = 24^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A263577.
a(n) = A115660(3*n + 1) = A192013(3*n + 1) = A128581(6*n + 2).
a(2*n) = A261115(n). a(2*n + 1) = A263548(n). a(4*n + 1) = a(n). a(4*n + 3) = 2 * A128582(n).
a(8*n + 4) = a(8*n + 6) = 0. a(8*n) = A113780(n). a(8*n + 2) = 2 * A260089(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/sqrt(6) = 1.282549... . - Amiram Eldar, Dec 28 2023

A046113 Coefficients in expansion of theta_3(q) * theta_3(q^6) in powers of q.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, May 18 2002

Keywords

Comments

Number of representations of n as a sum of six times a square and a square. - Ralf Stephan, May 14 2007
a(n) < 2 if and only if n is in A002480. a(n) > 0 if and only if n is in A002481. - Michael Somos, Mar 01 2011

Examples

			G.f. = 1 + 2*x + 2*x^4 + 2*x^6 + 4*x^7 + 2*x^9 + 4*x^10 + 4*x^15 + 2*x^16 + ...
		

References

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

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^6], {q, 0, n}]; (* Michael Somos, Apr 19 2015 *)
  • PARI
    {a(n) = my(G); if( n<0, 0, G = [ 1, 0; 0, 6]; polcoeff( 1 + 2 * x * Ser( qfrep( G, n)), n))}; /* Michael Somos, Mar 01 2011 */

Formula

G.f.: Sum_{ i, j = -oo..+oo } q^(i^2 + 6*j^2).
a(n) = A000377(n) + A115660(n). - Michael Somos, Mar 01 2011
a(0) = 1, a(n) = (1+(-1)^t)*b(n) for n > 0, where t is the number of prime factors of n, counting multiplicity, which are == 2,3,5,11 (mod 24), and b() is multiplicative with b(p^e) = (e+1) for primes p == 1,5,7,11 (mod 24) and b(p^e) = (1+(-1)^e)/2 for primes p == 13,17,19,23 (mod 24). (This formula is Corollary 4.2 in the Berkovich-Yesilyurt paper). - Jeremy Lovejoy, Nov 14 2024

A108563 Number of representations of n as sum of twice a square plus thrice a square.

Original entry on oeis.org

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

Views

Author

Ralf Stephan, May 13 2007

Keywords

Comments

Number of solutions to n = 2*a^2 + 3*b^2 in integers.
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
a(n) > 0 if and only if n is in A002480. a(n) < 2 if n is in A002481. - Michael Somos, Mar 01 2011

Examples

			G.f. = 1 + 2*x^2 + 2*x^3 + 4*x^5 + 2*x^8 + 4*x^11 + 2*x^12 + 4*x^14 + 2*x^18 + ...
a(0) = 1 since 0 = 2*0^2 + 3*0^2, a(5) = 4 since 5 = 2*1^2 + 3*1^2 = 2*(-1)^2 + 3*1^2 = 2*1^2 + 3*(-1)^2 = 2*(-1^2) + 3*(-1)^2.
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q^2] EllipticTheta[ 3, 0, q^3], {q, 0, n}]; (* Michael Somos, Apr 19 2015 *)
    a[n_] := Module[{a, b, r}, r = Reduce[n == 2a^2 + 3b^2, {a, b}, Integers]; Which[r === False, 0, r[[0]] === And, 1, r[[0]] === Or, Length[r]]];
    Table[a[n], {n, 0, 105}] (* Jean-François Alcover, Jan 09 2019 *)
  • PARI
    for(n=0,120,print1(if(n<1,n==0,qfrep([2,0;0,3],n)[n]*2),","))
    
  • PARI
    {a(n) = my(G); if( n<0, 0, G = [2, 0; 0, 3]; polcoeff( 1 + 2 * x * Ser( qfrep( G, n)), n))}; /* Michael Somos, Mar 01 2011 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^4 + A)^5 * eta(x^6 + A)^5 / (eta(x^2 + A)^2 * eta(x^3 + A)^2 * eta(x^8 + A)^2 * eta(x^12 + A)^2), n))}; /* Michael Somos, Jan 20 2017 */
    
  • Sage
    Q = DiagonalQuadraticForm(ZZ,[3, 2])
    Q.representation_number_list(102) # Peter Luschny, Jun 20 2014

Formula

G.f.: 1 + Sum_{k>0} x^k * (1 + x^(4*k)) * (1 + x^(6*k)) / (1 + x^(12*k)) - Sum_{k>0} Kronecker( k, 3) * x^k * (1 - x^(2*k)) / (1 + x^(4*k)).
G.f.: Sum_{i, j in Z} x^(2*i^2 + 3*j^2). - Michael Somos, Mar 01 2011
Expansion of phi(q^2) * phi(q^3) in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Mar 01 2011
A115660(n) = A000377(n) - a(n). - Michael Somos, Mar 01 2011
Euler transform of period 24 sequence [0, 2, 2, -3, 0, -1, 0, -1, 2, 2, 0, -4, 0, 2, 2, -1, 0, -1, 0, -3, 2, 2, 0, -2, ...]. - Michael Somos, Jan 20 2017
Expansion of eta(q^4)^5 * eta(q^6)^5 / (eta(q^2)^2 * eta(q^3)^2 * eta(q^8)^2 * eta(q^12)^2) in powers of q. - Michael Somos, Jan 20 2017
a(0) = 1, a(n) = (1-(-1)^t)*b(n) for n > 0, where t is the number of prime factors of n, counting multiplicity, which are == 2,3,5,11 (mod 24), and b() is multiplicative with b(p^e) = (e+1) for primes p == 1,5,7,11 (mod 24) and b(p^e) = (1+(-1)^e)/2 for primes p == 13,17,19,23 (mod 24). (This formula is Corollary 4.2 in the Berkovich-Yesilyurt paper). - Jeremy Lovejoy, Nov 14 2024

Extensions

Edited by Charles R Greathouse IV, Oct 28 2009
Edited by N. J. A. Sloane, Mar 04 2011

A257921 Expansion of f(x^2, -x^4) * f(-x, -x^5)^2 / f(-x^12, -x^12) in powers of x where f(, ) is Ramanujan's general theta functions.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jul 12 2015

Keywords

Comments

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

Examples

			G.f. = 1 - 2*x + 2*x^2 - 2*x^3 + x^6 - 2*x^7 + 4*x^8 - 4*x^13 + 2*x^14 + ...
G.f. = q^3 - 2*q^7 + 2*q^11 - 2*q^15 + q^27 - 2*q^31 + 4*q^35 - 4*q^55 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, With[ {m = 4 n + 3}, DivisorSum[ m, KroneckerSymbol[ 12, #] KroneckerSymbol[ -2, m/#] &]]];
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 2, Pi/4, x^(1/2)] EllipticTheta[ 2, 0, x^(3/2)])^2 / (2^(5/2) x^(3/4) EllipticTheta[ 2, Pi/4, x] EllipticTheta[ 4, 0, x^12]), {x, 0, n};
  • PARI
    {a(n) = if( n<0, 0, n = 4*n + 3; sumdiv(n, d, kronecker( 12, d) * kronecker( -2, n/d)))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^4 + A)^3 * eta(x^6 + A)^4 * eta(x^24 + A) / (eta(x^2 + A)^3 * eta(x^3 + A)^2 *eta(x^8 + A) * eta(x^12 + A)^2), n))};

Formula

Expansion of psi(-x)^2 * psi(x^3)^2 / (psi(-x^2) * phi(-x^12)) in powers of x where phi(), psi() are Ramanujan theta functions.
Expansion of q^(-3/4) * eta(q)^2 * eta(q^4)^3 * eta(q^6)^4 * eta(q^24) / (eta(q^2)^3 * eta(q^3)^2 *eta(q^8) * eta(q^12)^2) in powers of q.
Euler transform of period 24 sequence [ -2, 1, 0, -2, -2, -1, -2, -1, 0, 1, -2, -2, -2, 1, 0, -1, -2, -1, -2, -2, 0, 1, -2, -2, ...].
a(n) = b(4*n + 3) where b() is multiplicative with b(2^e) = 0^e, b(3^e) = 1, b(p^e) = e+1 if p == 1, 11 (mod 24), b(p^e) = (e+1) * (-1)^e if p == 5, 7 (mod 24), b(p^e) = (1 + (-1)^e) / 2 if p == 13, 17, 19, 23 (mod 24).
a(n) = (-1)^n * A261119(n) = A134177(2*n + 1) = - A128580(2*n + 1) = - A115660(4*n+3).
a(6*n + 4) = a(6*n + 5) = 0.
a(2*n) = A257920(n). a(2*n + 1) = -2 * A259896(n). a(3*n) = A259668(n). a(6*n + 2) = 2 * A128591(n).
Showing 1-10 of 11 results. Next