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

A033687 Theta series of hexagonal lattice A_2 with respect to deep hole divided by 3.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
a(n)=0 if and only if A000731(n)=0 (see the Han-Ono paper). - Emeric Deutsch, May 16 2008
Number of 3-core partitions of n (denoted c_3(n) in Granville and Ono, p. 340). - Brian Hopkins, May 13 2008
Denoted by g_1(q) in Cynk and Hulek in Remark 3.4 on page 12 (but not explicitly listed).
This is a member of an infinite family of integer weight modular forms. g_1 = A033687, g_2 = A030206, g_3 = A130539, g_4 = A000731. - Michael Somos, Aug 24 2012
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 + x + 2*x^2 + 2*x^4 + x^5 + 2*x^6 + x^8 + 2*x^9 + 2*x^10 + 2*x^12 + 2*x^14 + ...
G.f. = q + q^4 + 2*q^7 + 2*q^13 + q^16 + 2*q^19 + q^25 + 2*q^28 + 2*q^31 + 2*q^37 + ...
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 111.
  • N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 79, Eq. (32.35) and (32.351).

Crossrefs

Programs

  • Magma
    Basis( ModularForms( Gamma1(9), 1), 316) [2]; /* Michael Somos, May 06 2015 */
  • Mathematica
    a[ n_] := If[ n < 0, 0, DivisorSum[ 3 n + 1, KroneckerSymbol[ -3, #] &]]; (* Michael Somos, Sep 23 2013 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x^3]^3 / QPochhammer[ x], {x, 0, n}]; (* Michael Somos, Sep 01 2015 *)
  • PARI
    {a(n) = if( n<0, 0, sumdiv( 3*n + 1, d, kronecker( -3, d)))}; /* Michael Somos, Nov 03 2005 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^3 + A)^3 / eta(x + A), n))};
    
  • PARI
    {a(n) = my(A, p, e); if( n<0, 0, A = factor( 3*n + 1); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, 0, p%6==1, e+1, 1-e%2)))}; /* Michael Somos, May 06 2015 */
    

Formula

Euler transform of period 3 sequence [1, 1, -2, ...].
Expansion of q^(-1/3) * eta(q^3)^3 / eta(q) in powers of q.
a(4*n + 1) = a(n). - Michael Somos, Dec 06 2004
a(n) = b(3*n + 1) where b(n) is multiplicative and b(p^e) = 0^e if p = 3, b(p^e) = e+1 if p == 1 (mod 6), b(p^e) = (1+(-1)^e)/2 if p == 2, 5 (mod 6). - Michael Somos, May 20 2005
Given g.f. A(x), B(q) = q * A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = u^2*w - 2*u*w^2 - v^3. - Michael Somos, Dec 06 2004
Given g.f. A(x), B(q)= q * A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u1*u3^2 + u1*u6^2 - u1*u3*u6 - u2^2*u3. - Michael Somos, May 20 2005
Given g.f. A(x), B(q) = q * A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u2*u3^2 + 2*u2*u3*u6 + 4*u2*u6^2 - u1^2*u6. - Michael Somos, May 20 2005
G.f.: Product_{k>0} (1 - x^(3*k))^3 / (1 - x^k).
G.f.: Sum_{k in Z} x^k / (1 - x^(3*k + 1)) = Sum_{k in Z} x^k / (1 - x^(6*k + 2)). - Michael Somos, Nov 03 2005
Expansion of q^(-1) * c(q^3) / 3 = q^(-1) * (a(q) - b(q)) / 9 in powers of q^3 where a(), b(), c() are cubic AGM theta functions. - Michael Somos, Dec 25 2007
G.f. is a period 1 Fourier series which satisfies f(-1 / (27 t)) = 3^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g(t) is the g.f. for A005928.
a(n) = Sum_{d|3n+1} LegendreSymbol{d,3} - Brian Hopkins, May 13 2008
q-series for a(n): Sum_{n >= 0} q^(n^2+n)(1-q)(1-q^2)...(1-q^n)/((1-q^(n+1))(1-q^(n+2))...(1-q^(2n+1))). [From Jeremy Lovejoy, Jun 12 2009]
a(n) = A002324(3*n + 1). 3*a(n) = A005882(n) = A033685(3*n + 1). - Michael Somos, Apr 04 2003
G.f.: (2 * psi(x^2) * f(x^2, x^4) + phi(x) * f(x^1, x^5)) / 3 where phi(), psi() are Ramanujan theta functions and f(, ) is Ramanujan's general theta function. - Michael Somos, Sep 07 2018
Sum_{k=1..n} a(k) ~ 2*Pi*n/3^(3/2). - Vaclav Kotesovec, Dec 17 2022

A000727 Expansion of Product_{k >= 1} (1 - x^k)^4.

Original entry on oeis.org

1, -4, 2, 8, -5, -4, -10, 8, 9, 0, 14, -16, -10, -4, 0, -8, 14, 20, 2, 0, -11, 20, -32, -16, 0, -4, 14, 8, -9, 20, 26, 0, 2, -28, 0, -16, 16, -28, -22, 0, 14, 16, 0, 40, 0, -28, 26, 32, -17, 0, -32, -16, -22, 0, -10, 32, -34, -8, 14, 0, 45, -4, 38, 8, 0, 0, -34, -8, 38, 0, -22, -56, 2, -28, 0, 0, -10, 20, 64, -40, -20, 44
Offset: 0

Views

Author

Keywords

Comments

Number 51 of the 74 eta-quotients listed in Table I of Martin (1996).
Ramanujan (see the link, pp. 155 and 157 Nr. 23.) conjectured the expansion coefficients called Psi_4(n) of eta^4(6*z) in powers of q = exp(2*Pi*i*z), Im(z) > 0, where i is the imaginary unit. In the Finch link on p. 5, multiplicity is used and Psi_4(p^r), called f(p^r), is given (see also b(p^e) formula given by Michael Somos, Aug 23 2006). Mordell proved this conjecture on pp. 121-122 based on Klein-Fricke, Theorie der elliptischen Modulfunktionen, 1892, Band II, p. 374. The product formula for the Dirichlet series, Mordell, eq. (7) for a=2,is used to find Psi_4(n), called f_2(n), from f_2(p) for primes p. The primes p = 2 and 3 do not appear in the product. - Wolfdieter Lang, May 03 2016

Examples

			G.f. = 1 - 4*x + 2*x^2 + 8*x^3 - 5*x^4 - 4*x^5 - 10*x^6 + 8*x^7 + 9*x^8 + ...
G.f. = q - 4*q^7 + 2*q^13 + 8*q^19 - 5*q^25 - 4*q^31 - 10*q^37 + 8*q^43 + ...
		

References

  • Morris Newman, A table of the coefficients of the powers of eta(tau). Nederl. Akad. Wetensch. Proc. Ser. A. 59 = Indag. Math. 18 (1956), 204-216.
  • J. H. Silverman, A Friendly Introduction to Number Theory, 3rd ed., Pearson Education, Inc, 2006, p. 415. Exer. 47.2.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Julia
    # DedekindEta is defined in A000594.
    L000727List(len) = DedekindEta(len, 4)
    L000727List(82) |> println # Peter Luschny, Mar 09 2018
    
  • Magma
    qEigenform( EllipticCurve( [0, 0, 0, 0, 1]), 493); /* Michael Somos, Jun 12 2014 */
    
  • Magma
    A := Basis( ModularForms( Gamma0(36), 2), 493); A[2] - 4*A[8]; /* Michael Somos, Jun 12 2014 */
    
  • Magma
    Basis( CuspForms( Gamma0(36), 2), 493)[1]; /* Michael Somos, May 17 2015 */
    
  • Magma
    Coefficients(&*[(1-x^m)^4:m in [1..100]])[1..100] where x is PolynomialRing(Integers()).1; // Vincenzo Librandi, Mar 10 2018
  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d,j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:= etr(n-> -4): seq(a(n), n=0..81); # Alois P. Heinz, Sep 08 2008
  • Mathematica
    etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n - j], {j, 1, n}]/n]; b]; a = etr[-4&]; Table[a[n], {n, 0, 81}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x]^4, {x, 0, n}]; (* Michael Somos, Jun 12 2014 *)
    nmax = 80; CoefficientList[Series[Sum[Sum[(-1)^(k+m) * (2*k+1) * q^(k*(k+1)/2 + m*(3*m-1)/2), {k, 0, nmax}], {m, -nmax, nmax}], {q, 0, nmax}], q] (* Vaclav Kotesovec, Dec 06 2015 *)
  • PARI
    {a(n) = my(A, p, e, x, y, a0, a1); if( n<0, 0, n = 6*n + 1; A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p<5, 0, p%6==5, if(e%2, 0, (-1)^(e/2) * p^(e/2)), for( y=1, sqrtint(p\3), if( issquare( p - 3*y^2, &x), break)); a0=1; if( x%3!=1, x=-x); a1 = x = 2*x; for( i=2, e, y = x*a1 - p*a0; a0=a1; a1=y); a1)))}; /* Michael Somos, Aug 23 2006 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff(eta(x + x * O(x^n))^4, n))};
    
  • PARI
    {a(n) = if( n<0, 0, ellak( ellinit( [0, 0, 0, 0, 1], 1), 6*n + 1))}; /* Michael Somos, Jul 01 2004 */
    
  • Sage
    ModularForms( Gamma0(36), 2, prec=493).0; # Michael Somos, Jun 12 2014
    

Formula

Euler transform of period 1 sequence [-4, -4, ...]. - Michael Somos, Apr 02 2005
Given g.f. A(x), then B(q) = q * A(q^3)^2 satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = w*u^2 - v^3 + 16 * u*w^2. - Michael Somos, Apr 02 2005
a(n) = b(6*n + 1) where b() is multiplicative with b(2^e) = b(3^e) = 0^e, b(p^e) = b(p) * b(p^(e-1)) - p * b(p^(e-2)), b(p) = 0 if p == 5 (mod 6), b(p) = 2*x where p = x^2 + 3*y^2 == 1 (mod 6) and x == 1 (mod 3). - Michael Somos, Aug 23 2006
Coefficients of L-series for elliptic curve "36a1": y^2 = x^3 + 1. - Michael Somos, Jul 01 2004
a(n) = (-1)^n * A187076(n). a(2*n + 1) = -4 * A187150(n). a(25*n + 9) = a(25*n + 14) = a(25*n + 19) = a(25*n + 24) = 0. a(25*n + 4) = -5 * a(n). Convolution inverse of A023003. Convolution square of A002107. Convolution square is A000731.
a(0) = 1, a(n) = -(4/n)*Sum_{k=1..n} A000203(k)*a(n-k) for n > 0. - Seiichi Manyama, Mar 26 2017
G.f.: exp(-4*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 05 2018
Let M = p_1*...*p_k be a positive integer whose prime factors p_i (not necessarily distinct) are all congruent to 5 (mod 6). Then a( M^2*n + (M^2 - 1)/6 ) = (-1)^k*M*a(n). See Cooper et al., equation 4. - Peter Bala, Dec 01 2020
a(n) = b(6*n + 1) where b() is multiplicative with b(3^e) = 0^e, b(p^e) = (1 + (-1)^e)/2 * (-p)^(e/2) if p == 2 (mod 3), b(p^e) = (((x+sqrt(-3)*y)/2)^(e+1) - ((x-sqrt(-3)*y)/2)^(e+1))/x if p == 1 (mod 3) where p = x^2 + 3*y^2 and x == 1 (mod 3). - Jianing Song, Mar 19 2022

A002107 Expansion of Product_{k>=1} (1 - x^k)^2.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Number of partitions of n into an even number of distinct parts minus number of partitions of n into an odd number of distinct parts, with 2 types of each part. E.g., for n=4, we consider k and k* to be different versions of k and so we have 4, 4*, 31, 31*, 3*1, 3*1*, 22*, 211*, 2*11*. The even partitions number 5 and the odd partitions number 4, so a(4)=5-4=1. - Jon Perry, Apr 04 2004
Also, number of partitions of n into parts of -2 different kinds (based upon formal analogy). - Michele Dondi (blazar(AT)lcm.mi.infn.it), Jun 29 2004
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Number 68 of the 74 eta-quotients listed in Table I of Martin (1996).

Examples

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

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000712 (reciprocal of g.f.), A010815, A010816, A258406.
Powers of Euler's product: A000594, A000727 - A000731, A000735, A000739, A010815 - A010840.

Programs

  • Julia
    # DedekindEta is defined in A000594.
    A002107List(len) = DedekindEta(len, 2)
    A002107List(78) |> println # Peter Luschny, Mar 09 2018
  • Magma
    Basis( CuspForms( Gamma1(144), 1), 926) [1]; /* Michael Somos, May 17 2015 */
    
  • Maple
    A010816 := proc (n); if frac(sqrt(8*n+1)) = 0 then (-1)^((1/2)*isqrt(8*n+1)-1/2)*isqrt(8*n+1) else 0 end if; end proc:
    N := 10:
    a := proc (n) option remember; if n < 0 then 0 else A010816(n) + add( (-1)^(k+1)*a(n - (1/2)*k*(3*k-1) ), k = -N..-1) + add( (-1)^(k+1)*a(n - (1/2)*k*(3*k-1) ), k = 1..N) end if; end proc:
    seq(a(n), n = 0..100); # Peter Bala, Apr 06 2022
  • Mathematica
    terms = 78; Clear[s]; s[n_] := s[n] = Product[(1 - x^k)^2, {k, 1, n}] // Expand // CoefficientList[#, x]& // Take[#, terms]&; s[n = 10]; s[n = 2*n]; While[s[n] != s[n - 1], n = 2*n]; A002107 = s[n] (* Jean-François Alcover, Jan 17 2013 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x]^2, {x, 0, n}]; (* Michael Somos, Jan 31 2015 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 - x^k, {k, n}]^2, {x, 0, n}]; (* Michael Somos, Jan 31 2015 *)
  • PARI
    {a(n) = my(A, p, e, x); if( n<0, 0, n = 12*n + 1; A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p<5, 0, p%12>1, if( e%2, 0, (-1)^((p%12==5) * e/2)), for( i=1, sqrtint(p\9), if( issquare(p - 9*i^2), x=i; break)); (e + 1) * (-1)^(e*x))))}; /* Michael Somos, Aug 30 2006 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( eta(x + x * O(x^n))^2, n))}; /* Michael Somos, Aug 30 2006 */
    
  • PARI
    Vec(eta(x)^2) \\ Charles R Greathouse IV, Apr 22 2016
    

Formula

Expansion of q^(-1/12) * eta(q)^2 in powers of q. - Michael Somos, Mar 06 2012
Euler transform of period 1 sequence [ -2, ...]. - Michael Somos, Mar 06 2012
a(n) = b(12*n + 1) where b(n) is multiplicative and b(2^e) = b(3^e) = 0^e, b(p^e) = (1 + (-1)^e) / 2 if p == 7, 11 (mod 12), b(p^e) = (-1)^(e/2) * (1 + (-1)^e) / 2 if p == 5 (mod 12), b(p^e) = (e + 1) * (-1)^(e*x) if p == 1 (mod 12) where p = x^2 + 9*y^2. - Michael Somos, Sep 16 2006
Convolution inverse of A000712.
a(0) = 1, a(n) = -(2/n)*Sum{k = 0..n-1} a(k)*sigma_1(n-k). - Joerg Arndt, Feb 05 2011
Expansion of f(-x)^2 in powers of x where f() is a Ramanujan theta function. - Michael Somos, May 17 2015
G.f. is a period 1 Fourier series which satisfies f(-1 / (144 t)) = 12 (t/i) f(t) where q = exp(2 Pi i t). - Michael Somos, May 17 2015
a(n) = Sum_{k=0..n} A010815(k)*A010815(n-k); self convolution of A010815. - Gevorg Hmayakyan, Sep 18 2016
G.f.: Sum_{m, n in Z, n >= 2*|m|} (-1)^n * x^((3*(2*n + 1)^2 - (6*m + 1)^2)/24). - Seiichi Manyama, Oct 01 2016
G.f.: exp(-2*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 05 2018
From Peter Bala, Jan 02 2021: (Start)
For prime p congruent to 5, 7 or 11 (mod 12), a(n*p^2 + (p^2 - 1)/12) = e*a(n), where e = 1 if p == 7 or 11 (mod 12) and e = -1 if p == 5 (mod 12).
If n and p are coprime then a(n*p + (p^2 - 1)/12) = 0. See Cooper et al., Theorem 1. (End)
With the convention that a(n) = 0 for n < 0 we have the recurrence a(n) = A010816(n) + Sum_{k a nonzero integer} (-1)^(k+1)*a(n - k*(3*k-1)/2), where A010816(n) = (-1)^m*(2*m+1) if n = m*(m + 1)/2, with m positive, is a triangular number else equals 0. For example, n = 10 = (4*5)/2 is a triangular number, A010816(10) = 9, and so a(10) = 9 + a(9) + a(8) - a(5) - a(3) = 9 - 2 - 2 - 2 - 2 = 1. - Peter Bala, Apr 06 2022

A286354 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=1} (1 - x^j)^k.

Original entry on oeis.org

1, 1, 0, 1, -1, 0, 1, -2, -1, 0, 1, -3, -1, 0, 0, 1, -4, 0, 2, 0, 0, 1, -5, 2, 5, 1, 1, 0, 1, -6, 5, 8, 0, 2, 0, 0, 1, -7, 9, 10, -5, 0, -2, 1, 0, 1, -8, 14, 10, -15, -4, -7, 0, 0, 0, 1, -9, 20, 7, -30, -6, -10, 0, -2, 0, 0, 1, -10, 27, 0, -49, 0, -5, 8, 0, -2, 0, 0, 1, -11, 35, -12, -70, 21, 11, 25, 9, 0, 1, 0, 0
Offset: 0

Views

Author

Ilya Gutkovskiy, May 08 2017

Keywords

Comments

A(n,k) number of partitions of n into an even number of distinct parts minus number of partitions of n into an odd number of distinct parts with k types of each part.

Examples

			A(3,2) = 2 because we have [2, 1], [2', 1], [2, 1'], [2', 1'] (number of partitions of 3 into an even number of distinct parts with 2 types of each part), [3], [3'] (number of partitions of 3 into an odd number of distinct parts with 2 types of each part) and 4 - 2 = 2.
Square array begins:
1,  1,  1,  1,  1,   1,  ...
0, -1, -2, -3, -4,  -5,  ...
0, -1, -1,  0,  2,   5,  ...
0,  0,  2,  5,  8,  10,  ...
0,  0,  1,  0, -5, -15,  ...
0,  1,  2,  0, -4,  -6,  ...
		

Crossrefs

Main diagonal gives A008705.
Antidiagonal sums give A299105.

Programs

  • Maple
    A:= proc(n, k) option remember; `if`(n=0, 1, -k*
          add(numtheory[sigma](j)*A(n-j, k), j=1..n)/n)
        end:
    seq(seq(A(n, d-n), n=0..d), d=0..12);  # Alois P. Heinz, Jun 21 2018
  • Mathematica
    Table[Function[k, SeriesCoefficient[Product[(1 - x^i)^k , {i, Infinity}], {x, 0, n}]][j - n], {j, 0, 12}, {n, 0, j}] // Flatten
    Table[Function[k, SeriesCoefficient[QPochhammer[x, x, Infinity]^k, {x, 0, n}]][j - n], {j, 0, 12}, {n, 0, j}] // Flatten
    Table[Function[k, SeriesCoefficient[Sum[(-1)^i*x^(i*(3*i + 1)/2), {i, -Infinity, Infinity}]^k, {x, 0, n}]][j - n], {j, 0, 12}, {n, 0, j}] // Flatten

Formula

G.f. of column k: Product_{j>=1} (1 - x^j)^k.
G.f. of column k: (Sum_{j=-inf..inf} (-1)^j*x^(j*(3*j+1)/2))^k.
Column k is the Euler transform of period 1 sequence [-k, -k, -k, ...].

A000729 Expansion of Product_{k >= 1} (1 - x^k)^6.

Original entry on oeis.org

1, -6, 9, 10, -30, 0, 11, 42, 0, -70, 18, -54, 49, 90, 0, -22, -60, 0, -110, 0, 81, 180, -78, 0, 130, -198, 0, -182, -30, 90, 121, 84, 0, 0, 210, 0, -252, -102, -270, 170, 0, 0, -69, 330, 0, -38, 420, 0, -190, -390, 0, -108, 0, 0, 0, -300, 99, 442, 210, 0, 418, -294, 0, 0, -510, 378, -540, 138, 0
Offset: 0

Views

Author

Keywords

Comments

This is Glaisher's function lambda(m). It appears to be defined only for odd m, and lambda(4t-1) = 0 (t >= 1), lambda(4t+1) = a(t) (t >= 0). - N. J. A. Sloane, Nov 25 2018
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Number 36 of the 74 eta-quotients listed in Table I of Martin (1996).
Dickson, v.2, p. 295 briefly states a result of Glaisher, 1883, pp 212-215. This result is that a(n) is the sum over all solutions of 16*n + 4 = x^2 + y^2 + z^2 + w^2 in nonnegative odd integers of chi(x) and is also the sum over all solutions of 8*n + 2 = x^2 + y^2 in nonnegative odd integers of chi(x) * chi(y) where chi(x) = x if x == 1 (mod 4) and -x if x == 3 (mod 4). [Michael Somos, Jun 18 2012]
Denoted by g_3(q) in Cynk and Hulek on page 8 as the unique weight 3 Hecke eigenform of level 16 with complex multiplication by i. - Michael Somos, Aug 24 2012
This is a member of an infinite family of integer weight modular forms. g_1 = A008441, g_2 = A002171, g_3 = A000729, g_4 = A215601, g_5 = A215472. - Michael Somos, Aug 24 2012

Examples

			G.f. = 1 - 6*x + 9*x^2 + 10*x^3 - 30*x^4 + 11*x^6 + 42*x^7 - 70*x^9 + 18*x^10 + ...
G.f. = q - 6*q^5 + 9*q^9 + 10*q^13 - 30*q^17 + 11*q^25 + 42*q^29 - 70*q^37 + ...
		

References

  • 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. 2, p. 295, and vol. 3, p. 134.
  • J. W. L. Glaisher, On the representations of a number as a sum of four squares, and on some allied arithmetical functions, Quarterly Journal of Pure and Applied Mathematics, 36 (1905), 305-358. See page 340.
  • J. W. L. Glaisher, The arithmetical functions P(m), Q(m), Omega(m), Quart. J. Math, 37 (1906), 36-48.
  • Morris Newman, A table of the coefficients of the powers of eta(tau). Nederl. Akad. Wetensch. Proc. Ser. A. 59 = Indag. Math. 18 (1956), 204-216.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Powers of Euler's product: A000594, A000727 - A000731, A000735, A000739, A002107, A010815 - A010840.

Programs

  • Magma
    A := Basis( ModularForms( Gamma1(16), 3), 274); A[2] - 6*A[6] + 9*A[10] + 10*A[14] - 30*A[18]; /* Michael Somos, May 17 2015 */
    
  • Magma
    A := Basis( CuspForms( Gamma1(16), 3), 274); A[1] - 6*A[5]; /* Michael Somos, Jan 09 2017 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ 1/16 EllipticTheta[ 4, 0, q] EllipticTheta[ 2, 0, q]^4 EllipticTheta[ 3, 0, q], {q, 0, 4 n + 1}]; (* Michael Somos, Jun 18 2012 *)
    a[ n_] := If[ n < 0, 0, With[ {m = Sqrt[ 16 n + 4]}, SeriesCoefficient[ Sum[ Mod[k, 2] q^k^2, {k, m}]^3 Sum[ KroneckerSymbol[ -4, k] k q^k^2, {k, m}], {q, 0, 16 n + 4}]]]; (* Michael Somos, Jun 12 2012 *)
    a[ n_] := With[ {m = InverseEllipticNomeQ @ q}, SeriesCoefficient[ Sqrt[(1 - m) m ] (EllipticK[m] 2/Pi)^3 / (4 q^(1/2)), {q, 0, 2 n}]]; (* Michael Somos, Jun 22 2012 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 - x^k, {k, n}]^6, {x, 0, n}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ x]^6, {x, 0, n}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := SeriesCoefficient[ (-1/4) EllipticThetaPrime[ 1, -Pi/4, q] EllipticTheta[ 1, -Pi/4, q]^3, {q, 0, 4 n + 1}]; (* Michael Somos, May 17 2015 *)
    a[ n_] := SeriesCoefficient[ (-1/16) EllipticThetaPrime[ 1, 0, q] EllipticTheta[ 1, -Pi/2, q]^3, {q, 0, 4 n + 1}]; (* Michael Somos, May 17 2015 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^6, n))};
    
  • PARI
    {a(n) = my(A, p, e, x, y, a0, a1); if( n<0, 0, n = 4*n + 1; A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p==2, 0, p%4==3, if( e%2, 0, p^e), forstep( i=1, sqrtint(p), 2, if( issquare( p - i^2, &y), x=i; break)); a0=1; a1 = y = 2*(x^2 - y^2); for( i=2, e, x = y*a1 - p^2*a0; a0=a1; a1=x); a1)))}; /* Michael Somos, Aug 21 2006 */
    
  • PARI
    {a(n)=local(tn=(sqrtint(8*n+1)+1)\2);polcoeff(sum(m=0,tn,(1+2*m)^2*x^(m^2+m)+x*O(x^n)) + 2*sum(m=0,tn,sum(k=1,tn,(1+4*(m^2+m-k^2))*x^(m^2+m+k^2)+x*O(x^n))),n)} /* Paul D. Hanna, Mar 15 2010 */
    

Formula

Expansion of q^(-1/4)/16 * theta_2(q)^4 * theta_3(q) * theta_4(q) in powers of q. - [Dickson, v. 3, p. 134] from Stieltjes footnote 160. Michael Somos, Jun 18 2012
Expansion of q^(-1/2) / 4 * k * k' * (K / (pi/2))^3 in powers of q^2 where k, k', K are Jacobi elliptic functions. - Michael Somos, Jun 22 2012
G.f.: Product_{k>0}(1 - x^k)^6.
Given g.f. A(x), then A(q^4) = f(-q^4)^6 = phi(q) * phi(-q) * psi(q^2)^4 where phi(), psi(), f() are Ramanujan theta functions. - Michael Somos, Aug 23 2006
a(n) = b(4*n + 1) where b(n) is multiplicative with b(2^e) = 0^e, b(p^e) = p^e * (1 + (-1)^e) / 2 if p == 3 (mod 4), b(p^e) = b(p) * b(p^(e-1)) - b(p^(e-2)) * p^2 if p == 1 (mod 4) and b(p) = 2 * (x^2 - y^2) where p = x^2 + y^2 and y is even. - Michael Somos, Aug 23 2006
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 64 (t/i)^3 f(t) where q = exp(2 Pi i t). - Michael Somos, Aug 24 2012
G.f.: Sum_{k>=0} a(k) * x^(4*k + 1) = (1/2) * Sum_{u,v in Z} (u*u - 4*v*v) * x^(u*u + 4*v*v). - Michael Somos, Jun 14 2007
G.f.: eta(x)^6 = Sum_{n>=0} (1+2n)^2*x^(n^2+n) + 2*Sum_{n>=0,k>=1} (1 + 4(n^2+n-k^2))*x^(n^2+n+k^2) - from the Milne and Leininger reference. [Paul D. Hanna, Mar 15 2010]
a(0) = 1, a(n) = -(6/n)*Sum_{k=1..n} A000203(k)*a(n-k) for n > 0. - Seiichi Manyama, Mar 26 2017
G.f.: exp(-6*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 05 2018
Let M be a positive integer whose prime factors are all congruent to 3 (mod 4) - see A004614. Then a( M^2*n + (M^2 - 1)/4 ) = M^2*a(n). See Cooper et al., equation 5. - Peter Bala, Dec 01 2020
a(n) = b(4*n + 1) where b(n) is multiplicative with b(2^e) = 0^e, b(p^e) = p^e * (1 + (-1)^e) / 2 if p == 3 (mod 4), b(p^e) = ((x+y*i)^(2*e+2) - (x-y*i)^(2*e+2))/((x+y*i)^2 - (x-y*i)^2) if p == 1 (mod 4) where p = x^2 + y^2 and x is odd. - Jianing Song, Mar 19 2022

A030206 Expansion of q^(-1/3) * eta(q)^2 * eta(q^3)^2 in powers of q.

Original entry on oeis.org

1, -2, -1, 0, 5, 4, -7, 0, -5, 2, -4, 0, 11, 0, 8, 0, -6, -10, 0, 0, -1, -8, 5, 0, -7, 14, 17, 0, 0, 0, -5, 0, -19, 10, -13, 0, 2, -4, 0, 0, -11, 8, 20, 0, 7, 0, 23, 0, 0, -22, -19, 0, 14, 0, -25, 0, 12, -16, 5, 0, -7, 0, 0, 0, 23, 12, 11, 0, 0, 20, -13, 0, 4, 0, -28, 0, -22, 0, 0, 0, 17, 2, -35, 0, 0, 16, -11, 0, 0, -10
Offset: 0

Views

Author

Keywords

Comments

Number 44 of the 74 eta-quotients listed in Table I of Martin (1996).
Denoted by g_2(q) in Cynk and Hulek in Remark 3.4 on page 12 as the unique weight 2 newform of level 27.
This is a member of an infinite family of integer weight modular forms. g_1 = A033687, g_2 = A030206, g_3 = A130539, g_4 = A000731. - Michael Somos, Aug 24 2012
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

Examples

			G.f. = 1 - 2*x - x^2 + 5*x^4 + 4*x^5 - 7*x^6 - 5*x^8 + 2*x^9 - 4*x^10 + 11*x^12 + ...
G.f. = q - 2*q^4 - q^7 + 5*q^13 + 4*q^16 - 7*q^19 - 5*q^25 + 2*q^28 - 4*q^31 + ...
		

Crossrefs

Programs

  • Magma
    A := Basis( ModularForms( Gamma0(27), 2), 271); A[2] - 2*A[5]; /* Michael Somos, Jun 12 2014 */
    
  • Magma
    qEigenform( EllipticCurve( [0, 0, 1, 0, 0]), 271); /* Michael Somos, Jun 12 2014 */
    
  • Magma
    Basis( CuspForms( Gamma0(27), 2), 271)[1]; /* Michael Somos, Mar 24 2015 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ (QPochhammer[ x] QPochhammer[ x^3])^2, {x, 0, n}]; (* Michael Somos, Jun 12 2014 *)
  • PARI
    {a(n) = my(A, p, e, x, y, a0, a1); if( n<0, 0, n = 3*n + 1; A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k,]; if( p==3, 0, p%3==2, if( e%2, 0, (-1)^(e/2) * p^(e/2)), for( i=1, sqrtint(4*p\27), if( issquare(4*p - 27*i^2, &y), break)); a0=1; a1 = y*= (-1)^(y%3); for( i=2, e, x = y*a1 - p*a0; a0=a1; a1=x); a1)))}; /* Michael Somos, Aug 13 2006 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x + A)^2 * eta(x^3 + A)^2, n))}; /* Michael Somos, Feb 19 2007 */
    
  • PARI
    {a(n) = ellak( ellinit( [0, 0, 1, 0, 0], 1), 3*n + 1)}; /* Michael Somos, Jun 12 2014 */
    
  • Sage
    ModularForms( Gamma0(27), 2, prec=271).0; # Michael Somos, Jun 12 2014
    

Formula

Expansion of q^(-1/3) * b(q) * c(q) / 3 in powers of q where b(), c() are cubic AGM theta functions. - Michael Somos, Nov 01 2006
Coefficients of L-series for elliptic curve "27a3": y^2 + y = x^3. - Michael Somos, Aug 13 2006
Euler transform of period 3 sequence [-2, -2, -4, ...]. - Michael Somos, Dec 06 2004
G.f. is a period 1 Fourier series which satisfies f(-1 / (27 t)) = 27 (t/i)^2 f(t) where q = exp(2 Pi i t).
G.f.: Product_{k>0} (1 - x^k)^2 * (1 - x^(3*k))^2.
a(n) = b(3*n + 1) where b(n) is multiplicative with b(3^e) = 0^e, b(p^e) = (1 + (-1)^e) / 2 * (-1)^(e/2) * p^(e/2), if p == 2 (mod 3), otherwise b(p^e) = b(p) * b(p^(e-1)) - p * b(p^(e-2)). - Michael Somos, Aug 13 2006
Given g.f. A(x), then B(q)= q*A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = v^3 - u*w * (u + 4*w). - Michael Somos, Dec 06 2004
a(4*n + 3) = a(16*n + 13) = 0. - Michael Somos, Oct 19 2005
a(4*n + 1) = -2 * a(n). - Michael Somos, Dec 06 2004
a(25*n + 8) = -5 * a(n). Convolution square of A030203. - Michael Somos, Mar 13 2012

A322430 Numbers k such that the coefficient of x^k in the expansion of Product_{j>=1} (1-x^j)^8 is zero.

Original entry on oeis.org

3, 7, 11, 13, 15, 18, 19, 23, 27, 28, 29, 31, 35, 38, 39, 43, 45, 47, 48, 51, 53, 55, 59, 61, 62, 63, 67, 68, 71, 73, 75, 77, 78, 79, 83, 84, 87, 88, 91, 93, 95, 98, 99, 103, 106, 107, 109, 111, 113, 115, 117, 118, 119, 123, 125, 127, 128, 130, 131, 135, 138, 139, 141
Offset: 1

Views

Author

Seiichi Manyama, Dec 07 2018

Keywords

Comments

Indices of zero entries in A000731.
Complement of A267137. - Kemoneilwe Thabo Moseki, Dec 12 2019

Crossrefs

Numbers k such that the coefficient of x^k in the expansion of Product_{j>=1} (1 - x^j)^m is zero: A090864 (m=1), A213250 (m=2), A014132 (m=3), A302056 (m=4), A302057 (m=5), A020757 (m=6), this sequence (m=8), A322431 (m=10), A322432 (m=14), A322043 (m=15), A322433 (m=26).

Programs

  • PARI
    my(x='x+O('x^160)); Vec(select(x->(x==0), Vec(eta(x)^8 - 1), 1)) \\ Michel Marcus, Dec 08 2018

A339707 Dirichlet g.f.: Product_{k>=2} (1 - k^(-s))^8.

Original entry on oeis.org

1, -8, -8, 20, -8, 56, -8, 0, 20, 56, -8, -104, -8, 56, 56, -70, -8, -104, -8, -104, 56, 56, -8, -104, 20, 56, 0, -104, -8, -328, -8, 64, 56, 56, 56, 36, -8, 56, 56, -104, -8, -328, -8, -104, -104, 56, -8, 456, 20, -104, 56, -104, -8, -104, 56, -104, 56, 56, -8, 344, -8, 56, -104, 56, 56
Offset: 1

Views

Author

Ilya Gutkovskiy, Dec 13 2020

Keywords

Crossrefs

Formula

a(1) = 1; a(n) = -Sum_{d|n, d < n} A339323(n/d) * a(d).
a(p^k) = A000731(k) for prime p.

A010818 Expansion of Product (1 - x^k)^10 in powers of x.

Original entry on oeis.org

1, -10, 35, -30, -105, 238, 0, -260, -165, 140, 1054, -770, -595, 0, -715, 2162, 455, 0, -2380, -1820, 2401, -680, 1495, 3080, 1615, -6958, -1925, 0, 0, 5100, -1442, 8330, -5355, 1330, 0, -16790, 0, 8190, 8265, 0, 1918, 0, 8415, -10230, -7140, -9362
Offset: 0

Views

Author

Keywords

Examples

			G.f. = 1 - 10*x + 35*x^2 - 30*x^3 - 105*x^4 + 238*x^5 - 260*x^7 - 165*x^8 + ...
G.f. = q^5 - 10*q^17 + 35*q^29 - 30*q^41 - 105*q^53 + 238*q^65 - 260*q^89 + ...
		

References

  • Newman, Morris; A table of the coefficients of the powers of eta(tau). Nederl. Akad. Wetensch. Proc. Ser. A. 59 = Indag. Math. 18 (1956), 204-216.

Crossrefs

Powers of Euler's product: A000594, A000727 - A000731, A000735, A000739, A002107, A010815 - A010840.

Programs

  • Mathematica
    a[ n_] := SeriesCoefficient[ QPochhammer[ x]^10, {x, 0, n}]; (* Michael Somos, Jun 24 2013 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( eta(x + x * O(x^n))^10, n))}; /* Michael Somos, Jun 09 2011 */
    
  • PARI
    {a(n) = local(m, x, y, z); if( n<0, 0, m = 12*n + 5; z = 0; for( x = -sqrtint(m), sqrtint(m), if( x%6 != 1, next); if( issquare( m - x^2, &y), if( y%6 == 2, y = -y); if( y%6 == 4, z += x*y * (x*x - y*y) ))); z / 6)}; /* Michael Somos, Jun 09 2011 */
    
  • PARI
    {a(n) = local(A, p, e, i, x, y, a0, a1); if( n<0, 0, n = 12*n + 5; A = factor(n); 1 / 48 * prod( k=1, matsize(A)[1], if( p=A[k,1], e=A[k,2]; if( p<5, 0, if( p%12 > 6, if( e%2, 0, p^(2*e)), forstep( i = 1, sqrtint( p), 2, if( issquare( p - i^2, &y), x=i; break)); if( p%12 == 5, a1 = 8 * x*y * (x-y) * (x+y) * (-1)^((x%6==1) + (y%6==4)), a1 = 2 * (x^2-y^2+2*x*y) * (x^2-y^2-2*x*y) * (-1)^(x%6==3) ); a0 = 1; y = a1; for( i=2, e, x = y * a1 - p^4 * a0; a0=a1; a1=x); a1 )))))}; /* Michael Somos, Jun 24 2013 */

Formula

Expansion of f(-x)^10 in powers of x where f() is a Ramanujan theta function.
Expansion of q^(-5/12) * eta(q)^10 in powers of q. - Michael Somos, Jun 09 2011
a(n) = b(12*n + 5) / 48 where b() is multiplicative with b(2^e) = b(3^e) = 0^e, b(p^e) = (1 + (-1)^e) / 2 * p^(2*e) if p == 7 or 11 (mod 12), b(p^e) = b(p) * b(p^(e-1)) - p^4 * b(p^(e-2)) if p == 1 or 5 (mod 12). - Michael Somos, Jun 24 2013
G.f. is a period 1 Fourier series which satisfies f(-1 / (144 t)) = 12^5 (t/i)^5 f(t) where q = exp(2 Pi i t). - Michael Somos, Jan 06 2014
G.f.: Product_{k>0} (1 - x^k)^10. a(49*n + 20) = 2401 * a(n).
48 * a(n) = A234565(3*n + 1). a(7*n + 2) = 0 unless n == 2 (mod 7). - Michael Somos, Jul 18 2014
a(0) = 1, a(n) = -(10/n)*Sum_{k=1..n} A000203(k)*a(n-k) for n > 0. - Seiichi Manyama, Mar 27 2017
G.f.: exp(-10*Sum_{k>=1} x^k/(k*(1 - x^k))). - Ilya Gutkovskiy, Feb 05 2018
Let M = p_1*...*p_k be a positive integer whose prime factors p_i (not necessarily distinct) are congruent to 7 (mod 12) or 11 (mod 12). Then a( M^2*n + 10*(M^2 - 1)/24 ) = M^4*a(n). See Cooper et al., Theorem 1. - Peter Bala, Dec 01 2020

A033690 Theta series of A2[hole]^4.

Original entry on oeis.org

1, 4, 14, 28, 57, 84, 148, 196, 312, 364, 546, 624, 910, 988, 1352, 1456, 1974, 2072, 2710, 2800, 3705, 3724, 4816, 4788, 6188, 6076, 7658, 7644, 9620, 9352, 11536, 11284, 14183, 13468, 16542, 15996, 19864, 18928, 22820, 21904, 26880, 25284
Offset: 0

Views

Author

Keywords

Examples

			q^4 + 4*q^7 + 14*q^10 + 28*q^13 + 57*q^16 + 84*q^19 + 148*q^22 + ...
		

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 111, Eq (63)^4.

Crossrefs

Cf. A033685.

Programs

  • Mathematica
    s = (QPochhammer[q^3]^3/QPochhammer[q])^4 + O[q]^45; CoefficientList[s, q] (* Jean-François Alcover, Nov 04 2015 *)
  • PARI
    {a(n) = local(A); if(n<0, 0, A = x*O(x^n); polcoeff( (eta(x^3 +A)^3 / eta(x +A) )^4, n))} /* Michael Somos, Aug 22 2007 */

Formula

a(n) = A033685^4.
Expansion of q^(-4/3) * (eta(q^3)^3 / eta(q))^4 in powers of q. - Michael Somos, Aug 22 2007
Expansion of c(q)^4 / (81 * q^(4/3)) in powers of q where c() is a cubic AGM function. - Michael Somos, Aug 22 2007
Euler transform of period 3 sequence [ 4, 4, -8, ...]. - Michael Somos, Aug 22 2007
A092342(n) = A000731(n) + 81*a(n-1). - Michael Somos, Aug 22 2007
Showing 1-10 of 20 results. Next