[go: up one dir, main page]

login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A256550 Triangle read by rows, T(n,k) = EL(n,k)/(n-k+1)! and EL(n,k) the matrix-exponential of the unsigned Lah numbers scaled by exp(-1), for n>=0 and 0<=k<=n. 1

%I #8 Apr 19 2015 08:49:49

%S 1,0,1,0,1,1,0,2,3,1,0,5,12,6,1,0,15,50,40,10,1,0,52,225,250,100,15,1,

%T 0,203,1092,1575,875,210,21,1,0,877,5684,10192,7350,2450,392,28,1,0,

%U 4140,31572,68208,61152,26460,5880,672,36,1

%N Triangle read by rows, T(n,k) = EL(n,k)/(n-k+1)! and EL(n,k) the matrix-exponential of the unsigned Lah numbers scaled by exp(-1), for n>=0 and 0<=k<=n.

%F T(n+1,1) = Bell(n) = A000110(n).

%F T(n+2,2) = C(n+2,2)*Bell(n) = A105479(n+2).

%F T(n+1,n) = A000217(n).

%F T(n+2,n) = A008911(n+1).

%e Triangle starts:

%e 1;

%e 0, 1;

%e 0, 1, 1;

%e 0, 2, 3, 1;

%e 0, 5, 12, 6, 1;

%e 0, 15, 50, 40, 10, 1;

%e 0, 52, 225, 250, 100, 15, 1;

%e 0, 203, 1092, 1575, 875, 210, 21, 1;

%o (Sage)

%o def T(dim) :

%o M = matrix(ZZ, dim)

%o for n in range(dim) :

%o M[n, n] = 1

%o for k in range(n) :

%o M[n,k] = (k*n*gamma(n)^2)/(gamma(k+1)^2*gamma(n-k+1))

%o E = M.exp()/exp(1)

%o for n in range(dim) :

%o for k in range(n) :

%o M[n,k] = E[n,k]/factorial(n-k+1)

%o return M

%o T(8) # Computes the sequence as a lower triangular matrix.

%Y Cf. A000110, A000217, A008911, A105479, A256551 (matrix inverse).

%K nonn,tabl,easy

%O 0,8

%A _Peter Luschny_, Apr 01 2015

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 31 08:34 EDT 2024. Contains 375560 sequences. (Running on oeis4.)