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-10 of 12 results. Next

A049209 a(n) = -Product_{k=0..n} (7*k-1); sept-factorial numbers.

Original entry on oeis.org

1, 6, 78, 1560, 42120, 1432080, 58715280, 2818333440, 155008339200, 9610517030400, 663125675097600, 50397551307417600, 4182996758515660800, 376469708266409472000, 36517561701841718784000, 3797826416991538753536000, 421558732286060801642496000
Offset: 0

Views

Author

Keywords

Crossrefs

Row sums of triangle A051186 (scaled Stirling1 triangle).
Sequences of the form m^n*Pochhammer((m-1)/m, n): A000007 (m=1), A001147 (m=2), A008544 (m=3), A008545 (m=4), A008546 (m=5), A008543 (m=6), this sequence (m=7), A049210 (m=8), A049211 (m=9), A049212 (m=10), A254322 (m=11), A346896 (m=12).

Programs

  • Magma
    [ -&*[ (7*k-1): k in [0..n-1] ]: n in [1..15] ]; // Klaus Brockhaus, Nov 10 2008
    
  • Mathematica
    CoefficientList[Series[(1-7*x)^(-6/7),{x,0,20}],x] * Range[0,20]! (* Vaclav Kotesovec, Jan 28 2015 *)
    With[{m=7}, Table[m^n*Pochhammer[(m-1)/m, n], {n, 0, 30}]] (* G. C. Greubel, Feb 16 2022 *)
  • Sage
    m=7; [m^n*rising_factorial((m-1)/m, n) for n in (0..30)] # G. C. Greubel, Feb 16 2022

Formula

a(n) = 6*A034833(n) = (7*n-1)*(!^7), n >= 1, a(0) := 1.
a(n) = Product_{k=1..n} (7*k - 1). a(0) = 1; a(n) = (7*n - 1)*a(n-1) for n > 0. - Klaus Brockhaus, Nov 10 2008
G.f.: 1/(1-6*x/(1-7*x/(1-13*x/(1-14*x/(1-20*x/(1-21*x/(1-27*x/(1-28*x/(1-...(continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-1)^n*Sum_{k=0..n} 7^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
a(n) = 7^n * Gamma(n+6/7) / Gamma(6/7). - Vaclav Kotesovec, Jan 28 2015
E.g.f.: (1-7*x)^(-6/7). - Vaclav Kotesovec, Jan 28 2015
From Nikolaos Pantelidis, Dec 19 2020: (Start)
G.f.: 1/G(0) where G(k) = 1 - (14*k+6)*x - 7*(k+1)*(7*k+6)*x^2/G(k+1); (continued fraction).
which starts as 1/(1-6*x-42*x^2/(1-20*x-182*x^2/(1-34*x-420*x^2/(1-48*x-756*x^2/(1-62*x-1190*x^2/(1-... )))))) (Jacobi continued fraction).
G.f.: 1/Q(0) where Q(k) = 1 - (7*k+6)*x/(1 - (7*k+7)*x/Q(k+1) ); (continued fraction). (End)
Sum_{n>=0} 1/a(n) = 1 + (e/7)^(1/7)*(Gamma(6/7) - Gamma(6/7, 1/7)). - Amiram Eldar, Dec 19 2022

A049211 a(n) = Product_{k=1..n} (9*k - 1); 9-factorial numbers.

Original entry on oeis.org

1, 8, 136, 3536, 123760, 5445440, 288608320, 17893715840, 1270453824640, 101636305971200, 9045631231436800, 886471860680806400, 94852489092846284800, 11002888734770169036800, 1375361091846271129600000, 184298386307400331366400000, 26354669241958247385395200000
Offset: 0

Views

Author

Keywords

Crossrefs

Sequences of the form m^n*Pochhammer((m-1)/m, n): A000007 (m=1), A001147 (m=2), A008544 (m=3), A008545 (m=4), A008546 (m=5), A008543 (m=6), A049209 (m=7), A049210 (m=8), this sequence (m=9), A049212 (m=10), A254322 (m=11), A346896 (m=12).

Programs

  • Magma
    m:=9; [Round(m^n*Gamma(n +(m-1)/m)/Gamma((m-1)/m)): n in [0..20]]; // G. C. Greubel, Feb 08 2022
    
  • Mathematica
    CoefficientList[Series[(1-9*x)^(-8/9),{x,0,20}],x] * Range[0,20]! (* Vaclav Kotesovec, Jan 28 2015 *)
  • PARI
    a(n) = prod(k=1, n, 9*k-1); \\ Michel Marcus, Jan 08 2015
    
  • Sage
    m=9; [m^n*rising_factorial((m-1)/m, n) for n in (0..20)] # G. C. Greubel, Feb 08 2022

Formula

a(n) = 8*A035022(n) = (9*n-1)(!^9), n >= 1, a(0) = 1.
a(n) = (-1)^n*Sum_{k=0..n} 9^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
a(n) = 9^n * Gamma(n+8/9) / Gamma(8/9). - Vaclav Kotesovec, Jan 28 2015
E.g.f: (1-9*x)^(-8/9). - Vaclav Kotesovec, Jan 28 2015
From Nikolaos Pantelidis, Dec 09 2020: (Start)
G.f.: 1/(1-8*x-72*x^2/(1-26*x-306*x^2/(1-44*x-702*x^2/(1-62*x-1260*x^2/(1-80*x-1980*x^2/(1-...)))))) (Jacobi continued fraction).
G.f.: 1/(1-8*x/(1-9*x/(1-17*x/(1-18*x/(1-26*x/(1-27*x/(1-35*x/(1-36*x/(1-44*x/(1-45*x/(1-...))))))))))) (Stieltjes continued fraction). (End)
From Nikolaos Pantelidis, Dec 19 2020: (Start)
G.f.: 1/G(0) where G(k) = 1 - (18*k+8)*x - 9*(k+1)*(9*k+8)*x^2/G(k+1) (continued fraction).
G.f.: 1/Q(0) where Q(k) = 1 - x*(9*k+8)/(1 - x*(9*k+9)/Q(k+1) ) (continued fraction). (End)
G.f.: hypergeometric2F0([1, 8/9], [--], 9*x). - G. C. Greubel, Feb 08 2022
Sum_{n>=0} 1/a(n) = 1 + (e/9)^(1/9)*(Gamma(8/9) - Gamma(8/9, 1/9)). - Amiram Eldar, Dec 21 2022

Extensions

a(9) (originally given incorrectly as 1011636305971200) corrected by Peter Bala, Feb 20 2015
a(15)-a(16) from Vincenzo Librandi, Feb 20 2015
a(16) corrected and incorrect MAGMA program removed by Georg Fischer, May 10 2021

A049212 a(n) = -Product_{k=0..n} (10*k - 1); deca-factorial numbers.

Original entry on oeis.org

1, 9, 171, 4959, 193401, 9476649, 559122291, 38579438079, 3047775608241, 271252029133449, 26853950884211451, 2927080646379048159, 348322596919106730921, 44933615002564768288809, 6245772485356502792144451, 930620100318118916029523199, 147968595950580907648694188641
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [Round(10^n*Gamma(n+9/10)/Gamma(9/10)): n in [0..25]]; // G. C. Greubel, Feb 03 2022
    
  • Mathematica
    CoefficientList[Series[(1-10*x)^(-9/10),{x,0,20}],x] * Range[0,20]! (* Vaclav Kotesovec, Jan 28 2015 *)
  • PARI
    a(n) = {-prod(k=0, n, 10*k-1)} \\ Andrew Howroyd, Jan 02 2020
    
  • Sage
    [10^n*rising_factorial(9/10, n) for n in (0..25)] # G. C. Greubel, Feb 03 2022

Formula

a(n) = 9*A035278(n) = (10*n-1)(!^10), n >= 1, a(0) = 1.
a(n) = (-1)^n*Sum_{k=0..n} 10^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
a(n) = 10^n * Gamma(n+9/10) / Gamma(9/10). - Vaclav Kotesovec, Jan 28 2015
E.g.f.: (1-10*x)^(-9/10). - Vaclav Kotesovec, Jan 28 2015
From Nikolaos Pantelidis, Jan 17 2021: (Start)
G.f.: 1/G(0) where G(k) = 1 - (20*k+9)*x - 10*(k+1)*(10*k+9)*x^2/G(k+1) (continued fraction).
G.f.: 1/(1-9*x-90*x^2/(1-29*x-380*x^2/(1-49*x-870*x^2/(1-69*x-1560*x^2/(1-89*x-2450*x^2/(1-...)))))) (Jacobi continued fraction).
G.f.: 1/Q(0) where Q(k) = 1 - x*(10*k+9)/(1 - x*(10*k+10)/Q(k+1)) (continued fraction).
G.f.: 1/(1-9*x/(1-10*x/(1-19*x/(1-20*x/(1-29*x/(1-30*x/(1-39*x/(1-40*x/(1-49*x/(1-50*x/(1-...))))))))))) (Stieltjes continued fraction).
(End)
G.f.: Hypergeometric2F0([1, 9/10], --; 10*x). - G. C. Greubel, Feb 03 2022
Sum_{n>=0} 1/a(n) = 1 + (e/10)^(1/10)*(Gamma(9/10) - Gamma(9/10, 1/10)). - Amiram Eldar, Dec 22 2022

Extensions

Terms a(14) and beyond from Andrew Howroyd, Jan 02 2020

A254322 Expansion of e.g.f.: (1-11*x)^(-10/11).

Original entry on oeis.org

1, 10, 210, 6720, 288960, 15603840, 1014249600, 77082969600, 6706218355200, 657209398809600, 71635824470246400, 8596298936429568000, 1126115160672273408000, 159908352815462823936000, 24465977980765812062208000, 4012420388845593178202112000
Offset: 0

Views

Author

Vaclav Kotesovec, Jan 28 2015

Keywords

Comments

Generally, for k > 1, if e.g.f. = (1-k*x)^(-(k-1)/k) then a(n) ~ n! * k^n / (n^(1/k) * Gamma((k-1)/k)).

Crossrefs

Sequences of the form k^n*Pochhammer((k-1)/k, n): A000007 (k=1), A001147 (k=2), A008544 (k=3), A008545 (k=4), A008546 (k=5), A008543 (k=6), A049209 (k=7), A049210 (k=8), A049211 (k=9), A049212 (k=10), this sequence (k=11), A346896 (k=12).

Programs

  • Magma
    m=11; [Round(m^n*Gamma(n +(m-1)/m)/Gamma((m-1)/m)): n in [0..20]]; // G. C. Greubel, Feb 08 2022
    
  • Mathematica
    CoefficientList[Series[(1-11*x)^(-10/11), {x, 0, 20}], x] * Range[0, 20]!
    FullSimplify[Table[11^n * Gamma[n+10/11] / Gamma[10/11], {n, 0, 18}]]
  • Sage
    m=11; [m^n*rising_factorial((m-1)/m, n) for n in (0..20)] # G. C. Greubel, Feb 08 2022

Formula

D-finite with recurrence: a(0) = 1; a(n) = (11*n-1) * a(n-1) for n > 0. [corrected by Georg Fischer, Dec 23 2019]
a(n) = 11^n * Gamma(n+10/11) / Gamma(10/11).
a(n) ~ n! * 11^n / (n^(1/11) * Gamma(10/11)).
From Nikolaos Pantelidis, Jan 17 2021: (Start)
G.f.: 1/G(0) where G(k) = 1 - (22*k+10)*x - 11*(k+1)*(11*k+10)*x^2/G(k+1) (continued fraction).
G.f.: 1/(1-10*x-110*x^2/(1-32*x-462*x^2/(1-54*x-1056*x^2/(1-76*x-1892*x^2/(1-98*x-2970*x^2/(1-...)))))) (Jacobi continued fraction).
G.f.: 1/Q(0) where Q(k) = 1 - x*(11*k+10)/(1 - x*(11*k+11)/Q(k+1)) (continued fraction).
G.f.: 1/(1-10*x/(1-11*x/(1-21*x/(1-22*x/(1-32*x/(1-33*x/(1-43*x/(1-44*x/(1-54*x/(1-55*x/(1-...))))))))))) (Stieltjes continued fraction).
(End)
G.f.: hypergeometric2F0([1, 10/11], [--], 11*x). - G. C. Greubel, Feb 08 2022
Sum_{n>=0} 1/a(n) = 1 + (e/11)^(1/11)*(Gamma(10/11) - Gamma(10/11, 1/11)). - Amiram Eldar, Dec 22 2022

A034975 One seventh of octo-factorial numbers.

Original entry on oeis.org

1, 15, 345, 10695, 417105, 19603935, 1078216425, 67927634775, 4822862069025, 381006103452975, 33147531000408825, 3149015445038838375, 324348590839000352625, 36002693583129039141375, 4284320536392355657823625, 544108708121829168543600375, 73454675596446937753386050625
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 1 select 1 else (8*n-1)*Self(n-1): n in [1..40]]; // G. C. Greubel, Oct 21 2022
    
  • Mathematica
    Table[8^n*Pochhammer[7/8, n]/7, {n, 40}] (* G. C. Greubel, Oct 21 2022 *)
  • SageMath
    [8^n*rising_factorial(7/8,n)/7 for n in range(1,40)] # G. C. Greubel, Oct 21 2022

Formula

7*a(n) = (8*n-1)!^8 = Product_{j=1..n} (8*j-1) = (8*n)!/((2*n)!*2^(6*n)*3^2*5 * A045755(n)*A007696(n)*A034909(n)*A034911(n)*A034176(n)).
E.g.f.: (-1+(1-8*x)^(-7/8))/7.
G.f.: x/(1-15*x/(1-8*x/(1-23*x/(1-16*x/(1-31*x/(1-24*x/(1-39*x/(1-32*x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
a(n) = (1/7) * 8^n * Pochhammer(n, 7/8). - G. C. Greubel, Oct 21 2022
From Amiram Eldar, Dec 20 2022: (Start)
a(n) = A049210(n)/7.
Sum_{n>=1} 1/a(n) = 7*(e/8)^(1/8)*(Gamma(7/8) - Gamma(7/8, 1/8)). (End)

A346896 Expansion of e.g.f.: (1-12*x)^(-11/12).

Original entry on oeis.org

1, 11, 253, 8855, 416185, 24554915, 1743398965, 144702114095, 13746700839025, 1470896989775675, 175036741783305325, 22929813173612997575, 3278963283826658653225, 508239308993132091249875, 84875964601853059238729125, 15192797663731697603732513375
Offset: 0

Views

Author

Nikolaos Pantelidis, Aug 06 2021

Keywords

Crossrefs

Sequences of the form m^n*Pochhammer((m-1)/m, n): A000007 (m=1), A001147 (m=2), A008544 (m=3), A008545 (m=4), A008546 (m=5), A008543 (m=6), A049209 (m=7), A049210 (m=8), A049211 (m=9), A049212 (m=10), A254322 (m=11), this sequence (m=12).

Programs

  • Magma
    m:=12; [Round(m^n*Gamma(n +(m-1)/m)/Gamma((m-1)/m)): n in [0..20]]; // G. C. Greubel, Feb 16 2022
  • Mathematica
    CoefficientList[Series[(1-12*x)^(-11/12),{x,0,20}], x] * Range[0, 20]!
    FullSimplify[Table[12^n Gamma[n+11/12]/Gamma[11/12],{n,0,15}]] (* Stefano Spezia, Aug 07 2021 *)
  • Sage
    m=12; [m^n*rising_factorial((m-1)/m, n) for n in (0..20)] # G. C. Greubel, Feb 16 2022
    

Formula

G.f.: 1/(1-11*x/(1-12*x/(1-23*x/(1-24*x/(1-35*x/(1-36*x/(1-47*x/(1-48*x/(1-59*x/(1-60*x/(1-...))))))))))) (Stieltjes continued fraction).
G.f.: 1/Q(0) where Q(k) = 1 - x*(12*k+11)/(1 - x*(12*k+12)/Q(k+1) ) (continued fraction).
G.f.: 1/(1-11*x-132*x^2/(1-35*x-552*x^2/(1-59*x-1260*x^2/(1-83*x-2256*x^2/(1-107*x-3540*x^2/(1-...)))))) (Jacobi continued fraction).
G.f.: 1/G(0) where G(k) = 1 - x*(24*k+11) - 12*(k+1)*(12*k+11)*x^2/G(k+1) (continued fraction).
a(n) = 12^n*Gamma(n+11/12)/Gamma(11/12). - Stefano Spezia, Aug 07 2021
Sum_{n>=0} 1/a(n) = 1 + (e/12)^(1/12)*(Gamma(11/12) - Gamma(11/12, 1/12)). - Amiram Eldar, Dec 22 2022

A051187 Generalized Stirling number triangle of the first kind.

Original entry on oeis.org

1, -8, 1, 128, -24, 1, -3072, 704, -48, 1, 98304, -25600, 2240, -80, 1, -3932160, 1122304, -115200, 5440, -120, 1, 188743680, -57802752, 6651904, -376320, 11200, -168, 1, -10569646080, 3425697792, -430309376, 27725824, -1003520, 20608, -224, 1
Offset: 1

Views

Author

Keywords

Comments

T(n,m)= R_n^m(a=0, b=8) in the notation of the given 1962 reference.
T(n,m) is a Jabotinsky matrix, i.e. the monic row polynomials E(n,x) := Sum_{m=1..n} T(n,m)*x^m = Product_{j=0..n-1} (x - 8*j), n >= 1, and E(0,x) := 1 are exponential convolution polynomials (see A039692 for the definition and a Knuth reference).
From Petros Hadjicostas, Jun 07 2020: (Start)
For integers n, m >= 0 and complex numbers a, b (with b <> 0), the numbers R_n^m(a,b) were introduced by Mitrinovic (1961) and further examined by Mitrinovic and Mitrinovic (1962). Such numbers are related to the work of Nörlund (1924).
They are defined via Product_{r=0..n-1} (x - (a + b*r)) = Sum_{m=0..n} R_n^m(a,b)*x^m for n >= 0. As a result, R_n^m(a,b) = R_{n-1}^{m-1}(a,b) - (a + b*(n-1))*R_{n-1}^m(a,b) for n >= m >= 1 with R_1^0(a,b) = a, R_1^1(a,b) = 1, R_n^m(a,b) = 0 for n < m, and R_0^0(a,b) = 1.
With a = 0 and b = 1, we get the Stirling numbers of the first kind S1(n,m) = R_n^m(a=0, b=1) = A048994(n,m).
We have R_n^m(a,b) = Sum_{k=0}^{n-m} (-1)^k * a^k * b^(n-m-k) * binomial(m+k, k) * S1(n, m+k) for n >= m >= 0.
For the current array, T(n,m) = R_n^m(a=0, b=8) but with no zero row or column. (End)

Examples

			Triangle T(n,m) (with rows n >= 1 and columns m = 1..n) begins:
          1;
         -8,         1;
        128,       -24,       1;
      -3072,       704,     -48,       1;
      98304,    -25600,    2240,     -80,     1;
   -3932160,   1122304, -115200,    5440,  -120,    1;
  188743680, -57802752, 6651904, -376320, 11200, -168, 1;
  ...
3rd row o.g.f.: E(3,x) = Product_{j=0..2} (x - 8*j) = 128*x - 24*x^2 + x^3.
		

Crossrefs

First (m=1) column sequence is: A051189(n-1).
Row sums (signed triangle): A049210(n-1)*(-1)^(n-1).
Row sums (unsigned triangle): A045755(n).
The b=1..7 triangles are: A008275 (Stirling1 triangle), A039683, A051141, A051142, A051150, A051151, A051186.

Formula

T(n, m) = T(n-1, m-1) - 8*(n-1)*T(n-1, m) for n >= m >= 1; T(n, m) := 0 for n < m; T(n, 0) := 0 for n >= 1; T(0, 0) = 1.
E.g.f. for the m-th column of the signed triangle: (log(1 + 8*x)/8)^m/m!.
From Petros Hadjicostas, Jun 07 2020: (Start)
T(n,m) = 8^(n-m)*Stirling1(n,m) = 8^(n-m)*A048994(n,m) = 8^(n-m)*A008275(n,m) for n >= m >= 1.
Bivariate e.g.f.-o.g.f.: Sum_{n,m >= 1} T(n,m)*x^n*y^m/n! = exp((y/8)*log(1 + 8*x)) - 1 = (1 + 8*x)^(y/8) - 1. (End)

A303007 Expansion of (1-240*x)^(1/8).

Original entry on oeis.org

1, -30, -3150, -472500, -81506250, -15160162500, -2956231687500, -595469525625000, -122815589660156250, -25791273828632812500, -5493541325498789062500, -1183608449221102734375000, -257434837705589844726562500, -56437637496994696728515625000
Offset: 0

Views

Author

Seiichi Manyama, Jun 15 2018

Keywords

Crossrefs

(1-b*x)^(1/A003557(b)): A002420 (b=4), A004984 (b=8), A004990 (b=9), (-1)^n * A108735 (b=12), A301271 (b=16), (-1)^n * A108733 (b=18), A049393 (b=25), A004996 (b=36), this sequence (b=240), A303055 (b=504), A305886 (b=1728).

Programs

  • Mathematica
    CoefficientList[Series[Surd[1-240x,8],{x,0,20}],x] (* Harvey P. Dale, Aug 29 2024 *)
  • PARI
    N=20; x='x+O('x^N); Vec((1-240*x)^(1/8))

Formula

a(n) = 30^n/n! * Product_{k=0..n-1} (8*k - 1) for n > 0.
a(n) = 15^n * A301271(n).
a(n) ~ -2^(4*n - 3) * 15^n / (Gamma(7/8) * n^(9/8)). - Vaclav Kotesovec, Jun 16 2018
D-finite with recurrence: n*a(n) +30*(-8*n+9)*a(n-1)=0. - R. J. Mathar, Jan 20 2020

A088996 Triangle T(n, k) read by rows: T(n, k) = Sum_{j=0..n} binomial(j, n-k) * |Stirling1(n, n-j)|.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 2, 7, 6, 0, 6, 29, 46, 24, 0, 24, 146, 329, 326, 120, 0, 120, 874, 2521, 3604, 2556, 720, 0, 720, 6084, 21244, 39271, 40564, 22212, 5040, 0, 5040, 48348, 197380, 444849, 598116, 479996, 212976, 40320
Offset: 0

Views

Author

Philippe Deléham, Dec 01 2003, Aug 17 2007

Keywords

Examples

			Triangle begins:
  1;
  0,    1;
  0,    1,     2;
  0,    2,     7,      6;
  0,    6,    29,     46,     24;
  0,   24,   146,    329,    326,    120;
  0,  120,   874,   2521,   3604,   2556,    720;
  0,  720,  6084,  21244,  39271,  40564,  22212,   5040;
  0, 5040, 48348, 197380, 444849, 598116, 479996, 212976, 40320;
  ...
		

Crossrefs

Variant: A059364, diagonals give A000007, A000142, A067318.
Cf. A001147 (row sums), A048994, A084938.

Programs

  • Magma
    A088996:= func< n,k | (&+[(-1)^j*Binomial(j,n-k)*StirlingFirst(n,n-j): j in [0..n]]) >;
    [A088996(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Feb 23 2022
  • Maple
    A059364 := (n, k) -> add(abs(Stirling1(n, n - j))*binomial(j, n - k), j = 0..n);
    seq(seq(A059364(n, k), k = 0..n), n = 0..8);  # Peter Luschny, Aug 27 2025
  • Mathematica
    T[n_, k_]:= T[n, k]= Sum[(-1)^(n-i)*Binomial[i, k] StirlingS1[n+1, n+1-i], {i, 0, n}]; {{1}}~Join~Table[Abs@ T[n, k], {n,0,10}, {k,n+1,0,-1}] (* Michael De Vlieger, Jun 19 2018 *)
  • Sage
    def A088996(n,k): return add((-1)^(n-i)*binomial(i,k)*stirling_number1(n+1,n+1-i) for i in (0..n))
    for n in (0..10): [A088996(n,k) for k in (0..n)]  # Peter Luschny, May 12 2013
    

Formula

T(n, k) given by [0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, ...] where DELTA is the operator defined in A084938. [Original name.]
Sum_{k=0..n} (-1)^k*T(n,k) = (-1)^n.
From Vladeta Jovovic, Dec 15 2004: (Start)
E.g.f.: (1-y-y*x)^(-1/(1+x)).
Sum_{k=0..n} T(n, k)*x^k = Product_{k=1..n} (k*x+k-1). (End)
T(n, k) = n*T(n-1, k-1) + (n-1)*T(n-1, k); T(0, 0) = 1, T(0, k) = 0 if k > 0, T(n, k) = 0 if k < 0. - Philippe Deléham, May 22 2005
Sum_{k=0..n} T(n,k)*x^(n-k) = A019590(n+1), A000012(n), A000142(n), A001147(n), A007559(n), A007696(n), A008548(n), A008542(n), A045754(n), A045755(n) for x = -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, respectively. Sum_{k=0..n} T(n,k)*x^k = A033999(n), A000007(n), A001147(n), A008544(n), A008545(n), A008546(n), A008543(n), A049209(n), A049210(n), A049211(n), A049212(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, respectively. - Philippe Deléham, Aug 10 2007
T(n, k) = Sum_{j=0..n} (-1)^j*binomial(j, n-k)*StirlingS1(n, n-j). - G. C. Greubel, Feb 23 2022

Extensions

New name using a formula of G. C. Greubel by Peter Luschny, Aug 27 2025

A147626 Octo-factorial numbers (5).

Original entry on oeis.org

1, 6, 84, 1848, 55440, 2106720, 96909120, 5233092480, 324451733760, 22711621363200, 1771506466329600, 152349556104345600, 14320858273808486400, 1460727543928465612800, 160680029832131217408000, 18960243520191483654144000, 2388990683544126940422144000
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 1 select 1 else (8*n-10)*Self(n-1): n in [1..40]]; // G. C. Greubel, Oct 21 2022
    
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst,s],{n,5,2*5!,8}];lst
    Table[8^(n-1)*Pochhammer[3/4, n-1], {n,40}] (* G. C. Greubel, Oct 21 2022 *)
  • SageMath
    [8^(n-1)*rising_factorial(3/4, n-1) for n in range(1,40)] # G. C. Greubel, Oct 21 2022

Formula

a(n+1) = Sum_{k=0..n} A132393(n,k)*6^k*8^(n-k). - Philippe Deléham, Nov 09 2008
a(n) = (-2)^n*Sum_{k=0..n} 4^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 2*x/G(0), where G(k) = 1 + 1/(1 - 2*x*(8*k+6)/(2*x*(8*k+6) - 1 + 16*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
From G. C. Greubel, Oct 21 2022: (Start)
a(n) = 8^n * Pochhammer(n, 3/4) = -2^(3*n-1) * Pochhammer(n, -1/4).
a(n) = (8*n - 10)*a(n-1). (End)
Sum_{n>=1} 1/a(n) = 1 + (e/8^2)^(1/8)*(Gamma(3/4) - Gamma(3/4, 1/8)). - Amiram Eldar, Dec 20 2022
Showing 1-10 of 12 results. Next