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

A033715 Number of integer solutions (x, y) to the equation x^2 + 2y^2 = n.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Theta series of lattice C2 with Gram matrix [ 1, 0; 0, 2]. a(n) is nonzero if and only if n is in A002479. - Michael Somos, Dec 15 2011
Number 17 of the 74 eta-quotients listed in Table I of Martin (1996).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Denoted by |a_4(n)| in Kassel and Reutenauer 2015. - Michael Somos, Jun 16 2015

Examples

			G.f. = 1 + 2*q + 2*q^2 + 4*q^3 + 2*q^4 + 4*q^6 + 2*q^8 + 6*q^9 + 4*q^11 + 4*q^12 + ...
		

References

  • Bruce C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, 1991, see p. 114 Entry 8(iii).
  • J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, Springer-Verlag, 1999, p. 102, eq. 9.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 3, p. 19.
  • Nathan J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 78, Eq. (32.24).
  • J. W. L. Glaisher, Table of the excess of the number of (8k+1)- and (8k+3)-divisors of a number over the number of (8k+5)- and (8k+7)-divisors, Messenger Math., 31 (1901), 82-91.
  • J. V. Uspensky and M. A. Heaslet, Elementary Number Theory, McGraw-Hill, NY, 1939, p. 346.

Crossrefs

Number of integer solutions to f(x,y) = n where f(x,y) is the principal binary quadratic form with discriminant d: A004016 (d=-3), A004018 (d=-4), A002652 (d=-7), this sequence (d=-8), A028609 (d=-11), A028641 (d=-19), A138811 (d=-43).

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(8), 1), 105); A[1] + 2*A[2] + 2*A[3]; /* Michael Somos, Aug 29 2014 */
  • Maple
    d:=proc(r,m,n) local i,t1; t1:=0; for i from 1 to n do if n mod i = 0 and i-r mod m = 0 then t1:=t1+1; fi; od: t1; end; [seq(2*(d(1,8,n)+d(3,8,n)-d(5,8,n)-d(7,8,n)),n=1..120)];
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^2], {q, 0, n}]; (* Michael Somos, Sep 09 2012 *)
    a[ n_] := If[ n < 1, Boole[ n == 0], 2 DivisorSum[ n, KroneckerSymbol[ -2, #] &]]; (* Michael Somos, Aug 29 2014 *)
    a[ n_] := SeriesCoefficient[ (QPochhammer[ q^2] QPochhammer[ q^4])^3 / (QPochhammer[ q] QPochhammer[ q^8])^2, {q, 0, n}]; (* Michael Somos, Aug 29 2014 *)
  • PARI
    {a(n) = if( n<1, n==0, 2 * (issquare(n) - issquare(2*n) + 2 * sum( i=1, sqrtint(n\2), issquare(n - 2*i^2))))};
    
  • PARI
    {a(n) = if( n<1, n==0, 2 * sumdiv( n, d, kronecker( -2, d)))}; /* Michael Somos, Aug 23 2005 */
    
  • PARI
    {a(n) = if( n<1, n==0, 2 * qfrep([ 1, 0; 0, 2], n)[n])}; /* Michael Somos, Aug 23 2005 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^-2 * eta(x^2 + A)^3 * eta(x^4 + A)^3 * eta(x^8 + A)^-2, n))};
    
  • Sage
    Q = DiagonalQuadraticForm(ZZ, [1,2]); Q.representation_number_list(104); # Peter Luschny, Jun 20 2014
    

Formula

Fine gives an explicit formula for a(n) in terms of the divisors of n.
Euler transform of period 8 sequence [ 2, -1, 2, -4, 2, -1, 2, -2, ...].
Expansion of (eta(q^2) * eta(q^4))^3 / (eta(q) * eta(q^8))^2 in powers of q.
Coefficients in expansion of Sum_{i,j=-inf..inf} q^(i^2 + 2*j^2).
G.f. = s(2)^3*s(4)^3/(s(1)^2*s(8)^2), where s(k) := subs(q=q^k, eta(q)), where eta(q) is Dedekind's function, cf. A010815. [Fine]
G.f.: 1 + 2 * Sum_{k>0} Kronecker(-2, n) * x^k / (1 - x^k) = 1 + 2 * Sum_{k>0} (x^k + x^(3*k)) / (1 + x^(4*k)).
G.f.: theta_3(q) * theta_3(q^2) = Product_{k>0} (1 + x^(2*k)) * ((1 + x^k) * (1 - x^(2*k)) / (1 + x^(4*k)))^2.
From Michael Somos, Oct 23 2006: (Start)
Moebius transform is period 8 sequence [ 2, 0, 2, 0, -2, 0, -2, 0, ...].
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 - 3*u3) * (u1 - u2 - u3 + u6) - (u2 - 3*u6) * (u1 - 2*u2 - u3 + 2*u6). (End)
a(n) = 2 * A002325(n) unless n = 0.
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 8^(1/2) (t/i) f(t) where q = exp(2 Pi i t). - Michael Somos, Sep 09 2012
From Michael Somos, Aug 29 2014: (Start)
Expansion of phi(q) * phi(q^2) in powers of q where phi() is a Ramanujan theta function.
a(2*n) = a(n). a(2*n + 1) = 2 * A113411(n). (End)
From Michael Somos, May 17 2015: (Start)
a(n) = A028572(4*n) = A133692(2*n) = A139093(8*n) = A226225(8*n) = A226240(4*n) = A242609(4*n) = A245572(4*n) / 3 = (-1)^floor((n + 1)/2) * A082564(n).
a(8*n + 5) = a(8*n + 7) = 0. a(8*n + 1) = 2 * A112603(n). a(8*n + 3) = 4 * A033761(n). (End)
a(0) = 1, a(n) = 2 * b(n) for n > 0, where b() is multiplicative with b(2^e) = 1, b(p^e) = e + 1 if p == 1, 3 (mod 8), b(p^e) = (1 + (-1)^e)/2 if p == 5, 7 (mod 8). - Jianing Song, Sep 04 2018 [Corrected by Jeremy Lovejoy, Nov 12 2024]
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=0..m} a(k) = Pi/sqrt(2) = 2.221441... (A247719). - Amiram Eldar, Dec 16 2023

A121444 Expansion of f(x^3, x^9) * f(x, x^2) in powers of x where f(, ) is Ramanujan's general theta functions.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jul 30 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^4 + 2*x^5 + x^7 + x^8 + x^9 + 2*x^10 + x^11 + ...
G.f. = q^5 + q^17 + q^29 + q^41 + q^53 + 2*q^65 + q^89 + q^101 + q^113 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, Sum[ I^d, {d, Divisors[12 n + 5]}] / (2 I)]; (* Michael Somos, Jul 25 2015 *)
    a[ n_] := SeriesCoefficient[ 2 x^(3/8) QPochhammer[ x^6]^3 / (QPochhammer[ x, x^2] EllipticTheta[ 2, 0, x^(3/2)]), {x, 0, n}]; (* Michael Somos, Jan 31 2015 *)
    a[ n_] := Length @ FindInstance[ 24 n + 10 == (6 j + 3)^2 + (6 k + 1)^2 && j >= 0, {j, k}, Integers, 10^9]; (* Michael Somos, Jul 02 2015 *)
    a[ n_] := If[ n < 0, 0, DivisorSum[ 12 n + 5, KroneckerSymbol[ -4, #] &] / 2]; (* Michael Somos, Nov 11 2015 *)
    a[ n_] := If[ n < 0, 0, Sum[ Boole[ Mod[d, 4] == 1] - Boole[ Mod[d, 4] == 3], {d, Divisors[12 n + 5]}] / 2]; (* Michael Somos, Nov 11 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x] QPochhammer[ x^3] QPochhammer[ x^6], {x, 0, n}]; (* Michael Somos, Nov 11 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) * eta(x^3 + A) * eta(x^6 + A) / eta(x + A), n))};
    
  • PARI
    {a(n) = if( n<0, 0, n = 12*n + 5; sumdiv(n, d, (d%4==1) - (d%4==3)) / 2)};

Formula

Expansion of f(-x^3) * f(-x^6) / chi(-x) in powers of x where chi(), f() are Ramanujan theta functions.
Expansion of q^(-5/12) * eta(q^2) * eta(q^3) * eta(q^6) / eta(q) in powers of q.
Euler transform of period 6 sequence [ 1, 0, 0, 0, 1, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (72 t)) = (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A258210.
G.f.: Product_{k>0} (1 + x^k) * (1 - x^(3*k)) * (1 - x^(6*k)).
-2 * a(n) = A121363(3*n + 1).
Convolution square is A098098.
a(n) = (-1)^n * A258832(n) = A052343(3*n + 1). -a(n) = A258291(3*n + 1). 2 * a(n) = A008441(3*n + 1). - Michael Somos, Jul 02 2015
From Peter Bala, Jan 07 2021: (Start)
G.f. A(x) = Sum_{n = -oo..oo} x^n/(1 - x^(12*n + 5)). See Agarwal, p. 285, equation 6.19.
A(x^2) = Sum_{n = -oo..oo} x^(2*n)/(1 - x^(12*n + 5)). Cf. A033761. (End)

A045828 One fourth of theta series of cubic lattice with respect to face.

Original entry on oeis.org

1, 2, 2, 4, 3, 2, 6, 4, 4, 6, 4, 4, 7, 8, 2, 8, 8, 4, 10, 4, 4, 10, 10, 8, 9, 4, 6, 12, 8, 6, 10, 12, 4, 14, 8, 4, 16, 10, 8, 8, 9, 10, 12, 12, 8, 12, 12, 4, 20, 10, 6, 20, 8, 6, 10, 12, 8, 20, 18, 8, 11, 12, 12, 16, 8, 6, 20, 16, 12, 14, 8, 12, 20, 14, 6, 12, 20, 8, 26, 12, 8, 22, 8, 12, 15
Offset: 0

Views

Author

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Number of solutions to n = t1 + t2 + 2*t3 where t1, t2, t3 are triangular numbers. - Michael Somos, Jan 02 2006
The cubic lattice is the set of triples [a, b, c] where the entries are all integers. A face is centered at a triple where one entry is an integer and the other two are one half an odd integer. - Michael Somos, Jun 29 2012

Examples

			G.f. = 1 + 2*x + 2*x^2 + 4*x^3 + 3*x^4 + 2*x^5 + 6*x^6 + 4*x^7 + 4*x^8 + 6*x^9 + ...
G.f. = q + 2*q^3 + 2*q^5 + 4*q^7 + 3*q^9 + 2*q^11 + 6*q^13 + 4*q^15 + 4*q^17 + ...
		

References

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

Crossrefs

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ 1/4 EllipticTheta[ 3, 0, x] EllipticTheta[ 2, 0, x]^2, {x, 0, n + 1/2}]; (* Michael Somos, Jun 29 2012 *)
    a[ n_] := SeriesCoefficient[ 1/8 EllipticTheta[ 2, 0, x^2] EllipticTheta[ 2, 0, x]^2, {x, 0, 2 n + 1}]; (* Michael Somos, Jun 29 2012 *)
    QP = QPochhammer; s = (QP[q^2]^3*QP[q^4]^2)/QP[q]^2 + O[q]^90; CoefficientList[s, q] (* Jean-François Alcover, Nov 27 2015, adapted from PARI *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^3 * eta(x^4 + A)^2 / eta(x + A)^2, n))}; /* Michael Somos, Oct 25 2006 */

Formula

Expansion of q^(-1/2) * (eta(q^2)^3 * eta(q^4)^2) / eta(q)^2 in powers of q. - Michael Somos, Jan 02 2006
Expansion of phi(x) * psi(x^2)^2 = psi(x)^2 * psi(x^2) = psi(x)^4 / phi(x) in powers of x where phi(), psi() are Ramanujan theta functions. - Michael Somos, Jun 29 2012
Euler transform of period 4 sequence [2, -1, 2, -3, ...]. - Michael Somos, Mar 05 2003
Convolution of A033761 and A010054. - Michael Somos, Jun 29 2012
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = (1/2)^(1/2) (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A212885. - Michael Somos, Sep 08 2018

Extensions

Edited by Michael Somos, Mar 05 2003

A113411 Excess of number of divisors of 2n+1 of form 8k+1, 8k+3 over those of form 8k+5, 8k+7.

Original entry on oeis.org

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

Views

Author

Michael Somos, Oct 29 2005

Keywords

Comments

Ramanujan theta functions: f(q) := Prod_{k>=1} (1-(-q)^k) (see A121373), phi(q) := theta_3(q) := Sum_{k=-oo..oo} q^(k^2) (A000122), psi(q) := Sum_{k=0..oo} q^(k*(k+1)/2) (A010054), chi(q) := Prod_{k>=0} (1+q^(2k+1)) (A000700).
Bisection of A002325. Number of ways to write n as a sum of a square plus four times a triangular number [Hirschhorn]. - R. J. Mathar, Mar 23 2011

Examples

			1 + 2*x + 3*x^4 + 2*x^5 + 2*x^8 + 2*x^9 + x^12 + 4*x^13 + 4*x^16 + ...
q + 2*q^3 + 3*q^9 + 2*q^11 + 2*q^17 + 2*q^19 + q^25 + 4*q^27 + 4*q^33 + ...
		

References

  • Nathan J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 82, Eq. (32.55).

Crossrefs

Programs

  • Mathematica
    a[n_] := DivisorSum[2n+1, Switch[Mod[#, 8], 1|3, 1, 5|7, -1]&]; Table[a[n], {n, 0, 104}] (* Jean-François Alcover, Dec 04 2015 *)
  • PARI
    a(n) = if( n<0, 0, n = 2*n + 1; sumdiv(n, d, (-1)^(d%8>3)))
    
  • PARI
    a(n) = local(n1); if( n<0, 0, n1 = sqrtint(n); polcoeff( sum(k=1,n1, 2*x^k^2, 1 + x*O(x^n)) * sum(k=0,n1, x^(2*k^2 + 2*k)), n))
    
  • PARI
    a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^5 * eta(x^8 + A)^2 / (eta(x + A)^2 * eta(x^4 + A)^3), n))
    
  • PARI
    a(n) = local(A, p, e); if( n<0, 0, n = 2*n + 1; A = factor(n); prod(k=1, matsize(A)[1], if( p=A[k,1], e=A[k,2]; if( p==2, 0, if( abs(p%8-6)==1, (1+(-1)^e)/2, e+1)))))

Formula

Expansion of phi(q) * psi(q^4) in powers of q where psi(), phi() are Ramanujan theta functions.
Expansion of q^(-1) * (eta(q^4)^5 * eta(q^16)^2) / (eta(q^2)^2 * eta(q^8)^3) in powers of q^2.
a(n) = b(2*n + 1) where b(n) is multiplicative and b(2^e) = 0^e, b(p^e) = e+1 if p == 1, 3 (mod 8), b(p^e) = (1+(-1)^e)/2 if p == 5, 7 (mod 8).
Euler transform of period 8 sequence [ 2, -3, 2, 0, 2, -3, 2, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 2^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is g.f. for A133692. - Michael Somos, Mar 16 2011
G.f.: (Sum_{k} x^k^2) * (Sum_{k>=0} x^(2*k^2 + 2*k)).
G.f.: Sum_{k>=0} a(k) * x^(2*k + 1) = Sum_{k>=0} F(x^(2*k + 1), x^(3*(2*k + 1))) where F(x, y) = (x + y) / (1 + x*y).
a(4*n + 2) = a(4*n + 3) = 0. a(4*n) = A112603(n). a(4*n + 1) = 2 * A033761(n).
From Peter Bala, Jan 07 2021: (Start)
Conjectural g.f.s: A(x) = Sum_{n >= 0} (-1)^(n*(n-1)/2)*x^n/(1 - x^(2*n+1)).
A(x) = Sum_{n = -oo..oo} (-1)^n*x^(2*n)/(1 - x^(4*n+1)) = Sum_{n = -oo..oo} (-1)^n*x^(2*n+1)/(1 - x^(4*n+3)). (End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(2*sqrt(2)) = 1.1107207... (A093954). - Amiram Eldar, Dec 28 2023

A082564 Expansion of eta(q)^2 * eta(q^2) / eta(q^4) in powers of q.

Original entry on oeis.org

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

Views

Author

Benoit Cloitre, May 05 2003

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
a(n) is nonzero if and only if n is in A002479. - Michael Somos, Dec 15 2011
Absolute values appear to give A033715 = 2*A002325.
Denoted by a_4(n) in Kassel and Reutenauer 2015. - Michael Somos, Jun 04 2015

Examples

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

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(32), 1), 105); A[1] - 2*A[2] - 2*A[3] + 4*A[4] + 2*A[5] - 4*A[7] + 2*A[9] - 6*A[10] + 4*A[12] + 4*A[13] - 4*A[16]; /* Michael Somos, Aug 29 2014 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ q]^2 QPochhammer[ q^2] / QPochhammer[ q^4], {q, 0, n}]; (* Michael Somos, Aug 29 2014 *)
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 4, 0, q] EllipticTheta[ 4, 0, q^2], {q, 0, n}]; (* Michael Somos, Aug 29 2014 *)
    a[ n_] := If[ n < 1, Boole[ n == 0], 2 (-1)^Quotient[ n + 1, 2] DivisorSum[ n, KroneckerSymbol[ -2, #] &]]; (* Michael Somos, Aug 29 2014 *)
  • PARI
    {a(n) = if( n<1, n==0, 2 * (-1)^((n+1) \ 2) * sumdiv( n, d, kronecker( -2, d)))}; /* Michael Somos, Mar 30 2007 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^2 + A) / eta(x^4 + A), n))};
    

Formula

Expansion of phi(-q) * phi(-q^2) in powers of q where phi() is a Ramanujan theta function. - Michael Somos, Mar 30 2007
Euler transform of period 4 sequence [ -2, -3, -2, -2, ...]. - Michael Somos, Mar 30 2007
G.f. is a period 1 Fourier series which satisfies f(-1 / (32 t)) = 2^(11/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A033761. - Michael Somos, Aug 29 2014
G.f.: Product_{k>0} (1 - x^k)^2 / (1 + x^(2*k)). - Michael Somos, Mar 30 2007
a(n) = -2 * A129134(n) unless n=0. - Michael Somos, Mar 30 2007
a(n) = (-1)^floor( (n+1)/2 ) * A033715(n). - Michael Somos, Aug 29 2014
a(2*n) = A133692(n). a(2*n + 1) = -2 * A125095(n). - Michael Somos, Aug 29 2014
a(3*n + 1) = -2 * A258747(n). a(3*n + 2) = -2 * A258764(n). - Michael Somos, Jun 09 2015

A139093 Expansion of phi(q) * phi(-q^2) in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Apr 08 2008

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
a(n) is nonzero if and only if n is in A002479.

Examples

			G.f. = 1 + 2*q - 2*q^2 - 4*q^3 + 2*q^4 - 4*q^6 + 2*q^8 + 6*q^9 - 4*q^11 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(32), 1), 105); A[1] + 2*A[2] - 2*A[3] - 4*A[4] + 2*A[5] - 4*A[7] + 2*A[9] + 6*A[10] - 4*A[12] + 4*A[13] + 4*A[16]; /* Michael Somos, Aug 29 2014 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 4, 0, q^2], {q, 0, n}]; (* Michael Somos, Aug 29 2014 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ q^2]^7 / (QPochhammer[ q]^2 QPochhammer[ q^4]^3), {q, 0, n}]; (* Michael Somos, Feb 18 2015 *)
    a[ n_] := If[ n < 1, Boole[n == 0], 2 (-1)^Quotient[n, 2] Sum[ JacobiSymbol[ -2, d], {d, Divisors @ n}]]; (* Michael Somos, Feb 18 2015 *)
  • PARI
    {a(n) = if( n<1, n==0, 2 * (-1)^(n\2) * sumdiv(n, d, kronecker( -2, d)))};
    
  • PARI
    {a(n) = my(A); if ( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^7 / (eta(x + A)^2 * eta(x^4 + A)^3), n))};
    

Formula

Expansion of eta(q^2)^7 / (eta(q)^2 * eta(q^4)^3) in powers of q.
Euler transform of period 4 sequence [ 2, -5, 2, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (4 t)) = 8^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A112603.
G.f.: Product_{k>0} (1 - x^(2*k))^2 * (1 + x^(2*k-1))^2 / (1 + x^(2*k)).
a(8*n + 5) = a(8*n + 7) = 0.
a(n) = (-1)^n * A082564(n). a(2*n) = A133692(n). a(2*n + 1) = 2 * A125095(n). a(4*n) = a(8*n) = A033715(n). a(8*n + 1) = 2 * A112603(n). a(8*n + 3) = -4 * A033761(n).

A227395 Expansion of q^2 * phi(-q) * psi(q^16) in powers of q where phi(), psi() are Ramanujan theta functions.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jul 10 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 + 2*q^6 - 2*q^11 + 3*q^18 - 2*q^19 + 2*q^22 - 4*q^27 + 2*q^34 + ...
		

Crossrefs

Programs

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

Formula

Expansion of eta(q)^2 * eta(q^32)^2 / (eta(q^2) * eta(q^16)) in powers of q.
Euler transform of period 32 sequence [ -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, 0, -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, -1, -2, -2, ...].
G.f. is a period 1 Fourier series which satisfies f(-1 / (32 t)) = 32^(1/2) (t/i) f(t) where q = exp(2 Pi i t).
a(4*n) = a(4*n + 1) = a(8*n + 7) = 0. a(4*n + 2) = A113411(n). a(8*n + 3) = -2 * A033761(n).
G.f.: x^2 * Product_{k>0} (1 - x^k)^2 * (1 - x^(32*k))^2 / ((1 - x^(2*k)) * (1 - x^(16*k))).
a(n) = (-1)^n * A255258(n). - Michael Somos, Feb 20 2015

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

Original entry on oeis.org

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

Views

Author

Michael Somos, Nov 20 2006

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Number 45 of the 74 eta-quotients listed in Table I of Martin (1996). - Michael Somos, Mar 14 2012

Examples

			G.f. = 1 - 2*x + 3*x^4 - 2*x^5 + 2*x^8 - 2*x^9 + x^12 - 4*x^13 + 4*x^16 + ...
G.f. = q - 2*q^3 + 3*q^9 - 2*q^11 + 2*q^17 - 2*q^19 + q^25 - 4*q^27 + 4*q^33 + ...
		

Crossrefs

Programs

  • Mathematica
    a[ n_] := If[ n < 0, 0, (-1)^n DivisorSum[ 2 n + 1, If[ Mod[#, 8] > 3, -1, 1] &]]; (* Michael Somos, Jul 09 2015 *)
    a[ n_] := SeriesCoefficient[ (QPochhammer[ x] QPochhammer[ x^8])^2 / (QPochhammer[ x^2] QPochhammer[ x^4]), {x, 0, n}]; (* Michael Somos, Jul 09 2015 *)
  • PARI
    {a(n) = if( n<0, 0, (-1)^n * sumdiv( 2*n + 1, d, (-1)^(d%8>3)))};
    
  • PARI
    {a(n) = if( n<0, 0, n = 2*n + 1; qfrep( [1, 0; 0, 8], n)[n] - qfrep( [3, 1; 1, 3], n)[n])};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^8 + A)^2 / (eta(x^2 + A) * eta(x^4 + A)), n))}
    
  • PARI
    {a(n) = if( n<0, 0, n = 2*n + 1; sumdiv( n, d, kronecker(2, d) * kronecker( -4, n/d)))};

Formula

Expansion of q^(-1/2) * (eta(q)^2 * eta(q^8)^2) / (eta(q^2) * eta(q^4)) in powers of q.
Given g.f. A(x), then B(q) = q * A(q^2) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u1^2*u6 * (u1 + 3*u3) + 2 * u2^2*u3 * (u2 + 3*u6) - 3 * u3^2*u2 * (u1 + u3) - 6 * u6^2*u1 * (u2 + u6).
a(n) = b(2*n + 1) where b() is multiplicative with b(2^e) = 0^e, b(p^e) = (e+1) * (-1)^e if p == 1, 3 (mod 8), b(p^e) = (1 + (-1)^e) / 2 if p == 5, 7 (mod 8).
Euler transform of period 8 sequence [ -2, -1, -2, 0, -2, -1, -2, -2, ...].
G.f.: (Sum_{k in Z} (-1)^k * x^k^2) * (Sum_{k>=0} x^(2*k^2 + 2*k)).
a(4*n + 2) = a(4*n + 3) = 0. a(n) = (-1)^n * A113411(n). a(4*n) = A112603(n). a(4*n + 1) = -2 * A033761(n).

A244540 Expansion of phi(q) * (phi(q) + phi(q^2)) / 2 in powers of q where phi() is a Ramanujan theta function.

Original entry on oeis.org

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

Views

Author

Michael Somos, Jun 29 2014

Keywords

Comments

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

Examples

			G.f. = 1 + 3*q + 3*q^2 + 2*q^3 + 3*q^4 + 4*q^5 + 2*q^6 + 3*q^8 + 5*q^9 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(8), 1), 33); A[1] + 3*A[2] + 3*A[3];
  • Mathematica
    a[ n_] := If[ n < 1, Boole[n == 0], Sum[ {3, 0, -1, 0, 1, 0, -3, 0}[[ Mod[ d, 8, 1] ]], {d, Divisors @ n}]];
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] (EllipticTheta[ 3, 0, q] + EllipticTheta[ 3, 0, q^2]) / 2, {q, 0, n}];
  • PARI
    {a(n) = if( n<1, n==0, sumdiv(n, d, [0, 3, 0, -1, 0, 1, 0, -3][d%8 + 1]))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = sum(k=1, sqrtint(n), 2 * x^k^2, 1 + x * O(x^n)); polcoeff( A * (A + subst(A, x, x^2)) / 2, n))};
    
  • Sage
    A = ModularForms( Gamma1(8), 1, prec=33) . basis(); A[0] + 3*A[1] + 3*A[2];
    

Formula

Expansion of f(-q^3, -q^5)^2 * phi(q) / psi(-q) = f(-q^3, -q^5)^2 * chi(q)^3 in powers of q where phi(), psi(), chi(), f() are Ramanujan theta functions.
Euler transform of period 8 sequence [3, -3, 1, 0, 1, -3, 3, -2, ...].
Moebius transform is period 8 sequence [3, 0, -1, 0, 1, 0, -3, 0, ...].
Convolution product of A244526 and A107635. Convolution product of A000122 and A093709.
a(n) = (A004018(n) + A033715(n)) / 2 = A244543(2*n).
a(2*n) = a(n). a(8*n + 3) = 2*A033761(n). a(8*n + 5) = 4*A053692(n). a(8*n + 7) = 0.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=0..m} a(k) = Pi*(1 + 1/sqrt(2))/2 = 2.681517... . - Amiram Eldar, Jun 08 2025

A125096 Expansion of -1 + (phi(q) * phi(q^2) + phi(-q^2) * phi(q^4)) / 2 in powers of q.

Original entry on oeis.org

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

Views

Author

Michael Somos, Nov 20 2006

Keywords

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[MemberQ[{1, 3}, Mod[p, 8]], e + 1, (1 + (-1)^e)/2]; f[2, e_] := If[e > 1, 2, 0]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Oct 13 2022 *)
  • PARI
    {a(n) = if( n<1, 0, qfrep([1, 0; 0, 8], n)[n] + qfrep([3, 1; 1, 3], n)[n])}

Formula

a(n) is multiplicative with a(2) = 0, a(2^e) = 2 if e>1, a(p^e) = e+1 if p == 1, 3 (mod 8), a(p^e) = (1+(-1)^e)/2 if p == 5, 7 (mod 8).
a(4*n + 2) = a(8*n + 5) = a(8*n + 7) = 0. a(4*n) = 2 * A002325(n). a(8*n + 1) = A112603(n). a(8*n + 3) = A033761(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/(2*sqrt(2)) = 1.110720... (A093954). - Amiram Eldar, Oct 13 2022
Showing 1-10 of 21 results. Next