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!)
A238385 Shifted lower triangular matrix A238363 with a main diagonal of ones. 16

%I #42 Jul 27 2017 09:09:48

%S 1,1,1,-1,2,1,2,-3,3,1,-6,8,-6,4,1,24,-30,20,-10,5,1,-120,144,-90,40,

%T -15,6,1,720,-840,504,-210,70,-21,7,1,-5040,5760,-3360,1344,-420,112,

%U -28,8,1,40320,-45360,25920,-10080,3024,-756,168,-36,9,1,-362880,403200,-226800,86400,-25200,6048,-1260,240,-45,10,1

%N Shifted lower triangular matrix A238363 with a main diagonal of ones.

%C Shift A238363 and add a main diagonal of ones to obtain this array. The row polynomials form a special Sheffer sequence of polynomials, an Appell sequence.

%F a(n,k) = (-1)^(n+k-1)*n!/((n-k)*k!) for k<n and a(n,n)=1.

%F Along the n-th diagonal (n>0) Diag(n,k) = a(n+k,k) = (-1)^(n-1)(n-1)! * A007318(n+k,k).

%F E.g.f.: (log(1+t)+1)*exp(x*t).

%F E.g.f. for inverse: exp(x*t)/(log(1+t)+1).

%F The lowering/annihilation and raising/creation operators for the row polynomials are L=D=d/dx and R=x+1/[(1+D)(1+log(1+D))], i.e., L p(n,x)= n*p(n-1,x) and R p(n,x)= p(n+1,x).

%F E.g.f. of row sums: (log(1+t)+1)*exp(t). Cf. |row sums-1|=|A002741|.

%F E.g.f. of unsigned row sums: (-log(1-t)+1)*exp(t). Cf. A002104 + 1.

%F Let dP = A132440, the infinitesimal generator for the Pascal matrix, I, the identity matrix, and T, this entry's lower triangular matrix, then exp(T-I)=I+dP, i.e., T=I+log(I+dP). Also, ((T-I)_n)^n=0, where (T-I)_n denotes the n X n submatrix, i.e., (T-I)_n is nilpotent of order n. - _Tom Copeland_, Mar 02 2014

%F Dividing each subdiagonal by its first element (-1)^(n-1)*(n-1)! yields Pascal's triangle A007318. This is equivalent to multiplying the e.g.f. by exp(t)/(log(1+t)+1). - _Tom Copeland_, Apr 16 2014

%F From _Tom Copeland_, Apr 25 2014: (Start)

%F A) T = [St1]*[dP]*[St2] + I = [padded A008275]*A132440*A048993 + I

%F B) = [St1]*[dP]*[St1]^(-1) + I

%F C) = [St2]^(-1)*[dP]*[St2] + I

%F D) = [St2]^(-1)*[dP]*[St1]^(-1) + I,

%F where [St1]=padded A008275 just as [St2]=A048993=padded A008277 and I=identity matrix. Cf. A074909. (End)

%F From _Tom Copeland_, Jul 26 2017: (Start)

%F p_n(x) = (1 + log(1+D)) x^n = (1 + D - D^2/2 + D^3/3- ...) x^n = x^n + n! * Sum_(k=1,..,n) (-1)^(k+1) (1/k) x^(n-k)/(n-k)!.

%F Unsigned T with the first two diagonals nulled gives an exponential infinitesimal generator M (infinigen) for the rencontres numbers A008290, and negated M gives the infinigen for A055137; i.e., with M = |T| - I - dP = -log(I-dP)-dP, then e^M = e^(-dP) / (I-dP) = lower triangular A008290, and e^(-M) = e^dP (I-dP) = A007318 * (I-dP) = lower triangular A055137. The matrix formulation is consistent with the operator relations e^(-D) / (1-D) x^n = n-th row polynomial of A008290 and e^D (1-D) x^n = n-th row polynomial of A055137. (End)

%e The triangle a(n,k) begins:

%e n\k 0 1 2 3 4 5 6 7 8 9 10 ...

%e 0: 1

%e 1: 1 1

%e 2: -1 2 1

%e 3: 2 -3 3 1

%e 4: -6 8 -6 4 1

%e 5: 24 -30 20 -10 5 1

%e 6: -120 144 -90 40 -15 6 1

%e 7: 720 -840 504 -210 70 -21 7 1

%e 8: -5040 5760 -3360 1344 -420 112 -28 8 1

%e 9: 40320 -45360 25920 -10080 3024 -756 168 -36 9 1

%e 10: -362880 403200 -226800 86400 -25200 6048 -1260 240 -45 10 1

%e ... formatted by _Wolfdieter Lang_, Mar 09 2014

%e ----------------------------------------------------------------------------

%Y Cf. A002741, A007318, A008275, A008277, A008290, A048993, A055137, A074909, A132440, A238363.

%K sign,tabl,easy

%O 0,5

%A _Tom Copeland_, Feb 25 2014

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 March 29 08:59 EDT 2024. Contains 371268 sequences. (Running on oeis4.)