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.

A094418 Generalized ordered Bell numbers Bo(5,n).

Original entry on oeis.org

1, 5, 55, 905, 19855, 544505, 17919055, 687978905, 30187495855, 1490155456505, 81732269223055, 4931150091426905, 324557348772511855, 23141780973332248505, 1776997406800302687055, 146197529083891406394905, 12829862285488250150167855, 1196280147496701351115120505
Offset: 0

Views

Author

Ralf Stephan, May 02 2004

Keywords

Comments

Fifth row of array A094416, which has more information.

Crossrefs

Programs

  • Magma
    A094416:= func< n,k | (&+[Factorial(j)*n^j*StirlingSecond(k,j): j in [0..k]]) >;
    A094418:= func< k | A094416(5,k) >;
    [A094418(n): n in [0..30]]; // G. C. Greubel, Jan 12 2024
    
  • Mathematica
    t = 30; Range[0, t]! CoefficientList[Series[1/(6 - 5 Exp[x]), {x, 0, t}], x] (* Vincenzo Librandi, Mar 16 2014 *)
  • PARI
    my(N=25,x='x+O('x^N)); Vec(serlaplace(1/(6 - 5*exp(x)))) \\ Joerg Arndt, Jan 15 2024
  • SageMath
    def A094416(n,k): return sum(factorial(j)*n^j*stirling_number2(k,j) for j in range(k+1)) # array
    def A094418(k): return A094416(5,k)
    [A094418(n) for n in range(31)] # G. C. Greubel, Jan 12 2024
    

Formula

E.g.f.: 1/(6 - 5*exp(x)).
a(n) = Sum_{k=0..n} A131689(n,k) * 5^k. - Philippe Deléham, Nov 03 2008
a(n) ~ n! / (6*(log(6/5))^(n+1)). - Vaclav Kotesovec, Mar 14 2014
a(0) = 1; a(n) = 5 * Sum_{k=1..n} binomial(n,k) * a(n-k). - Ilya Gutkovskiy, Jan 17 2020
a(0) = 1; a(n) = 5 * a(n-1) - 6 * Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023
From Seiichi Manyama, Jun 01 2025: (Start)
a(n) = (-1)^(n+1)/6 * Li_{-n}(6/5), where Li_{n}(x) is the polylogarithm function.
a(n) = (1/6) * Sum_{k>=0} k^n * (5/6)^k.
a(n) = (5/6) * Sum_{k=0..n} 6^k * (-1)^(n-k) * A131689(n,k) for n > 0. (End)

A365569 Expansion of e.g.f. 1 / (6 - 5 * exp(x))^(3/5).

Original entry on oeis.org

1, 3, 27, 387, 7659, 193491, 5948091, 215446563, 8984708235, 423944899443, 22328393101659, 1298429924941251, 82625791930962219, 5711012035686681363, 426058604580805219323, 34121803137713388036963, 2919847869159667841599947, 265868538017899566748612275
Offset: 0

Views

Author

Seiichi Manyama, Sep 09 2023

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := Sum[Product[5*j + 3, {j, 0, k - 1}] * StirlingS2[n, k], {k, 0, n}]; Array[a, 18, 0] (* Amiram Eldar, Sep 11 2023 *)
    With[{nn=20},CoefficientList[Series[1/(6-5*Exp[x])^(3/5),{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Nov 03 2024 *)
  • PARI
    a(n) = sum(k=0, n, prod(j=0, k-1, 5*j+3)*stirling(n, k, 2));

Formula

a(n) = Sum_{k=0..n} (Product_{j=0..k-1} (5*j+3)) * Stirling2(n,k).
a(0) = 1; a(n) = Sum_{k=1..n} (5 - 2*k/n) * binomial(n,k) * a(n-k).
a(n) ~ sqrt(2*Pi) * n^(n + 1/10) / (6^(3/5) * Gamma(3/5) * exp(n) * log(6/5)^(n + 3/5)). - Vaclav Kotesovec, Nov 11 2023
a(0) = 1; a(n) = 3*a(n-1) - 6*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023

A365570 Expansion of e.g.f. 1 / (6 - 5 * exp(x))^(4/5).

Original entry on oeis.org

1, 4, 40, 616, 12856, 338728, 10781176, 402250216, 17213590840, 831013114792, 44675458306168, 2646758624166760, 171319908334752184, 12028779733435667752, 910538645035885918456, 73918475291961325824232, 6406179168820339231897144
Offset: 0

Views

Author

Seiichi Manyama, Sep 09 2023

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := Sum[Product[5*j + 4, {j, 0, k - 1}] * StirlingS2[n, k], {k, 0, n}]; Array[a, 17, 0] (* Amiram Eldar, Sep 11 2023 *)
  • PARI
    a(n) = sum(k=0, n, prod(j=0, k-1, 5*j+4)*stirling(n, k, 2));

Formula

a(n) = Sum_{k=0..n} (Product_{j=0..k-1} (5*j+4)) * Stirling2(n,k).
a(0) = 1; a(n) = Sum_{k=1..n} (5 - k/n) * binomial(n,k) * a(n-k).
a(n) ~ sqrt(2*Pi) * n^(n + 3/10) / (6^(4/5) * Gamma(4/5) * exp(n) * log(6/5)^(n + 4/5)). - Vaclav Kotesovec, Nov 11 2023
a(0) = 1; a(n) = 4*a(n-1) - 6*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023

A365585 Expansion of e.g.f. 1 / (1 + 5 * log(1-x))^(2/5).

Original entry on oeis.org

1, 2, 16, 214, 4030, 98020, 2923580, 103306320, 4219788720, 195631761360, 10148327972160, 582405469831920, 36635844203963760, 2506613821744700640, 185327181909308762400, 14724431257109269113600, 1251088847268683450630400, 113202071235423519573369600
Offset: 0

Views

Author

Seiichi Manyama, Sep 10 2023

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := Sum[Product[5*j + 2, {j, 0, k - 1}] * Abs[StirlingS1[n, k]], {k, 0, n}]; Array[a, 18, 0] (* Amiram Eldar, Sep 10 2023 *)
  • PARI
    a(n) = sum(k=0, n, prod(j=0, k-1, 5*j+2)*abs(stirling(n, k, 1)));

Formula

a(n) = Sum_{k=0..n} (Product_{j=0..k-1} (5*j+2)) * |Stirling1(n,k)|.
a(0) = 1; a(n) = Sum_{k=1..n} (5 - 3*k/n) * (k-1)! * binomial(n,k) * a(n-k).

A367374 Expansion of the e.g.f. (exp(x) / (5 - 4*exp(x)))^(2/5).

Original entry on oeis.org

1, 2, 12, 128, 1944, 38264, 924936, 26507672, 878565000, 33058419032, 1392125985864, 64864749910424, 3313075222410504, 184071465908101592, 11051901784679926728, 713107430713993422872, 49208366812318404125832, 3616200105869781814285400
Offset: 0

Views

Author

Seiichi Manyama, Nov 15 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n, (-1)^(n-k)*prod(j=0, k-1, 5*j+2)*stirling(n, k, 2));

Formula

a(n) = Sum_{k=0..n} (-1)^(n-k) * (Product_{j=0..k-1} (5*j+2)) * Stirling2(n,k).
a(0) = 1; a(n) = Sum_{k=1..n} (-1)^k * (3*k/n - 5) * binomial(n,k) * a(n-k).
a(0) = 1; a(n) = 2*a(n-1) + 4*Sum_{k=1..n-1} binomial(n-1,k) * a(n-k).
Showing 1-5 of 5 results.