# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a225466 Showing 1-1 of 1 %I A225466 #61 Mar 14 2024 15:43:32 %S A225466 1,2,3,4,21,9,8,117,135,27,16,609,1431,702,81,32,3093,13275,12015, %T A225466 3240,243,64,15561,115479,171990,81405,13851,729,128,77997,970515, %U A225466 2238327,1655640,479682,56133,2187,256,390369,7998111,27533142,29893941,13121514,2561706 %N A225466 Triangle read by rows, 3^k*S_3(n, k) where S_m(n, k) are the Stirling-Frobenius subset numbers of order m; n >= 0, k >= 0. %C A225466 The definition of the Stirling-Frobenius subset numbers of order m is in A225468. %C A225466 From _Wolfdieter Lang_, Apr 09 2017: (Start) %C A225466 This is the Sheffer triangle (exp(2*x), exp(3*x) - 1), denoted by S2[3,2]. See also A282629 for S2[3,1]. The stirling2 triangle A048993 is in this notation denoted by S2[1,0]. %C A225466 The a-sequence for this Sheffer triangle has e.g.f. 3*x/log(1+x) and is 3*A006232(n)/A006233(n) (Cauchy numbers of the first kind). For a- and z-sequences for Sheffer triangles see the W. Lang link under A006232, also with references). %C A225466 The z-sequence has e.g.f. (3/(log(1+x)))*(1 - 1/(1+x)^(2/3)) and gives 2*A284862/A284863. %C A225466 The first column k sequences divided by 3^k are A000079, A016127, A016297, A025999. For the e.g.f.s and o.g.f.s see below. %C A225466 The row sums give A284864. The alternating row sums give A284865. %C A225466 This triangle appears in the o.g.f. G(n, x) of the sequence {(2 + 3*m)^n}_{m>=0}, as G(n, x) = Sum_{k=0..n} T(n, k)*k!*x^k/(1-x)^(k+1), n >= 0. Hence the corresponding e.g.f. is, by the linear inverse Laplace transform, E(n, t) = Sum_{m >=0} (2 + 3*m)^n t^m/m! = exp(t)*Sum_{k=0..n} T(n, k)*t^k. %C A225466 The corresponding Eulerian number triangle is A225117(n, k) = Sum_{m=0..k} (-1)^(k-m)*binomial(n-m, k-m)*T(n, m)*m!, 0 <= k <= n. (End) %H A225466 Vincenzo Librandi, Rows n = 0..50, flattened %H A225466 Paweł Hitczenko, A class of polynomial recurrences resulting in (n/log n, n/log^2 n)-asymptotic normality, arXiv:2403.03422 [math.CO], 2024. See p. 9. %H A225466 Peter Luschny, Eulerian polynomials. %H A225466 Peter Luschny, The Stirling-Frobenius numbers. %H A225466 Shi-Mei Ma, Toufik Mansour, and Matthias Schork, Normal ordering problem and the extensions of the Stirling grammar, Russian Journal of Mathematical Physics, 2014, 21(2), arXiv:1308.0169 [math.CO], 2013, p. 12. %F A225466 T(n, k) = (1/k!)*Sum_{j=0..n} binomial(j, n-k)*A_3(n, j) where A_m(n, j) are the generalized Eulerian numbers A225117. %F A225466 For a recurrence see the Maple program. %F A225466 T(n, 0) ~ A000079; T(n, 1) ~ A005057; T(n, n) ~ A000244. %F A225466 From _Wolfdieter Lang_, Apr 09 2017: (Start) %F A225466 T(n, k) = Sum_{j=0..k} binomial(k,j)*(-1)^(j-k)*(2 + 3*j)^n/k!, 0 <= k <= n. %F A225466 E.g.f. of triangle: exp(2*z)*exp(x*(exp(3*z)-1)) (Sheffer type). %F A225466 E.g.f. for sequence of column k is exp(2*x)*((exp(3*x) - 1)^k)/k! (Sheffer property). %F A225466 O.g.f. for sequence of column k is 3^k*x^k/Product_{j=0..k} (1 - (2+3*j)*x). %F A225466 A nontrivial recurrence for the column m=0 entries T(n, 0) = 2^n from the z-sequence given above: T(n,0) = n*Sum_{k=0..n-1} z(k)*T(n-1,k), n >= 1, T(0, 0) = 1. %F A225466 The corresponding recurrence for columns k >= 1 from the a-sequence is T(n, k) = (n/k)* Sum_{j=0..n-k} binomial(k-1+j, k-1)*a(j)*T(n-1, k-1+j). %F A225466 Recurrence for row polynomials R(n, x) (Meixner type): R(n, x) = ((3*x+2) + 3*x*d_x)*R(n-1, x), with differentiation d_x, for n >= 1, with input R(0, x) = 1. %F A225466 (End) %F A225466 Boas-Buck recurrence for column sequence m: T(n, k) = (1/(n - m))*[(n/2)*(4 + 3*m)*T(n-1, k) + m* Sum_{p=m..n-2} binomial(n, p)(-3)^(n-p)*Bernoulli(n-p)*T(p, k)], for n > k >= 0, with input T(k, k) = 3^k. See a comment and references in A282629, An example is given below. - _Wolfdieter Lang_, Aug 11 2017 %e A225466 [n\k][ 0, 1, 2, 3, 4, 5, 6, 7] %e A225466 [0] 1, %e A225466 [1] 2, 3, %e A225466 [2] 4, 21, 9, %e A225466 [3] 8, 117, 135, 27, %e A225466 [4] 16, 609, 1431, 702, 81, %e A225466 [5] 32, 3093, 13275, 12015, 3240, 243, %e A225466 [6] 64, 15561, 115479, 171990, 81405, 13851, 729, %e A225466 [7] 128, 77997, 970515, 2238327, 1655640, 479682, 56133, 2187. %e A225466 ... %e A225466 From _Wolfdieter Lang_, Aug 11 2017: (Start) %e A225466 Recurrence (see the Maple program): T(4, 2) = 3*T(3, 1) + (3*2+2)*T(3, 2) = 3*117 + 8*135 = 1431. %e A225466 Boas-Buck recurrence for column m = 2, and n = 4: T(4,2) = (1/2)*[2*(4 + 3*2)*T(3, 2) + 2*6*(-3)^2*Bernoulli(2)*T(2, 2))] = (1/2)*(20*135 + 12*9*(1/6)*9) = 1431. (End) %p A225466 SF_SS := proc(n, k, m) option remember; %p A225466 if n = 0 and k = 0 then return(1) fi; %p A225466 if k > n or k < 0 then return(0) fi; %p A225466 m*SF_SS(n-1, k-1, m) + (m*(k+1)-1)*SF_SS(n-1, k, m) end: %p A225466 seq(print(seq(SF_SS(n, k, 3), k=0..n)), n=0..5); %t A225466 EulerianNumber[n_, k_, m_] := EulerianNumber[n, k, m] = (If[ n == 0, Return[If[k == 0, 1, 0]]]; Return[(m*(n-k)+m-1)*EulerianNumber[n-1, k-1, m] + (m*k+1)*EulerianNumber[n-1, k, m]]); SFSS[n_, k_, m_] := Sum[ EulerianNumber[n, j, m]*Binomial[j, n-k], {j, 0, n}]/k!; Table[ SFSS[n, k, 3], {n, 0, 8}, {k, 0, n}] // Flatten (* _Jean-François Alcover_, May 29 2013, translated from Sage *) %o A225466 (Sage) %o A225466 @CachedFunction %o A225466 def EulerianNumber(n, k, m) : %o A225466 if n == 0: return 1 if k == 0 else 0 %o A225466 return (m*(n-k)+m-1)*EulerianNumber(n-1,k-1,m) + (m*k+1)*EulerianNumber(n-1,k,m) %o A225466 def SF_SS(n, k, m): %o A225466 return add(EulerianNumber(n,j,m)*binomial(j,n-k) for j in (0..n))/ factorial(k) %o A225466 def A225466(n): return SF_SS(n, k, 3) %o A225466 (PARI) T(n, k) = sum(j=0, k, binomial(k, j)*(-1)^(j - k)*(2 + 3*j)^n/k!); %o A225466 for(n=0, 10, for(k=0, n, print1(T(n, k),", ");); print();) \\ _Indranil Ghosh_, Apr 10 2017 %o A225466 (Python) %o A225466 from sympy import binomial, factorial %o A225466 def T(n, k): return sum(binomial(k, j)*(-1)**(j - k)*(2 + 3*j)**n//factorial(k) for j in range(k + 1)) %o A225466 for n in range(11): print([T(n, k) for k in range(n + 1)]) # _Indranil Ghosh_, Apr 10 2017 %Y A225466 Cf. A048993 (m=1), A154537 (m=2), A225467 (m=4), A225468. %Y A225466 Cf. A000079, A000244, A005057, A016127, A016297, A025999, A006232/A006233, A225117, A225472, A225468, A282629, A284862/A284863, A284864, A284865. %K A225466 nonn,easy,tabl %O A225466 0,2 %A A225466 _Peter Luschny_, May 08 2013 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE