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

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

Original entry on oeis.org

1, 1, 0, 4, 4, 9, 15, 22, 52, 65, 129, 190, 335, 534, 814, 1399, 2074, 3462, 5135, 8303, 12658, 19562, 30182, 45542, 70620, 105034, 161223, 239532, 362929, 539252, 805320, 1197589, 1769483, 2624604, 3847755, 5681787, 8291848, 12165978, 17696362, 25796820
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 08 2017

Keywords

Comments

In general, if g.f. = Product_{k>=1} (1 + x^(2*k-1))^(c2*k^2 + c1*k + c0) and c2>0, then a(n) ~ exp(Pi*sqrt(2)/3 * (7*c2/15)^(1/4) * n^(3/4) + 3*(c1+c2) * Zeta(3) / (2*Pi^2) * sqrt(15*n/(7*c2)) + (Pi*(4*c0 + 2*c1 + c2) * (15/(7*c2))^(1/4) / (24*sqrt(2)) - 9*(c1+c2)^2 * Zeta(3)^2 * (15/(7*c2))^(5/4) / (2^(3/2) * Pi^5)) * n^(1/4) + 2025*(c1+c2)^3 * Zeta(3)^3 / (49 * c2^2 * Pi^8) - 15*(c1+c2) * (4*c0 + 2*c1 + c2) * Zeta(3) / (112 * c2 * Pi^2)) * (7/15)^(1/8) * 2^((c1+c2)/24 - 9/4) * c2^(1/8) / n^(5/8).

Crossrefs

Programs

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

Formula

a(n) ~ exp(Pi/3 * (7/15)^(1/4) * sqrt(2) * n^(3/4) + 3*Zeta(3) * sqrt(15*n/7) / (2*Pi^2) + (Pi * (15/7)^(1/4) / (24*sqrt(2)) - 9*Zeta(3)^2 * (15/7)^(5/4) / (2^(3/2) * Pi^5)) * n^(1/4) + 2025*Zeta(3)^3 / (49*Pi^8) - 15*Zeta(3) / (112*Pi^2)) * (7/15)^(1/8) / (2^(53/24) * n^(5/8)).

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

Original entry on oeis.org

1, -1, 5, -14, 36, -97, 246, -593, 1423, -3351, 7699, -17432, 38901, -85545, 185862, -399220, 848080, -1783682, 3716584, -7675916, 15722127, -31951330, 64452707, -129102947, 256876062, -507854808, 997954125, -1949631802, 3787674152, -7319306458, 14071371173
Offset: 0

Views

Author

Seiichi Manyama, Apr 09 2019

Keywords

Comments

This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = (-1)^(n+1) * n^2, g(n) = -1.

Crossrefs

Product_{k>=1} (1+x^k)^((-1)^k*k^b): A083365 (b=0), A284474 (b=1), this sequence (b=2).

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[Product[(1 + x^k)^((-1)^k*k^2), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 09 2019 *)
    nmax = 40; CoefficientList[Series[Product[(1 + x^(2*k))^(4*k^2) / (1 + x^(2*k - 1))^((2*k - 1)^2), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 09 2019 *)
  • PARI
    N=66; x='x+O('x^N); Vec(prod(k=1, N, (1+x^k)^((-1)^k*k^2)))

Formula

a(n) ~ (-1)^n * exp(2*Pi*n^(3/4)/3 + 3*Zeta(3)/(4*Pi^2)) / (4*n^(5/8)). - Vaclav Kotesovec, Apr 09 2019

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

Original entry on oeis.org

1, 2, 2, 10, 18, 36, 86, 150, 326, 608, 1164, 2230, 4046, 7632, 13622, 24868, 44222, 78304, 138312, 240138, 418648, 718292, 1233494, 2097350, 3552370, 5987642, 10026088, 16745600, 27779030, 45970868, 75650248, 124100970, 202720814, 329909400, 535132036
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 08 2017

Keywords

Comments

Convolution of A294749 and A294750.
In general, if g.f. = Product_{k>=1} ((1 + x^(2*k-1))/(1 - x^(2*k-1)))^(c2*k^2 + c1*k + c0) and c2 > 0, then a(n) ~ exp(Pi*sqrt(2) * c2^(1/4) * n^(3/4) / 3 + 7*(c1+c2) * Zeta(3) * sqrt(n) / (2*sqrt(c2) * Pi^2) + (Pi*(4*c0 + 2*c1 + c2) / (8*sqrt(2) * c2^(1/4)) - 49*(c1+c2)^2 * Zeta(3)^2 / (2^(3/2) * c2^(5/4) * Pi^5)) * n^(1/4) - (7*c0 + 21*c1/4 + c2 + 7*c0*c1/c2 + 7*c1^2/(2*c2)) * Zeta(3) / (4*Pi^2) + 22411*(c1+c2)^3 * Zeta(3)^3 / (196 * c2^2 * Pi^8) - (c1+c2)/24) * A^((c1+c2)/2) * (n^((c1+c2)/96 - 5/8) / (2^(c0/2 + (11*c1 + 5*c2)/48 + 9/4) * Pi^((c1+c2)/24) * c2^((c1+c2)/96 - 1/8))), where A is the Glaisher-Kinkelin constant A074962.

Crossrefs

Programs

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

Formula

a(n) ~ exp(sqrt(2)*Pi * n^(3/4)/3 + 7*Zeta(3) * sqrt(n) / (2*Pi^2) + (Pi / (8*sqrt(2)) - 49*Zeta(3)^2 / (2^(3/2) * Pi^5)) * n^(1/4) + 22411*Zeta(3)^3 / (196*Pi^8) - Zeta(3)/(4*Pi^2) - 1/24) * sqrt(A) / (2^(113/48) * Pi^(1/24) * n^(59/96)), where A is the Glaisher-Kinkelin constant A074962.

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

Original entry on oeis.org

1, 0, 0, 1, 0, 3, 1, 6, 3, 11, 12, 18, 29, 33, 69, 67, 138, 141, 275, 306, 516, 656, 972, 1353, 1828, 2712, 3477, 5280, 6654, 10038, 12756, 18789, 24369, 34796, 46167, 63990, 86629, 117189, 160698, 213984, 295092, 389517, 536683, 706590, 968289, 1276310
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 08 2017

Keywords

Crossrefs

Programs

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

Formula

a(n) ~ exp(2*Pi * n^(3/4) / (3^(5/4) * 5^(1/4)) - 5^(1/4) * Pi * n^(1/4) / (16*3^(3/4)) + 3*Zeta(3) / (32*Pi^2)) / (2^(31/16) * 15^(1/8) * n^(5/8)).
Showing 1-4 of 4 results.