Bell polynomials

In combinatorial mathematics, the Bell polynomials, named in honor of Eric Temple Bell, are used in the study of set partitions. They are related to Stirling and Bell numbers. They also occur in many applications, such as in the Faà di Bruno's formula.

Bell polynomials

Exponential Bell polynomials

The partial or incomplete exponential Bell polynomials are a triangular array of polynomials given by

B_{n,k}(x_{1},x_{2},\dots ,x_{n-k+1})=\sum {n! \over j_{1}!j_{2}!\cdots j_{n-k+1}!}\left({x_{1} \over 1!}\right)^{j_{1}}\left({x_{2} \over 2!}\right)^{j_{2}}\cdots \left({x_{n-k+1} \over (n-k+1)!}\right)^{j_{n-k+1}},

where the sum is taken over all sequences j₁, j₂, j₃, ..., j_n−k+1 of non-negative integers such that these two conditions are satisfied:

j_{1}+j_{2}+\cdots +j_{n-k+1}=k,

j_{1}+2j_{2}+3j_{3}+\cdots +(n-k+1)j_{n-k+1}=n.

The sum

B_{n}(x_{1},\dots ,x_{n})=\sum _{{k=1}}^{n}B_{{n,k}}(x_{1},x_{2},\dots ,x_{{n-k+1}})

is called the nth complete exponential Bell polynomial.

Ordinary Bell polynomials

Likewise, the partial ordinary Bell polynomial, in contrast to the usual exponential Bell polynomial defined above, is given by

{\hat {B}}_{n,k}(x_{1},x_{2},\ldots ,x_{n-k+1})=\sum {\frac {k!}{j_{1}!j_{2}!\cdots j_{n-k+1}!}}x_{1}^{j_{1}}x_{2}^{j_{2}}\cdots x_{n-k+1}^{j_{n-k+1}},

where the sum runs over all sequences j₁, j₂, j₃, ..., j_n−k+1 of non-negative integers such that

j_{1}+j_{2}+\cdots +j_{n-k+1}=k,

j_{1}+2j_{2}+\cdots +(n-k+1)j_{n-k+1}=n.

The ordinary Bell polynomials can be expressed in the terms of exponential Bell polynomials:

{\hat {B}}_{n,k}(x_{1},x_{2},\ldots ,x_{n-k+1})={\frac {k!}{n!}}B_{n,k}(1!\cdot x_{1},2!\cdot x_{2},\ldots ,(n-k+1)!\cdot x_{n-k+1}).

In general, Bell polynomial refers to the exponential Bell polynomial, unless otherwise explicitly stated.

Combinatorial meaning

The exponential Bell polynomial encodes the information related to the ways a set can be partitioned. For example, if we consider a set {A, B, C}, it can be partitioned into two non-empty, non-overlapping subsets, which is also referred to as parts or blocks, in 3 different ways:

{{A}, {B, C}}

{{B}, {A, C}}

{{C}, {B, A}}

Thus, we can encode the information regarding these partitions as

B_{3,2}(x_{1},x_{2})=3x_{1}x_{2}.

Here, the subscripts of B_3,2 tells us that we are considering the partitioning of set with 3 elements into 2 blocks. The subscript of each x_i indicates the presence of block with i elements (or block of size i) in a given partition. So here, x₂ indicates the presence of a block with two elements. Similarly, x₁ indicates the presence of a block with a single element. The exponent of x_i^j indicates that there are j such blocks of size i in a single partition. Here, since both x₁ and x₂ has exponent 1, it indicates that there is only one such block in a given partition. The coefficient of the monomial indicates how many such partitions there are. For our case, there are 3 partitions of a set with 3 elements into 2 blocks, where in each partition the elements are divided into two blocks of sizes 1 and 2.

Since any set can be divided into a single block in only one way, the above interpretation would mean that B_n,1 = x_n. Similarly, since there is only one way that a set with n elements be divided into n singletons, B_n,n = x₁ⁿ.

As a more complicated example, consider

B_{6,2}(x_{1},x_{2},x_{3},x_{4},x_{5})=6x_{5}x_{1}+15x_{4}x_{2}+10x_{3}^{2}.

This tells us that if a set with 6 elements is divided into 2 blocks, then we can have 6 partitions with blocks of size 1 and 5, 15 partitions with blocks of size 4 and 2, and 10 partitions with 2 blocks of size 3.

Note that the sum of the subscripts in a monomials is equal to the total number of elements. Thus, the number of monomials that appear in the partial Bell polynomial is equal to the number of ways the integer n can be expressed as a summation of k positive integers. This is the same as the integer partition of n into k parts. For instance, in the above examples, the integer 3 can be partitioned into two parts as 2+1 only. Thus, there is only one monomial in B_3,2. However, the integer 6 can be partitioned into two parts as 5+1, 4+2, and 3+3. Thus, there are three monomials in B_6,2. Indeed, the subscripts of the variables in a monomial are the same as those given by the integer partition, indicating the sizes of the different blocks. The total number of monomials appearing in a complete Bell polynomial B_n is thus equal to the total number of integer partitions of n.

Also note that the degree of each monomial, which is the sum of the exponents of each variable in the monomial, is equal to the number of blocks the set is divided into. That is, j₁ + j₂ + ... = k . Thus, given a complete Bell polynomial B_n, we can separate the partial Bell polynomial B_n,k by collecting all those monomials with degree k.

Finally, if we disregard the sizes of the blocks and put all x_i = x, then the summation of the coefficients of the partial Bell polynomial B_n,k will give the total number of ways that a set with n elements can be partitioned into k blocks, which is the same as the Stirling numbers of the second kind. Also, the summation of all the coefficients of the complete Bell polynomial B_n will give us the total number of ways a set with n elements can be partitioned into non-overlapping subsets, which is the same as the Bell number.

In general, if the integer n is partitioned into a sum in which "1" appears j₁ times, "2" appears j₂ times, and so on, then the number of partitions of a set of size n that collapse to that partition of the integer n when the members of the set become indistinguishable is the corresponding coefficient in the polynomial.

Examples

For example, we have

B_{{6,2}}(x_{1},x_{2},x_{3},x_{4},x_{5})=6x_{5}x_{1}+15x_{4}x_{2}+10x_{3}^{2}

because there are

6 ways to partition a set of 6 as 5 + 1,

15 ways to partition a set of 6 as 4 + 2, and

10 ways to partition a set of 6 as 3 + 3.

Similarly,

B_{{6,3}}(x_{1},x_{2},x_{3},x_{4})=15x_{4}x_{1}^{2}+60x_{3}x_{2}x_{1}+15x_{2}^{3}

because there are

15 ways to partition a set of 6 as 4 + 1 + 1,

60 ways to partition a set of 6 as 3 + 2 + 1, and

15 ways to partition a set of 6 as 2 + 2 + 2.

Properties

Generating function

The exponential partial Bell polynomials can be defined by the double series expansion of its generating function:

{\begin{aligned}\Phi (t,u)&=\exp \left(u\sum _{j=1}^{\infty }x_{j}{\frac {t^{j}}{j!}}\right)=\sum _{n,k\geq 0}B_{n,k}(x_{1},\ldots ,x_{n-k+1}){\frac {t^{n}}{n!}}u^{k}\\&=1+\sum _{n=1}^{\infty }{\frac {t^{n}}{n!}}\left\{\sum _{k=1}^{n}u^{k}B_{n,k}(x_{1},\ldots ,x_{n-k+1})\right\}\end{aligned}}

In other words, by what amounts to the same, by the series expansion of the exponential:

{\frac {1}{k!}}\left(\sum _{j=1}^{\infty }x_{j}{\frac {t^{j}}{j!}}\right)^{k}=\sum _{n=k}^{\infty }B_{n,k}(x_{1},\ldots ,x_{n-k+1}){\frac {t^{n}}{n!}},\qquad k=0,1,2,\ldots

The complete exponential Bell polynomial is defined by $\Phi (t,1)$ , or in other words:

\Phi (t,1)=\exp \left(\sum _{j=1}^{\infty }x_{j}{\frac {t^{j}}{j!}}\right)=\sum _{n=0}^{\infty }B_{n}(x_{1},\ldots ,x_{n}){\frac {t^{n}}{n!}}.

Thus, the n-th complete Bell polynomial is given by

B_{n}(x_{1},\ldots ,x_{n})=\left.{\frac {\mathrm {d} ^{n}}{\mathrm {d} t^{n}}}\exp \left(\sum _{j=1}^{\infty }x_{j}{\frac {t^{j}}{j!}}\right)\right|_{t=0}.

Likewise, the ordinary partial Bell polynomial can be defined by the generating function

{\hat {\Phi }}(t,u)=\exp(u\sum _{j=1}^{\infty }x_{j}t^{j})=\sum _{k\leq n}{\hat {B}}_{n,k}(x_{1},\ldots ,x_{n-k+1})t^{n}{\frac {u^{k}}{k!}}.

Or, equivalently, by series expansion of the exponential

\left(\sum _{j=1}^{\infty }x_{j}t^{j}\right)^{k}=\sum _{n=k}^{\infty }{\hat {B}}_{n,k}(x_{1},\ldots ,x_{n-k+1})t^{n}.

Recurrence relations

The complete Bell polynomials can be recursively defined as,

B_{n+1}(x_{1},\ldots ,x_{n+1})=\sum _{i=0}^{n}{n \choose i}B_{n-i}(x_{1},\ldots ,x_{n-i})x_{i+1},

with the initialization $B_{0}=1$ .

The partial Bell polynomials can also be computed efficiently by a recursion relation

B_{n,k}=\sum _{i=1}^{n-k+1}{\binom {n-1}{i-1}}x_{i}B_{n-i,k-1}

where $B_{0,0}=1;$ $B_{n,0}=0\;\mathrm {for} \;n\geq 1;$ $B_{0,k}=0\;\mathrm {for} \;k\geq 1.$

The complete Bell polynomials also satisfy the binomial type relation

B_{n}(x_{1}+y_{1},\ldots ,x_{n}+y_{n})=\sum _{i=0}^{n}{n \choose i}B_{n-i}(x_{1},\ldots ,x_{n-i})B_{i}(y_{1},\ldots ,y_{i}).

Determinant form

The complete Bell polynomial can be expressed as a determinant, as given by

B_{n}(x_{1},\dots ,x_{n})=\det {\begin{bmatrix}x_{1}&{n-1 \choose 1}x_{2}&{n-1 \choose 2}x_{3}&{n-1 \choose 3}x_{4}&{n-1 \choose 4}x_{5}&\cdots &\cdots &x_{n}\\\\-1&x_{1}&{n-2 \choose 1}x_{2}&{n-2 \choose 2}x_{3}&{n-2 \choose 3}x_{4}&\cdots &\cdots &x_{{n-1}}\\\\0&-1&x_{1}&{n-3 \choose 1}x_{2}&{n-3 \choose 2}x_{3}&\cdots &\cdots &x_{{n-2}}\\\\0&0&-1&x_{1}&{n-4 \choose 1}x_{2}&\cdots &\cdots &x_{{n-3}}\\\\0&0&0&-1&x_{1}&\cdots &\cdots &x_{{n-4}}\\\\0&0&0&0&-1&\cdots &\cdots &x_{{n-5}}\\\\\vdots &\vdots &\vdots &\vdots &\vdots &\ddots &\ddots &\vdots \\\\0&0&0&0&0&\cdots &-1&x_{1}\end{bmatrix}}.

Stirling numbers and Bell numbers

The value of the Bell polynomial B_n,k(x₁,x₂,...) on the sequence of factorials equals an unsigned Stirling number of the first kind:

B_{n,k}(0!,1!,\dots ,(n-k)!)=c(n,k)=|s(n,k)|=\left[{n \atop k}\right].

The value of the Bell polynomial B_n,k(x₁,x₂,...) on the sequence of ones equals a Stirling number of the second kind:

B_{{n,k}}(1,1,\dots ,1)=S(n,k)=\left\{{n \atop k}\right\}.

The sum of these values gives the value of the complete Bell polynomial on the sequence of ones:

B_{n}(1,1,\dots ,1)=\sum _{{k=1}}^{n}B_{{n,k}}(1,1,\dots ,1)=\sum _{{k=1}}^{n}\left\{{n \atop k}\right\},

which is the nth Bell number.

Inverse relations

If we define

y_{n}=\sum _{k=1}^{n}B_{n,k}(x_{1},\ldots ,x_{n-k+1}),

then we have the inverse relationship

x_{n}=\sum _{k=1}^{n}(-1)^{k-1}(k-1)!B_{n,k}(y_{1},\ldots ,y_{n-k+1}).

Touchard polynomials

Touchard polynomial $T_{n}(x)=\sum _{{k=0}}^{n}\left\{{n \atop k}\right\}\cdot x^{k}$ can be expressed as the value of the complete Bell polynomial on all arguments being x:

T_{n}(x)=B_{n}(x,x,\dots ,x).

Convolution identity

For sequences x_n, y_n, n = 1, 2, ..., define a sort of convolution by:

(x{\mathbin {\diamondsuit }}y)_{n}=\sum _{j=1}^{n-1}{n \choose j}x_{j}y_{n-j}.

Note that the bounds of summation are 1 and n − 1, not 0 and n .

Let $x_{n}^{{k\diamondsuit }}\,$ be the nth term of the sequence

\displaystyle \underbrace {x{\mathbin {\diamondsuit }}\cdots {\mathbin {\diamondsuit }}x} _{k\ \mathrm {factors} }.\,

Then

B_{n,k}(x_{1},\dots ,x_{n-k+1})={x_{n}^{k\diamondsuit } \over k!}.\,

For example, let us compute $B_{{4,3}}(x_{1},x_{2})$ . We have

x=(x_{1}\ ,\ x_{2}\ ,\ x_{3}\ ,\ x_{4}\ ,\dots )

x{\mathbin {\diamondsuit }}x=(0,\ 2x_{1}^{2}\ ,\ 6x_{1}x_{2}\ ,\ 8x_{1}x_{3}+6x_{2}^{2}\ ,\dots )

x{\mathbin {\diamondsuit }}x{\mathbin {\diamondsuit }}x=(0\ ,\ 0\ ,\ 6x_{1}^{3}\ ,\ 36x_{1}^{2}x_{2}\ ,\dots )

and thus,

B_{4,3}(x_{1},x_{2})={\frac {(x{\mathbin {\diamondsuit }}x{\mathbin {\diamondsuit }}x)_{4}}{3!}}=6x_{1}^{2}x_{2}.

Other identities

$B_{n,k}(1!,2!,\ldots ,(n-k+1)!)={\binom {n-1}{k-1}}{\frac {n!}{k!}}=L(n,k)$ which gives the Lah number
$B_{n,k}(1,2,3,\ldots ,n-k+1)={\binom {n}{k}}k^{n-k}$ which gives the idempotent number

Examples

The first few complete Bell polynomials are:

{\begin{aligned}B_{0}={}&1,\\[8pt]B_{1}(x_{1})={}&x_{1},\\[8pt]B_{2}(x_{1},x_{2})={}&x_{1}^{2}+x_{2},\\[8pt]B_{3}(x_{1},x_{2},x_{3})={}&x_{1}^{3}+3x_{1}x_{2}+x_{3},\\[8pt]B_{4}(x_{1},x_{2},x_{3},x_{4})={}&x_{1}^{4}+6x_{1}^{2}x_{2}+4x_{1}x_{3}+3x_{2}^{2}+x_{4},\\[8pt]B_{5}(x_{1},x_{2},x_{3},x_{4},x_{5})={}&x_{1}^{5}+10x_{2}x_{1}^{3}+15x_{2}^{2}x_{1}+10x_{3}x_{1}^{2}+10x_{3}x_{2}+5x_{4}x_{1}+x_{5}\\[8pt]B_{6}(x_{1},x_{2},x_{3},x_{4},x_{5},x_{6})={}&x_{1}^{6}+15x_{2}x_{1}^{4}+20x_{3}x_{1}^{3}+45x_{2}^{2}x_{1}^{2}+15x_{2}^{3}+60x_{3}x_{2}x_{1}\\&{}+15x_{4}x_{1}^{2}+10x_{3}^{2}+15x_{4}x_{2}+6x_{5}x_{1}+x_{6},\\[8pt]B_{7}(x_{1},x_{2},x_{3},x_{4},x_{5},x_{6},x_{7})={}&x_{1}^{7}+21x_{1}^{5}x_{2}+35x_{1}^{4}x_{3}+105x_{1}^{3}x_{2}^{2}+35x_{1}^{3}x_{4}\\&{}+210x_{1}^{2}x_{2}x_{3}+105x_{1}x_{2}^{3}+21x_{1}^{2}x_{5}+105x_{1}x_{2}x_{4}\\&{}+70x_{1}x_{3}^{2}+105x_{2}^{2}x_{3}+7x_{1}x_{6}+21x_{2}x_{5}+35x_{3}x_{4}+x_{7}.\end{aligned}}

Applications

Faà di Bruno's formula

Faà di Bruno's formula may be stated in terms of Bell polynomials as follows:

{d^{n} \over dx^{n}}f(g(x))=\sum _{{k=1}}^{n}f^{{(k)}}(g(x))B_{{n,k}}\left(g'(x),g''(x),\dots ,g^{{(n-k+1)}}(x)\right).

Similarly, a power-series version of Faà di Bruno's formula may be stated using Bell polynomials as follows. Suppose

f(x)=\sum _{{n=1}}^{\infty }{a_{n} \over n!}x^{n}\qquad {\mathrm {and}}\qquad g(x)=\sum _{{n=1}}^{\infty }{b_{n} \over n!}x^{n}.

Then

g(f(x))=\sum _{{n=1}}^{\infty }{\sum _{{k=1}}^{{n}}b_{k}B_{{n,k}}(a_{1},\dots ,a_{{n-k+1}}) \over n!}x^{n}.

In particular, the complete Bell polynomials appear in the exponential of a formal power series:

\exp \left(\sum _{i=1}^{\infty }{a_{i} \over i!}x^{i}\right)=\sum _{n=0}^{\infty }{B_{n}(a_{1},\dots ,a_{n}) \over n!}x^{n},

which also represents the exponential generating function of the complete Bell polynomials on a fixed sequence of arguments $a_{1},a_{2},\dots$ .

Reversion of series

Let two functions f and g be expressed in formal power series as

f(w)=\sum _{k=0}^{\infty }f_{k}{\frac {w^{k}}{k!}},\qquad \mathrm {and} \qquad g(z)=\sum _{k=0}^{\infty }g_{k}{\frac {z^{k}}{k!}},

such that g is the compositional inverse of f defined by g(f(w)) = w or f(g(z)) = z. If f₀ = 0 and f₁ ≠ 0, then an explicit form of the coefficients of the inverse can be given in term of Bell polynomials as ^[1]

g_{n}={\frac {1}{f_{1}^{n}}}\sum _{k=1}^{n-1}(-1)^{k}n^{\bar {k}}B_{n-1,k}({\hat {f}}_{1},{\hat {f}}_{2},\ldots ,{\hat {f}}_{n-k}),\qquad n\geq 2,

with ${\hat {f}}_{k}={\frac {f_{k+1}}{(k+1)f_{1}}},$ and $n^{\bar {k}}=n(n+1)\cdots (n+k-1)$ is the rising factorial, and $g_{1}={\frac {1}{f_{1}}}.$

Asymptotic expansion of Laplace-type integrals

Consider the integral of the form

I(\lambda )=\int _{a}^{b}e^{-\lambda f(x)}g(x)\mathrm {d} x,

where (a,b) is a real (finite or infinite) interval, λ is a large positive parameter and the functions f and g are continuous. Let f have a single minimum in [a,b] which occurs at x=a. Assume that as x→a⁺,

f(x)\sim f(a)+\sum _{k=0}^{\infty }a_{k}(x-a)^{k+\alpha },

g(x)\sim \sum _{k=0}^{\infty }b_{k}(x-a)^{k+\beta -1},

with α > 0, Re(β) > 0; and that the expansion of f can be term wise differentiated. Then, Laplace-Erdelyi theorem states that the asymptotic expansion of the integral I(λ) is given by

I(\lambda )\sim e^{-\lambda f(a)}\sum _{n=0}^{\infty }\Gamma {\Big (}{\frac {n+\beta }{\alpha }}{\Big )}{\frac {c_{n}}{\lambda ^{(n+\beta )/\alpha }}}\qquad \mathrm {as} \quad \lambda \rightarrow \infty ,

where the coefficients c_n are expressible in terms of a_n and b_n using partial ordinary Bell polynomials, as given by Campbell-Froman-Walles-Wojdylo formula:

c_{n}={\frac {1}{\alpha a_{0}^{(n+\beta )/\alpha }}}\sum _{k=0}^{n}b_{n-k}\sum _{j=0}^{k}{\binom {-{\frac {n+\beta }{\alpha }}}{j}}{\frac {1}{a_{0}^{j}}}{\hat {B}}_{k,j}(a_{1},a_{2},\ldots ,a_{k-j+1}).

Symmetric polynomials

The elementary symmetric polynomial $e_{n}$ and the power sum symmetric polynomial $p_{n}$ can be related to each other using Bell polynomials as:

e_{n}={\frac {(-1)^{n}}{n!}}B_{n}(p_{1},-1!p_{2},2!p_{3},-3!p_{4}\ldots ,(-1)^{n-1}(n-1)!p_{n}),

p_{n}={\frac {(-1)^{n-1}}{(n-1)!}}\sum _{k=1}^{n}(-1)^{k-1}(k-1)!B_{n,k}(e_{1},2!e_{2},3!e_{3},\ldots ,(n-k+1)!e_{n-k+1}).

This fact allows us to express the coefficients of monic polynomials in term the Bell polynomials of its roots. For instance, this can be applied with Cayley-Hamilton theorem to obtain the determinant of a n × n square matrix A in terms of its trace

\det(A)={\frac {(-1)^{n}}{n!}}B_{n}(s_{1},s_{2},\ldots ,s_{n}),~\qquad \mathrm {where} \qquad s_{k}=(-1)^{k-1}(k-1)!\;\mathrm {tr} (A^{k}).

Cycle index of symmetric groups

The cycle index of the symmetric group $S_{n}$ can be expressed in terms of complete Bell polynomials as follows:

Z(S_{n})={\frac {B_{n}(0!\,a_{1},1!\,a_{2},\dots ,(n-1)!\,a_{n})}{n!}}.

Moments and cumulants

The sum

\mu _{n}'=B_{n}(\kappa _{1},\dots ,\kappa _{n})=\sum _{k=1}^{n}B_{n,k}(\kappa _{1},\dots ,\kappa _{n-k+1})

is the nth raw moment of a probability distribution whose first n cumulants are κ₁, ..., κ_n. In other words, the nth moment is the nth complete Bell polynomial evaluated at the first n cumulants. Likewise, the nth cumulant can be given in terms of the moments as

\kappa _{n}=\sum _{k=1}^{n}(-1)^{k-1}(k-1)!B_{n,k}(\mu '_{1},\ldots ,\mu '_{n-k+1}).

Hermite polynomials

The probabilists' Hermite polynomials can be expressed in terms of Bell polynomials as

He_{n}(x)=B_{n}(x,-1,0,\ldots ,0),

where x_i = 0 for all i > 2; thus allowing for a combinatorial interpretation of the coefficients of the Hermite polynomials. This can be seen by comparing the generating function of the Hermite polynomials

\exp \left(xt-{\frac {t^{2}}{2}}\right)=\sum _{n=0}^{\infty }{\mathit {He}}_{n}(x){\frac {t^{n}}{n!}}

with that of Bell polynomials.

Representation of polynomial sequences of binomial type

For any sequence a₁, a₂, …, a_n of scalars, let

p_{n}(x)=\sum _{{k=1}}^{n}B_{{n,k}}(a_{1},\dots ,a_{{n-k+1}})x^{k}.

Then this polynomial sequence is of binomial type, i.e. it satisfies the binomial identity

p_{n}(x+y)=\sum _{{k=0}}^{n}{n \choose k}p_{k}(x)p_{{n-k}}(y).

Example: For a₁ = … = a_n = 1, the polynomials

p_{n}(x)

represent Touchard polynomials.

More generally, we have this result:

Theorem: All polynomial sequences of binomial type are of this form.

If we define a formal power series

h(x)=\sum _{{k=1}}^{\infty }{a_{k} \over k!}x^{k},

then for all n,

h^{{-1}}\left({d \over dx}\right)p_{n}(x)=np_{{n-1}}(x).

Software

Bell polynomials are implemented in:

Notes

↑ Eqn (11.43), p. 437, C.A. Charalambides, Enumerative Combinatorics, Chapman & Hall / CRC, 2002

References

Eric Temple Bell (1927–1928). "Partition Polynomials". Annals of Mathematics. 29 (1/4): 38–46. JSTOR 1967979. MR 1502817. doi:10.2307/1967979.
Louis Comtet (1974). Advanced Combinatorics: The Art of Finite and Infinite Expansions. Dordrecht, Holland / Boston, U.S.: Reidel Publishing Company.
Steven Roman. The Umbral Calculus. Dover Publications.
Vassily G. Voinov, Mikhail S. Nikulin (1994). "On power series, Bell polynomials, Hardy-Ramanujan-Rademacher problem and its statistical applications". Kybernetika. 30 (3): 343–358. ISSN 0023-5954.
Andrews, George E. (1998). The Theory of Partitions. Cambridge Mathematical Library (1st pbk ed.). Cambridge University Press. pp. 204–211. ISBN 0-521-63766-X.
Silvia Noschese, Paolo E. Ricci (2003). "Differentiation of Multivariable Composite Functions and Bell Polynomials". Journal of Computational Analysis and Applications. 5 (3): 333–340. doi:10.1023/A:1023227705558.
Abbas, Moncef; Bouroubi, Sadek (2005). "On new identities for Bell's polynomial". Disc. Math (293): 5–10. MR 2136048. doi:10.1016/j.disc.2004.08.023.
Khristo N. Boyadzhiev (2009). "Exponential Polynomials, Stirling Numbers, and Evaluation of Some Gamma Integrals". Abstract and Applied Analysis. 2009: Article ID 168672. doi:10.1155/2009/168672. (contains also elementary review of the concept Bell-polynomials)
V. V. Kruchinin (2011). "Derivation of Bell Polynomials of the Second Kind". arXiv:1104.5065 .
Griffiths, Martin (2012). "Families of sequences from a class of multinomial sums". Journal of Integer Sequences. 15. MR 2872465.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

Bell polynomials