|
COMMENT
|
The matrix operation b = T*a can be characterized in several ways in terms of the coefficients a(n) and b(n) , their o.g.f.'s A(x) and B(x), or e.g.f's EA(x) and EB(x).
1) b(0) = a(0) , b(n) = a(n) - n a(n-1) for n > 0.
2) b(n) = n! Lag{n,(.)!*Lag[.,a(.),0],-1} , umbrally,
3) b(n) = n! sum(j=0,min(1,n)) (-1)^j * binomial(n,j)*a(n-j)/(n-j)!
4) b(n) = (-1)^n n! Lag(n,a(.),1-n)
5) B(x) = (1-xDx) A(x) = [1-x*Lag(1,-xD:,0)] A(x)
6) EB(x) = (1-x) EA(x)
where D is the derivative w.r.t. x and Lag(n,x,m) is the associated Laguerre polynomial of order m. These formulae are easily generalized for repeated applications of the operator.
c = (1,-1,0,0,0,...) is the sequence associated to T under the list partition transform and the associated operations described in A133314. The reciprocal sequence is d = (0!,1!,2!,3!,4!,...).
Consequently, the inverse of T is TI(n,k) = binomial(n,k)*d(n-k) = A094587 , which has the property that the terms at and below TI(m,m) are the associated sequence under the list partition transform for the inverse for T^(m+1) for m=0,1,2,3,... .
Row sums of T = [formula 3 with all a(n) = 1] = [binomial transform of c] = [coefficients of B(x) with A(x) = 1/(1-x)] = A024000 = (1,0,-1,-2,-3,...) , with e.g.f. = [EB(x) with EA(x) = exp(x)] = (1-x) * exp(x) = exp(x)*exp(c(.)*x) = exp[(1+c(.))*x].
Alternating row sums of T = [formula 3 with all a(n) = (-1)^n] = [finite differences of c] = [coefficients of B(x) with A(x) = 1/(1+x)] = (1,-2,3,-4,...) , with e.g.f. = [EB(x) with EA(x) = exp(-x)] = (1-x) * exp(-x) = exp(-x)*exp(c(.)*x) = exp[-(1-c(.))*x].
An e.g.f. for the o.g.f's for repeated applications of T on A(x) is given by
exp[t*(1-xDx)] A(x) = e^t * sum(n=0,1,...) (-t*x)^n * Lag(n,-:xD:,0) A(x)
= e^t * exp{[-t*u/(1+t*u)]*:xD:} / (1+t*u) A(x) (eval. at u=x)
= e^t * A[x/(1+t*x)]/(1+t*x) .
See A132014 for more notes on repeated applications.
|
