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.

A349331 G.f. A(x) satisfies A(x) = 1 + x * A(x)^4 / (1 - x).

Original entry on oeis.org

1, 1, 5, 31, 219, 1678, 13570, 114014, 985542, 8708099, 78298727, 714105907, 6590200215, 61427125994, 577456943614, 5468604044500, 52122539760992, 499613409224137, 4813105582181533, 46576519080852235, 452545041339982871, 4413071971740021275, 43177663974461532959
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 15 2021

Keywords

Crossrefs

Programs

  • Maple
    a:= n-> coeff(series(RootOf(1+x*A^4/(1-x)-A, A), x, n+1), x, n):
    seq(a(n), n=0..22);  # Alois P. Heinz, Nov 15 2021
  • Mathematica
    nmax = 22; A[] = 0; Do[A[x] = 1 + x A[x]^4/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n - 1, k - 1] Binomial[4 k, k]/(3 k + 1), {k, 0, n}], {n, 0, 22}]
  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
    A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^4, k)) )); A[n+1]}
    for(n=0, 30, print1(a(n), ", ")) \\ Vaclav Kotesovec, Nov 23 2024, after Paul D. Hanna

Formula

a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(4*k,k) / (3*k+1).
a(n) ~ 283^(n + 1/2) / (2^(7/2) * sqrt(Pi) * n^(3/2) * 3^(3*n + 3/2)). - Vaclav Kotesovec, Nov 15 2021

A349332 G.f. A(x) satisfies A(x) = 1 + x * A(x)^5 / (1 - x).

Original entry on oeis.org

1, 1, 6, 46, 406, 3901, 39627, 418592, 4551672, 50610692, 572807157, 6577068383, 76426719408, 897078662538, 10620634999318, 126676885170703, 1520759193166329, 18361269213121164, 222814883564042704, 2716125963857227904, 33244557641365865109
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 15 2021

Keywords

Crossrefs

Programs

  • Maple
    a:= n-> coeff(series(RootOf(1+x*A^5/(1-x)-A, A), x, n+1), x, n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Nov 15 2021
  • Mathematica
    nmax = 20; A[] = 0; Do[A[x] = 1 + x A[x]^5/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n - 1, k - 1] Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 20}]
  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
    A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^5, k)) )); A[n+1]}
    for(n=0, 30, print1(a(n), ", ")) \\ Vaclav Kotesovec, Nov 23 2024, after Paul D. Hanna

Formula

a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(5*k,k) / (4*k+1).
a(n) ~ 3381^(n + 1/2) / (25 * sqrt(Pi) * n^(3/2) * 2^(8*n + 7/2)). - Vaclav Kotesovec, Nov 15 2021
Recurrence: 8*n*(2*n - 1)*(4*n - 1)*(4*n + 1)*a(n) = (4405*n^4 - 10346*n^3 + 9575*n^2 - 4354*n + 840)*a(n-1) - 12*(n-2)*(1255*n^3 - 3957*n^2 + 4492*n - 1820)*a(n-2) + 2*(n-3)*(n-2)*(10655*n^2 - 32733*n + 26908)*a(n-3) - 4*(n-4)*(n-3)*(n-2)*(3445*n - 6986)*a(n-4) + 3381*(n-5)*(n-4)*(n-3)*(n-2)*a(n-5). - Vaclav Kotesovec, Nov 17 2021

A349333 G.f. A(x) satisfies A(x) = 1 + x * A(x)^6 / (1 - x).

Original entry on oeis.org

1, 1, 7, 64, 678, 7836, 95838, 1219527, 15979551, 214151601, 2921712145, 40444378948, 566634504256, 8019501351103, 114484746457075, 1646614155398872, 23837794992712680, 347081039681365623, 5079306905986689309, 74670702678690897079, 1102218694940440851877
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 15 2021

Keywords

Crossrefs

Programs

  • Maple
    a:= n-> coeff(series(RootOf(1+x*A^6/(1-x)-A, A), x, n+1), x, n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Nov 15 2021
  • Mathematica
    nmax = 20; A[] = 0; Do[A[x] = 1 + x A[x]^6/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n - 1, k - 1] Binomial[6 k, k]/(5 k + 1), {k, 0, n}], {n, 0, 20}]
  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
    A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^6, k)) )); A[n+1]}
    for(n=0, 30, print1(a(n), ", ")) \\ Vaclav Kotesovec, Nov 23 2024, after Paul D. Hanna

Formula

a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(6*k,k) / (5*k+1).
a(n) ~ 49781^(n + 1/2) / (72 * sqrt(3*Pi) * n^(3/2) * 5^(5*n + 3/2)). - Vaclav Kotesovec, Nov 15 2021

A349335 G.f. A(x) satisfies A(x) = 1 + x * A(x)^8 / (1 - x).

Original entry on oeis.org

1, 1, 9, 109, 1541, 23823, 390135, 6651051, 116798643, 2098313686, 38382509118, 712447023590, 13385500614902, 254065657922154, 4864482597112186, 93840443376075810, 1822169236520766546, 35586928273002974487, 698572561837366684479, 13775697096997873764647
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 15 2021

Keywords

Comments

In general, for m > 1, Sum_{k=0..n} binomial(n-1,k-1) * binomial(m*k,k) / ((m-1)*k+1) ~ (m-1)^(m/2 - 2) * (1 + m^m/(m-1)^(m-1))^(n + 1/2) / (sqrt(2*Pi) * m^((m-1)/2) * n^(3/2)). - Vaclav Kotesovec, Nov 15 2021

Crossrefs

Programs

  • Maple
    a:= n-> coeff(series(RootOf(1+x*A^8/(1-x)-A, A), x, n+1), x, n):
    seq(a(n), n=0..19);  # Alois P. Heinz, Nov 15 2021
  • Mathematica
    nmax = 19; A[] = 0; Do[A[x] = 1 + x A[x]^8/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[Binomial[n - 1, k - 1] Binomial[8 k, k]/(7 k + 1), {k, 0, n}], {n, 0, 19}]
  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
    A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^8, k)) )); A[n+1]}
    for(n=0, 30, print1(a(n), ", ")) \\ Vaclav Kotesovec, Nov 23 2024, after Paul D. Hanna

Formula

a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(8*k,k) / (7*k+1).
a(n) ~ 17600759^(n + 1/2) / (2048 * sqrt(Pi) * n^(3/2) * 7^(7*n + 3/2)). - Vaclav Kotesovec, Nov 15 2021

A349363 G.f. A(x) satisfies: A(x) = 1 + x * A(x)^7 / (1 + x).

Original entry on oeis.org

1, 1, 6, 57, 629, 7589, 96942, 1288729, 17643920, 247089010, 3522891561, 50964747400, 746241617226, 11038241689188, 164696773030055, 2475832560808858, 37462189433509758, 570112127356828846, 8720472842436039280, 133997057207982607092, 2067402314984991892461
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 15 2021

Keywords

Crossrefs

Programs

  • Maple
    a:= n-> coeff(series(RootOf(1+x*A^7/(1+x)-A, A), x, n+1), x, n):
    seq(a(n), n=0..20);  # Alois P. Heinz, Nov 15 2021
  • Mathematica
    nmax = 20; A[] = 0; Do[A[x] = 1 + x A[x]^7/(1 + x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    Table[Sum[(-1)^(n - k) Binomial[n - 1, k - 1] Binomial[7 k, k]/(6 k + 1), {k, 0, n}], {n, 0, 20}]

Formula

a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(n-1,k-1) * binomial(7*k,k) / (6*k+1).
a(n) = (-1)^(n+1)* F([8/7, 9/7, 10/7, 11/7, 12/7, 13/7, 1-n], [4/3, 3/2, 5/3, 11/6, 2, 13/6], 7^7/6^6), where F is the generalized hypergeometric function. - Stefano Spezia, Nov 15 2021
a(n) ~ 776887^(n + 1/2) / (343 * sqrt(Pi) * n^(3/2) * 2^(6*n + 2) * 3^(6*n + 3/2)). - Vaclav Kotesovec, Nov 17 2021

A378326 a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(n*k,k) / ((n-1)*k+1).

Original entry on oeis.org

1, 1, 3, 19, 219, 3901, 95838, 3022909, 116798643, 5350403737, 283728025998, 17104314563843, 1155635807408096, 86513627563199279, 7109252862969177287, 636268582522962837475, 61610670571434193189443, 6418044336586421956746033, 715718717341021991299583730
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 23 2024

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[n-1, k-1]*Binomial[n*k, k]/((n-1)*k+1), {k, 0, n}], {n, 0, 20}]

Formula

a(n) ~ exp(n + exp(-1) - 1/2) * n^(n - 5/2) / sqrt(2*Pi).

A378325 G.f. A(x) = Sum_{n>=0} a(n)*x^n, where a(n) = Sum_{k=0..n-1} [x^k] A(x)^k for n >= 1 with a(0) = 1.

Original entry on oeis.org

1, 1, 2, 7, 41, 338, 3499, 42969, 606351, 9633640, 169888025, 3290314970, 69409429043, 1584105116525, 38894316619948, 1022411500472240, 28653072049382809, 852911635849385778, 26876978490909421289, 893929164892155754432, 31296785296935394097351, 1150551256823546563078988
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 23 2024

Keywords

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
    A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^k, k)) )); A[n+1]}
    for(n=0, 30, print1(a(n), ", ")) \\ after Paul D. Hanna

Formula

a(n) ~ c * n! / (n^alpha * LambertW(1)^n), where alpha = 2 - 2*LambertW(1) - 1/(1 + LambertW(1)) = 0.22760967581532... and c = 0.323194722450152336...
Showing 1-7 of 7 results.