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 14 results. Next

A062199 Second (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 12, 126, 1344, 15120, 181440, 2328480, 31933440, 467026560, 7264857600, 119870150400, 2092278988800, 38532804710400, 746943599001600, 15205637551104000, 324386934423552000, 7237883474325504000, 168600109166641152000, 4093235983656787968000
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

a(n) is the total number of ascending runs of length 5 over all permutations of {1,2,...,n+5}. a(1) = 12 because we have: [1,2,3,4,6,5], [1,2,3,5,6,4], [1,2,4,5,6,3], [1,3,4,5,6,2], [2,1,3,4,5,6], [2,3,4,5,6,1], [3,1,2,4,5,6], [4,1,2,3,5,6], [5,1,2,3,4,6], [6,1,2,3,4,5], and [1,2,3,4,5,6] which has two runs of length 5. - Geoffrey Critzer, Feb 21 2014

Crossrefs

Cf. A001720 (first column of A062140), A264781.

Programs

  • Magma
    [Binomial(n, 5)*Factorial(n-4): n in [5..25]]; // Vincenzo Librandi, Feb 23 2014
    
  • Mathematica
    Table[Sum[n!/5!, {i, 5, n}], {n, 5, 21}] (* Zerinvary Lajos, Jul 12 2009 *)
    With[{nn=20},CoefficientList[Series[(1+5x)/(1-x)^7,{x,0,nn}],x] Range[ 0,nn]!] (* Harvey P. Dale, Nov 10 2016 *)
  • PARI
    x='x+O('x^30); Vec(serlaplace((1+5*x)/(1-x)^7)) \\ G. C. Greubel, Feb 07 2018
  • Sage
    [binomial(n,5)*factorial (n-4) for n in range(5, 22)] # Zerinvary Lajos, Jul 07 2009
    

Formula

E.g.f.: (1+5*x)/(1-x)^7.
a(n) = A062140(n+1, 1) = (n+1)!*binomial(n+5, 5).
If we define f(n,i,x)= Sum(Sum(binomial(k,j)*Stirling1(n,k)*Stirling2(j,i)*x^(k-j),j=i..k),k=i..n) then a(n-1)=(-1)^(n-1)*f(n,1,-6), (n>=1). [Milan Janjic, Mar 01 2009]
a(n) = Sum_{k>0} k * A264781(n+5,k). - Alois P. Heinz, Nov 24 2015
Assuming offset 1: a(n) = -n!*binomial(-n,5). - Peter Luschny, Apr 29 2016
From Amiram Eldar, Sep 24 2022: (Start)
Sum_{n>=0} 1/a(n) = 1565/12 - 50*e - 5*gamma + 5*Ei(1), where gamma is Euler's constant (A001620) and Ei(1) is the exponential integral at 1 (A091725).
Sum_{n>=0} (-1)^n/a(n) = -125/12 + 20/e + 5*gamma - 5*Ei(-1), where -Ei(-1) is the negated exponential integral at -1 (A099285). (End)

Extensions

More terms from Vincenzo Librandi, Feb 23 2014

A062260 Third (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 21, 336, 5040, 75600, 1164240, 18627840, 311351040, 5448643200, 99891792000, 1917922406400, 38532804710400, 809188898918400, 17739910476288000, 405483668029440000, 9650511299100672000
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Crossrefs

Programs

  • Magma
    [Factorial(n+2)*Binomial(n+6, 6)/2: n in [0..30]]; // G. C. Greubel, Feb 06 2018
  • Maple
    a:=n->sum((n-j)*n!/6!, j=5..n): seq(a(n), n=6..21); # Zerinvary Lajos, Apr 29 2007
  • Mathematica
    Table[(n + 2)! Binomial[n + 6, 6]/2, {n, 0, 20}] (* Wesley Ivan Hurt, Jan 23 2017 *)
  • PARI
    { f=1; for (n=0, 100, f*=n + 2; write("b062260.txt", n, " ", f*binomial(n + 6, 6)/2) ) } \\ Harry J. Smith, Aug 03 2009
    
  • Sage
    [binomial(n,6)*factorial (n-4)/2 for n in range(6, 22)] # Zerinvary Lajos, Jul 07 2009
    

Formula

E.g.f.: (1+12*x+15*x^2)/(1-x)^9.
a(n) = A062140(n+2, 2) = (n+2)!*binomial(n+6, 6)/2!.
If we define f(n,i,x) = Sum_{k=1..n} Sum_{j=i..k} binomial(k,j) * Stirling1(n,k) * Stirling2(j,i) * x^(k-j) then a(n-2) = (-1)^n * f(n,2,-7), (n>=2). - Milan Janjic, Mar 01 2009
a(n) = binomial(n,6)*(n-4)!/2, n >= 6. - Zerinvary Lajos, Jul 07 2009

A062266 Row sums of unsigned triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 6, 43, 358, 3393, 36046, 424051, 5470158, 76751233, 1163391958, 18941512731, 329604456886, 6103575192193, 119823200043678, 2485452283923043, 54309931242376606, 1246803623807490561, 29999359707124127398, 754865494585690965643, 19824604328577866107398
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Crossrefs

Programs

  • Magma
    [(&+[Factorial(n)*Binomial(n+4,n-m)/Factorial(m): m in [0..n]]): n in [0..20]]; // G. C. Greubel, Feb 06 2018
    
  • Maple
    A062266:= n -> n!*simplify(LaguerreL(n,4,-1), 'LaguerreL');
    seq(A062266(n), n = 0 .. 30); # G. C. Greubel, Mar 10 2021
  • Mathematica
    Table[n!*SeriesCoefficient[E^(x/(1-x))/(1-x)^5,{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 11 2012 *)
    With[{nn=20},CoefficientList[Series[Exp[x/(1-x)]/(1-x)^5,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Jul 11 2025 *)
  • PARI
    my(x='x+O('x^66)); Vec(serlaplace(exp(x/(1-x))/(1-x)^5)) \\ Joerg Arndt, May 06 2013
    
  • PARI
    a(n) = vecsum(apply(abs,Vec(n!*pollaguerre(n, 4)))); \\ Michel Marcus, Feb 06 2021
    
  • Sage
    [factorial(n)*gen_laguerre(n, 4, -1) for n in (0..30)] # G. C. Greubel, Mar 10 2021

Formula

E.g.f.: exp(x/(1-x))/(1-x)^5.
a(n) = Sum_{m=0..n} n!*binomial(n+4, n-m)/m!.
a(n) = 2*(n+2)*a(n-1) - (n-1)*(n+3)*a(n-2). - Vaclav Kotesovec, Oct 11 2012
a(n) ~ exp(2*sqrt(n)-n-1/2)*n^(n+9/4)/sqrt(2). - Vaclav Kotesovec, Oct 11 2012
a(n) = n!*LaguerreL(n, 4, -1). - G. C. Greubel, Mar 10 2021

A062261 Fourth (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 32, 720, 14400, 277200, 5322240, 103783680, 2075673600, 42810768000, 913296384000, 20183850086400, 462393656524800, 10981849342464000, 270322445352960000, 6893222356500480000, 181981070211612672000
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Crossrefs

Programs

  • Magma
    [Factorial(n+3)*Binomial(n+7,7)/6: n in [0..30]]; // G. C. Greubel, May 13 2018
  • Mathematica
    Table[(n+3)!*Binomial[n+7, 7]/3!, {n, 0, 30}] (* G. C. Greubel, May 13 2018 *)
  • PARI
    { f=2; for (n=0, 100, f*=n + 3; write("b062261.txt", n, " ", f*binomial(n + 7, 7)/6) ) } \\ Harry J. Smith, Aug 03 2009
    

Formula

E.g.f.: (1+21*x+63*x^2+35*x^3)/(1-x)^11.
a(n) = A062140(n+3, 3).
a(n) = (n+3)!*binomial(n+7, 7)/3!.
If we define f(n,i,x) = Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)* Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-3)=(-1)^(n-1)*f(n,3,-8), (n>=3). - Milan Janjic, Mar 01 2009

A062262 Fifth (unsigned) column sequence of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 45, 1350, 34650, 831600, 19459440, 454053600, 10702692000, 256864608000, 6307453152000, 158947819430400, 4118193503424000, 109818493424640000, 3015784780968960000, 85303626661693440000
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Crossrefs

Programs

  • Magma
    [Factorial(n+4)*Binomial(n+8,8)/24: n in [0..30]]; // G. C. Greubel, May 13 2018
  • Mathematica
    Table[(n+4)!*Binomial[n+8, 8]/4!, {n, 0, 30}] (* G. C. Greubel, May 13 2018 *)
    With[{nn=20},CoefficientList[Series[(1+32x+168x^2+224x^3+70x^4)/(1-x)^13,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Mar 18 2025 *)
  • PARI
    { f=6; for (n=0, 100, f*=n + 4; write("b062262.txt", n, " ", f*binomial(n + 8, 8)/24) ) } \\ Harry J. Smith, Aug 03 2009
    

Formula

E.g.f.: (1+32*x+168*x^2+224*x^3+70*x^4)/(1-x)^13.
a(n) = A062140(n+4, 4).
a(n) = (n+4)!*binomial(n+8, 8)/4!.
If we define f(n,i,x)= Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)* Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-4) = (-1)^n*f(n,4,-9), (n>=4). - Milan Janjic, Mar 01 2009

A062263 Sixth (unsigned) column of triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 60, 2310, 73920, 2162160, 60540480, 1664863200, 45664819200, 1261490630400, 35321737651200, 1006669523059200, 29284931579904000, 871226714502144000, 26538906072526848000, 828392996692445184000
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Crossrefs

Programs

  • Magma
    [Factorial(n+5)*Binomial(n+9, 9)/Factorial(5): n in [0..20]]; // G. C. Greubel, May 12 2018
  • Mathematica
    Table[(n+5)!*Binomial[n+9,9]/5!, {n, 0, 20}] (* G. c. Greubel, May 12 2018 *)
  • PARI
    { f=24; for (n=0, 100, f*=n + 5; write("b062263.txt", n, " ", f*binomial(n + 9, 9)/120) ) } \\ Harry J. Smith, Aug 03 2009
    

Formula

E.g.f.: N(4;5, x)/(1-x)^15, with N(4;5, x) := Sum_{k=0..5} A062264(5, k)* x^k = 1 + 45*x + 360*x^2 + 840*x^3 + 630*x^4 + 226*x^5.
a(n) = A062140(n+5, 5).
a(n) = (n+5)!*binomial(n+9, 9)/5!.
If we define f(n,i,x) = Sum_{k=i..n} Sum_{j=i..k} binomial(k,j) * Stirling1(n,k) * Stirling2(j,i) * x^(k-j) then a(n-9) = (-1)^(n-1)*f(n,9,-6), (n>=9). - Milan Janjic, Mar 01 2009

A062265 Row sums of signed triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 4, 19, 104, 641, 4364, 32251, 254176, 2091841, 17435924, 138844931, 891248984, 263059969, -163754125796, -4970760027029, -117798281164336, -2588474951884159, -55489648295242204, -1184521077396558989, -25406942370946446776, -549455868757454486399, -11980725887273702949076
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Crossrefs

Cf. A062140.

Programs

  • Magma
    [Factorial(n)*(&+[(-1)^k*Binomial(n+4,n-k)/Factorial(k): k in [0..n]]): n in [0..30]]; // G. C. Greubel, May 13 2018
  • Mathematica
    Table[n!*LaguerreL[n, 4, 1],{n,0,20}] (* Vaclav Kotesovec, Aug 01 2013 *)
  • PARI
    for(n=0,30, print1(n!*sum(k=0,n, (-1)^k*binomial(n+4,n-k)/k!), ", ")) \\ G. C. Greubel, May 13 2018
    
  • PARI
    a(n) = vecsum(Vec(n!*pollaguerre(n, 4))); \\ Michel Marcus, Feb 06 2021
    

Formula

E.g.f.: exp(-x/(1-x))/(1-x)^5.
a(n) = n! * Sum_{m=0..n} (-1)^m * binomial(n+4, n-m)/m!.
a(n) = 2*(n+1)*a(n-1) - (n-1)*(n+3)*a(n-2). - Vaclav Kotesovec, Aug 01 2013

A021009 Triangle of coefficients of Laguerre polynomials n!*L_n(x) (rising powers of x).

Original entry on oeis.org

1, 1, -1, 2, -4, 1, 6, -18, 9, -1, 24, -96, 72, -16, 1, 120, -600, 600, -200, 25, -1, 720, -4320, 5400, -2400, 450, -36, 1, 5040, -35280, 52920, -29400, 7350, -882, 49, -1, 40320, -322560, 564480, -376320, 117600, -18816, 1568, -64, 1, 362880, -3265920
Offset: 0

Views

Author

Keywords

Comments

In absolute values, this sequence also gives the lower triangular readout of the exponential of a matrix whose entry {j+1,j} equals (j-1)^2 (and all other entries are zero). - Joseph Biberstine (jrbibers(AT)indiana.edu), May 26 2006
A partial permutation on a set X is a bijection between two subsets of X. |T(n,n-k)| equals the numbers of partial permutations of an n-set having domain cardinality equal to k. Let E denote the operator D*x*D, where D is the derivative operator d/dx. Then E^n = Sum_{k = 0..n} |T(n,k)|*x^k*D^(n+k). - Peter Bala, Oct 28 2008
The unsigned triangle is the generalized Riordan array (exp(x), x) with respect to the sequence n!^2 as defined by Wang and Wang (the generalized Riordan array (exp(x), x) with respect to the sequence n! is Pascal's triangle A007318, and with respect to the sequence n!*(n+1)! is A105278). - Peter Bala, Aug 15 2013
The unsigned triangle appears on page 83 of Ser (1933). - N. J. A. Sloane, Jan 16 2020

Examples

			The triangle a(n,m) starts:
n\m   0       1      2       3      4      5    6  7  8
0:    1
1:    1      -1
2:    2      -4      1
3:    6     -18      9      -1
4:   24     -96     72     -16      1
5:  120    -600    600    -200     25     -1
6:  720   -4320   5400   -2400    450    -36    1
7: 5040  -35280  52920  -29400   7350   -882   49  -1
8:40320 -322560 564480 -376320 117600 -18816 1568 -64 1
...
From _Wolfdieter Lang_, Jan 31 2013 (Start)
Recurrence (usual one): a(4,1) = 7*(-18) - 6 - 3^2*(-4) = -96.
Recurrence (simplified version): a(4,1) = 5*(-18) - 6 = -96.
Recurrence (Sage program): |a(4,1)| = 6 + 3*18 + 4*9 = 96. (End)
Embedded recurrence (Maple program): a(4,1) = -4!*(1 + 3) = -96.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 799.
  • G. Rota, Finite Operator Calculus, Academic Press, New York, 1975.
  • J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, p. 83.

Crossrefs

Row sums give A009940, alternating row sums are A002720.
Column sequences (unsigned): A000142, A001563, A001809-A001812 for m=0..5.
Central terms: A295383.
For generators and generalizations see A132440.

Programs

  • Magma
    /* As triangle: */ [[((-1)^k)*Factorial(n)*Binomial(n, k)/Factorial(k): k in [0..n]]: n in [0.. 10]]; // Vincenzo Librandi, Jan 18 2020
  • Maple
    A021009 := proc(n,k) local S; S := proc(n,k) option remember; `if`(k = 0, 1, `if`( k > n, 0, S(n-1,k-1)/k + S(n-1,k))) end: (-1)^k*n!*S(n,k) end: seq(seq(A021009(n,k), k=0..n), n=0..8); # Peter Luschny, Jun 21 2017
    # Alternative for the unsigned case (function RiordanSquare defined in A321620):
    RiordanSquare(add(x^m, m=0..10), 10, true); # Peter Luschny, Dec 06 2018
  • Mathematica
    Flatten[ Table[ CoefficientList[ n!*LaguerreL[n, x], x], {n, 0, 9}]] (* Jean-François Alcover, Dec 13 2011 *)
  • PARI
    p(n) = denominator(bestapprPade(Ser(vector(2*n, k, (k-1)!))));
    concat(1, concat(vector(9, n, Vec(-p(n)))))  \\ Gheorghe Coserea, Dec 01 2016
    
  • PARI
    {T(n, k) = if( n<0, 0, n! * polcoeff( sum(i=0, n, binomial(n, n-i) * (-x)^i / i!), k))}; /* Michael Somos, Dec 01 2016 */
    
  • PARI
    row(n) = Vecrev(n!*pollaguerre(n)); \\ Michel Marcus, Feb 06 2021
    
  • Sage
    def A021009_triangle(dim): # computes unsigned T(n,k).
        M = matrix(ZZ,dim,dim)
        for n in (0..dim-1): M[n,n] = 1
        for n in (1..dim-1):
            for k in (0..n-1):
                M[n,k] = M[n-1,k-1]+(2*k+1)*M[n-1,k]+(k+1)^2*M[n-1,k+1]
        return M
    A021009_triangle(9) # Peter Luschny, Sep 19 2012
    

Formula

a(n, m) = ((-1)^m)*n!*binomial(n, m)/m! = ((-1)^m)*((n!/m!)^2)/(n-m)! if n >= m, otherwise 0.
E.g.f. for m-th column: (-x/(1-x))^m /((1-x)*m!), m >= 0.
Representation (of unsigned a(n, m)) as special values of Gauss hypergeometric function 2F1, in Maple notation: n!*(-1)^m*hypergeom([ -m, n+1 ], [ 1 ], 1)/m!. - Karol A. Penson, Oct 02 2003
Sum_{m>=0} (-1)^m*a(n, m) = A002720(n). - Philippe Deléham, Mar 10 2004
E.g.f.: (1/(1-x))*exp(x*y/(x-1)). - Vladeta Jovovic, Apr 07 2005
Sum_{n>=0, m>=0} a(n, m)*(x^n/n!^2)*y^m = exp(x)*BesselJ(0, 2*sqrt(x*y)). - Vladeta Jovovic, Apr 07 2005
Matrix square yields the identity matrix: L^2 = I. - Paul D. Hanna, Nov 22 2008
From Tom Copeland, Oct 20 2012: (Start)
Symbolically, with D=d/dx and LN(n,x)=n!L_n(x), define :Dx:^j = D^j x^j, :xD:^j = x^j D^j, and LN(.,x)^j = LN(j,x) = row polynomials of A021009.
Then some useful relations are
1) (:Dx:)^n = LN(n,-:xD:) [Rodriguez formula]
2) (xDx)^n = x^n D^n x^n = x^n LN(n,-:xD:) [See Al-Salam ref./A132440]
3) (DxD)^n = D^n x^n D^n = LN(n,-:xD:) D^n [See ref. in A132440]
4) umbral composition LN(n,LN(.,x))= x^n [See Rota ref.]
5) umbral comp. LN(n,-:Dx:) = LN(n,-LN(.,-:xD:)) = 2^n LN(n,-:xD:/2)= n! * (n-th row e.g.f.(x) of A038207 with x replaced by :xD:).
An example for 2) is the operator (xDx)^2 = (xDx)(xDx) = xD(x^2 + x^3D)= 2x^2 + 4x^3 D + x^4 D^2 = x^2 (2 + 4x D + x^2 D^2) = x^2 (2 + 4 :xD: + :xD:^2) = x^2 LN(2,-:xD:) = x^2 2! L_2(-:xD:).
An example of the umbral composition in 5) is given in A038207.
The op. xDx is related to the Euler/binomial transformation for power series/o.g.f.s. through exp(t*xDx) f(x) = f[x/(1-t*x)]/(1-t*x) and to the special Moebius/linear fractional/projective transformation z exp(-t*zDz)(1/z)f(z) = f(z/(1+t*z)).
For a general discussion of umbral calculus see the Gessel link. (End)
From Wolfdieter Lang, Jan 31 2013: (Start)
Standard recurrence derived from the three term recurrence of the orthogonal polynomials system {n!*L(n,x)}: L(n,x) = (2*n - 1 - x)*L(n-1,x) - (n-1)^2*L(n-2,x), n>=1, L(-1,x) = 0, L(0,x) = 1.
a(n,m) = (2*n-1)*a(n-1,m) - a(n-1,m-1) - (n-1)^2*a(n-2,m),
n >=1, with a(n,-1) = 0, a(0,0) = 1, a(n,m) = 0 if n < m. (compare this with Peter Luschny's program for the unsigned case |a(n,m)| = (-1)^m*a(n,m)).
Simplified recurrence (using column recurrence from explicit form for a(n,m) given above):
a(n,m) = (n+m)*a(n-1,m) - a(n-1,m-1), n >= 1, a(0,0) = 1, a(n,-1) = 0, a(n,m) = 0 if n < m. (End)
|T(n,k)| = [x^k] (-1)^n*U(-n,1,-x), where U(a,b,x) is Kummer's hypergeometric U function. - Peter Luschny, Apr 11 2015
T(n,k) = (-1)^k*n!*S(n,k) where S(n,k) is recursively defined by: "if k = 0 then 1 else if k > n then 0 else S(n-1,k-1)/k + S(n-1,k)". - Peter Luschny, Jun 21 2017
The unsigned case is the exponential Riordan square (see A321620) of the factorial numbers. - Peter Luschny, Dec 06 2018
Omitting the diagonal and signs, this array is generated by the commutator [D^n,x^n] = D^n x^n - x^n D^n = Sum_{i=0..n-1} ((n!/i!)^2/(n-i)!) x^i D^i on p. 9 of both papers by Belov-Kanel and Kontsevich. - Tom Copeland, Jan 23 2020

Extensions

Name changed and table given by Wolfdieter Lang, Nov 28 2011

A066667 Coefficient triangle of generalized Laguerre polynomials (a=1).

Original entry on oeis.org

1, 2, -1, 6, -6, 1, 24, -36, 12, -1, 120, -240, 120, -20, 1, 720, -1800, 1200, -300, 30, -1, 5040, -15120, 12600, -4200, 630, -42, 1, 40320, -141120, 141120, -58800, 11760, -1176, 56, -1, 362880, -1451520, 1693440, -846720, 211680, -28224, 2016
Offset: 0

Views

Author

Christian G. Bower, Dec 17 2001

Keywords

Comments

Same as A008297 (Lah triangle) except for signs.
Row sums give A066668. Unsigned row sums give A000262.
The Laguerre polynomials L(n;x;a=1) under discussion are connected with Hermite-Bell polynomials p(n;x) for power -1 (see also A215216) defined by the following relation: p(n;x) := x^(2n)*exp(x^(-1))*(d^n exp(-x^(-1))/dx^n). We have L(n;x;a=1)=(-x)^(n-1)*p(n;1/x), p(n+1;x)=x^2(dp(n;x)/dx)+(1-2*n*x)p(n;x), and p(1;x)=1, p(2;x)=1-2*x, p(3;x)=1-6*x+6*x^2, p(4;x)=1-12*x+36*x^2-24*x^3, p(5;x)=1-20*x+120*x^2-240*x^3+120*x^4. Note that if we set w(n;x):=x^(2n)*p(n;1/x) then w(n+1;x)=(w(n;x)-(dw(n;x)/dx))*x^2, which is almost analogous to the recurrence formula for Bell polynomials B(n+1;x)=(B(n;x)+(dB(n;x)/dx))*x. - Roman Witula, Aug 06 2012.

Examples

			Triangle a(n,m) begins
n\m     0        1       2       3      4      5    6   7  8
0:      1
1:      2       -1
2:      6       -6       1
3:     24      -36      12      -1
4:    120     -240     120     -20      1
5:    720    -1800    1200    -300     30     -1
6:   5040   -15120   12600   -4200    630    -42    1
7:  40320  -141120  141120  -58800  11760  -1176   56  -1
8: 362880 -1451520 1693440 -846720 211680 -28224 2016 -72  1
9: 3628800, -16329600, 21772800, -12700800, 3810240, -635040, 60480, -3240, 90, -1.
Reformatted and extended by _Wolfdieter Lang_, Jan 31 2013.
From _Wolfdieter Lang_, Jan 31 2013 (Start)
Recurrence (standard): a(4,2) = 2*4*12 - (-36) - 4*3*1 = 120.
Recurrence (simple): a(4,2) = 7*12 - (-36) = 120. (End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 778 (22.5.17).
  • F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge, 1998, p. 95 (4.1.62)
  • R. Witula, E. Hetmaniok, and D. Slota, The Hermite-Bell polynomials for negative powers, (submitted, 2012)

Crossrefs

Programs

  • Maple
    A066667 := (n, k) -> (-1)^k*binomial(n, k)*(n + 1)!/(k + 1)!:
    for n from 0 to 9 do seq(A066667(n,k), k = 0..n) od; # Peter Luschny, Jun 22 2022
  • Mathematica
    Table[(-1)^m*Binomial[n, m]*(n + 1)!/(m + 1)!, {n, 0, 8}, {m, 0, n}] // Flatten (* Michael De Vlieger, Sep 04 2019 *)
  • PARI
    row(n) = Vecrev(n!*pollaguerre(n, 1)); \\ Michel Marcus, Feb 06 2021

Formula

E.g.f. (relative to x, keep y fixed): A(x)=(1/(1-x))^2*exp(x*y/(x-1)).
From Wolfdieter Lang, Jan 31 2013: (Start)
a(n,m) = (-1)^m*binomial(n,m)*(n+1)!/(m+1)!, n >= m >= 0. [corrected by Georg Fischer, Oct 26 2022]
Recurrence from standard three term recurrence for orthogonal generalized Laguerre polynomials {P(n,x):=n!*L(1,n,x)}:
P(n,x) = (2*n-x)*P(n-1,x) - n*(n-1)*P(n-2), n>=1, P(-1,x) = 0, P(0,x) = 1.
a(n,m) = 2*n*a(n-1,m) - a(n-1,m-1) - n*(n-1)*a(n-2,m), n>=1, a(0,0) =1, a(n,-1) = 0, a(n,m) = 0 if n < m.
Simplified recurrence from explicit form of a(n,m):
a(n,m) = (n+m+1)*a(n-1,m) - a(n-1,m-1), n >= 1, a(0,0) =1, a(n,-1) = 0, a(n,m) = 0 if n < m.
(End)

A062264 Coefficient triangle of certain polynomials N(4; m,x).

Original entry on oeis.org

1, 1, 5, 1, 12, 15, 1, 21, 63, 35, 1, 32, 168, 224, 70, 1, 45, 360, 840, 630, 126, 1, 60, 675, 2400, 3150, 1512, 210, 1, 77, 1155, 5775, 11550, 9702, 3234, 330, 1, 96, 1848, 12320, 34650, 44352, 25872, 6336, 495, 1, 117, 2808, 24024, 90090, 162162, 144144, 61776, 11583, 715
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

The e.g.f. of the m-th (unsigned) column sequence without leading zeros of the generalized (a=4) Laguerre triangle L(4; n+m,m) = A062140(n+m,m), n >= 0, is N(4; m,x)/(1-x)^(5+2*m), with the row polynomials N(4; m,x) := Sum_{k=0..m} T(m,k)*x^k.

Examples

			Triangle begins as:
  1;
  1,   5;
  1,  12,   15;
  1,  21,   63,    35;
  1,  32,  168,   224,     70;
  1,  45,  360,   840,    630,    126;
  1,  60,  675,  2400,   3150,   1512,    210;
  1,  77, 1155,  5775,  11550,   9702,   3234,    330;
  1,  96, 1848, 12320,  34650,  44352,  25872,   6336,    495;
  1, 117, 2808, 24024,  90090, 162162, 144144,  61776,  11583,   715;
  1, 140, 4095, 43680, 210210, 504504, 630630, 411840, 135135, 20020, 1001;
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), A132813 (c=2), A062196 (c=3), A062145 (c=4), this sequence (c=5), A062190 (c=6).
Columns: A028347 (k=2), A104473 (k=3), A104474 (k=4), A104475 (k=5), A027814 (k=6), A103604 (k=7), A104476 (k=8), A104478 (k=9).
Diagonals: A000332 (k=n), A027810 (k=n-1), A105249 (k=n-2), A105250 (k=n-3), A105251 (k=n-4), A105252 (k=n-5), A105253 (k=n-6), A105254 (k=n-7).
Sums: A002694 (row).

Programs

  • Magma
    A062264:= func< n,k | Binomial(n,k)*Binomial(n+4,k) >;
    [A062264(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Mar 03 2025
    
  • Mathematica
    A062264[n_, k_]:= Binomial[n,k]*Binomial[n+4,k];
    Table[A062264[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Mar 03 2025 *)
  • SageMath
    def A062264(n,k): return binomial(n,k)*binomial(n+4,k)
    print(flatten([[A062264(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 03 2025

Formula

T(m, k) = [x^k] N(4; m, x), with N(4; m, x) = ((1-x)^(2*m+5))*(d^m/dx^m)((x^m)/(m!*(1-x)^(m+5))).
N(4; m, x) = Sum_{j=0..m} (binomial(m, j)*(2*m+4-j)!/((m+4)!*(m-j)!)*(x^(m-j))*(1-x)^j).
From G. C. Greubel, Mar 03 2025: (Start)
T(n, k) = binomial(n,k)*binomial(n+4,k).
Sum_{k=0..n} (-1)^k*T(n, k) = (1/4)*( (1+(-1)^n)*(-1)^((n+2)/2)*(n^2 + 5*n - 2)*Catalan((n+2)/2)/(n+1) + 8*(1-(-1)^n)*(-1)^((n+1)/2)*Catalan((n+1)/2) ). (End)
Showing 1-10 of 14 results. Next