login
This site is supported by donations to The OEIS Foundation.

 

Logo

Invitation: celebrating 50 years of OEIS, 250000 sequences, and Sloane's 75th, there will be a conference at DIMACS, Rutgers, Oct 9-10 2014.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A133314 Coefficients of list partition transform: reciprocal of an exponential generating function (e.g.f.) 40
1, -1, -1, 2, -1, 6, -6, -1, 8, 6, -36, 24, -1, 10, 20, -60, -90, 240, -120, -1, 12, 30, -90, 20, -360, 480, -90, 1080, -1800, 720, -1, 14, 42, -126, 70, -630, 840, -420, -630, 5040, -4200, 2520, -12600, 15120, -5040, -1, 16, 56, -168, 112 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

COMMENTS

The list partition transform of a sequence a(n) for which a(0)=1 is illustrated by:

b(0) = 1

b(1) = -a(1)

b(2) = -a(2) + 2 a(1)^2

b(3) = -a(3) + 6 a(2)a(1) - 6 a(1)^3

b(4) = -a(4) + 8 a(3)a(1) + 6 a(2)^2 - 36 a(2)*a(1)^2 + 24 a(1)^4

The unsigned coefficients are A049019 with a leading 1. The sign is dependent on the partition as evident from inspection (replace a(n)'s by -1).

Expressed umbrally,

exp(a(.)*x)*exp(b(.)*x) = 1, i.e., (a(.)+b(.))^n = 1 for n=0 and 0 for all other values of n.

Expressed recursively,

b(0) = 1, b(n) = -sum(j=1,...,n) binomial(n,j)*a(j)*b(n-j); which is conditionally self-inverse, i.e., the roles of a(k) and b(k) may be reversed with a(0) = b(0) = 1.

Expressed in matrix form, b(n) form the first column of B = matrix inverse of A .

A = Pascal matrix diagonally multipied by a(n), i.e., A(n,k) = binomial(n,k)*a(n-k).

Some examples of reciprocal pairs of sequences under these operations are:

1) A084358 and -A000262 with the first term set to 1.

2) (1,-1,0,0,...) and (0!,1!,2!,3!,...) with the unsigned associated matrices A128229 and A094587.

3) (1,-1,-1,-1,...) and A000670.

4) A000110 and A000587.

5) (1,-2,-2,0,0,0,...) and (0! c(1),1! c(2),2! c(3),3! c(4),...) where c(n) = A000129(n) with the associated matrices A110327 and A110330.

6) (1,-2,2,0,0,0,...) and (1!,2!,3!,4!,...).

7) Sequences of rising and signed lowering factorials form reciprocal pairs where a(n) = (-1)^n m!/(m-n)! and b(n) = (m-1+n)!/(m-1)! for m=0,1,2,... .

Denote the action of the list partition transform on the sequence a(.) or an invertible matrix M by LPT(a) = b or LPT(M)= M^(-1).

If the matrix equation M = exp(T) also holds, then exp[a(.)*T] * exp[b(.)*T] = exp[(a+b)*T] = Identity matrix because (a(.)+b(.))^n = delta(n), the Kronecker delta.

Therefore, [exp(a*T)]^(-1) = exp[b*T] = exp[LPT(a)*T] = LPT[exp(a*T)].

The fundamental Pascal (A007318), unsigned Lah (A105278) and associated Laguerre matrices can be generated by exponentiation of special infinitesimal matrices (see A132440, A132710 and A132681) such that finding LPT(a) amounts to diagonally multiplying any of the fundamental matrices by a(.), followed by matrix inversion and then extraction of the b(.) factors from the first column (simplest for the Pascal formulae above).

Conversely, the inverses of matrices formed by diagonally multiplying the three fundamental matrices by a(.) are given by diagonally multiplying the fundamental matrices by b(.).

If LPT(M) is defined differently as application of the top formula to a(n) = M^n, then b(n) = (-M)^n and the formalism could even be applied to more general sequences of matrices M(.), providing the reciprocal of exp[t*M(.)].

The group of fundamental lower triangular matrices M = exp(T) such that LPT[exp(a*T)] = exp[LPT(a)*T] = [exp[(a)*T]]^(-1) are obtained by infinitesimal generator matrices of the form T =

0;

t(0), 0;

0, t(1), 0;

0, 0, t(2), 0;

0, 0, 0, t(3), 0;

...

T^m has trivially vanishing terms except along the m'th subdiagonal, which is a sequence of generalized factorials:

[ t(0)*t(1)...t(m-2)*t(m-1), t(1)*t(2)...t(m-1)*t(m), t(2)*t(3)...t(m)*t(m+1), ... ].

Therefore the principal submatrices of T (given by setting t(j) = 0 for j>n-1) are nilpotent with at least [Tsub_n]^(n+1) = 0.

The general group of matrices GM[a(.)] = exp[a(.)*T] can also be obtained through diagonal multiplication of M = exp(T) by the sequence a(.), as in the Pascal matrix example above and their inverses by diagonal multiplication by LPT(a) = b.

Weighted-mappings interpretation for the top partition equation:

Given n pre-nodes (Pre) and k post-nodes (Post), each Pre is connected to only one Post and each Post has at least one Pre connected to it (surjections or onto functions/maps). Weight each Post by -a(m) where m is the number of connections to the Post.

Weight each map by the product of the Post weights and multiply by the number of maps that share the same connectivity. Sum over the possible mappings for n Pre. The result is b(n).

E.g., b(3) = [ 3 Pre to 1 Post ] + [ 3 Pre to 2 Post ] + [ 3 Pre to 3 Post ]

= [1 map with 1 Post with 3 connections] + [ 6 maps with 1 Post with 2 connections and 1 Post with 1 connection] + [6 maps with 3 Post with 1 connection each]

= -a(3) + 6 * [-a(2)*-a(1)] + 6 * [-a(1)*-a(1)*-a(1)].

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 0..250

FORMULA

b(n-1) = (1/n)(d/da(1))p_n[a(1),a(2),...,a(n)] where p_n are the row partition polynomials of the cumulant generator A127671. - Tom Copeland, Oct 13 2012

(E.g.f. of matrix B) = (e.g.f. of b)·exp(x·t) = exp(b(.)·t)·exp(x·t) = exp(x·t)/exp(a(·)·t) = (e.g.f. of A^(-1)) and (e.g.f. of matrix A) = exp(a(·)·t)·exp(x·t) = exp(x·t)/exp(b(·)·t) = (e.g.f. of B^(-1)). These e.g.f. define Appell sequences of polynomials. - Tom Copeland, Mar 22 2014

MATHEMATICA

Clear[a, b]; b[0] = 1; b[n_] := b[n] = -Sum[Binomial[n, j]*a[j]*b[n-j], {j, 1, n}]; row[0] = {1}; row[n_] := Coefficient[b[n], #]& /@ (Times @@ (a /@ #)&) /@ IntegerPartitions[n]; Table[row[n], {n, 0, 8}] // Flatten (* Jean-François Alcover, Apr 23 2014 *)

CROSSREFS

Row lengths in A000041.

Sequence in context: A048998 A213615 A049019 * A208909 A229565 A208919

Adjacent sequences:  A133311 A133312 A133313 * A133315 A133316 A133317

KEYWORD

sign,tabf

AUTHOR

Tom Copeland, Oct 18 2007, Oct 29 2007, Nov 16 2007

EXTENSIONS

More terms from Jean-François Alcover, Apr 23 2014

STATUS

approved

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

Content is available under The OEIS End-User License Agreement .

Last modified September 21 03:47 EDT 2014. Contains 247022 sequences.