A136656 Coefficients for rewriting generalized falling factorials into ordinary falling factorials.
1, 0, -2, 0, 6, 4, 0, -24, -36, -8, 0, 120, 300, 144, 16, 0, -720, -2640, -2040, -480, -32, 0, 5040, 25200, 27720, 10320, 1440, 64, 0, -40320, -262080, -383040, -199920, -43680, -4032, -128, 0, 362880, 2963520, 5503680, 3764880, 1142400, 163968, 10752, 256, 0, -3628800, -36288000, -82978560
Offset: 0
Examples
Triangle starts: [1] [0, -2] [0, 6, 4] [0, -24, -36, -8] [0, 120, 300, 144, 16] ... Recurrence: a(4,2) = -7*a(3,2)-2*a(3,1) = -7*(-36) -2*(-24) = 300. a(4,2)=300 as sum over the M3 numbers A036040 for the 2 parts partitions of 4: 4*fallfac(-2,1)^1*fallfac(-2,3)^1 + 3*fallfac(-2,2)^2 = 4*(-2)*(-24)+3*6^2 = 300. Row n=3: [0,-24,-36,-8] for the coefficients in rewriting fallfac(-2*t,3)=((-1)^3)*risefac(2*t,3) = ((-1)^3)*(2*t)*(2*t+1)*(2*t+2) = 0*1 -24*t -36*t*(t-1) -8*t*(t-1)*(t-2).
References
- Ch. A. Charalambides, Enumerative Combinatorics, Chapman & Hall/CRC, Boca Raton, Florida, 2002, ch. 8.4 p. 301 ff. Table 8.3 (with row n=0 and column k=0 and s=-2).
Links
- M. Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3
- W. Lang, First ten rows and more.
- Peter Luschny, The Bell Transform
Crossrefs
Programs
-
Maple
# The function BellMatrix is defined in A264428. BellMatrix(n -> (-1)^(n+1)*(n+2)!, 9); # Peter Luschny, Jan 27 2016
-
Mathematica
fallfac[n_, k_] := Pochhammer[n - k + 1, k]; a[n_, k_] := Sum[(-1)^(k - r)*Binomial[k, r]*fallfac[-2*r, n], {r, 0, k}]/k!; Table[a[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 09 2013 *) BellMatrix[f_Function, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]]; rows = 10; M = BellMatrix[(-1)^(# + 1)*(# + 2)!&, rows]; Table[M[[n, k]], {n, 1, rows}, {k, 1, n}] // Flatten (*Jean-François Alcover, Jun 23 2018, after Peter Luschny *)
-
Sage
@CachedFunction def T(n, k): # unsigned case if k == 0: return 1 if n == 0 else 0 return sum(T(n-j-1,k-1)*(j+1)*(j+2)*gamma(n)/gamma(n-j) for j in (0..n-k)) for n in range(7): [T(n,k) for k in (0..n)] # Peter Luschny, Mar 31 2015
Formula
Recurrence: a(n,k) = 0 if n
E.g.f. column k: (1/(1+x)^2 - 1)^k/k!, k>=0. From the Charalambides reference Theorem 8.14, p. 305 for s=-2.(hence a Sheffer triangle of Jabotinsky type).
a(n,k) = sum(((-1)^(k-r))*binomial(k,r)*fallfac(-2*r,n),r=0..k)/k!, n>=k>=0. From the Charalambides reference Theorem 8.15, p. 306 for s=-2.
a(n,k) = sum(S1(n,r)*S2(r,k)*(-2)^r,r=k..n) with the Stirling numbers S1(n,r)= A048993(n,r) and S2(r,k)= A048993(r,k). From the Charalambides reference Theorem 8.13, p.304 for s=-2.
a(n,k) = sum(M_3(n,k,q)*product(fallfac(-2,j)^e(n,m,q,j),j=1..n),q=1..p(n,k)) if n>=k>=1, else 0. Here p(n,k)=A008284(n,m), the number of k parts partitions of n, the M_3 partition numbers are given in A036040 and e(n,m,q,j) is the exponent of j in the q-th k parts partition of n. Rewritten eq. (8.50), Theorem 8.16, p. 307, from the Charalambides reference for s=-2.
Recurrence for the unsigned case: a(n,k) = Sum_{j=0..n-k} a(n-j-1,k-1)*C(n-1,j)*(j+2)! if k<>0 else k^n. - Peter Luschny, Mar 31 2015
Comments