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 10 results.

A080054 G.f.: Product_{n >= 0} (1+x^(2n+1))/(1-x^(2n+1)).

Original entry on oeis.org

1, 2, 2, 4, 6, 8, 12, 16, 22, 30, 40, 52, 68, 88, 112, 144, 182, 228, 286, 356, 440, 544, 668, 816, 996, 1210, 1464, 1768, 2128, 2552, 3056, 3648, 4342, 5160, 6116, 7232, 8538, 10056, 11820, 13872, 16248, 18996, 22176, 25844, 30068, 34936, 40528
Offset: 0

Views

Author

Michael Somos, Jan 26 2003

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
G.f. for pairs of partitions of type R.
G.f. for the number of partitions of 2n in which all odd parts occur with multiplicity 2 and the even parts occur with multiplicity 1. Also g.f. for the number of partitions of 2n free of multiples of 4. All odd parts occur with even multiplicities. The even parts occur with multiplicity 1. - Noureddine Chair, Feb 10 2005
This is also the number of overpartitions of an integer into odd parts. - James Sellers, Feb 18 2008
The Higher Algebra reference on page 517 has an unnumbered example between 251 and 252: "If u^6-v^6+5u^2v^2(u^2-v^2)+4uv(1-u^4v^4)=0, prove that (u^2-v^2)^6=16u^2v^2(1-u^8)(1-v^8). [PEMB. COLL. CAMB.]". It turns out that this is two forms of the modular equation of degree 5. - Michael Somos, May 12 2011
Convolution of A000009 and A000700. - Vaclav Kotesovec, Aug 23 2015
Let F(x) = Product_{n >= 0} (1 + x^(2*n+1))/(1 - x^(2*n+1)). Let m be a nonzero integer. The simple continued fractions expansions of the real numbers F(1/m) may be predictable - given a positive integer n, the sequence of the n-th partial denominators of the continued fraction expansion of F(1/m) may be polynomial or quasi-polynomial in m for sufficiently large m. An example is given below. - Peter Bala, Nov 03 2019

Examples

			G.f. = 1 + 2*q + 2*q^2 + 4*q^3 + 6*q^4 + 8*q^5 + 12*q^6 + 16*q^7 + 22*q^8 + 30*q^9 + ...
From _Peter Bala_, Nov 03 2019: (Start)
F(x) := Product_{n >= 0} (1 + x^(2*n+1))/(1 - x^(2*n+1)).
Simple continued fraction expansions of F(1/(2*m)):
  m=2 [1; 1, 2, 1, 1, 1, 1, 2, 1, 2,   33, 1, 3,  7, 4,  33, 1, 8,  4,    2, 1,...]
  m=3 [1; 2, 2, 2, 1, 1, 2, 2, 2, 2,  110, 1, 2, 46, 3, 110, 1, 3, 12,    1, 7,...]
  m=4 [1; 3, 2, 3, 1, 1, 3, 2, 3, 2,  259, 1, 1,  1, 2,  15, 2, 1,  2,  259, 1,...]
  m=5 [1; 4, 2, 4, 1, 1, 4, 2, 4, 2,  504, 1, 1,  1, 1,  78, 1, 1,  2,  504, 1,...]
  m=6 [1; 5, 2, 5, 1, 1, 5, 2, 5, 2,  869, 1, 1,  2, 2,  23, 2, 2,  2,  869, 1,...]
  m=7 [1; 6, 2, 6, 1, 1, 6, 2, 6, 2, 1378, 1, 1,  2, 1, 110, 1, 2,  2, 1378, 1,...]
  m=8 [1; 7, 2, 7, 1, 1, 7, 2, 7, 2, 2055, 1, 1,  3, 2,  31, 2, 3,  2, 2055, 1,...]
  m=9 [1; 8, 2, 8, 1, 1, 8, 2, 8, 2, 2924, 1, 1,  3, 1, 142, 1, 3,  2, 2924, 1,...]
The sequence of the 10th partial denominators [33,110,259,504,...], starting at m = 2, appears to be given by the polynomial 4*m^3 + m - 1.
The sequence of the 15th partial denominators [15,78,23,110,31,142,...], starting at m = 4, appears to be quasi-polynomial in m, with constituent polynomials 4*m - 1 and 16*m - 2. (End)
		

References

  • B. C. Berndt, Ramanujan's theory of theta-functions, Theta functions: from the classical to the modern, Amer. Math. Soc., Providence, RI, 1993, pp. 1-63. MR 94m:11054.
  • A. Cayley, An Elementary Treatise on Elliptic Functions, 2nd ed., G. Bell and Sons, 1895, p. 245, Art. 333.
  • J. W. L. Glaisher, Identities, Messenger of Mathematics, 5 (1876), pp. 111-112. see Eq. VI
  • J. W. L. Glaisher, On Some Continued Fractions, Messenger of Mathematics, 7 (1878), pp. 67-68, see p. 68
  • H. S. Hall and S. R. Knight, Higher Algebra, Macmillan, 1957, p. 517.

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-2) +add(2*b(n-i*j, i-2), j=1..n/i)))
        end:
    a:= n-> b(n, n-1+irem(n, 2)):
    seq(a(n), n=0..50);  # Alois P. Heinz, Feb 10 2014
    # alternative program using expansion of f(x, x^3) / f(-x, -x^3):
    with(gfun): series( add(x^(n*(2*n-1)), n = -8..8)/add((-1)^n*x^(n*(2*n-1)), n = -8..8), x, 100): seriestolist(%); # Peter Bala, Feb 05 2021
  • Mathematica
    a[ n_] := With[ {m = InverseEllipticNomeQ @ q}, SeriesCoefficient[ (1 - m )^(-1/8), {q, 0, n}]]; (* Michael Somos, Aug 03 2011 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q] / EllipticTheta[ 4, 0, q])^(1/2), {q, 0, n}]; (* Michael Somos, Aug 03 2011 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -q] / QPochhammer[ q], {q, 0, n}]; (* Michael Somos, May 10 2014 *)
    a[ n_] := SeriesCoefficient[ QHypergeometricPFQ[ {-1}, {}, q^2, q], {q, 0, n}]; (* Michael Somos, May 10 2014 *)
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, b[n, i - 2] + Sum[2*b[n - i*j, i - 2], {j, 1, n/i}]]];
    a[n_] := b[n, n - 1 + Mod[n, 2]];
    Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Nov 05 2017, after Alois P. Heinz *)
  • PARI
    {a(n) = my(A, m); if( n<0, 0, m=1; A = 1 + 2*x + O(x^2); while( m
    				
  • PARI
    a(n)=polcoeff(exp(2*sum(k=0,n\2,sigma(2*k+1)/(2*k+1)*x^(2*k+1))),n) /* Paul D. Hanna */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^3 / (eta(x + A)^2 * eta(x^4 + A)), n))}; /* Michael Somos, Jul 07 2005 */

Formula

Expansion of f(q) / f(-q) in powers of q where f() is a Ramanujan theta function.
Expansion of (1 - k^2)^(-1/8) = k'^(-1/4) in powers of the nome q = exp(-Pi K'/K).
Expansion of eta(q^2)^3 / (eta(q^4) * eta(q)^2) in powers of q.
Euler transform of period 4 sequence [ 2, -1, 2, 0, ...].
(theta_3(q) / theta_4(q))^(1/2) = (phi(q) / phi(-q))^(1/2) = chi(q) / chi(-q) = psi(q) / psi(-q) = f(q) / f(-q) where phi{}, chi(), psi(), f() are Ramanujan theta functions.
G.f.: A(x) = exp( 2*sum_{n>=0} sigma(2*n+1)/(2*n+1)*x^(2*n+1) ). - Paul D. Hanna, Mar 01 2004
G.f. satisfies: A(-x) = 1/A(x), (A(x)+A(-x))/2 = A(x^2)*A(x^4)^2, A(x) = sqrt((A(x^2)^4+1)/2) + sqrt((A(x^2)^4-1)/2). - Paul D. Hanna, Mar 27 2004
Another g.f.: 1/product_{ k>= 1 } (1+x^(2*k))*(1-x^(2*k-1))^2. - Vladeta Jovovic, Mar 29 2004
G.f. A(x) satisfies 0 = f(A(x), A(x^3)) where f(u, v) = (u - v^3) * (v + 2*u^3) - u * (u^3 - v). - Michael Somos, Aug 03 2011
G.f. A(x) satisfies 0 = f(A(x), A(x^5)) where f(u, v) = (u^2 - v^2)^6 - 16 * u^2 * v^2 * (1 - u^8) * (1 - v^8). - Michael Somos, May 12 2011
G.f. A(x) satisfies 0 = f(A(x), A(x^7)) where f(u, v) = (1 - u^8) * (1 - v^8) - (1 - u*v)^8. - Michael Somos, Jan 01 2006
G.f. is a period 1 Fourier series which satisfies f(-1 / (32 t)) = 2^(-1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A029838. - Michael Somos, Aug 03 2011
G.f.: (theta_3/theta_4)^(1/2) = ((Sum_{k in Z} x^(k^2))/(Sum_{k in Z} (-x)^(k^2)))^(1/2) = Product_{k>0} (1 - x^(4k-2))/((1 - x^(4k-1))(1 - x^(4k-3)))^2.
G.f.: Product_{ k >= 1 } (1 + x^(2*k-1))*(1 + x^k) = product_{ k >= 1 } (1 + x^(2*k-1))/(1 - x^(2*k-1)).
G.f.: 1 + 2*x / (1 - x) + 2*x^3 * (1 + x) / ((1 - x)*(1 - x^2)) + 2*x^6 * (1 + x)*(1 + x^2) / ((1 - x)*(1 - x^2)*(1 - x^3)) + ... [Glaisher 1876] - Michael Somos, Jun 20 2012
G.f.: 1 / (1 - 2*x / (1 + x - (x^2 - x^4) / (1 + x^3 - (x^3 - x^7) / (1 + x^5 - (x^4 - x^10) / (1 + x^7 - ...))))) [Glaisher 1878] - Michael Somos, Jun 24 2012
a(n) = (-1)^floor(n/2) * A080015(n) = (-1)^n * A108494(n). Convolution inverse is A108494. Convolution square is A007096.
Empirical : Sum_{n>=0} exp(-Pi)^n * a(n) = 2^(1/8). - Simon Plouffe, Feb 20 2011
Empirical : Sum_{n>=0} (-exp(-Pi))^n * a(n) = 1/2^(1/8). - Simon Plouffe, Feb 20 2011
a(n) ~ Pi * BesselI(1, Pi*sqrt(n/2)) / (4*sqrt(n)) ~ exp(Pi*sqrt(n/2)) / (2^(9/4) * n^(3/4)) * (1 - 3/(4*Pi*(sqrt(2*n))) - 15/(64*Pi^2*n)). - Vaclav Kotesovec, Aug 23 2015, extended Jan 09 2017
Simon Plouffe's empirical observations are true. Furthermore, for every positive rational p, Sum_{n>=0} exp(-Pi*sqrt(p))^n * a(n) = 1/(Sum_{n>=0} (-exp(-Pi*sqrt(p)))^n * a(n)) is an algebraic number (see the MathOverflow link). - Vladimir Reshetnikov, Nov 23 2016
G.f.: f(x,x^3)/f(-x,-x^3) = ( Sum_{n = -oo..oo} x^(n*(2*n-1)) )/( Sum_{n = -oo..oo} (-1)^n*x^(n*(2*n-1)) ), where f(a,b) = Sum_{n = -oo..oo} a^(n*(n+1)/2)*b^(n*(n-1)/2) is Ramanujan's 2-variable theta function. - Peter Bala, Feb 05 2021
G.f. A(q) = (-lambda(-q)/lambda(q))^(1/8), where lambda(q) = 16*q - 128*q^2 + 704*q^3 - 3072*q^4 + ... is the elliptic modular function in powers of the nome q = exp(i*Pi*t), the g.f. of A115977; lambda(q) = k(q)^2, where k(q) = (theta_2(q) / theta_3(q))^2 is the elliptic modulus. - Peter Bala, Sep 26 2023
Recurrence: a(n) = c(n) + Sum_{k = 1..floor((-1 + sqrt(1 + 8*n))/2)} (-1)^(1 + k*(k+1)/2) * a(n - k*(k+1)/2), where c(n) = 1 if n is a triangular number, otherwise c(n) = 0. See A010054. - Peter Bala, Jun 08 2025

Extensions

Definition simplified by N. J. A. Sloane, Apr 24 2014

A169975 Expansion of Product_{i>=0} (1 + x^(4*i+1)).

Original entry on oeis.org

1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 2, 1, 0, 1, 2, 1, 0, 1, 3, 2, 0, 1, 3, 3, 1, 1, 4, 4, 1, 1, 4, 5, 2, 1, 5, 7, 3, 1, 5, 8, 5, 2, 6, 10, 6, 2, 6, 12, 9, 3, 7, 14, 11, 4, 7, 16, 15, 6, 8, 19, 18, 8, 9, 21, 23, 11, 10, 24, 27, 14, 11, 27, 34, 19, 13, 30, 39, 24, 15, 33, 47, 31, 18, 37, 54, 38
Offset: 0

Views

Author

N. J. A. Sloane, Aug 29 2010

Keywords

Comments

Number of partitions into distinct parts of the form 4*k+1.
In general, if a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} (1 + x^(a*k + b)), then a(n) ~ exp(Pi*sqrt(n/(3*a))) / (2^(1 + b/a) * (3*a)^(1/4) * n^(3/4)) [Meinardus, 1954]. - Vaclav Kotesovec, Aug 26 2015
Convolution of A147599 and A169975 is A000700. - Vaclav Kotesovec, Jan 18 2017

Crossrefs

Programs

  • Mathematica
    nmax = 200; CoefficientList[Series[Product[(1 + x^(4*k+1)), {k, 0, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 26 2015 *)
    nmax = 200; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[If[Mod[k, 4] == 1, Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}]; ], {k, 2, nmax}]; poly (* Vaclav Kotesovec, Jan 18 2017 *)

Formula

G.f.: Sum_{n>=0} (x^(2*n^2 - n) / Product_{k=1..n} (1 - x^(4*k))). - Joerg Arndt, Mar 10 2011
G.f.: G(0)/x where G(k) = 1 - 1/(1 - 1/(1 - 1/(1+(x)^(4*k+1))/G(k+1) )); (recursively defined continued fraction, see A006950). - Sergei N. Gladkovskii, Jan 28 2013
a(n) ~ exp(Pi*sqrt(n)/(2*sqrt(3))) / (2^(7/4) * 3^(1/4) * n^(3/4)) * (1 - (3*sqrt(3)/(4*Pi) + Pi/(192*sqrt(3))) / sqrt(n)). - Vaclav Kotesovec, Aug 26 2015, extended Jan 18 2017

A261612 Expansion of Product_{k>=0} (1 + x^(3*k+1)).

Original entry on oeis.org

1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 2, 1, 1, 2, 1, 1, 3, 2, 1, 3, 3, 2, 4, 4, 2, 4, 5, 3, 5, 7, 4, 5, 8, 6, 7, 10, 7, 7, 12, 10, 9, 14, 12, 10, 16, 16, 13, 19, 19, 15, 22, 24, 19, 25, 28, 22, 29, 35, 28, 33, 40, 33, 38, 48, 41, 44, 55, 48, 51, 66, 59, 58, 74, 69
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 26 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 100; CoefficientList[Series[Product[(1 + x^(3*k+1)), {k, 0, nmax}], {x, 0, nmax}], x]
    nmax = 100; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[If[Mod[k, 3] == 1, Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}];], {k, 2, nmax}]; poly (* Vaclav Kotesovec, Jan 13 2017 *)

Formula

a(n) ~ exp(Pi*sqrt(n)/3) / (2^(4/3) * sqrt(3) * n^(3/4)) * (1 - (Pi/144 + 9/(8*Pi)) / sqrt(n)). - Vaclav Kotesovec, Aug 26 2015, extended Jan 16 2017
G.f.: Sum_{k>=0} x^(k*(3*k - 1)/2) / Product_{j=1..k} (1 - x^(3*j)). - Ilya Gutkovskiy, Nov 24 2020

A261611 Expansion of Product_{k>=0} (1 + x^(4*k+1))/(1 - x^(4*k+1)).

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 6, 6, 6, 8, 12, 14, 14, 16, 22, 28, 30, 32, 40, 50, 56, 60, 70, 86, 98, 106, 120, 144, 166, 180, 200, 234, 270, 296, 324, 372, 428, 472, 514, 580, 664, 736, 800, 890, 1010, 1124, 1222, 1346, 1514, 1684, 1834, 2008, 2240, 2488, 2712, 2956
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 26 2015

Keywords

Comments

In general, if a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} (1 + x^(a*k+b))/(1 - x^(a*k+b)), then a(n) ~ Gamma(b/a) * a^(b/(2*a) - 1/2) * Pi^(b/a - 1) * exp(Pi*sqrt(n/a)) / (2^(2*b/a + 1) * n^(b/(2*a) + 1/2)).

Crossrefs

Programs

  • Mathematica
    nmax = 60; CoefficientList[Series[Product[(1 + x^(4*k+1))/(1 - x^(4*k+1)), {k, 0, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(Pi*sqrt(n)/2) * Gamma(1/4) / (2^(9/4) * Pi^(3/4) * n^(5/8)).

A261616 Expansion of Product_{k>=0} 1/(1 - x^(3*k+1))^2.

Original entry on oeis.org

1, 2, 3, 4, 7, 10, 13, 18, 26, 34, 44, 58, 76, 96, 123, 156, 196, 244, 304, 374, 461, 566, 690, 836, 1015, 1224, 1470, 1762, 2110, 2512, 2987, 3542, 4191, 4944, 5825, 6842, 8025, 9392, 10971, 12788, 14891, 17300, 20068, 23242, 26883, 31034, 35787, 41204
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 26 2015

Keywords

Comments

Self-convolution of A035382.
In general, if a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} 1/(1 - x^(a*k+b))^2, then a(n) ~ Gamma(b/a)^2 * a^(b/a - 3/4) * exp(2*Pi*sqrt(n/(3*a))) * Pi^(2*b/a - 2) / (4 * 3^(b/a - 1/4) * n^(b/a + 1/4)).

Crossrefs

Programs

  • Mathematica
    nmax = 60; CoefficientList[Series[Product[1/(1 - x^(3*k+1))^2, {k, 0, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(2*Pi*sqrt(n)/3) * Gamma(1/3)^2 / (4 * sqrt(3) * Pi^(4/3) * n^(7/12)).

A261615 Expansion of Product_{k>=0} (1 + x^(3*k+1))^2.

Original entry on oeis.org

1, 2, 1, 0, 2, 4, 2, 2, 5, 4, 3, 8, 10, 6, 9, 14, 11, 14, 22, 18, 17, 30, 32, 28, 41, 46, 39, 54, 68, 60, 73, 94, 85, 96, 131, 128, 130, 170, 175, 176, 229, 246, 237, 294, 330, 320, 386, 446, 430, 492, 582, 578, 642, 762, 763, 818, 977, 1008, 1061, 1254, 1311
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 26 2015

Keywords

Comments

Self-convolution of A261612.
In general, if a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} (1 + x^(a*k+b))^2, then a(n) ~ exp(Pi*sqrt(2*n/(3*a))) / (2^(2*b/a + 1/4) * 3^(1/4) * a^(1/4) * n^(3/4)).

Crossrefs

Programs

  • Mathematica
    nmax = 60; CoefficientList[Series[Product[(1 + x^(3*k+1))^2, {k, 0, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(Pi*sqrt(2*n)/3) / (2^(11/12) * sqrt(3) * n^(3/4)).

A261647 Expansion of Product_{k>=0} ((1+x^(2*k+1))/(1-x^(2*k+1)))^3.

Original entry on oeis.org

1, 6, 18, 44, 102, 216, 428, 816, 1494, 2650, 4584, 7740, 12804, 20808, 33264, 52400, 81462, 125100, 189966, 285516, 425016, 627040, 917436, 1331856, 1919332, 2746926, 3905784, 5519352, 7754064, 10833192, 15055216, 20817600, 28647414, 39241336, 53517060
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 28 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax=60; CoefficientList[Series[Product[((1+x^(2*k+1))/(1-x^(2*k+1)))^3,{k,0,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ exp(Pi*sqrt(3*n/2)) * 3^(1/4) / (8 * 2^(1/4) * n^(3/4)).

A261649 Expansion of Product_{k>=0} ((1+x^(3*k+1))/(1-x^(3*k+1)))^2.

Original entry on oeis.org

1, 4, 8, 12, 20, 36, 56, 80, 120, 180, 252, 348, 492, 680, 912, 1228, 1652, 2180, 2856, 3744, 4860, 6256, 8044, 10284, 13048, 16520, 20848, 26140, 32672, 40756, 50596, 62576, 77256, 95060, 116540, 142592, 174036, 211736, 257056, 311448, 376332, 453764, 546160
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 28 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax=60; CoefficientList[Series[Product[((1+x^(3*k+1))/(1-x^(3*k+1)))^2,{k,0,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ exp(Pi*sqrt(2*n/3)) * Gamma(1/3)^2 / (2^(7/4) * 3^(5/12) * Pi^(4/3) * n^(7/12)).

A261651 Expansion of Product_{k>=0} ((1+x^(3*k+1))/(1-x^(3*k+1)))^3.

Original entry on oeis.org

1, 6, 18, 38, 72, 138, 254, 432, 708, 1154, 1836, 2826, 4288, 6456, 9552, 13902, 20070, 28722, 40614, 56916, 79242, 109448, 149988, 204318, 276672, 372288, 498264, 663602, 879252, 1159470, 1522564, 1990788, 2592162, 3362638, 4346244, 5597100, 7183792, 9191004
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 28 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax=60; CoefficientList[Series[Product[((1+x^(3*k+1))/(1-x^(3*k+1)))^3,{k,0,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ exp(Pi*sqrt(n)) * Gamma(1/3)^3 / (4 * Pi^2 * sqrt(3*n)).

A261652 Expansion of Product_{k>=0} ((1+x^(4*k+1))/(1-x^(4*k+1)))^3.

Original entry on oeis.org

1, 6, 18, 38, 66, 108, 182, 306, 486, 728, 1068, 1578, 2318, 3312, 4614, 6388, 8862, 12192, 16488, 22038, 29400, 39156, 51702, 67554, 87810, 113982, 147384, 189200, 241446, 307356, 390408, 493662, 621006, 778712, 974628, 1216284, 1511756, 1872840, 2315538
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 28 2015

Keywords

Comments

In general, if j > 0, a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} ((1 + x^(a*k+b))/(1 - x^(a*k+b)))^j, then a(n) ~ Gamma(b/a)^j * 2^(j/2 - 3/2 - 2*b*j/a) * a^(-j/4 - 1/4 + b*j/(2*a)) * exp(Pi*sqrt(j*n/a)) * j^(1/4 - j/4 + b*j/(2*a)) * Pi^(b*j/a - j) * n^(j/4 - 3/4 - b*j/(2*a)).

Crossrefs

Cf. A015128 (a=1, b=1, j=1), A156616.
Cf. A080054 (a=2, b=1, j=1), A007096 (a=2, b=1, j=2), A261647 (a=2, b=1, j=3), A014969 (a=2, b=1, j=4), A261648 (a=2, b=1, j=5), A014970 (a=2, b=1, j=6), A014972 (a=2, b=1, j=8), A103261 (a=2, b=1, j=10).
Cf. A261610 (a=3, b=1, j=1), A261649 (a=3, b=1, j=2), A261651 (a=3, b=1, j=3).
Cf. A261611 (a=4, b=1, j=1), A261650 (a=4, b=1, j=2), A261652 (a=4, b=1, j=3).

Programs

  • Mathematica
    nmax=60; CoefficientList[Series[Product[((1+x^(4*k+1))/(1-x^(4*k+1)))^3,{k,0,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ exp(Pi*sqrt(3*n)/2) * 2^(1/4) * Gamma(1/4)^3 / (8 * 3^(1/8) * Pi^(9/4) * n^(3/8)).
Showing 1-10 of 10 results.