The Bell number ,
for a nonnegative integer
, counts the ways a set of
elements can be partitioned
into nonempty subsets
(not to be confused with the Bernoulli number,
which is also commonly denoted
).
Equivalently, in the terminology of the twelvefold way,
is the number of ways of distributing
distinct balls into an unrestricted number of nonempty identical
boxes.
For example, there are five ways the numbers can be partitioned:
,
,
,
, and
, so
.
, and the first few Bell numbers
for
, 2, ... are 1, 2, 5, 15, 52, 203,
877, 4140, 21147, 115975, ... (OEIS A000110).
The numbers of digits in
for
, 1, ... are given by 1, 6, 116, 1928, 27665, ... (OEIS A113015).
Bell numbers are implemented in the Wolfram Language as BellB[n].
Though Bell numbers have traditionally been attributed to E. T. Bell as a result of the general theory he developed in his 1934 paper (Bell 1934), the first systematic study of Bell numbers was made by Ramanujan in chapter 3 of his second notebook approximately 25-30 years prior to Bell's work (B. C. Berndt, pers. comm., Jan. 4 and 13, 2010).
The first few prime Bell numbers occur at indices , 3, 7, 13, 42, 55, 2841, ... (OEIS A051130),
with no others less than
(Weisstein, Apr. 23, 2006). These correspond to the numbers 2, 5, 877, 27644437,
... (OEIS A051131).
was proved prime by I. Larrosa Canestro in 2004
after 17 months of computation using the elliptic
curve primality proving program PRIMO.
Bell numbers are closely related to Catalan numbers. The diagram above shows the constructions giving and
, with line segments representing elements in the same
subset and dots representing subsets
containing a single element (Dickau). The integers
can be defined by the sum
|
(1)
|
where
is a Stirling number of the second
kind, i.e., as the Stirling transform of
the sequence 1, 1, 1, ....
For , the Bell numbers are given in
terms of generalized hypergeometric
functions by
|
(2)
|
(K. A. Penson, pers. comm., Jan. 14, 2007).
The Bell numbers can also be generated using the sum and recurrence relation
|
(3)
|
where
is a binomial coefficient, using the formula
of Comtet (1974)
|
(4)
|
for , where
denotes the ceiling function.
Dobiński's formula gives the
th Bell number for
,
|
(5)
|
For , a variation of Dobiński's
formula gives
|
(6)
| |||
|
(7)
|
where
is a subfactorial (Pitman 1997).
For , a double sum
is given by
|
(8)
|
The ordinary generating function is a formal
power series with radius of convergence
zero. For
,
the function
|
(9)
| |||
|
(10)
| |||
|
(11)
| |||
|
(12)
| |||
|
(13)
|
has the Bell numbers as its asymptotic coefficients as . The exponential
generating function is
|
(14)
|
An amazing integral representation for was given by Cesàro (1885),
|
(15)
| |||
|
(16)
|
(Becker and Browne 1941, Callan 2005) for , where
denotes the imaginary part
of
.
The Bell number
is also equal to
,
where
is a Bell polynomial.
de Bruijn (1981) gave the asymptotic formula
|
(17)
|
Lovász (1993) showed that this formula gives the asymptotic limit
|
(18)
|
where
and
is the Lambert
W-function (Graham et al. 1994, p. 493). Odlyzko (1995) gave
|
(19)
|
Touchard's congruence states
|
(20)
|
when is prime
and
is a nonnegative
integer. This gives as a special case for
the congruence
|
(21)
|
for prime.
For every prime number
, the congruence
|
(22)
|
holds for all nonnegative integers (Lunnon et al. 1979, Montgomery et al. 2010).
It is conjectured that
is always the least period modulo
. The sequence of Bell numbers
is periodic modulo every
positive integer
(Levine and Dalton 1962, Lunnon et al. 1979). The least periods for moduli
, 2, ... are 1, 3, 13, 12, 781, 39, 137257, 24, 39, 2343,
28531167061, 156, ... (OEIS A054767).
The Bell numbers also have the curious property that
|
(23)
| |||
|
(24)
|
(Lenard 1992; Radoux 1992, p. 6; Ehrenborg 2000, pp. 557-558), where the product is simply a superfactorial
and is a Barnes
G-function, the first few of which for
, 1, 2, ... are 1, 1, 2, 12, 288, 34560, 24883200, ... (OEIS
A000178).