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

A047241 Numbers that are congruent to {1, 3} mod 6.

Original entry on oeis.org

1, 3, 7, 9, 13, 15, 19, 21, 25, 27, 31, 33, 37, 39, 43, 45, 49, 51, 55, 57, 61, 63, 67, 69, 73, 75, 79, 81, 85, 87, 91, 93, 97, 99, 103, 105, 109, 111, 115, 117, 121, 123, 127, 129, 133, 135, 139, 141, 145, 147, 151, 153, 157, 159, 163, 165, 169, 171, 175, 177, 181, 183
Offset: 1

Views

Author

Keywords

Comments

Also the numbers k such that 10^p+k could possibly be prime. - Roderick MacPhee, Nov 20 2011 This statement can be written as follows. If 10^m + k = prime, for any m >= 1, then k is in this sequence. See the pink box comments by Roderick MacPhee from Dec 09 2014. - Wolfdieter Lang, Dec 09 2014
The odd-indexed terms are one more than the arithmetic mean of their neighbors; the even-indexed terms are one less than the arithmetic mean of their neighbors. - Amarnath Murthy, Jul 29 2003
Partial sums are A212959. - Philippe Deléham, Mar 16 2014
12*a(n) is conjectured to be the length of the boundary after n iterations of the hexagon and square expansion shown in the link. The squares and hexagons have side length 1 in some units. The pattern is supposed to become the planar Archimedean net 4.6.12 when n -> infinity. - Kival Ngaokrajang, Nov 30 2014
Positive numbers k for which 1/2 + k/3 + k^2/6 is an integer. - Bruno Berselli, Apr 12 2018

References

  • L. Lovasz, J. Pelikan, K. Vesztergombi, Discrete Mathematics, Springer (2003); 14.4, p. 225.

Crossrefs

Subsequence of A186422.
Union of A016921 and A016945. - Wesley Ivan Hurt, Sep 28 2013

Programs

  • Haskell
    a047241 n = a047241_list !! (n-1)
    a047241_list = 1 : 3 : map (+ 6) a047241_list
    -- Reinhard Zumkeller, Feb 19 2013
    
  • Maple
    seq(3*k-2-((k+1) mod 2), k=1..100); # Wesley Ivan Hurt, Sep 28 2013
  • Mathematica
    Table[{2, 4}, {30}] // Flatten // Prepend[#, 1]& // Accumulate (* Jean-François Alcover, Jun 10 2013 *)
    Select[Range[200], MemberQ[{1, 3}, Mod[#, 6]]&] (* or *) LinearRecurrence[{1, 1, -1}, {1, 3, 7}, 70] (* Harvey P. Dale, Oct 01 2013 *)
  • PARI
    a(n)=bitor(3*n-3,1) \\ Charles R Greathouse IV, Sep 28 2013
    
  • Python
    for n in range(1,10**5):print(3*n-2-((n+1)%2)) # Soumil Mandal, Apr 14 2016

Formula

From Paul Barry, Sep 04 2003: (Start)
O.g.f.: (1 + 2*x + 3*x^2)/((1 + x)*(1 - x)^2) = (1 + 2*x + 3*x^2)/((1 - x)*(1 - x^2)).
E.g.f.: (6*x + 1)*exp(x)/2 + exp(-x)/2;
a(n) = 3*n - 5/2 - (-1)^n/2. (End)
a(n) = 2*floor((n-1)/2) + 2*n - 1. - Gary Detlefs, Mar 18 2010
a(n) = 6*n - a(n-1) - 8 with n > 1, a(1)=1. - Vincenzo Librandi, Aug 05 2010
a(n) = 3*n - 2 - ((n+1) mod 2). - Wesley Ivan Hurt, Jun 29 2013
a(1)=1, a(2)=3, a(3)=7; for n>3, a(n) = a(n-1) + a(n-2) - a(n-3). - Harvey P. Dale, Oct 01 2013
From Benedict W. J. Irwin, Apr 13 2016: (Start)
A005408(a(n)+1) = A016813(A001651(n)),
A007310(a(n)) = A005408(A087444(n)-1),
A007310(A005408(a(n)+1)) = A017533(A001651(n)). (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(4*sqrt(3)) + log(3)/4. - Amiram Eldar, Dec 11 2021

Extensions

Formula corrected by Bruno Berselli, Jun 24 2010

A098884 Number of partitions of n into distinct parts in which each part is congruent to 1 or 5 mod 6.

Original entry on oeis.org

1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 2, 2, 1, 0, 1, 2, 3, 3, 2, 1, 1, 3, 5, 5, 3, 1, 2, 5, 7, 7, 5, 3, 3, 7, 11, 11, 7, 4, 6, 11, 15, 15, 11, 7, 8, 15, 22, 22, 15, 10, 13, 22, 30, 30, 23, 16, 18, 30, 42, 42, 31, 22, 27, 43, 56, 56, 44, 33, 37, 57, 77, 77, 59, 45, 53, 79, 101, 101, 82, 64, 71
Offset: 0

Views

Author

Noureddine Chair, Oct 14 2004

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Convolution of A281244 and A280456. - Vaclav Kotesovec, Jan 18 2017

Examples

			E.g. a(25)=5 because 25=19+5+1=17+7+1=13+7+5=13+11+1.
G.f. = 1 + x + x^5 + x^6 + x^7 + x^8 + x^11 + 2*x^12 + 2*x^13 + x^14 + x^16 + ...
G.f. = q + q^13 + q^61 + q^73 + q^85 + q^97 + q^133 + 2*q^145 + 2*q^157 + q^169 + ...
		

Crossrefs

Programs

  • Haskell
    a098884 = p a007310_list where
       p _  0     = 1
       p (k:ks) m = if k > m then 0 else p ks (m - k) + p ks m
    -- Reinhard Zumkeller, Feb 19 2013
  • Maple
    series(product((1+x^(6*k-1))*(1+x^(6*k-5)),k=1..100),x=0,100);
  • Mathematica
    a[ n_] := SeriesCoefficient[ Product[ 1 - (-x)^k + x^(2 k), {k, n}], {x, 0, n}]; (* Michael Somos, Sep 20 2013 *)
    a[ n_] := SeriesCoefficient[ 1 / Product[ 1 - x^k + x^(2 k), {k, 1, n, 2}], {x, 0, n}]; (* Michael Somos, Sep 20 2013 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 + x^k, {k, 1, n, 2}] / Product[ 1 + x^k, {k, 3, n, 6}], {x, 0, n}]; (* Michael Somos, Sep 20 2013 *)
    a[ n_] := SeriesCoefficient[ Product[ 1 + x^k, {k, 1, n, 6}] Product[ 1 + x^k, {k, 5, n, 6}], {x, 0, n}]; (* Michael Somos, Sep 20 2013 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^6] QPochhammer[ -x^5, x^6], {x, 0, n}]; (* Michael Somos, Sep 20 2013 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^2] / QPochhammer[ -x^3, x^6], {x, 0, n}]; (* Michael Somos, Sep 20 2013 *)
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^2 * eta(x^3 + A) * eta(x^12 + A) / (eta(x + A) * eta(x^4 + A) * eta(x^6 + A)^2), n))}; /* Michael Somos, Jun 26 2005 */
    
  • PARI
    {a(n) = my(A, m); if( n<0, 0, A = x * O(x^n); m = sqrtint(3*n + 1); polcoeff( sum(k= -((m-1)\3), (m+1)\3, x^(k * (3*k - 2)), A) / eta(x^6 + A), n))}; /* Michael Somos, Sep 20 2013 */
    

Formula

Expansion of chi(x) / chi(x^3) in powers of x where chi() is a Ramanujan theta function. - Michael Somos, Sep 20 2013
Expansion of f(x^1, x^5) / f(-x^6) in powers of x where f(,) is a Ramanujan theta function. - Michael Somos, Sep 20 2013
Expansion of G(x^6) * H(-x) + x * G(-x) * H(x^6) where G() (A003114), H() (A003106) are Rogers-Ramanujan functions.
Expansion of q^(-1/12) * eta(q^2)^2 * eta(q^3) * eta(q^12) / (eta(q) * eta(q^4) * eta(q^6)^2) in powers of q.
Euler transform of period 12 sequence [ 1, -1, 0, 0, 1, 0, 1, 0, 0, -1, 1, 0, ...]. - Michael Somos, Jun 26 2005
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = g(t) where q = exp(2 Pi i t) and g() is the g.f. for A227398. - Michael Somos, Sep 20 2013
G.f.: Product_{k>0} (1 - (-x)^k + x^(2*k)).
G.f.: 1 / Product_{k>0} (1 - x^(2*k - 1) + x^(4*k - 2)).
G.f.: 1 / Product_{k>0} ((1 + x^(6*k - 3)) / (1 + x^(2*k - 1))).
G.f.: Product_{k>0} ((1 + x^(6*k - 1)) * (1 + x^(6*k - 5))).
G.f.: 1 / Product_{k>0} (1 + (-x)^(3*k - 1)) * (1 + (-x)^(3*k - 2)).
G.f.: (Sum_{k in Z} x^(k * (3*k - 2))) / (Sum_{k in Z} (-1)^k * x^(3*k * (3*k-1))).
A109389(n) = (-1)^n * a(n). Convolution inverse of A227398.
a(n) ~ exp(sqrt(n)*Pi/3)/ (2*sqrt(6)*n^(3/4)) * (1 + (Pi/72 - 9/(8*Pi)) / sqrt(n)). - Vaclav Kotesovec, Aug 30 2015, extended Jan 18 2017

Extensions

Typo in Maple program fixed by Vaclav Kotesovec, Nov 15 2016

A047261 Numbers that are congruent to {2, 4, 5} mod 6.

Original entry on oeis.org

2, 4, 5, 8, 10, 11, 14, 16, 17, 20, 22, 23, 26, 28, 29, 32, 34, 35, 38, 40, 41, 44, 46, 47, 50, 52, 53, 56, 58, 59, 62, 64, 65, 68, 70, 71, 74, 76, 77, 80, 82, 83, 86, 88, 89, 92, 94, 95, 98, 100, 101, 104, 106, 107, 110, 112, 113, 116, 118, 119, 122, 124
Offset: 1

Views

Author

Keywords

Comments

If B and C are terms in the sequence then 2*B*C is a term. B (resp. C) is a term iff B (resp. C) mod 6 = 2, 4 or 5. It follows that (2*B*C) mod 6 = (2*(B mod 6)*(C mod 6)) mod 6 = 2 or 4 and therefore 2*B*C is a term. Examples: for B=16 and C=29, 2*16*29 = 928 is a term: (2*B*C) mod 6 = (2*16*29) mod 6 = 4; (2*2*2) mod 6 = 2. - Jerzy R Borysowicz, May 24 2018

Crossrefs

Cf. A047242 (complement).

Programs

  • Haskell
    a047261 n = a047261_list !! n
    a047261_list = 2 : 4 : 5 : map (+ 6) a047261_list
    -- Reinhard Zumkeller, Feb 19 2013, Jul 06 2012
    
  • Magma
    [n : n in [0..150] | n mod 6 in [2, 4, 5]]; // Wesley Ivan Hurt, Jun 14 2016
  • Maple
    A047261:=n->(6*n-1-2*cos(2*n*Pi/3))/3: seq(A047261(n), n=1..100); # Wesley Ivan Hurt, Jun 14 2016
  • Mathematica
    CoefficientList[Series[(1 + x)*(x^2 + 2)/((1 + x + x^2)*(x - 1)^2), {x, 0, 50}], x] (* Wesley Ivan Hurt, Aug 16 2014 *)
    Select[ Range@ 125, MemberQ[{2, 4, 5}, Mod[#, 6]] &] (* or *)
    LinearRecurrence[{1, 0, 1, -1}, {2, 4, 5, 8}, 62] (* Robert G. Wilson v, Jun 13 2018 *)

Formula

G.f.: x*(1+x)*(x^2+2) / ((1+x+x^2)*(x-1)^2). - R. J. Mathar, Oct 08 2011
A214090(a(n)) = 1. - Reinhard Zumkeller, Jul 06 2012
From Wesley Ivan Hurt, Jun 14 2016: (Start)
a(n) = a(n-1) + a(n-3) - a(n-4) for n>4.
a(n) = (6*n - 1 - 2*cos(2*n*Pi/3))/3.
a(3k) = 6k-1, a(3k-1) = 6k-2, a(3k-2) = 6k-4. (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/6 - log(2+sqrt(3))/(2*sqrt(3)) + log(2)/3. - Amiram Eldar, Dec 16 2021
E.g.f.: (3 + exp(x)*(6*x - 1) - 2*exp(-x/2)*cos(sqrt(3)*x/2))/3. - Stefano Spezia, Jul 26 2024

Extensions

More terms from Wesley Ivan Hurt, Aug 16 2014

A097451 Number of partitions of n into parts congruent to {2, 3, 4} mod 6.

Original entry on oeis.org

1, 0, 1, 1, 2, 1, 3, 2, 5, 4, 7, 6, 11, 9, 15, 14, 22, 20, 31, 29, 43, 41, 58, 57, 80, 78, 106, 107, 142, 143, 188, 191, 247, 253, 321, 332, 418, 432, 537, 561, 690, 721, 880, 924, 1118, 1178, 1412, 1493, 1781, 1884, 2231, 2370, 2789, 2965, 3472, 3698, 4309, 4596
Offset: 0

Views

Author

Vladeta Jovovic, Aug 23 2004

Keywords

Comments

Number of partitions of n in which no part is 1, no part appears more than twice and no two parts differ by 1. Example: a(6)=3 because we have [6],[4,2] and [3,3]. - Emeric Deutsch, Feb 16 2006
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			a(8)=5 because we have [8],[44],[422],[332] and [2222].
G.f. = 1 + x^2 + x^3 + 2*x^4 + x^5 + 3*x^6 + 2*x^7 + 5*x^8 + 4*x^9 + ...
G.f. = q^7 + q^55 + q^79 + 2*q^103 + q^127 + 3*q^151 + 2*q^175 + 5*q^199 + ...
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, Exercise 7.9.

Crossrefs

Programs

  • Haskell
    a097451 n = p a047228_list n where
       p _  0         = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Nov 16 2012
    
  • Maple
    g:=1/product((1-x^(2+6*j))*(1-x^(3+6*j))*(1-x^(4+6*j)),j=0..15): gser:=series(g,x=0,75): seq(coeff(gser,x,n),n=0..67); # Emeric Deutsch, Feb 16 2006
  • Mathematica
    a[ n_] := SeriesCoefficient[ 1 / Product[ 1 - Boole[ OddQ[ Quotient[ k + 1, 3]]] x^k, {k, n}], {x, 0, n}]; (* Michael Somos, Sep 24 2013 *)
    a[ n_] := SeriesCoefficient[ QPochhammer[ -x^3, x^3] QPochhammer[ x^6] / QPochhammer[ x^2], {x, 0, n}]; (* Michael Somos, Sep 24 2013 *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 1 / prod(k=1, n, 1 - ( (k+1)\3 % 2) * x^k, 1 + x * O(x^n)), n))}; /* Michael Somos, Sep 24 2013 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^6 + A)^2 / (eta(x^2 + A) * eta(x^3 + A)), n))}; /* Michael Somos, Sep 24 2013 */

Formula

Euler transform of period 6 sequence [ 0, 1, 1, 1, 0, 0, ...].
G.f.: 1/Product_{j>=0} ((1-x^(2+6j))(1-x^(3+6j))(1-x^(4+6j))). - Emeric Deutsch, Feb 16 2006
Expansion of psi(x^3) / f(-x^2) in powers of x where psi(), f() are Ramanujan theta functions. - Michael Somos, Sep 24 2013
Expansion of q^(-7/24) * eta(q^6)^2 / (eta(q^2) * eta(q^3)) in powers of q. - Michael Somos, Sep 24 2013
a(n) ~ exp(Pi*sqrt(n/3)) / (4*3^(3/4)*n^(3/4)). - Vaclav Kotesovec, Aug 30 2015
Expansion of f(-x, -x^5) / f(-x, -x^2) in powers of x where f(, ) is Ramanujan's general theta function. - Michael Somos, Oct 06 2015

Extensions

More terms from Emeric Deutsch, Feb 16 2006

A096981 Number of partitions of n into parts congruent to {0, 1, 3, 5} mod 6.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 5, 6, 7, 10, 12, 15, 21, 25, 30, 39, 46, 56, 72, 85, 101, 125, 147, 175, 215, 252, 296, 356, 415, 487, 582, 676, 786, 927, 1072, 1244, 1460, 1682, 1939, 2255, 2588, 2976, 3446, 3942, 4510, 5189, 5916, 6751, 7739, 8797, 9999, 11406, 12927, 14657
Offset: 0

Views

Author

Noureddine Chair, Aug 19 2004

Keywords

Comments

Also, number of partitions of n in which the distinct parts are prime to 3 and the unrestricted parts are multiples of 3.
The inverted graded parafermionic partition function. This g.f. is a generalization of A003105, A006950 and A096938

Examples

			a(11) = 15 because we can write 11 = 10+1 = 8+2+1 = 7+4 = 5+4+2 (parts do not contain multiple of 3) = 9+2 = 8+3 = 7+3+1 = 6+5 = 6+4+1 = 6+3+2 = 5+3+3 = 5+3+2+1 = 4+3+3+1 = 3+3+3+2.
1 + x + x^2 + 2*x^3 + 2*x^4 + 3*x^5 + 5*x^6 + 6*x^7 + 7*x^8 + 10*x^9 + ...
q^-5 + q^19 + q^43 + 2*q^67 + 2*q^91 + 3*q^115 + 5*q^139 + 6*q^163 + 7*q^187 + ...
		

References

  • T. M. Apostol, An Introduction to Analytic Number Theory, Springer-Verlag, NY, 1976

Crossrefs

Programs

  • Haskell
    a096981 = p $ tail a047273_list where
       p _  0         = 1
       p ks'@(k:ks) m = if k > m then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Feb 19 2013
  • Maple
    series(product(1/(1-x^k+x^(2*k)-x^(3*k)+x^(4*k)-x^(5*k)), k=1..150), x=0,100);
  • Mathematica
    CoefficientList[ Series[ Product[ 1/(1 - x^k + x^(2k) - x^(3k) + x^(4k) - x^(5k)), {k, 55}], {x, 0, 53}], x] (* Robert G. Wilson v, Aug 21 2004 *)
    nmax = 100; CoefficientList[Series[x^3*QPochhammer[-1/x^2, x^3] * QPochhammer[-1/x, x^3]/((1 + x)*(1 + x^2) * QPochhammer[x^3, x^3]), {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 31 2015 *)
  • PARI
    {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A) / (eta(x + A) * eta(x^6 + A)), n))} /* Michael Somos, Jun 08 2012 */
    

Formula

Expansion of q^(5/24) * eta(q^2) / (eta(q) * eta(q^6)) in powers of q. - Michael Somos, Jun 08 2012
Euler transform of period 6 sequence [1, 0, 1, 0, 1, 1, ...]. - Vladeta Jovovic, Aug 20 2004
G.f.: 1/product_{k>=1}(1-x^k+x^(2*k)-x^(3*k)+x^(4*k)-x^(5*k)) = Product_{k>=1}(1+x^(3*k-1))(1+x^(3*k-2))/(1-x^(3*k)).
a(n) ~ exp(2*Pi*sqrt(n)/3) / (2*sqrt(6)*n). - Vaclav Kotesovec, Aug 31 2015

Extensions

Better definition from Vladeta Jovovic, Aug 20 2004
More terms from Robert G. Wilson v, Aug 21 2004
Incorrect b-file replaced by Vaclav Kotesovec, Aug 31 2015

A253144 Number of partitions of n into distinct parts congruent to 1, 2, or 4 modulo 6.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 5, 5, 5, 6, 6, 7, 9, 10, 11, 12, 13, 14, 16, 18, 19, 21, 23, 25, 28, 31, 34, 37, 40, 43, 47, 52, 56, 61, 66, 71, 78, 85, 92, 99, 107, 115, 124, 135, 145, 156, 168, 180, 194, 210, 226, 242, 260, 278, 297, 320, 343, 367, 393, 420, 449, 481, 516, 550, 587, 626, 666, 712, 760, 810, 863, 919, 978, 1041, 1110, 1180, 1254, 1333, 1414, 1503, 1598, 1697, 1801
Offset: 0

Views

Author

Jeremy Lovejoy, Mar 23 2015

Keywords

Comments

a(n) is also the number of partitions of n into parts congruent to 1, 7, or 10 modulo 12.
a(n) is also the number of partitions of n into parts that differ by at least 6, where the inequality is strict if the larger part is 0, 3, or 5 modulo 6, with the exception that 6+1 may appear.

Examples

			a(14) = 5, the valid partitions being 14, 13+1, 10+4, 8+4+2, and 7+4+2+1.
		

Crossrefs

Cf. A056970.

Programs

  • Maple
    series(mul((1+x^(6*k+1))*(1+x^(6*k+2))*(1+x^(6*k+4)),k=0..100),x=0,100)

Formula

a(n) ~ exp(sqrt(n/6)*Pi) / (2^(17/12) * 3^(1/4) * n^(3/4)). - Vaclav Kotesovec, May 24 2018

A253136 The number of overpartitions of n into parts congruent to 2, 4, or 5 modulo 6.

Original entry on oeis.org

1, 0, 2, 0, 4, 2, 6, 4, 10, 8, 18, 14, 28, 24, 44, 42, 68, 66, 102, 104, 154, 160, 226, 238, 330, 354, 476, 516, 676, 742, 958, 1056, 1342, 1486, 1862, 2076, 2568, 2872, 3516, 3940, 4782, 5370, 6464, 7268, 8686, 9774, 11606, 13070, 15428, 17380, 20408, 22986
Offset: 0

Views

Author

Jeremy Lovejoy, Mar 23 2015

Keywords

Crossrefs

Cf. A056970.

Programs

  • Maple
    series(mul((1+x^(6*k+2))*(1+x^(6*k+4))*(1+x^(6*k+5))/((1-x^(6*k+2))*(1-x^(6*k+4))*(1-x^(6*k+5))), k=0..100), x=0, 100);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+
         `if`(irem(i, 6) in [2, 4, 5], add(2*b(n-i*j, i-1), j=1..n/i), 0)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..60);  # Alois P. Heinz, Jan 04 2019
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, b[n, i - 1] + If[ MemberQ[{2, 4, 5}, Mod[i, 6]], Sum[2b[n - i j, i-1], {j, 1, n/i}], 0]]];
    a[n_] := b[n, n];
    a /@ Range[0, 60] (* Jean-François Alcover, Nov 23 2020, after Alois P. Heinz *)
    nmax = 60; CoefficientList[Series[Product[(1 + x^(6*k+2)) * (1 + x^(6*k+4)) * (1 + x^(6*k+5)) / ((1 - x^(6*k+2)) * (1 - x^(6*k+4)) * (1 - x^(6*k+5))), {k, 0, nmax/6}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 14 2021 *)

Formula

a(n) ~ Pi^(5/6) * exp(Pi*sqrt(n/2)) / (2^(7/4) * 3^(1/6) * Gamma(1/6) * n^(11/12)). - Vaclav Kotesovec, Jan 14 2021
Showing 1-7 of 7 results.