|
%I #45 Feb 23 2021 16:00:33
%S 1,3,1,8,9,1,42,59,18,1,264,450,215,30,1,2160,4114,2475,565,45,1,
%T 20880,43512,30814,9345,1225,63,1,236880,528492,420756,154609,27720,
%U 2338,84,1,3064320,7235568,6316316,2673972,594489,69552,4074,108,1
%N Jabotinsky-triangle related to A039647.
%C Triangle gives the nonvanishing entries of the Jabotinsky matrix for F(z)= A(z)/z = 1/(1-z-z^2) where A(z) is the g.f. of the Fibonacci numbers A000045. (Notation of F(z) as in Knuth's paper.)
%C E(n,x) := Sum_{m=1..n} a(n,m)*x^m, E(0,x)=1, are exponential convolution polynomials: E(n,x+y) = Sum_{k=0..n} binomial(n,k)*E(k,x)*E(n-k,y) (cf. Knuth's paper with E(n,x)= n!*F(n,x)).
%C E.g.f. for E(n,x): (1 - z - z^2)^(-x).
%C Explicit a(n,m) formula: see Knuth's paper for f(n,m) formula with f(k)= A039647(n).
%C E.g.f. for the m-th column sequence: ((-log(1 - z - z^2))^m)/m!.
%C Also the Bell transform of n!*(F(n)+F(n+2)), F(n) the Fibonacci numbers. For the definition of the Bell transform see A264428 and the link. - _Peter Luschny_, Jan 16 2016
%H Vincenzo Librandi, <a href="/A039692/b039692.txt">Rows n = 1..50, flattened</a>
%H D. E. Knuth, <a href="http://arxiv.org/abs/math/9207221">Convolution polynomials</a>, Mathematica J. 2.1 (1992), no. 4, 67-78.
%H Peter Luschny, <a href="https://oeis.org/wiki/User:Peter_Luschny/BellTransform">The Bell transform</a>
%F a(n, 1)= A039647(n)=(n-1)!*L(n), L(n) := A000032(n) (Lucas); a(n, m) = Sum_{j=1..n-m+1} binomial(n-1, j-1)*A039647(j)*a(n-j, m-1), n >= m >= 2.
%F Conjectured row sums: sum_{m=1..n} a(n,m) = A005442(n). - _R. J. Mathar_, Jun 01 2009
%F T(n,m) = n! * Sum_{k=m..n} stirling1(k,m)*binomial(k,n-k)*(-1)^(k+m)/k!. - _Vladimir Kruchinin_, Mar 26 2013
%e 1;
%e 3, 1;
%e 8, 9, 1;
%e 42, 59, 18, 1;
%e 264, 450, 215, 30, 1;
%p A000032 := proc(n) option remember; coeftayl( (2-x)/(1-x-x^2),x=0,n) ; end: A039647 := proc(n) (n-1)!*A000032(n) ; end: A039692 := proc(n,m) option remember ; if m = 1 then A039647(n) ; else add( binomial(n-1,j-1)*A039647(j)*procname(n-j,m-1),j=1..n-m+1) ; fi; end: # _R. J. Mathar_, Jun 01 2009
%t t[n_, m_] := n!*Sum[StirlingS1[k, m]*Binomial[k, n-k]*(-1)^(k+m)/k!, {k, m, n}]; Table[t[n, m], {n, 1, 9}, {m, 1, n}] // Flatten (* _Jean-François Alcover_, Jun 21 2013, after _Vladimir Kruchinin_ *)
%o (Maxima) T(n,m) := n!*sum((stirling1(k,m)*binomial(k,n-k))*(-1)^(k+m)/k!,k,m,n); \\ _Vladimir Kruchinin_, Mar 26 2013
%o (PARI)
%o T(n,m) = n!*sum(k=m,n, (stirling(k,m,1)*binomial(k,n-k))*(-1)^(k+m)/k!);
%o for(n=1,10,for(k=1,n,print1(T(n,k),", "));print());
%o /* _Joerg Arndt_, Mar 27 2013 */
%o (Sage) # uses[bell_matrix from A264428]
%o # Adds 1,0,0,0, ... as column 0 to the left side of the triangle.
%o bell_matrix(lambda n: factorial(n)*(fibonacci(n)+fibonacci(n+2)), 8) # _Peter Luschny_, Jan 16 2016
%Y Cf. A039647, A000032, A000045. Another version of this triangle is in A194938.
%K nonn,tabl
%O 1,2
%A _Wolfdieter Lang_
|
