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

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

Original entry on oeis.org

1, 1, 3, 8, 18, 40, 88, 184, 384, 783, 1573, 3110, 6087, 11745, 22450, 42466, 79597, 147890, 272632, 498696, 905846, 1634270, 2929804, 5220581, 9249440, 16297659, 28567571, 49825296, 86487331, 149438681, 257077485, 440378787, 751313413, 1276765557, 2161511352
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 05 2017

Keywords

Comments

Euler transform of the generalized pentagonal numbers (A001318).

Crossrefs

Programs

  • Mathematica
    nmax = 34; CoefficientList[Series[Product[1/((1 - x^(2 k - 1))^(k (3 k - 1)/2) (1 - x^(2 k))^(k (3 k + 1)/2)), {k, 1, nmax}], {x, 0, nmax}], x]
    a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d Ceiling[d/2] Ceiling[(3 d + 1)/2]/2, {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 34}]

Formula

G.f.: Product_{k>=1} 1/(1 - x^k)^A001318(k).
a(n) ~ exp(Pi * 2^(5/4) / (3*5^(1/4)) * n^(3/4) + 3*Zeta(3) * sqrt(5*n) / (2^(3/2) * Pi^2) + (Pi/48 - 45*Zeta(3)^2 / (8*Pi^5)) * (5*n/2)^(1/4) + 225*Zeta(3)^3 / (8*Pi^8) - 11*Zeta(3) / (64*Pi^2))/ (2^(95/48) * 5^(1/8) * n^(5/8)). - Vaclav Kotesovec, Nov 07 2017

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

Original entry on oeis.org

1, 1, 9, 30, 106, 339, 1106, 3355, 10102, 29358, 83908, 234394, 644286, 1739933, 4631675, 12153197, 31485413, 80576160, 203902261, 510490213, 1265353568, 3106771717, 7559844833, 18239351931, 43650061720, 103657177941, 244346681972, 571930478187, 1329655624297, 3071230379625, 7049750442386, 16085170634548, 36489192684910
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 30 2016

Keywords

Comments

Euler transform of the octagonal numbers (A000567).

Crossrefs

Programs

  • Maple
    with(numtheory):
    a:= proc(n) option remember; `if`(n=0, 1, add(add(
          d^2*(3*d-2), d=divisors(j))*a(n-j), j=1..n)/n)
        end:
    seq(a(n), n=0..35);  # Alois P. Heinz, Dec 02 2016
  • Mathematica
    nmax=32; CoefficientList[Series[Product[1/(1 - x^k)^(k (3 k - 2)), {k, 1, nmax}], {x, 0, nmax}], x]
  • Python
    from sympy import divisors
    from sympy.core.cache import cacheit
    @cacheit
    def a(n): return 1 if n==0 else sum(sum(d**2*(3*d - 2) for d in divisors(j))*a(n - j) for j in range(1, n + 1))//n
    print([a(n) for n in range(51)]) # Indranil Ghosh, Aug 06 2017, after Maple code

Formula

G.f.: Product_{k>=1} 1/(1 - x^k)^(k*(3*k-2)).
a(n) ~ exp(4*Pi*n^(3/4) / (3*5^(1/4)) - 2*Zeta(3) * sqrt(5*n) / Pi^2 - 10*Zeta(3)^2 * (5*n)^(1/4) / Pi^5 - 200*Zeta(3)^3 / (3*Pi^8) - 3*Zeta(3) / (4*Pi^2) - 1/6) * A^2 / (2^(3/2) * 5^(1/12) * Pi^(1/6) * n^(7/12)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Nov 08 2017

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

Original entry on oeis.org

1, 5, 31, 148, 667, 2754, 10823, 40393, 145085, 502780, 1690603, 5530649, 17658430, 55141520, 168751779, 506933980, 1496999360, 4350994324, 12460305177, 35192973824, 98116587875, 270220568883, 735668636567, 1981082952258, 5279879097853, 13933764841202
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 07 2017

Keywords

Crossrefs

Programs

  • Maple
    N:= 50:
    S:= series(mul(1/(1-x^k)^(k*(3*k+2)), k=1..N),x,N+1):
    seq(coeff(S,x,n),n=0..N); # Robert Israel, Nov 07 2017
  • Mathematica
    nmax = 30; CoefficientList[Series[Product[1/(1-x^k)^(k*(3*k+2)), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(4*Pi*n^(3/4) / (3*5^(1/4)) + 2*Zeta(3) * sqrt(5*n) / Pi^2 - 2*5^(5/4) * Zeta(3)^2 * n^(1/4) / Pi^5 + 200*Zeta(3)^3 / (3*Pi^8) - 3*Zeta(3) / (4*Pi^2) + 1/6) * Pi^(1/6) / (A^2 * 2^(3/2) * 5^(1/6) * n^(2/3)), where A is the Glaisher-Kinkelin constant A074962.

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

Original entry on oeis.org

1, 2, 3, 11, 19, 42, 93, 170, 352, 658, 1266, 2351, 4316, 7926, 14146, 25458, 44748, 78687, 136747, 235988, 405139, 689108, 1168260, 1963940, 3289950, 5474700, 9070976, 14954802, 24537752, 40099905, 65225553, 105713691, 170600344, 274367688, 439568770, 701867457
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 06 2017

Keywords

Crossrefs

Programs

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

Formula

a(n) ~ exp(2*Pi * n^(3/4) / (3*5^(1/4)) + 2*Zeta(3) * sqrt(5*n) / Pi^2 + 5^(1/4)*(5*Pi/48 - 20*Zeta(3)^2 / Pi^5) * n^(1/4) + 800 * Zeta(3)^3 / (3*Pi^8) - 73*Zeta(3) / (96*Pi^2) - 1/12) * A / (2^(115/48) * 5^(5/48) * Pi^(1/12) * n^(29/48)), where A is the Glaisher-Kinkelin constant A074962.

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

Original entry on oeis.org

1, 2, 12, 46, 175, 610, 2107, 6918, 22256, 69498, 212649, 636910, 1874470, 5423332, 15457223, 43433088, 120467606, 330077358, 894193347, 2396636236, 6359325300, 16714566278, 43539016461, 112449776138, 288083439729, 732356943548, 1848098069644, 4630892393996
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 08 2018

Keywords

Comments

Euler transform of A006002.

Crossrefs

Programs

  • Mathematica
    nmax = 27; CoefficientList[Series[Product[1/(1 - x^k)^(k (k + 1)^2/2), {k, 1, nmax}], {x, 0, nmax}], x]
    a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d^2 (d + 1)^2/2, {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 27}]

Formula

G.f.: Product_{k>=1} 1/(1 - x^k)^A006002(k).
a(n) ~ exp(5 * (3*Zeta(5))^(1/5) * n^(4/5) / 2^(8/5) + Pi^4 * n^(3/5) / (90 * 2^(1/5) * (3*Zeta(5))^(3/5)) + (Zeta(3) / 2^(9/5) - Pi^8 / (27000 * 2^(4/5) * Zeta(5))) * n^(2/5) / (3*Zeta(5))^(2/5) + (Pi^8 / (12150000 * Zeta(5)) - Zeta(3) / 900) * Pi^4 * n^(1/5) / (2^(2/5) * 3^(1/5) * Zeta(5)^(6/5)) + 1/24 - Zeta(3) / (4*Pi^2) - Pi^16 / (5248800000 * Zeta(5)^3) + Pi^8 * Zeta(3) / (324000 * Zeta(5)^2) - Zeta(3)^2 / (120 * Zeta(5)) + Zeta'(-3)/2) * (3*Zeta(5))^(43/400) / (2^(57/200) * sqrt(5*A*Pi) * n^(243/400)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Apr 08 2018
Showing 1-5 of 5 results.