cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 76 results. Next

A029706 Sum C(n,k)*b(k), k=1..n, where b(k) is given by A001861.

Original entry on oeis.org

2, 10, 46, 226, 1214, 7106, 44958, 305090, 2206398, 16913986, 136823262, 1163490498, 10366252030, 96491364674, 935976996126, 9440144423874, 98800604237118, 1071092025420866, 12008090971866206, 139014305916844738
Offset: 0

Views

Author

Andre Poenitz (andre.poenitz(AT)mathematik.tu-chemnitz.de)

Keywords

Crossrefs

Arises if one of the two kinds of boxes mentioned in A001861 may 'fail'.

Programs

  • Mathematica
    a[n_] := Sum[Binomial[n, k]*BellB[k, 2], {k, 1, n}]; Table[a[n], {n, 1, 20}] (* Jean-François Alcover, Apr 26 2013 *)
  • PARI
    x='x+O('x^66); Vec(serlaplace(1/2-exp(x)+1/2*exp(2*exp(x)-2))) \\ Joerg Arndt, Apr 21 2013
    
  • PARI
    a(n) = sum(k=1, n, binomial(n, k)*sum(j=1, k, stirling(k, j, 2)*2^j)); \\ \\ Michel Marcus, Apr 26 2013

Formula

E.g.f.: 1/2 - exp(x) + 1/2*exp(2*exp(x)-2). a(n-2) = A035009(n)-1. - Ralf Stephan, Jan 26 2004

Extensions

Added missing digit in the last term from Jean-François Alcover, Apr 26 2013

A191643 Incorrect version of A001861.

Original entry on oeis.org

2, 8, 22, 94, 454, 2430, 14214, 89918, 610182, 4412958
Offset: 1

Views

Author

Keywords

References

  • Carlier, Jacques; and Lucet, Corinne; A decomposition algorithm for network reliability evaluation. In First International Colloquium on Graphs and Optimization (GOI), 1992 (Grimentz). Discrete Appl. Math. 65 (1996), 141-156.

Crossrefs

Cf. A191644.

Programs

  • PARI
    a(n) = sum(j=1, k, stirling(k,j,2)*2^(j)) \\ Michel Marcus, Apr 25 2013

Formula

There is an explicit formula in the reference (see p. 152).

A000110 Bell or exponential numbers: number of ways to partition a set of n labeled elements.

Original entry on oeis.org

1, 1, 2, 5, 15, 52, 203, 877, 4140, 21147, 115975, 678570, 4213597, 27644437, 190899322, 1382958545, 10480142147, 82864869804, 682076806159, 5832742205057, 51724158235372, 474869816156751, 4506715738447323, 44152005855084346, 445958869294805289, 4638590332229999353, 49631246523618756274
Offset: 0

Views

Author

Keywords

Comments

The leading diagonal of its difference table is the sequence shifted, see Bernstein and Sloane (1995). - N. J. A. Sloane, Jul 04 2015
Also the number of equivalence relations that can be defined on a set of n elements. - Federico Arboleda (federico.arboleda(AT)gmail.com), Mar 09 2005
a(n) = number of nonisomorphic colorings of a map consisting of a row of n+1 adjacent regions. Adjacent regions cannot have the same color. - David W. Wilson, Feb 22 2005
If an integer is squarefree and has n distinct prime factors then a(n) is the number of ways of writing it as a product of its divisors. - Amarnath Murthy, Apr 23 2001
Consider rooted trees of height at most 2. Letting each tree 'grow' into the next generation of n means we produce a new tree for every node which is either the root or at height 1, which gives the Bell numbers. - Jon Perry, Jul 23 2003
Begin with [1,1] and follow the rule that [1,k] -> [1,k+1] and [1,k] k times, e.g., [1,3] is transformed to [1,4], [1,3], [1,3], [1,3]. Then a(n) is the sum of all components: [1,1] = 2; [1,2], [1,1] = 5; [1,3], [1,2], [1,2], [1,2], [1,1] = 15; etc. - Jon Perry, Mar 05 2004
Number of distinct rhyme schemes for a poem of n lines: a rhyme scheme is a string of letters (e.g., 'abba') such that the leftmost letter is always 'a' and no letter may be greater than one more than the greatest letter to its left. Thus 'aac' is not valid since 'c' is more than one greater than 'a'. For example, a(3)=5 because there are 5 rhyme schemes: aaa, aab, aba, abb, abc; also see example by Neven Juric. - Bill Blewett, Mar 23 2004
In other words, number of length-n restricted growth strings (RGS) [s(0),s(1),...,s(n-1)] where s(0)=0 and s(k) <= 1 + max(prefix) for k >= 1, see example (cf. A080337 and A189845). - Joerg Arndt, Apr 30 2011
Number of partitions of {1, ..., n+1} into subsets of nonconsecutive integers, including the partition 1|2|...|n+1. E.g., a(3)=5: there are 5 partitions of {1,2,3,4} into subsets of nonconsecutive integers, namely, 13|24, 13|2|4, 14|2|3, 1|24|3, 1|2|3|4. - Augustine O. Munagi, Mar 20 2005
Triangle (addition) scheme to produce terms, derived from the recurrence, from Oscar Arevalo (loarevalo(AT)sbcglobal.net), May 11 2005:
1
1 2
2 3 5
5 7 10 15
15 20 27 37 52
... [This is Aitken's array A011971]
With P(n) = the number of integer partitions of n, p(i) = the number of parts of the i-th partition of n, d(i) = the number of different parts of the i-th partition of n, p(j,i) = the j-th part of the i-th partition of n, m(i,j) = multiplicity of the j-th part of the i-th partition of n, one has: a(n) = Sum_{i=1..P(n)} (n!/(Product_{j=1..p(i)} p(i,j)!)) * (1/(Product_{j=1..d(i)} m(i,j)!)). - Thomas Wieder, May 18 2005
a(n+1) is the number of binary relations on an n-element set that are both symmetric and transitive. - Justin Witt (justinmwitt(AT)gmail.com), Jul 12 2005
If the rule from Jon Perry, Mar 05 2004, is used, then a(n-1) = [number of components used to form a(n)] / 2. - Daniel Kuan (dkcm(AT)yahoo.com), Feb 19 2006
a(n) is the number of functions f from {1,...,n} to {1,...,n,n+1} that satisfy the following two conditions for all x in the domain: (1) f(x) > x; (2) f(x)=n+1 or f(f(x))=n+1. E.g., a(3)=5 because there are exactly five functions that satisfy the two conditions: f1={(1,4),(2,4),(3,4)}, f2={(1,4),(2,3),(3,4)}, f3={(1,3),(2,4),(3,4)}, f4={(1,2),(2,4),(3,4)} and f5={(1,3),(2,3),(3,4)}. - Dennis P. Walsh, Feb 20 2006
Number of asynchronic siteswap patterns of length n which have no zero-throws (i.e., contain no 0's) and whose number of orbits (in the sense given by Allen Knutson) is equal to the number of balls. E.g., for n=4, the condition is satisfied by the following 15 siteswaps: 4444, 4413, 4242, 4134, 4112, 3441, 2424, 1344, 2411, 1313, 1241, 2222, 3131, 1124, 1111. Also number of ways to choose n permutations from identity and cyclic permutations (1 2), (1 2 3), ..., (1 2 3 ... n) so that their composition is identity. For n=3 we get the following five: id o id o id, id o (1 2) o (1 2), (1 2) o id o (1 2), (1 2) o (1 2) o id, (1 2 3) o (1 2 3) o (1 2 3). (To see the bijection, look at Ehrenborg and Readdy paper.) - Antti Karttunen, May 01 2006
a(n) is the number of permutations on [n] in which a 3-2-1 (scattered) pattern occurs only as part of a 3-2-4-1 pattern. Example: a(3) = 5 counts all permutations on [3] except 321. See "Eigensequence for Composition" reference a(n) = number of permutation tableaux of size n (A000142) whose first row contains no 0's. Example: a(3)=5 counts {{}, {}, {}}, {{1}, {}}, {{1}, {0}}, {{1}, {1}}, {{1, 1}}. - David Callan, Oct 07 2006
From Gottfried Helms, Mar 30 2007: (Start)
This sequence is also the first column in the matrix-exponential of the (lower triangular) Pascal-matrix, scaled by exp(-1): PE = exp(P) / exp(1) =
1
1 1
2 2 1
5 6 3 1
15 20 12 4 1
52 75 50 20 5 1
203 312 225 100 30 6 1
877 1421 1092 525 175 42 7 1
First 4 columns are A000110, A033306, A105479, A105480. The general case is mentioned in the two latter entries. PE is also the Hadamard-product Toeplitz(A000110) (X) P:
1
1 1
2 1 1
5 2 1 1
15 5 2 1 1 (X) P
52 15 5 2 1 1
203 52 15 5 2 1 1
877 203 52 15 5 2 1 1
(End)
The terms can also be computed with finite steps and precise integer arithmetic. Instead of exp(P)/exp(1) one can compute A = exp(P - I) where I is the identity-matrix of appropriate dimension since (P-I) is nilpotent to the order of its dimension. Then a(n)=A[n,1] where n is the row-index starting at 1. - Gottfried Helms, Apr 10 2007
When n is prime, a(n) == 2 (mod n), but the converse is not always true. Define a Bell pseudoprime to be a composite number n such that a(n) == 2 (mod n). W. F. Lunnon recently found the Bell pseudoprimes 21361 = 41*521 and C46 = 3*23*16218646893090134590535390526854205539989357 and conjectured that Bell pseudoprimes are extremely scarce. So the second Bell pseudoprime is unlikely to be known with certainty in the near future. I confirmed that 21361 is the first. - David W. Wilson, Aug 04 2007 and Sep 24 2007
This sequence and A000587 form a reciprocal pair under the list partition transform described in A133314. - Tom Copeland, Oct 21 2007
Starting (1, 2, 5, 15, 52, ...), equals row sums and right border of triangle A136789. Also row sums of triangle A136790. - Gary W. Adamson, Jan 21 2008
This is the exponential transform of A000012. - Thomas Wieder, Sep 09 2008
From Abdullahi Umar, Oct 12 2008: (Start)
a(n) is also the number of idempotent order-decreasing full transformations (of an n-chain).
a(n) is also the number of nilpotent partial one-one order-decreasing transformations (of an n-chain).
a(n+1) is also the number of partial one-one order-decreasing transformations (of an n-chain). (End)
From Peter Bala, Oct 19 2008: (Start)
Bell(n) is the number of n-pattern sequences [Cooper & Kennedy]. An n-pattern sequence is a sequence of integers (a_1,...,a_n) such that a_i = i or a_i = a_j for some j < i. For example, Bell(3) = 5 since the 3-pattern sequences are (1,1,1), (1,1,3), (1,2,1), (1,2,2) and (1,2,3).
Bell(n) is the number of sequences of positive integers (N_1,...,N_n) of length n such that N_1 = 1 and N_(i+1) <= 1 + max{j = 1..i} N_j for i >= 1 (see the comment by B. Blewett above). It is interesting to note that if we strengthen the latter condition to N_(i+1) <= 1 + N_i we get the Catalan numbers A000108 instead of the Bell numbers.
(End)
Equals the eigensequence of Pascal's triangle, A007318; and starting with offset 1, = row sums of triangles A074664 and A152431. - Gary W. Adamson, Dec 04 2008
The entries f(i, j) in the exponential of the infinite lower-triangular matrix of binomial coefficients b(i, j) are f(i, j) = b(i, j) e a(i - j). - David Pasino, Dec 04 2008
Equals lim_{k->oo} A071919^k. - Gary W. Adamson, Jan 02 2009
Equals A154107 convolved with A014182, where A014182 = expansion of exp(1-x-exp(-x)), the eigensequence of A007318^(-1). Starting with offset 1 = A154108 convolved with (1,2,3,...) = row sums of triangle A154109. - Gary W. Adamson, Jan 04 2009
Repeated iterates of (binomial transform of [1,0,0,0,...]) will converge upon (1, 2, 5, 15, 52, ...) when each result is prefaced with a "1"; such that the final result is the fixed limit: (binomial transform of [1,1,2,5,15,...]) = (1,2,5,15,52,...). - Gary W. Adamson, Jan 14 2009
From Karol A. Penson, May 03 2009: (Start)
Relation between the Bell numbers B(n) and the n-th derivative of 1/Gamma(1+x) evaluated at x=1:
a) produce a number of such derivatives through seq(subs(x=1, simplify((d^n/dx^n)GAMMA(1+x)^(-1))), n=1..5);
b) leave them expressed in terms of digamma (Psi(k)) and polygamma (Psi(k,n)) functions and unevaluated;
Examples of such expressions, for n=1..5, are:
n=1: -Psi(1),
n=2: -(-Psi(1)^2 + Psi(1,1)),
n=3: -Psi(1)^3 + 3*Psi(1)*Psi(1,1) - Psi(2,1),
n=4: -(-Psi(1)^4 + 6*Psi(1)^2*Psi(1,1) - 3*Psi(1,1)^2 - 4*Psi(1)*Psi(2,1) + Psi(3, 1)),
n=5: -Psi(1)^5 + 10*Psi(1)^3*Psi(1,1) - 15*Psi(1)*Psi(1,1)^2 - 10*Psi(1)^2*Psi(2,1) + 10*Psi(1,1)*Psi(2,1) + 5*Psi(1)*Psi(3,1) - Psi(4,1);
c) for a given n, read off the sum of absolute values of coefficients of every term involving digamma or polygamma functions.
This sum is equal to B(n). Examples: B(1)=1, B(2)=1+1=2, B(3)=1+3+1=5, B(4)=1+6+3+4+1=15, B(5)=1+10+15+10+10+5+1=52;
d) Observe that this decomposition of the Bell number B(n) apparently does not involve the Stirling numbers of the second kind explicitly.
(End)
The numbers given above by Penson lead to the multinomial coefficients A036040. - Johannes W. Meijer, Aug 14 2009
Column 1 of A162663. - Franklin T. Adams-Watters, Jul 09 2009
Asymptotic expansions (0!+1!+2!+...+(n-1)!)/(n-1)! = a(0) + a(1)/n + a(2)/n^2 + ... and (0!+1!+2!+...+n!)/n! = 1 + a(0)/n + a(1)/n^2 + a(2)/n^3 + .... - Michael Somos, Jun 28 2009
Starting with offset 1 = row sums of triangle A165194. - Gary W. Adamson, Sep 06 2009
a(n+1) = A165196(2^n); where A165196 begins: (1, 2, 4, 5, 7, 12, 14, 15, ...). such that A165196(2^3) = 15 = A000110(4). - Gary W. Adamson, Sep 06 2009
The divergent series g(x=1,m) = 1^m*1! - 2^m*2! + 3^m*3! - 4^m*4! + ..., m >= -1, which for m=-1 dates back to Euler, is related to the Bell numbers. We discovered that g(x=1,m) = (-1)^m * (A040027(m) - A000110(m+1) * A073003). We observe that A073003 is Gompertz's constant and that A040027 was published by Gould, see for more information A163940. - Johannes W. Meijer, Oct 16 2009
a(n) = E(X^n), i.e., the n-th moment about the origin of a random variable X that has a Poisson distribution with (rate) parameter, lambda = 1. - Geoffrey Critzer, Nov 30 2009
Let A000110 = S(x), then S(x) = A(x)/A(x^2) when A(x) = A173110; or (1, 1, 2, 5, 15, 52, ...) = (1, 1, 3, 6, 20, 60, ...) / (1, 0, 1, 0, 3, 0, 6, 0, 20, ...). - Gary W. Adamson, Feb 09 2010
The Bell numbers serve as the upper limit for the number of distinct homomorphic images from any given finite universal algebra. Every algebra homomorphism is determined by its kernel, which must be a congruence relation. As the number of possible congruence relations with respect to a finite universal algebra must be a subset of its possible equivalence classes (given by the Bell numbers), it follows naturally. - Max Sills, Jun 01 2010
For a proof of the o.g.f. given in the R. Stephan comment see, e.g., the W. Lang link under A071919. - Wolfdieter Lang, Jun 23 2010
Let B(x) = (1 + x + 2x^2 + 5x^3 + ...). Then B(x) is satisfied by A(x)/A(x^2) where A(x) = polcoeff A173110: (1 + x + 3x^2 + 6x^3 + 20x^4 + 60x^5 + ...) = B(x) * B(x^2) * B(x^4) * B(x^8) * .... - Gary W. Adamson, Jul 08 2010
Consider a set of A000217(n) balls of n colors in which, for each integer k = 1 to n, exactly one color appears in the set a total of k times. (Each ball has exactly one color and is indistinguishable from other balls of the same color.) a(n+1) equals the number of ways to choose 0 or more balls of each color without choosing any two colors the same positive number of times. (See related comments for A000108, A008277, A016098.) - Matthew Vandermast, Nov 22 2010
A binary counter with faulty bits starts at value 0 and attempts to increment by 1 at each step. Each bit that should toggle may or may not do so. a(n) is the number of ways that the counter can have the value 0 after n steps. E.g., for n=3, the 5 trajectories are 0,0,0,0; 0,1,0,0; 0,1,1,0; 0,0,1,0; 0,1,3,0. - David Scambler, Jan 24 2011
No Bell number is divisible by 8, and no Bell number is congruent to 6 modulo 8; see Theorem 6.4 and Table 1.7 in Lunnon, Pleasants and Stephens. - Jon Perry, Feb 07 2011, clarified by Eric Rowland, Mar 26 2014
a(n+1) is the number of (symmetric) positive semidefinite n X n 0-1 matrices. These correspond to equivalence relations on {1,...,n+1}, where matrix element M[i,j] = 1 if and only if i and j are equivalent to each other but not to n+1. - Robert Israel, Mar 16 2011
a(n) is the number of monotonic-labeled forests on n vertices with rooted trees of height less than 2. We note that a labeled rooted tree is monotonic-labeled if the label of any parent vertex is greater than the label of any offspring vertex. See link "Counting forests with Stirling and Bell numbers". - Dennis P. Walsh, Nov 11 2011
a(n) = D^n(exp(x)) evaluated at x = 0, where D is the operator (1+x)*d/dx. Cf. A000772 and A094198. - Peter Bala, Nov 25 2011
B(n) counts the length n+1 rhyme schemes without repetitions. E.g., for n=2 there are 5 rhyme schemes of length 3 (aaa, aab, aba, abb, abc), and the 2 without repetitions are aba, abc. This is basically O. Munagi's result that the Bell numbers count partitions into subsets of nonconsecutive integers (see comment above dated Mar 20 2005). - Eric Bach, Jan 13 2012
Right and left borders and row sums of A212431 = A000110 or a shifted variant. - Gary W. Adamson, Jun 21 2012
Number of maps f: [n] -> [n] where f(x) <= x and f(f(x)) = f(x) (projections). - Joerg Arndt, Jan 04 2013
Permutations of [n] avoiding any given one of the 8 dashed patterns in the equivalence classes (i) 1-23, 3-21, 12-3, 32-1, and (ii) 1-32, 3-12, 21-3, 23-1. (See Claesson 2001 reference.) - David Callan, Oct 03 2013
Conjecture: No a(n) has the form x^m with m > 1 and x > 1. - Zhi-Wei Sun, Dec 02 2013
Sum_{n>=0} a(n)/n! = e^(e-1) = 5.57494152476..., see A234473. - Richard R. Forberg, Dec 26 2013 (This is the e.g.f. for x=1. - Wolfdieter Lang, Feb 02 2015)
Sum_{j=0..n} binomial(n,j)*a(j) = (1/e)*Sum_{k>=0} (k+1)^n/k! = (1/e) Sum_{k=1..oo} k^(n+1)/k! = a(n+1), n >= 0, using the Dobinski formula. See the comment by Gary W. Adamson, Dec 04 2008 on the Pascal eigensequence. - Wolfdieter Lang, Feb 02 2015
In fact it is not really an eigensequence of the Pascal matrix; rather the Pascal matrix acts on the sequence as a shift. It is an eigensequence (the unique eigensequence with eigenvalue 1) of the matrix derived from the Pascal matrix by adding at the top the row [1, 0, 0, 0 ...]. The binomial sum formula may be derived from the definition in terms of partitions: label any element X of a set S of N elements, and let X(k) be the number of subsets of S containing X with k elements. Since each subset has a unique coset, the number of partitions p(N) of S is given by p(N) = Sum_{k=1..N} (X(k) p(N-k)); trivially X(k) = N-1 choose k-1. - Mason Bogue, Mar 20 2015
a(n) is the number of ways to nest n matryoshkas (Russian nesting dolls): we may identify {1, 2, ..., n} with dolls of ascending sizes and the sets of a set partition with stacks of dolls. - Carlo Sanna, Oct 17 2015
Number of permutations of [n] where the initial elements of consecutive runs of increasing elements are in decreasing order. a(4) = 15: `1234, `2`134, `23`14, `234`1, `24`13, `3`124, `3`2`14, `3`24`1, `34`12, `34`2`1, `4`123, `4`2`13, `4`23`1, `4`3`12, `4`3`2`1. - Alois P. Heinz, Apr 27 2016
Taking with alternating signs, the Bell numbers are the coefficients in the asymptotic expansion (Ramanujan): (-1)^n*(A000166(n) - n!/exp(1)) ~ 1/n - 2/n^2 + 5/n^3 - 15/n^4 + 52/n^5 - 203/n^6 + O(1/n^7). - Vladimir Reshetnikov, Nov 10 2016
Number of treeshelves avoiding pattern T231. See A278677 for definitions and examples. - Sergey Kirgizov, Dec 24 2016
Presumably this satisfies Benford's law, although the results in Hürlimann (2009) do not make this clear. - N. J. A. Sloane, Feb 09 2017
a(n) = Sum(# of standard immaculate tableaux of shape m, m is a composition of n), where this sum is over all integer compositions m of n > 0. This formula is easily seen to hold by identifying standard immaculate tableaux of size n with set partitions of { 1, 2, ..., n }. For example, if we sum over integer compositions of 4 lexicographically, we see that 1+1+2+1+3+3+3+1 = 15 = A000110(4). - John M. Campbell, Jul 17 2017
a(n) is also the number of independent vertex sets (and vertex covers) in the (n-1)-triangular honeycomb bishop graph. - Eric W. Weisstein, Aug 10 2017
Even-numbered entries represent the numbers of configurations of identity and non-identity for alleles of a gene in n diploid individuals with distinguishable maternal and paternal alleles. - Noah A Rosenberg, Jan 28 2019
Number of partial equivalence relations (PERs) on a set with n elements (offset=1), i.e., number of symmetric, transitive (not necessarily reflexive) relations. The idea is to add a dummy element D to the set, and then take equivalence relations on the result; anything equivalent to D is then removed for the partial equivalence relation. - David Spivak, Feb 06 2019
Number of words of length n+1 with no repeated letters, when letters are unlabeled. - Thomas Anton, Mar 14 2019
Named by Becker and Riordan (1948) after the Scottish-American mathematician and writer Eric Temple Bell (1883 - 1960). - Amiram Eldar, Dec 04 2020
Also the number of partitions of {1,2,...,n+1} with at most one n+1 singleton. E.g., a(3)=5: {13|24, 12|34, 123|4, 14|23, 1234}. - Yuchun Ji, Dec 21 2020
a(n) is the number of sigma algebras on the set of n elements. Note that each sigma algebra is generated by a partition of the set. For example, the sigma algebra generated by the partition {{1}, {2}, {3,4}} is {{}, {1}, {2}, {1,2}, {3,4}, {1,3,4}, {2,3,4}, {1,2,3,4}}. - Jianing Song, Apr 01 2021
a(n) is the number of P_3-free graphs on n labeled nodes. - M. Eren Kesim, Jun 04 2021
a(n) is the number of functions X:([n] choose 2) -> {+,-} such that for any ordered 3-tuple abc we have X(ab)X(ac)X(bc) not in {+-+,++-,-++}. - Robert Lauff, Dec 09 2022
From Manfred Boergens, Mar 11 2025: (Start)
The partitions in the definition can be described as disjoint covers of the set. "Covers" in general give rise to the following amendments:
For disjoint covers which may include one empty set see A186021.
For arbitrary (including non-disjoint) covers see A003465.
For arbitrary (including non-disjoint) covers which may include one empty set see A000371. (End)

Examples

			G.f. = 1 + x + 2*x^2 + 5*x^3 + 15*x^4 + 52*x^5 + 203*x^6 + 877*x^7 + 4140*x^8 + ...
From Neven Juric, Oct 19 2009: (Start)
The a(4)=15 rhyme schemes for n=4 are
  aaaa, aaab, aaba, aabb, aabc, abaa, abab, abac, abba, abbb, abbc, abca, abcb, abcc, abcd
The a(5)=52 rhyme schemes for n=5 are
  aaaaa, aaaab, aaaba, aaabb, aaabc, aabaa, aabab, aabac, aabba, aabbb, aabbc, aabca, aabcb, aabcc, aabcd, abaaa, abaab, abaac, ababa, ababb, ababc, abaca, abacb, abacc, abacd, abbaa, abbab, abbac, abbba, abbbb, abbbc, abbca, abbcb, abbcc, abbcd, abcaa, abcab, abcac, abcad, abcba, abcbb, abcbc, abcbd, abcca, abccb, abccc, abccd, abcda, abcdb, abcdc, abcdd, abcde
(End)
From _Joerg Arndt_, Apr 30 2011: (Start)
Restricted growth strings (RGS):
For n=0 there is one empty string;
for n=1 there is one string [0];
for n=2 there are 2 strings [00], [01];
for n=3 there are a(3)=5 strings [000], [001], [010], [011], and [012];
for n=4 there are a(4)=15 strings
1: [0000], 2: [0001], 3: [0010], 4: [0011], 5: [0012], 6: [0100], 7: [0101], 8: [0102], 9: [0110], 10: [0111], 11: [0112], 12: [0120], 13: [0121], 14: [0122], 15: [0123].
These are one-to-one with the rhyme schemes (identify a=0, b=1, c=2, etc.).
(End)
Consider the set S = {1, 2, 3, 4}. The a(4) = 1 + 3 + 6 + 4 + 1 = 15 partitions are: P1 = {{1}, {2}, {3}, {4}}; P21 .. P23 = {{a,4}, S\{a,4}} with a = 1, 2, 3; P24 .. P29 = {{a}, {b}, S\{a,b}} with 1 <= a < b <= 4;  P31 .. P34 = {S\{a}, {a}} with a = 1 .. 4; P4 = {S}. See the Bottomley link for a graphical illustration. - _M. F. Hasler_, Oct 26 2017
		

References

  • Stefano Aguzzoli, Brunella Gerla and Corrado Manara, Poset Representation for Goedel and Nilpotent Minimum Logics, in Symbolic and Quantitative Approaches to Reasoning with Uncertainty, Lecture Notes in Computer Science, Volume 3571/2005, Springer-Verlag. [Added by N. J. A. Sloane, Jul 08 2009]
  • S. Ainley, Problem 19, QARCH, No. IV, Nov 03, 1980.
  • J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003, p. 205.
  • W. Asakly, A. Blecher, C. Brennan, A. Knopfmacher, T. Mansour, S. Wagner, Set partition asymptotics and a conjecture of Gould and Quaintance, Journal of Mathematical Analysis and Applications, Volume 416, Issue 2, Aug 15 2014, Pages 672-682.
  • J. Balogh, B. Bollobas and D. Weinreich, A jump to the Bell numbers for hereditary graph properties, J. Combin. Theory Ser. B 95 (2005), no. 1, 29-48.
  • R. E. Beard, On the coefficients in the expansion of exp(exp(t)) and exp(-exp(t)), J. Institute Actuaries, 76 (1951), 152-163.
  • H. W. Becker, Abstracts of two papers related to Bell numbers, Bull. Amer. Math. Soc., 52 (1946), p. 415.
  • E. T. Bell, The iterated exponential numbers, Ann. Math., 39 (1938), 539-557.
  • C. M. Bender, D. C. Brody and B. K. Meister, Quantum Field Theory of Partitions, J. Math. Phys., 40,7 (1999), 3239-45.
  • E. A. Bender and J. R. Goldman, Enumerative uses of generating functions, Indiana Univ. Math. J., 20 (1971), 753-765.
  • G. Birkhoff, Lattice Theory, Amer. Math. Soc., Revised Ed., 1961, p. 108, Ex. 1.
  • M. T. L. Bizley, On the coefficients in the expansion of exp(lambda exp(t)), J. Inst. Actuaries, 77 (1952), p. 122.
  • J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 41.
  • Carlier, Jacques; and Lucet, Corinne; A decomposition algorithm for network reliability evaluation. In First International Colloquium on Graphs and Optimization (GOI), 1992 (Grimentz). Discrete Appl. Math. 65 (1996), 141-156.
  • Anders Claesson, Generalized Pattern Avoidance, European Journal of Combinatorics, 22 (2001) 961-971.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 210.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 92-93.
  • John H. Conway et al., The Symmetries of Things, Peters, 2008, p. 207.
  • Colin Defant, Highly sorted permutations and Bell numbers, ECA 1:1 (2021) Article S2R6.
  • De Angelis, Valerio, and Dominic Marcello. "Wilf's Conjecture." The American Mathematical Monthly 123.6 (2016): 557-573.
  • N. G. de Bruijn, Asymptotic Methods in Analysis, Dover, 1981, Sections 3.3. Case b and 6.1-6.3.
  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 52, p. 19, Ellipses, Paris 2008.
  • G. Dobinski, Summierung der Reihe Sum(n^m/n!) für m = 1, 2, 3, 4, 5, ..., Grunert Archiv (Arch. f. Math. und Physik), 61 (1877) 333-336.
  • L. F. Epstein, A function related to the series for exp(exp(z)), J. Math. and Phys., 18 (1939), 153-173.
  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 5.8, p. 321.
  • Flajolet, Philippe and Schott, Rene, Nonoverlapping partitions, continued fractions, Bessel functions and a divergent series, European J. Combin. 11 (1990), no. 5, 421-432.
  • Martin Gardner, Fractal Music, Hypercards and More (Freeman, 1992), Chapter 2.
  • H. W. Gould, Research bibliography of two special number sequences, Mathematica Monongaliae, Vol. 12, 1971.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, 2nd ed., p. 493.
  • Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
  • M. Kauers and P. Paule, The Concrete Tetrahedron, Springer 2011, p. 26.
  • D. E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, Section 7.2.1.5 (p. 418).
  • Christian Kramp, Der polynomische Lehrsatz (Leipzig: 1796), 113.
  • Lehmer, D. H. Some recursive sequences. Proceedings of the Manitoba Conference on Numerical Mathematics (Univ. Manitoba, Winnipeg, Man., 1971), pp. 15--30. Dept. Comput. Sci., Univ. Manitoba, Winnipeg, Man., 1971. MR0335426 (49 #208)
  • J. Levine and R. E. Dalton, Minimum periods, modulo p, of first-order Bell exponential integers, Math. Comp., 16 (1962), 416-423.
  • Levinson, H.; Silverman, R. Topologies on finite sets. II. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 699--712, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561090 (81c:54006)
  • S. Linusson, The number of M-sequences and f-vectors, Combinatorica, 19 (1999), 255-266.
  • L. Lovasz, Combinatorial Problems and Exercises, North-Holland, 1993, pp. 14-15.
  • M. Meier, On the number of partitions of a given set, Amer. Math. Monthly, 114 (2007), p. 450.
  • Merris, Russell, and Stephen Pierce. "The Bell numbers and r-fold transitivity." Journal of Combinatorial Theory, Series A 12.1 (1972): 155-157.
  • Moser, Leo, and Max Wyman. An asymptotic formula for the Bell numbers. Trans. Royal Soc. Canada, 49 (1955), 49-53.
  • A. Murthy, Generalization of partition function, introducing Smarandache factor partition, Smarandache Notions Journal, Vol. 11, No. 1-2-3, Spring 2000.
  • Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 1.4,1.8.
  • P. Peart, Hankel determinants via Stieltjes matrices. Proceedings of the Thirty-first Southeastern International Conference on Combinatorics, Graph Theory and Computing (Boca Raton, FL, 2000). Congr. Numer. 144 (2000), 153-159.
  • A. M. Robert, A Course in p-adic Analysis, Springer-Verlag, 2000; p. 212.
  • G.-C. Rota, Finite Operator Calculus.
  • Frank Ruskey, Jennifer Woodcock and Yuji Yamauchi, Counting and computing the Rand and block distances of pairs of set partitions, Journal of Discrete Algorithms, Volume 16, October 2012, Pages 236-248.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge; see Section 1.4 and Example 5.2.4.
  • Abdullahi Umar, On the semigroups of order-decreasing finite full transformations, Proc. Roy. Soc. Edinburgh 120A (1992), 129-142.
  • Abdullahi Umar, On the semigroups of partial one-to-one order-decreasing finite transformations, Proc. Roy. Soc. Edinburgh 123A (1993), 355-363.

Crossrefs

Equals row sums of triangle A008277 (Stirling subset numbers).
Partial sums give A005001. a(n) = A123158(n, 0).
See A061462 for powers of 2 dividing a(n).
Rightmost diagonal of triangle A121207. A144293 gives largest prime factor.
Equals row sums of triangle A152432.
Row sums, right and left borders of A212431.
A diagonal of A011971. - N. J. A. Sloane, Jul 31 2012
Diagonal of A102661. - Manfred Boergens, Mar 11 2025
Cf. A054767 (period of this sequence mod n).
Row sums are A048993. - Wolfdieter Lang, Oct 16 2014
Sequences in the Erné (1974) paper: A000110, A000798, A001035, A001927, A001929, A006056, A006057, A006058, A006059.
Bell polynomials B(n,x): A001861 (x=2), A027710 (x=3), A078944 (x=4), A144180 (x=5), A144223 (x=6), A144263 (x=7), A221159 (x=8).
Cf. A243991 (sum of reciprocals), A085686 (inv. Euler Transf.).

Programs

  • Haskell
    type N = Integer
    n_partitioned_k :: N -> N -> N
    1 `n_partitioned_k` 1 = 1
    1 `n_partitioned_k` _ = 0
    n `n_partitioned_k` k = k * (pred n `n_partitioned_k` k) + (pred n `n_partitioned_k` pred k)
    n_partitioned :: N -> N
    n_partitioned 0 = 1
    n_partitioned n = sum $ map (\k -> n `n_partitioned_k` k) $ [1 .. n]
    -- Felix Denis, Oct 16 2012
    
  • Haskell
    a000110 = sum . a048993_row -- Reinhard Zumkeller, Jun 30 2013
    
  • Julia
    function a(n)
        t = [zeros(BigInt, n+1) for _ in 1:n+1]
        t[1][1] = 1
        for i in 2:n+1
            t[i][1] = t[i-1][i-1]
            for j in 2:i
                t[i][j] = t[i-1][j-1] + t[i][j-1]
            end
        end
        return [t[i][1] for i in 1:n+1]
    end
    print(a(28)) # Paul Muljadi, May 07 2024
    
  • Magma
    [Bell(n): n in [0..40]]; // Vincenzo Librandi, Feb 07 2011
    
  • Maple
    A000110 := proc(n) option remember; if n <= 1 then 1 else add( binomial(n-1,i)*A000110(n-1-i),i=0..n-1); fi; end: # version 1
    A := series(exp(exp(x)-1),x,60): A000110 := n->n!*coeff(A,x,n): # version 2
    A000110:= n-> add(Stirling2(n, k), k=0..n): seq(A000110(n), n=0..22); # version 3, from Zerinvary Lajos, Jun 28 2007
    A000110 := n -> combinat[bell](n): # version 4, from Peter Luschny, Mar 30 2011
    spec:= [S, {S=Set(U, card >= 1), U=Set(Z, card >= 1)}, labeled]: G:={P=Set(Set(Atom, card>0))}: combstruct[gfsolve](G, unlabeled, x): seq(combstruct[count]([P, G, labeled], size=i), i=0..22);  # version 5, Zerinvary Lajos, Dec 16 2007
    BellList := proc(m) local A, P, n; A := [1, 1]; P := [1]; for n from 1 to m - 2 do
    P := ListTools:-PartialSums([A[-1], op(P)]); A := [op(A), P[-1]] od; A end: BellList(29); # Peter Luschny, Mar 24 2022
  • Mathematica
    f[n_] := Sum[ StirlingS2[n, k], {k, 0, n}]; Table[ f[n], {n, 0, 40}] (* Robert G. Wilson v *)
    Table[BellB[n], {n, 0, 40}] (* Harvey P. Dale, Mar 01 2011 *)
    B[0] = 1; B[n_] := 1/E Sum[k^(n - 1)/(k-1)!, {k, 1, Infinity}] (* Dimitri Papadopoulos, Mar 10 2015, edited by M. F. Hasler, Nov 30 2018 *)
    BellB[Range[0,40]] (* Eric W. Weisstein, Aug 10 2017 *)
    b[1] = 1; k = 1; Flatten[{1, Table[Do[j = k; k += b[m]; b[m] = j;, {m, 1, n-1}]; b[n] = k, {n, 1, 40}]}] (* Vaclav Kotesovec, Sep 07 2019 *)
    Table[j! Coefficient[Series[Exp[Exp[x] - 1], {x, 0, 20}], x, j], {j, 0, 20}] (* Nikolaos Pantelidis, Feb 01 2023 *)
    Table[(D[Exp[Exp[x]], {x, n}] /. x -> 0)/E, {n, 0, 20}] (* Joan Ludevid, Nov 05 2024 *)
  • Maxima
    makelist(belln(n),n,0,40); /* Emanuele Munarini, Jul 04 2011 */
    
  • PARI
    {a(n) = my(m); if( n<0, 0, m = contfracpnqn( matrix(2, n\2, i, k, if( i==1, -k*x^2, 1 - (k+1)*x))); polcoeff(1 / (1 - x + m[2,1] / m[1,1]) + x * O(x^n), n))}; /* Michael Somos */
    
  • PARI
    {a(n) = polcoeff( sum( k=0, n, prod( i=1, k, x / (1 - i*x)), x^n * O(x)), n)}; /* Michael Somos, Aug 22 2004 */
    
  • PARI
    a(n)=round(exp(-1)*suminf(k=0,1.0*k^n/k!)) \\ Gottfried Helms, Mar 30 2007 - WARNING! For illustration only: Gives silently a wrong result for n = 42 and an error for n > 42, with standard precision of 38 digits. - M. F. Hasler, Nov 30 2018
    
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( exp( exp( x + x * O(x^n)) - 1), n))}; /* Michael Somos, Jun 28 2009 */
    
  • PARI
    Vec(serlaplace(exp(exp('x+O('x^66))-1))) \\ Joerg Arndt, May 26 2012
    
  • PARI
    A000110(n)=sum(k=0,n,stirling(n,k,2)) \\ M. F. Hasler, Nov 30 2018
    
  • Perl
    use bignum;sub a{my($n)=@;my@t=map{[(0)x($n+1)]}0..$n;$t[0][0]=1;for my$i(1..$n){$t[$i][0]=$t[$i-1][$i-1];for my$j(1..$i){$t[$i][$j]=$t[$i-1][$j-1]+$t[$i][$j-1]}}return map{$t[$][0]}0..$n-1}print join(", ",a(28)),"\n" # Paul Muljadi, May 08 2024
  • Python
    # The objective of this implementation is efficiency.
    # m -> [a(0), a(1), ..., a(m)] for m > 0.
    def A000110_list(m):
        A = [0 for i in range(m)]
        A[0] = 1
        R = [1, 1]
        for n in range(1, m):
            A[n] = A[0]
            for k in range(n, 0, -1):
                A[k-1] += A[k]
            R.append(A[0])
        return R
    A000110_list(40) # Peter Luschny, Jan 18 2011
    
  • Python
    # requires python 3.2 or higher. Otherwise use def'n of accumulate in python docs.
    from itertools import accumulate
    A000110, blist, b = [1,1], [1], 1
    for _ in range(20):
        blist = list(accumulate([b]+blist))
        b = blist[-1]
        A000110.append(b) # Chai Wah Wu, Sep 02 2014, updated Chai Wah Wu, Sep 19 2014
    
  • Python
    from sympy import bell
    print([bell(n) for n in range(27)]) # Michael S. Branicky, Dec 15 2021
    
  • Python
    from functools import cache
    @cache
    def a(n, k=0): return int(n < 1) or k*a(n-1, k) + a(n-1, k+1)
    print([a(n) for n in range(27)])  # Peter Luschny, Jun 14 2022
    
  • Sage
    from sage.combinat.expnums import expnums2; expnums2(30, 1) # Zerinvary Lajos, Jun 26 2008
    
  • Sage
    [bell_number(n) for n in (0..40)] # G. C. Greubel, Jun 13 2019
    

Formula

E.g.f.: exp(exp(x) - 1).
Recurrence: a(n+1) = Sum_{k=0..n} a(k)*binomial(n, k).
a(n) = Sum_{k=0..n} Stirling2(n, k).
a(n) = Sum_{j=0..n-1} (1/(n-1)!)*A000166(j)*binomial(n-1, j)*(n-j)^(n-1). - André F. Labossière, Dec 01 2004
G.f.: (Sum_{k>=0} 1/((1-k*x)*k!))/exp(1) = hypergeom([-1/x], [(x-1)/x], 1)/exp(1) = ((1-2*x)+LaguerreL(1/x, (x-1)/x, 1)+x*LaguerreL(1/x, (2*x-1)/x, 1))*Pi/(x^2*sin(Pi*(2*x-1)/x)), where LaguerreL(mu, nu, z) = (gamma(mu+nu+1)/(gamma(mu+1)*gamma(nu+1)))* hypergeom([-mu], [nu+1], z) is the Laguerre function, the analytic extension of the Laguerre polynomials, for mu not equal to a nonnegative integer. This generating function has an infinite number of poles accumulating in the neighborhood of x=0. - Karol A. Penson, Mar 25 2002
a(n) = exp(-1)*Sum_{k >= 0} k^n/k! [Dobinski]. - Benoit Cloitre, May 19 2002
a(n) is asymptotic to n!*(2 Pi r^2 exp(r))^(-1/2) exp(exp(r)-1) / r^n, where r is the positive root of r exp(r) = n. See, e.g., the Odlyzko reference.
a(n) is asymptotic to b^n*exp(b-n-1/2)*sqrt(b/(b+n)) where b satisfies b*log(b) = n - 1/2 (see Graham, Knuth and Patashnik, Concrete Mathematics, 2nd ed., p. 493). - Benoit Cloitre, Oct 23 2002, corrected by Vaclav Kotesovec, Jan 06 2013
Lovasz (Combinatorial Problems and Exercises, North-Holland, 1993, Section 1.14, Problem 9) gives another asymptotic formula, quoted by Mezo and Baricz. - N. J. A. Sloane, Mar 26 2015
G.f.: Sum_{k>=0} x^k/(Product_{j=1..k} (1-j*x)) (see Klazar for a proof). - Ralf Stephan, Apr 18 2004
a(n+1) = exp(-1)*Sum_{k>=0} (k+1)^(n)/k!. - Gerald McGarvey, Jun 03 2004
For n>0, a(n) = Aitken(n-1, n-1) [i.e., a(n-1, n-1) of Aitken's array (A011971)]. - Gerald McGarvey, Jun 26 2004
a(n) = Sum_{k=1..n} (1/k!)*(Sum_{i=1..k} (-1)^(k-i)*binomial(k, i)*i^n + 0^n). - Paul Barry, Apr 18 2005
a(n) = A032347(n) + A040027(n+1). - Jon Perry, Apr 26 2005
a(n) = (2*n!/(Pi*e))*Im( Integral_{x=0..Pi} e^(e^(e^(ix))) sin(nx) dx ) where Im denotes imaginary part [Cesaro]. - David Callan, Sep 03 2005
O.g.f.: 1/(1-x-x^2/(1-2*x-2*x^2/(1-3*x-3*x^2/(.../(1-n*x-n*x^2/(...)))))) (continued fraction due to Ph. Flajolet). - Paul D. Hanna, Jan 17 2006
From Karol A. Penson, Jan 14 2007: (Start)
Representation of Bell numbers B(n), n=1,2,..., as special values of hypergeometric function of type (n-1)F(n-1), in Maple notation: B(n)=exp(-1)*hypergeom([2,2,...,2],[1,1,...,1],1), n=1,2,..., i.e., having n-1 parameters all equal to 2 in the numerator, having n-1 parameters all equal to 1 in the denominator and the value of the argument equal to 1.
Examples:
B(1)=exp(-1)*hypergeom([],[],1)=1
B(2)=exp(-1)*hypergeom([2],[1],1)=2
B(3)=exp(-1)*hypergeom([2,2],[1,1],1)=5
B(4)=exp(-1)*hypergeom([2,2,2],[1,1,1],1)=15
B(5)=exp(-1)*hypergeom([2,2,2,2],[1,1,1,1],1)=52
(Warning: this formula is correct but its application by a computer may not yield exact results, especially with a large number of parameters.)
(End)
a(n+1) = 1 + Sum_{k=0..n-1} Sum_{i=0..k} binomial(k,i)*(2^(k-i))*a(i). - Yalcin Aktar, Feb 27 2007
a(n) = [1,0,0,...,0] T^(n-1) [1,1,1,...,1]', where T is the n X n matrix with main diagonal {1,2,3,...,n}, 1's on the diagonal immediately above and 0's elsewhere. [Meier]
a(n) = ((2*n!)/(Pi * e)) * ImaginaryPart(Integral[from 0 to Pi](e^e^e^(i*theta))*sin(n*theta) dtheta). - Jonathan Vos Post, Aug 27 2007
From Tom Copeland, Oct 10 2007: (Start)
a(n) = T(n,1) = Sum_{j=0..n} S2(n,j) = Sum_{j=0..n} E(n,j) * Lag(n,-1,j-n) = Sum_{j=0..n} [ E(n,j)/n! ] * [ n!*Lag(n,-1, j-n) ] where T(n,x) are the Bell / Touchard / exponential polynomials; S2(n,j), the Stirling numbers of the second kind; E(n,j), the Eulerian numbers; and Lag(n,x,m), the associated Laguerre polynomials of order m. Note that E(n,j)/n! = E(n,j) / (Sum_{k=0..n} E(n,k)).
The Eulerian numbers count the permutation ascents and the expression [n!*Lag(n,-1, j-n)] is A086885 with a simple combinatorial interpretation in terms of seating arrangements, giving a combinatorial interpretation to n!*a(n) = Sum_{j=0..n} E(n,j) * [n!*Lag(n,-1, j-n)].
(End)
Define f_1(x), f_2(x), ... such that f_1(x)=e^x and for n=2,3,... f_{n+1}(x) = (d/dx)(x*f_n(x)). Then for Bell numbers B_n we have B_n=1/e*f_n(1). - Milan Janjic, May 30 2008
a(n) = (n-1)! Sum_{k=1..n} a(n-k)/((n-k)! (k-1)!) where a(0)=1. - Thomas Wieder, Sep 09 2008
a(n+k) = Sum_{m=0..n} Stirling2(n,m) Sum_{r=0..k} binomial(k,r) m^r a(k-r). - David Pasino (davepasino(AT)yahoo.com), Jan 25 2009. (Umbrally, this may be written as a(n+k) = Sum_{m=0..n} Stirling2(n,m) (a+m)^k. - N. J. A. Sloane, Feb 07 2009)
Sum_{k=1..n-1} a(n)*binomial(n,k) = Sum_{j=1..n}(j+1)*Stirling2(n,j+1). - [Zhao] - R. J. Mathar, Jun 24 2024
From Thomas Wieder, Feb 25 2009: (Start)
a(n) = Sum_{k_1=0..n+1} Sum_{k_2=0..n}...Sum_{k_i=0..n-i}...Sum_{k_n=0..1}
delta(k_1,k_2,...,k_i,...,k_n)
where delta(k_1,k_2,...,k_i,...,k_n) = 0 if any k_i > k_(i+1) and k_(i+1) <> 0
and delta(k_1,k_2,...,k_i,...,k_n) = 1 otherwise.
(End)
Let A be the upper Hessenberg matrix of order n defined by: A[i,i-1]=-1, A[i,j]:=binomial(j-1,i-1), (i<=j), and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det(A). - Milan Janjic, Jul 08 2010
G.f. satisfies A(x) = (x/(1-x))*A(x/(1-x)) + 1. - Vladimir Kruchinin, Nov 28 2011
G.f.: 1 / (1 - x / (1 - 1*x / (1 - x / (1 - 2*x / (1 - x / (1 - 3*x / ... )))))). - Michael Somos, May 12 2012
a(n+1) = Sum_{m=0..n} Stirling2(n, m)*(m+1), n >= 0. Compare with the third formula for a(n) above. Here Stirling2 = A048993. - Wolfdieter Lang, Feb 03 2015
G.f.: (-1)^(1/x)*((-1/x)!/e + (!(-1-1/x))/x) where z! and !z are factorial and subfactorial generalized to complex arguments. - Vladimir Reshetnikov, Apr 24 2013
The following formulas were proposed during the period Dec 2011 - Oct 2013 by Sergei N. Gladkovskii: (Start)
E.g.f.: exp(exp(x)-1) = 1 + x/(G(0)-x); G(k) = (k+1)*Bell(k) + x*Bell(k+1) - x*(k+1)*Bell(k)*Bell(k+2)/G(k+1) (continued fraction).
G.f.: W(x) = (1-1/(G(0)+1))/exp(1); G(k) = x*k^2 + (3*x-1)*k - 2 + x - (k+1)*(x*k+x-1)^2/G(k+1); (continued fraction Euler's kind, 1-step).
G.f.: W(x) = (1 + G(0)/(x^2-3*x+2))/exp(1); G(k) = 1 - (x*k+x-1)/( ((k+1)!) - (((k+1)!)^2)*(1-x-k*x+(k+1)!)/( ((k+1)!)*(1-x-k*x+(k+1)!) - (x*k+2*x-1)*(1-2*x-k*x+(k+2)!)/G(k+1))); (continued fraction).
G.f.: A(x) = 1/(1 - x/(1-x/(1 + x/G(0)))); G(k) = x - 1 + x*k + x*(x-1+x*k)/G(k+1); (continued fraction, 1-step).
G.f.: -1/U(0) where U(k) = x*k - 1 + x - x^2*(k+1)/U(k+1); (continued fraction, 1-step).
G.f.: 1 + x/U(0) where U(k) = 1 - x*(k+2) - x^2*(k+1)/U(k+1); (continued fraction, 1-step).
G.f.: 1 + 1/(U(0) - x) where U(k) = 1 + x - x*(k+1)/(1 - x/U(k+1)); (continued fraction, 2-step).
G.f.: 1 + x/(U(0)-x) where U(k) = 1 - x*(k+1)/(1 - x/U(k+1)); (continued fraction, 2-step).
G.f.: 1/G(0) where G(k) = 1 - x/(1 - x*(2*k+1)/(1 - x/(1 - x*(2*k+2)/G(k+1) ))); (continued fraction).
G.f.: G(0)/(1+x) where G(k) = 1 - 2*x*(k+1)/((2*k+1)*(2*x*k-1) - x*(2*k+1)*(2*k+3)*(2*x*k-1)/(x*(2*k+3) - 2*(k+1)*(2*x*k+x-1)/G(k+1) )); (continued fraction).
G.f.: -(1+2*x) * Sum_{k >= 0} x^(2*k)*(4*x*k^2-2*k-2*x-1) / ((2*k+1) * (2*x*k-1)) * A(k) / B(k) where A(k) = Product_{p=0..k} (2*p+1), B(k) = Product_{p=0..k} (2*p-1) * (2*x*p-x-1) * (2*x*p-2*x-1).
G.f.: (G(0) - 1)/(x-1) where G(k) = 1 - 1/(1-k*x)/(1-x/(x-1/G(k+1) )); (continued fraction).
G.f.: 1 + x*(S-1) where S = Sum_{k>=0} ( 1 + (1-x)/(1-x-x*k) )*(x/(1-x))^k/Product_{i=0..k-1} (1-x-x*i)/(1-x).
G.f.: (G(0) - 2)/(2*x-1) where G(k) = 2 - 1/(1-k*x)/(1-x/(x-1/G(k+1) )); (continued fraction).
G.f.: -G(0) where G(k) = 1 - (x*k - 2)/(x*k - 1 - x*(x*k - 1)/(x + (x*k - 2)/G(k+1) )); (continued fraction).
G.f.: G(0) where G(k) = 2 - (2*x*k - 1)/(x*k - 1 - x*(x*k - 1)/(x + (2*x*k - 1)/G(k+1) )); (continued fraction).
G.f.: (G(0) - 1)/(1+x) where G(k) = 1 + 1/(1-k*x)/(1-x/(x+1/G(k+1) )); (continued fraction).
G.f.: 1/(x*(1-x)*G(0)) - 1/x where G(k) = 1 - x/(x - 1/(1 + 1/(x*k-1)/G(k+1) )); (continued fraction).
G.f.: 1 + x/( Q(0) - x ) where Q(k) = 1 + x/( x*k - 1 )/Q(k+1); (continued fraction).
G.f.: 1+x/Q(0), where Q(k) = 1 - x - x/(1 - x*(k+1)/Q(k+1)); (continued fraction).
G.f.: 1/(1-x*Q(0)), where Q(k) = 1 + x/(1 - x + x*(k+1)/(x - 1/Q(k+1))); (continued fraction).
G.f.: Q(0)/(1-x), where Q(k) = 1 - x^2*(k+1)/( x^2*(k+1) - (1-x*(k+1))*(1-x*(k+2))/Q(k+1) ); (continued fraction).
(End)
a(n) ~ exp(exp(W(n))-n-1)*n^n/W(n)^(n+1/2), where W(x) is the Lambert W-function. - Vladimir Reshetnikov, Nov 01 2015
a(n) ~ n^n * exp(n/W(n)-1-n) / (sqrt(1+W(n)) * W(n)^n). - Vaclav Kotesovec, Nov 13 2015
From Natalia L. Skirrow, Apr 13 2025: (Start)
By taking logarithmic derivatives of the equivalent to Kotesovec's asymptotic for Bell polynomials at x=1, we obtain properties of the nth row of A008277 as a statistical distribution (where W=W(n),X=W(n)+1)
a(n+1)/a(n) ~ n/W + W/(2*(W+1)^2) is 1 more than the expectation.
(2*a(n+1)+a(n+2))/a(n) - (a(n+1)/a(n))^2 - a(n+2)/a(n+1) ~ n/(W*X)+1/(2*X^2)-3/(2*X^3)+1/X^4 is 1 more than the variance.
(This is a complete asymptotic characterisation, since they converge to normal distributions; see Harper, 1967)
(End)
a(n) are the coefficients in the asymptotic expansion of -exp(-1)*(-1)^x*x*Gamma(-x,0,-1), where Gamma(a,z0,z1) is the generalized incomplete Gamma function. - Vladimir Reshetnikov, Nov 12 2015
a(n) = 1 + floor(exp(-1) * Sum_{k=1..2*n} k^n/k!). - Vladimir Reshetnikov, Nov 13 2015
a(p^m) == m+1 (mod p) when p is prime and m >= 1 (see Lemma 3.1 in the Hurst/Schultz reference). - Seiichi Manyama, Jun 01 2016
a(n) = Sum_{k=0..n} hypergeom([1, -k], [], 1)*Stirling2(n+1, k+1) = Sum_{k=0..n} A182386(k)*Stirling2(n+1, k+1). - Mélika Tebni, Jul 02 2022
For n >= 1, a(n) = Sum_{i=0..n-1} a(i)*A074664(n-i). - Davide Rotondo, Apr 21 2024
a(n) is the n-th derivative of e^e^x divided by e at point x=0. - Joan Ludevid, Nov 05 2024

Extensions

Edited by M. F. Hasler, Nov 30 2018

A048993 Triangle of Stirling numbers of 2nd kind, S(n,k), n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 7, 6, 1, 0, 1, 15, 25, 10, 1, 0, 1, 31, 90, 65, 15, 1, 0, 1, 63, 301, 350, 140, 21, 1, 0, 1, 127, 966, 1701, 1050, 266, 28, 1, 0, 1, 255, 3025, 7770, 6951, 2646, 462, 36, 1, 0, 1, 511, 9330, 34105, 42525, 22827, 5880, 750, 45, 1
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Also known as Stirling set numbers.
S(n,k) enumerates partitions of an n-set into k nonempty subsets.
The o.g.f. for the sequence of diagonal k (k=0 for the main diagonal) is G(k,x) = ((x^k)/(1-x)^(2*k+1))*Sum_{m=0..k-1} A008517(k,m+1)*x^m. A008517 is the second-order Eulerian triangle. - Wolfdieter Lang, Oct 14 2005
From Philippe Deléham, Nov 14 2007: (Start)
Sum_{k=0..n} S(n,k)*x^k = B_n(x), where B_n(x) = Bell polynomials.
The first few Bell polynomials are:
B_0(x) = 1;
B_1(x) = 0 + x;
B_2(x) = 0 + x + x^2;
B_3(x) = 0 + x + 3x^2 + x^3;
B_4(x) = 0 + x + 7x^2 + 6x^3 + x^4;
B_5(x) = 0 + x + 15x^2 + 25x^3 + 10x^4 + x^5;
B_6(x) = 0 + x + 31x^2 + 90x^3 + 65x^4 + 15x^5 + x^6;
(End)
This is the Sheffer triangle (1, exp(x) - 1), an exponential (binomial) convolution triangle. The a-sequence is given by A006232/A006233 (Cauchy sequence). The z-sequence is the zero sequence. See the link under A006232 for the definition and use of these sequences. The row sums give A000110 (Bell), and the alternating row sums give A000587 (see the Philippe Deléham formulas and crossreferences below). - Wolfdieter Lang, Oct 16 2014
Also the inverse Bell transform of the factorial numbers (A000142). For the definition of the Bell transform see A264428 and for cross-references A265604. - Peter Luschny, Dec 31 2015
From Wolfdieter Lang, Feb 21 2017: (Start)
The transposed (trans) of this lower triagonal Sheffer matrix of the associated type S = (1, exp(x) - 1) (taken as N X N matrix for arbitrarily large N) provides the transition matrix from the basis {x^n/n!}, n >= 0, to the basis {y^n/n!}, n >= 0, with y^n/n! = Sum_{m>=n} S^{trans}(n, m) x^m/m! = Sum_{m>=0} x^m/m!*S(m, n).
The Sheffer transform with S = (g, f) of a sequence {a_n} to {b_n} for n >= 0, in matrix notation vec(b) = S vec(a), satisfies, with e.g.f.s A and B, B(x) = g(x)*A(f(x)) and B(x) = A(y(x)) identically, with vec(xhat) = S^{trans,-1} vec(yhat) in symbolic notation with vec(xhat)_n = x^n/n! (similarly for vec(yhat)).
(End)
Number of partitions of {1, 2, ..., n+1} into k+1 nonempty subsets such that no subset contains two adjacent numbers. - Thomas Anton, Sep 26 2022

Examples

			The triangle S(n,k) begins:
  n\k 0 1    2     3      4       5       6      7      8     9   10 11 12
  0:  1
  1:  0 1
  2:  0 1    1
  3:  0 1    3     1
  4:  0 1    7     6      1
  5:  0 1   15    25     10       1
  6:  0 1   31    90     65      15       1
  7:  0 1   63   301    350     140      21      1
  8:  0 1  127   966   1701    1050     266     28      1
  9:  0 1  255  3025   7770    6951    2646    462     36     1
 10:  0 1  511  9330  34105   42525   22827   5880    750    45    1
 11:  0 1 1023 28501 145750  246730  179487  63987  11880  1155   55  1
 12:  0 1 2047 86526 611501 1379400 1323652 627396 159027 22275 1705 66  1
 ... reformatted and extended - _Wolfdieter Lang_, Oct 16 2014
Completely symmetric function S(4, 2) = h^{(2)}_2 = 1^2 + 2^2 + 1^1*2^1 = 7; S(5, 2) = h^{(2)}_3 = 1^3 + 2^3 + 1^2*2^1 + 1^1*2^2 = 15. - _Wolfdieter Lang_, May 26 2017
From _Wolfdieter Lang_, Aug 11 2017: (Start)
Recurrence: S(5, 3) = S(4, 2) + 2*S(4, 3) = 7 + 3*6 = 25.
Boas-Buck recurrence for column m = 3, and n = 5: S(5, 3) = (3/2)*((5/2)*S(4, 3) + 10*Bernoulli(2)*S(3, 3)) = (3/2)*(15 + 10*(1/6)*1) = 25. (End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 835.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 310.
  • J. H. Conway and R. K. Guy, The Book of Numbers, Springer, p. 92.
  • F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 223.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 244.
  • J. Riordan, An Introduction to Combinatorial Analysis, p. 48.

Crossrefs

See especially A008277 which is the main entry for this triangle.
A000110(n) = sum(S(n, k)) k=0..n, n >= 0. Cf. A085693.
Cf. A084938, A106800 (mirror image), A138378, A213061 (mod 2).

Programs

  • Haskell
    a048993 n k = a048993_tabl !! n !! k
    a048993_row n = a048993_tabl !! n
    a048993_tabl = iterate (\row ->
       [0] ++ (zipWith (+) row $ zipWith (*) [1..] $ tail row) ++ [1]) [1]
    -- Reinhard Zumkeller, Mar 26 2012
  • Maple
    for n from 0 to 10 do seq(Stirling2(n,k),k=0..n) od; # yields sequence in triangular form # Emeric Deutsch, Nov 01 2006
  • Mathematica
    t[n_, k_] := StirlingS2[n, k]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Robert G. Wilson v *)
  • Maxima
    create_list(stirling2(n,k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • PARI
    for(n=0, 22, for(k=0, n, print1(stirling(n, k, 2), ", ")); print()); \\ Joerg Arndt, Apr 21 2013
    

Formula

S(n, k) = k*S(n-1, k) + S(n-1, k-1), n > 0; S(0, k) = 0, k > 0; S(0, 0) = 1.
Equals [0, 1, 0, 2, 0, 3, 0, 4, 0, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] where DELTA is Deléham's operator defined in A084938.
Sum_{k = 0..n} x^k*S(n, k) = A213170(n), A000587(n), A000007(n), A000110(n), A001861(n), A027710(n), A078944(n), A144180(n), A144223(n), A144263(n) respectively for x = -2, -1, 0, 1, 2, 3, 4, 5, 6, 7. - Philippe Deléham, May 09 2004, Feb 16 2013
S(n, k) = Sum_{i=0..k} (-1)^(k+i)binomial(k, i)i^n/k!. - Paul Barry, Aug 05 2004
Sum_{k=0..n} k*S(n,k) = B(n+1)-B(n), where B(q) are the Bell numbers (A000110). - Emeric Deutsch, Nov 01 2006
Equals the inverse binomial transform of A008277. - Gary W. Adamson, Jan 29 2008
G.f.: 1/(1-xy/(1-x/(1-xy/(1-2x/(1-xy/1-3x/(1-xy/(1-4x/(1-xy/(1-5x/(1-... (continued fraction equivalent to Deléham DELTA construction). - Paul Barry, Dec 06 2009
G.f.: 1/Q(0), where Q(k) = 1 - (y+k)*x - (k+1)*y*x^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 09 2013
Inverse of padded A008275 (padded just as A048993 = padded A008277). - Tom Copeland, Apr 25 2014
E.g.f. for the row polynomials s(n,x) = Sum_{k=0..n} S(n,k)*x^k is exp(x*(exp(z)-1)) (Sheffer property). E.g.f. for the k-th column sequence with k leading zeros is ((exp(x)-1)^k)/k! (Sheffer property). - Wolfdieter Lang, Oct 16 2014
G.f. for column k: x^k/Product_{j=1..k} (1-j*x), k >= 0 (with the empty product for k = 0 put to 1). See Abramowitz-Stegun, p. 824, 24.1.4 B. - Wolfdieter Lang, May 26 2017
Boas-Buck recurrence for column sequence m: S(n, k) = (k/(n - k))*(n*S(n-1, k)/2 + Sum_{p=k..n-2} (-1)^(n-p)*binomial(n,p)*Bernoulli(n-p)*S(p, k)), for n > k >= 0, with input T(k,k) = 1. See a comment and references in A282629. An example is given below. - Wolfdieter Lang, Aug 11 2017
The n-th row polynomial has the form x o x o ... o x (n factors), where o denotes the white diamond multiplication operator defined in Bala - see Example E4. - Peter Bala, Jan 07 2018
Sum_{k=1..n} k*S(n,k) = A138378(n). - Alois P. Heinz, Jan 07 2022
S(n,k) = Sum_{j=k..n} (-1)^(j-k)*A059297(n,j)*A354794(j,k). - Mélika Tebni, Jan 27 2023

A000248 Expansion of e.g.f. exp(x*exp(x)).

Original entry on oeis.org

1, 1, 3, 10, 41, 196, 1057, 6322, 41393, 293608, 2237921, 18210094, 157329097, 1436630092, 13810863809, 139305550066, 1469959371233, 16184586405328, 185504221191745, 2208841954063318, 27272621155678841, 348586218389733556, 4605223387997411873
Offset: 0

Views

Author

Keywords

Comments

Number of forests with n nodes and height at most 1.
Equivalently, number of idempotent mappings f from a set of n elements into itself (i.e., satisfying f o f = f). - Robert FERREOL, Oct 11 2007
In other words, a(n) = number of idempotents in the full semigroup of maps from [1..n] to itself. [Tainiter]
a(n) is the number of ways to select a set partition of {1,2,...,n} and then designate one element in each block (cell) of the partition.
Let set B have cardinality n. Then a(n) is the number of functions f:D->C over all partitions {D,C} of B. See the example in the Example Section below. We note that f:empty set->B is designated as the null function, whereas f:B->empty set is undefined unless B itself is empty. - Dennis P. Walsh, Dec 05 2013
In physics, a(n) would be interpreted as the number of projection operators P on S_n, i.e., ones satisfying P^2 = P. Example: a particle with a half-integer spin s has a spin space with 2s+1 base states which admits a(2s+1) linear projection operators (including the identity). These are important because they satisfy the operator identity exp(zU) = 1+(exp(z)-1)*U, valid for any complex z. - Stanislav Sykora, Nov 03 2016

Examples

			a(3)=10 since, for B={1,2,3}, we have 10 functions: 1 function of the type f:empty set->B; 6 functions of the type f:{x}->B\{x}; and 3 functions of the type f:{x,y}->B\{x,y}. - _Dennis P. Walsh_, Dec 05 2013
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 91.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.32(d).

Crossrefs

First row of array A098697.
Row sums of A133399.
Column k=1 of A210725, A279636.
Column k=2 of A245501.

Programs

  • Magma
    m:=25; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Exp(x*Exp(x)))); [Factorial(n-1)*b[n]: n in [1..m]]; // Vincenzo Librandi, Feb 01 2020
  • Maple
    A000248 := proc(n) local k; add(k^(n-k)*binomial(n, k), k=0..n); end; # Robert FERREOL, Oct 11 2007
    a:= proc(n) option remember; if n=0 then 1 else add(binomial(n-1, j) *(j+1) *a(n-1-j), j=0..n-1) fi end: seq(a(n), n=0..20); # Zerinvary Lajos, Mar 28 2009
  • Mathematica
    CoefficientList[Series[Exp[x Exp[x]],{x,0,20}],x]*Table[n!,{n,0,20}]
    a[0] = 1; a[1] = 1; a[n_] := a[n] = a[n - 1] + Sum[(Binomial[n - 1, j] + (n - 1) Binomial[n - 2, j]) a[j], {j, 0, n - 2}]; Table[a[n], {n, 0, 20}] (* David Callan, Oct 04 2013 *)
    Flatten[{1,Table[Sum[Binomial[n,k]*(n-k)^k,{k,0,n}],{n,1,20}]}] (* Vaclav Kotesovec, Jul 13 2014 *)
    Table[Sum[BellY[n, k, Range[n]], {k, 0, n}], {n, 0, 20}] (* Vladimir Reshetnikov, Nov 09 2016 *)
  • PARI
    a(n)=sum(k=0,n,binomial(n,k)*(n-k)^k); \\ Paul D. Hanna, Jun 26 2009
    
  • PARI
    x='x+O('x^66); Vec(serlaplace(exp(x*exp(x)))) \\ Joerg Arndt, Oct 06 2013
    
  • Sage
    # uses[bell_matrix from A264428]
    B = bell_matrix(lambda k: k+1, 20)
    print([sum(B.row(n)) for n in range(20)]) # Peter Luschny, Sep 03 2019
    

Formula

G.f.: Sum_{k>=0} x^k/(1-k*x)^(k+1). - Vladeta Jovovic, Oct 25 2003
a(n) = Sum_{k=0..n} C(n,k)*(n-k)^k. - Paul D. Hanna, Jun 26 2009
G.f.: G(0) where G(k) = 1 - x*(-1+2*k*x)^(2*k+1)/((x-1+2*k*x)^(2*k+2) - x*(x-1+2*k*x)^(4*k+4)/(x*(x-1+2*k*x)^(2*k+2) - (2*x-1+2*k*x)^(2*k+3)/G(k+1))) (continued fraction). - Sergei N. Gladkovskii, Jan 26 2013
E.g.f.: 1 + x/(1+x)*(G(0) - 1) where G(k) = 1 + exp(x)/(k+1)/(1-x/(x+(1)/G(k+1))) (continued fraction). - Sergei N. Gladkovskii, Feb 04 2013
Recurrence: a(0)=1, a(n) = Sum_{k=1..n} binomial(n-1,k-1)*k*a(n-k). - James East, Mar 30 2014
Asymptotics (Harris and Schoenfeld, 1968): a(n) ~ sqrt((r+1)/(2*Pi*(n+1)*(r^2+3*r+1))) * n! * exp((n+1)/(r+1)) / r^n, where r is the root of the equation r*(r+1)*exp(r) = n+1. - Vaclav Kotesovec, Jul 13 2014
a(n) = Sum_{k=0..n} A005727(k)*Stirling2(n, k). - Mélika Tebni, Jun 12 2022
More precise asymptotics: a(n) ~ n^(n + 1/2) / (sqrt(1 + 3*r + r^2) * exp(n - r*exp(r) + r/2) * r^(n + 1/2)), where r = 2*w - 1/(2*w) + 5/(8*w^2) - 19/(24*w^3) + 209/(192*w^4) - 763/(480*w^5) + 4657/(1920*w^6) - 6855/(1792*w^7) + 199613/(32256*w^8) + ... and w = LambertW(sqrt(n)/2). - Vaclav Kotesovec, Feb 20 2023

Extensions

In view of the multiple appearances of this sequence, I replaced the definition with the simple exponential generating function. - N. J. A. Sloane, Apr 16 2018

A038207 Triangle whose (i,j)-th entry is binomial(i,j)*2^(i-j).

Original entry on oeis.org

1, 2, 1, 4, 4, 1, 8, 12, 6, 1, 16, 32, 24, 8, 1, 32, 80, 80, 40, 10, 1, 64, 192, 240, 160, 60, 12, 1, 128, 448, 672, 560, 280, 84, 14, 1, 256, 1024, 1792, 1792, 1120, 448, 112, 16, 1, 512, 2304, 4608, 5376, 4032, 2016, 672, 144, 18, 1, 1024, 5120, 11520, 15360, 13440, 8064, 3360, 960, 180, 20, 1
Offset: 0

Views

Author

Keywords

Comments

This infinite matrix is the square of the Pascal matrix (A007318) whose rows are [ 1,0,... ], [ 1,1,0,... ], [ 1,2,1,0,... ], ...
As an upper right triangle, table rows give number of points, edges, faces, cubes,
4D hypercubes etc. in hypercubes of increasing dimension by column. - Henry Bottomley, Apr 14 2000. More precisely, the (i,j)-th entry is the number of j-dimensional subspaces of an i-dimensional hypercube (see the Coxeter reference). - Christof Weber, May 08 2009
Number of different partial sums of 1+[1,1,2]+[2,2,3]+[3,3,4]+[4,4,5]+... with entries that are zero removed. - Jon Perry, Jan 01 2004
Row sums are powers of 3 (A000244), antidiagonal sums are Pell numbers (A000129). - Gerald McGarvey, May 17 2005
Riordan array (1/(1-2x), x/(1-2x)). - Paul Barry, Jul 28 2005
T(n,k) is the number of elements of the Coxeter group B_n with descent set contained in {s_k}, 0<=k<=n-1. For T(n,n), we interpret this as the number of elements of B_n with empty descent set (since s_n does not exist). - Elizabeth Morris (epmorris(AT)math.washington.edu), Mar 01 2006
Let S be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xSy if x is a subset of y. Then T(n,k) = the number of elements (x,y) of S for which y has exactly k more elements than x. - Ross La Haye, Oct 12 2007
T(n,k) is number of paths in the first quadrant going from (0,0) to (n,k) using only steps B=(1,0) colored blue, R=(1,0) colored red and U=(1,1). Example: T(3,2)=6 because we have BUU, RUU, UBU, URU, UUB and UUR. - Emeric Deutsch, Nov 04 2007
T(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (0,1), and two kinds of step (1,0). - Joerg Arndt, Jul 01 2011
T(i,j) is the number of i-permutations of {1,2,3} containing j 1's. Example: T(2,1)=4 because we have 12, 13, 21 and 31; T(3,2)=6 because we have 112, 113, 121, 131, 211 and 311. - Zerinvary Lajos, Dec 21 2007
Triangle of coefficients in expansion of (2+x)^n. - N-E. Fahssi, Apr 13 2008
Sum of diagonals are Jacobsthal-numbers: A001045. - Mark Dols, Aug 31 2009
Triangle T(n,k), read by rows, given by [2,0,0,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 15 2009
Eigensequence of the triangle = A004211: (1, 3, 11, 49, 257, 1539, ...). - Gary W. Adamson, Feb 07 2010
f-vectors ("face"-vectors) for n-dimensional cubes [see e.g., Hoare]. (This is a restatement of Bottomley's above.) - Tom Copeland, Oct 19 2012
With P = Pascal matrix, the sequence of matrices I, A007318, A038207, A027465, A038231, A038243, A038255, A027466 ... = P^0, P^1, P^2, ... are related by Copeland's formula below to the evolution at integral time steps n= 0, 1, 2, ... of an exponential distribution exp(-x*z) governed by the Fokker-Planck equation as given in the Dattoli et al. ref. below. - Tom Copeland, Oct 26 2012
The matrix elements of the inverse are T^(-1)(n,k) = (-1)^(n+k)*T(n,k). - R. J. Mathar, Mar 12 2013
Unsigned diagonals of A133156 are rows of this array. - Tom Copeland, Oct 11 2014
Omitting the first row, this is the production matrix for A039683, where an equivalent differential operator can be found. - Tom Copeland, Oct 11 2016
T(n,k) is the number of functions f:[n]->[3] with exactly k elements mapped to 3. Note that there are C(n,k) ways to choose the k elements mapped to 3, and there are 2^(n-k) ways to map the other (n-k) elements to {1,2}. Hence, by summing T(n,k) as k runs from 0 to n, we obtain 3^n = Sum_{k=0..n} T(n,k). - Dennis P. Walsh, Sep 26 2017
Since this array is the square of the Pascal lower triangular matrix, the row polynomials of this array are obtained as the umbral composition of the row polynomials P_n(x) of the Pascal matrix with themselves. E.g., P_3(P.(x)) = 1 P_3(x) + 3 P_2(x) + 3 P_1(x) + 1 = (x^3 + 3 x^2 + 3 x + 1) + 3 (x^2 + 2 x + 1) + 3 (x + 1) + 1 = x^3 + 6 x^2 + 12 x + 8. - Tom Copeland, Nov 12 2018
T(n,k) is the number of 2-compositions of n+1 with some zeros allowed that have k zeros; see the Hopkins & Ouvry reference. - Brian Hopkins, Aug 16 2020
Also the convolution triangle of A000079. - Peter Luschny, Oct 09 2022

Examples

			Triangle begins with T(0,0):
   1;
   2,  1;
   4,  4,  1;
   8, 12,  6,  1;
  16, 32, 24,  8,  1;
  32, 80, 80, 40, 10,  1;
  ... -  corrected by _Clark Kimberling_, Aug 05 2011
Seen as an array read by descending antidiagonals:
[0] 1, 2,  4,   8,    16,    32,    64,     128,     256, ...     [A000079]
[1] 1, 4,  12,  32,   80,    192,   448,    1024,    2304, ...    [A001787]
[2] 1, 6,  24,  80,   240,   672,   1792,   4608,    11520, ...   [A001788]
[3] 1, 8,  40,  160,  560,   1792,  5376,   15360,   42240, ...   [A001789]
[4] 1, 10, 60,  280,  1120,  4032,  13440,  42240,   126720, ...  [A003472]
[5] 1, 12, 84,  448,  2016,  8064,  29568,  101376,  329472, ...  [A054849]
[6] 1, 14, 112, 672,  3360,  14784, 59136,  219648,  768768, ...  [A002409]
[7] 1, 16, 144, 960,  5280,  25344, 109824, 439296,  1647360, ... [A054851]
[8] 1, 18, 180, 1320, 7920,  41184, 192192, 823680,  3294720, ... [A140325]
[9] 1, 20, 220, 1760, 11440, 64064, 320320, 1464320, 6223360, ... [A140354]
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 155.
  • H. S. M. Coxeter, Regular Polytopes, Dover Publications, New York (1973), p. 122.

Crossrefs

Programs

  • GAP
    Flat(List([0..15], n->List([0..n], k->Binomial(n, k)*2^(n-k)))); # Stefano Spezia, Nov 21 2018
  • Haskell
    a038207 n = a038207_list !! n
    a038207_list = concat $ iterate ([2,1] *) [1]
    instance Num a => Num [a] where
       fromInteger k = [fromInteger k]
       (p:ps) + (q:qs) = p + q : ps + qs
       ps + qs         = ps ++ qs
       (p:ps) * qs'@(q:qs) = p * q : ps * qs' + [p] * qs
        *                = []
    -- Reinhard Zumkeller, Apr 02 2011
    
  • Haskell
    a038207' n k = a038207_tabl !! n !! k
    a038207_row n = a038207_tabl !! n
    a038207_tabl = iterate f [1] where
       f row = zipWith (+) ([0] ++ row) (map (* 2) row ++ [0])
    -- Reinhard Zumkeller, Feb 27 2013
    
  • Magma
    /* As triangle */ [[(&+[Binomial(n,i)*Binomial(i,k): i in [k..n]]): k in [0..n]]: n in [0..15]]; // Vincenzo Librandi, Nov 16 2018
    
  • Maple
    for i from 0 to 12 do seq(binomial(i, j)*2^(i-j), j = 0 .. i) end do; # yields sequence in triangular form - Emeric Deutsch, Nov 04 2007
    # Uses function PMatrix from A357368. Adds column 1, 0, 0, ... to the left.
    PMatrix(10, n -> 2^(n-1)); # Peter Luschny, Oct 09 2022
  • Mathematica
    Table[CoefficientList[Expand[(y + x + x^2)^n], y] /. x -> 1, {n, 0,10}] // TableForm (* Geoffrey Critzer, Nov 20 2011 *)
    Table[Binomial[n,k]2^(n-k),{n,0,10},{k,0,n}]//Flatten (* Harvey P. Dale, May 22 2020 *)
  • PARI
    {T(n, k) = polcoeff((x+2)^n, k)}; /* Michael Somos, Apr 27 2000 */
    
  • Sage
    def A038207_triangle(dim):
        M = matrix(ZZ,dim,dim)
        for n in range(dim): M[n,n] = 1
        for n in (1..dim-1):
            for k in (0..n-1):
                M[n,k] = M[n-1,k-1]+2*M[n-1,k]
        return M
    A038207_triangle(9)  # Peter Luschny, Sep 20 2012
    

Formula

T(n, k) = Sum_{i=0..n} binomial(n,i)*binomial(i,k).
T(n, k) = (-1)^k*A065109(n,k).
G.f.: 1/(1-2*z-t*z). - Emeric Deutsch, Nov 04 2007
Rows of the triangle are generated by taking successive iterates of (A135387)^n * [1, 0, 0, 0, ...]. - Gary W. Adamson, Dec 09 2007
From the formalism of A133314, the e.g.f. for the row polynomials of A038207 is exp(x*t)*exp(2x). The e.g.f. for the row polynomials of the inverse matrix is exp(x*t)*exp(-2x). p iterates of the matrix give the matrix with e.g.f. exp(x*t)*exp(p*2x). The results generalize for 2 replaced by any number. - Tom Copeland, Aug 18 2008
Sum_{k=0..n} T(n,k)*x^k = (2+x)^n. - Philippe Deléham, Dec 15 2009
n-th row is obtained by taking pairwise sums of triangle A112857 terms starting from the right. - Gary W. Adamson, Feb 06 2012
T(n,n) = 1 and T(n,k) = T(n-1,k-1) + 2*T(n-1,k) for kJon Perry, Oct 11 2012
The e.g.f. for the n-th row is given by umbral composition of the normalized Laguerre polynomials A021009 as p(n,x) = L(n, -L(.,-x))/n! = 2^n L(n, -x/2)/n!. E.g., L(2,x) = 2 -4*x +x^2, so p(2,x)= (1/2)*L(2, -L(.,-x)) = (1/2)*(2*L(0,-x) + 4*L(1,-x) + L(2,-x)) = (1/2)*(2 + 4*(1+x) + (2+4*x+x^2)) = 4 + 4*x + x^2/2. - Tom Copeland, Oct 20 2012
From Tom Copeland, Oct 26 2012: (Start)
From the formalism of A132440 and A218272:
Let P and P^T be the Pascal matrix and its transpose and H= P^2= A038207.
Then with D the derivative operator,
exp(x*z/(1-2*z))/(1-2*z)= exp(2*z D_z z) e^(x*z)= exp(2*D_x (x D_x)) e^(z*x)
= (1 z z^2 z^3 ...) H (1 x x^2/2! x^3/3! ...)^T
= (1 x x^2/2! x^3/3! ...) H^T (1 z z^2 z^3 ...)^T
= Sum_{n>=0} z^n * 2^n Lag_n(-x/2)= exp[z*EF(.,x)], an o.g.f. for the f-vectors (rows) of A038207 where EF(n,x) is an e.g.f. for the n-th f-vector. (Lag_n(x) are the un-normalized Laguerre polynomials.)
Conversely,
exp(z*(2+x))= exp(2D_x) exp(x*z)= exp(2x) exp(x*z)
= (1 x x^2 x^3 ...) H^T (1 z z^2/2! z^3/3! ...)^T
= (1 z z^2/2! z^3/3! ...) H (1 x x^2 x^3 ...)^T
= exp(z*OF(.,x)), an e.g.f for the f-vectors of A038207 where
OF(n,x)= (2+x)^n is an o.g.f. for the n-th f-vector.
(End)
G.f.: R(0)/2, where R(k) = 1 + 1/(1 - (2*k+1+ (1+y))*x/((2*k+2+ (1+y))*x + 1/R(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Nov 09 2013
A038207 = exp[M*B(.,2)] where M = A238385-I and (B(.,x))^n = B(n,x) are the Bell polynomials (cf. A008277). B(n,2) = A001861(n). - Tom Copeland, Apr 17 2014
T = (A007318)^2 = A112857*|A167374| = |A118801|*|A167374| = |A118801*A167374| = |P*A167374*P^(-1)*A167374| = |P*NpdP*A167374|. Cf. A118801. - Tom Copeland, Nov 17 2016
E.g.f. for the n-th subdiagonal, n = 0,1,2,..., equals exp(x)*P(n,x), where P(n,x) is the polynomial 2^n*Sum_{k = 0..n} binomial(n,k)*x^k/k!. For example, the e.g.f. for the third subdiagonal is exp(x)*(8 + 24*x + 12*x^2 + 4*x^3/3) = 8 + 32*x + 80*x^2/2! + 160*x^3/3! + .... - Peter Bala, Mar 05 2017
T(3*k+2,k) = T(3*k+2,k+1), T(2*k+1,k) = 2*T(2*k+1,k+1). - Yuchun Ji, May 26 2020
From Robert A. Russell, Aug 05 2020: (Start)
G.f. for column k: x^k / (1-2*x)^(k+1).
E.g.f. for column k: exp(2*x) * x^k / k!. (End)
Also the array A(n, k) read by descending antidiagonals, where A(n, k) = (-1)^n*Sum_{j= 0..n+k} binomial(n + k, j)*hypergeom([-n, j+1], [1], 1). - Peter Luschny, Nov 09 2021

A005493 2-Bell numbers: a(n) = number of partitions of [n+1] with a distinguished block.

Original entry on oeis.org

1, 3, 10, 37, 151, 674, 3263, 17007, 94828, 562595, 3535027, 23430840, 163254885, 1192059223, 9097183602, 72384727657, 599211936355, 5150665398898, 45891416030315, 423145657921379, 4031845922290572, 39645290116637023, 401806863439720943, 4192631462935194064
Offset: 0

Views

Author

Keywords

Comments

Number of Boolean sublattices of the Boolean lattice of subsets of {1..n}.
a(n) = p(n+1) where p(x) is the unique degree n polynomial such that p(k) = A000110(k+1) for k = 0, 1, ..., n. - Michael Somos, Oct 07 2003
With offset 1, number of permutations beginning with 12 and avoiding 21-3.
Rows sums of Bell's triangle (A011971). - Jorge Coveiro, Dec 26 2004
Number of blocks in all set partitions of an (n+1)-set. Example: a(2)=10 because the set partitions of {1,2,3} are 1|2|3, 1|23, 12|3, 13|2 and 123, with a total of 10 blocks. - Emeric Deutsch, Nov 13 2006
Number of partitions of n+3 with at least one singleton and with the smallest element in a singleton equal to 2. - Olivier Gérard, Oct 29 2007
See page 29, Theorem 5.6 of my paper on the arXiv: These numbers are the dimensions of the homogeneous components of the operad called ComTrip associated with commutative triplicial algebras. (Triplicial algebras are related to even trees and also to L-algebras, see A006013.) - Philippe Leroux, Nov 17 2007
Number of set partitions of (n+2) elements where two specific elements are clustered separately. Example: a(1)=3 because 1/2/3, 1/23, 13/2 are the 3 set partitions with 1, 2 clustered separately. - Andrey Goder (andy.goder(AT)gmail.com), Dec 17 2007
Equals A008277 * [1,2,3,...], i.e., the product of the Stirling number of the second kind triangle and the natural number vector. a(n+1) = row sums of triangle A137650. - Gary W. Adamson, Jan 31 2008
Prefaced with a "1" = row sums of triangle A152433. - Gary W. Adamson, Dec 04 2008
Equals row sums of triangle A159573. - Gary W. Adamson, Apr 16 2009
Number of embedded coalitions in an (n+1)-person game. - David Yeung (wkyeung(AT)hkbu.edu.hk), May 08 2008
If prefixed with 0, gives first differences of Bell numbers A000110 (cf. A106436). - N. J. A. Sloane, Aug 29 2013
Sum_{n>=0} a(n)/n! = e^(e+1) = 41.19355567... (see A235214). Contrast with e^(e-1) = Sum_{n>=0} A000110(n)/n!. - Richard R. Forberg, Jan 05 2014

Examples

			For example, a(1) counts (12), (1)-2, 1-(2) where dashes separate blocks and the distinguished block is parenthesized.
		

References

  • Olivier Gérard and Karol A. Penson, A budget of set partition statistics, in preparation. Unpublished as of 2017.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

A row or column of the array A108087.
Row sums of triangle A143494. - Wolfdieter Lang, Sep 29 2011. And also of triangle A362924. - N. J. A. Sloane, Aug 10 2023

Programs

  • Maple
    with(combinat): seq(bell(n+2)-bell(n+1),n=0..22); # Emeric Deutsch, Nov 13 2006
    seq(add(binomial(n, k)*(bell(n-k)), k=1..n), n=1..23); # Zerinvary Lajos, Dec 01 2006
    A005493  := proc(n) local a,b,i;
    a := [seq(3,i=1..n)]; b := [seq(2,i=1..n)];
    2^n*exp(-x)*hypergeom(a,b,x); round(evalf(subs(x=1,%),66)) end:
    seq(A005493(n),n=0..22); # Peter Luschny, Mar 30 2011
    BT := proc(n,k) option remember; if n = 0 and k = 0 then 1
    elif k = n then BT(n-1,0) else BT(n,k+1)+BT(n-1,k) fi end:
    A005493 := n -> add(BT(n,k),k=0..n):
    seq(A005493(i),i=0..22); # Peter Luschny, Aug 04 2011
    # For Maple code for r-Bell numbers, etc., see A232472. - N. J. A. Sloane, Nov 27 2013
  • Mathematica
    a=Exp[x]-1; Rest[CoefficientList[Series[a Exp[a],{x,0,20}],x] * Table[n!,{n,0,20}]]
    a[ n_] := If[ n<0, 0, With[ {m = n+1}, m! SeriesCoefficient[ # Exp@# &[ Exp@x - 1], {x, 0, m}]]]; (* Michael Somos, Nov 16 2011 *)
    Differences[BellB[Range[30]]] (* Harvey P. Dale, Oct 16 2014 *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( exp( exp( x + x * O(x^n)) + 2*x - 1), n))}; /* Michael Somos, Oct 09 2002 */
    
  • PARI
    {a(n) = if( n<0, 0, n+=2; subst( polinterpolate( Vec( serlaplace( exp( exp( x + O(x^n)) - 1) - 1))), x, n))}; /* Michael Somos, Oct 07 2003 */
    
  • Python
    # requires python 3.2 or higher. Otherwise use def'n of accumulate in python docs.
    from itertools import accumulate
    A005493_list, blist, b = [], [1], 1
    for _ in range(1001):
        blist = list(accumulate([b]+blist))
        b = blist[-1]
        A005493_list.append(blist[-2])
    # Chai Wah Wu, Sep 02 2014, updated Chai Wah Wu, Sep 20 2014

Formula

a(n-1) = Sum_{k=1..n} k*Stirling2(n, k) for n>=1.
E.g.f.: exp(exp(x) + 2*x - 1). First differences of Bell numbers (if offset 1). - Michael Somos, Oct 09 2002
G.f.: Sum_{k>=0} (x^k/Product_{l=1..k} (1-(l+1)x)). - Ralf Stephan, Apr 18 2004
a(n) = Sum_{i=0..n} 2^(n-i)*B(i)*binomial(n,i) where B(n) = Bell numbers A000110(n). - Fred Lunnon, Aug 04 2007 [Written umbrally, a(n) = (B+2)^n. - N. J. A. Sloane, Feb 07 2009]
Representation as an infinite series: a(n-1) = Sum_{k>=2} (k^n*(k-1)/k!)/exp(1), n=1, 2, ... This is a Dobinski-type summation formula. - Karol A. Penson, Mar 14 2002
Row sums of A011971 (Aitken's array, also called Bell triangle). - Philippe Deléham, Nov 15 2003
a(n) = exp(-1)*Sum_{k>=0} ((k+2)^n)/k!. - Gerald McGarvey, Jun 03 2004
Recurrence: a(n+1) = 1 + Sum_{j=1..n} (1+binomial(n, j))*a(j). - Jon Perry, Apr 25 2005
a(n) = A000296(n+3) - A000296(n+1). - Philippe Deléham, Jul 31 2005
a(n) = B(n+2) - B(n+1), where B(n) are Bell numbers (A000110). - Franklin T. Adams-Watters, Jul 13 2006
a(n) = A123158(n,2). - Philippe Deléham, Oct 06 2006
Binomial transform of Bell numbers 1, 2, 5, 15, 52, 203, 877, 4140, ... (see A000110).
Define f_1(x), f_2(x), ... such that f_1(x)=x*e^x, f_{n+1}(x) = (d/dx)(x*f_n(x)), for n=2,3,.... Then a(n-1) = e^(-1)*f_n(1). - Milan Janjic, May 30 2008
Representation of numbers a(n), n=0,1..., as special values of hypergeometric function of type (n)F(n), in Maple notation: a(n)=exp(-1)*2^n*hypergeom([3,3...3],[2.2...2],1), n=0,1..., i.e., having n parameters all equal to 3 in the numerator, having n parameters all equal to 2 in the denominator and the value of the argument equal to 1. Examples: a(0)= 2^0*evalf(hypergeom([],[],1)/exp(1))=1 a(1)= 2^1*evalf(hypergeom([3],[2],1)/exp(1))=3 a(2)= 2^2*evalf(hypergeom([3,3],[2,2],1)/exp(1))=10 a(3)= 2^3*evalf(hypergeom([3,3,3],[2,2,2],1)/exp(1))=37 a(4)= 2^4*evalf(hypergeom([3,3,3,3],[2,2,2,2],1)/exp(1))=151 a(5)= 2^5*evalf(hypergeom([3,3,3,3,3],[2,2,2,2,2],1)/exp(1)) = 674. - Karol A. Penson, Sep 28 2007
Let A be the upper Hessenberg matrix of order n defined by: A[i,i-1]=-1, A[i,j]=binomial(j-1,i-1), (i <= j), and A[i,j]=0 otherwise. Then, for n >= 1, a(n) = (-1)^(n)charpoly(A,-2). - Milan Janjic, Jul 08 2010
a(n) = D^(n+1)(x*exp(x)) evaluated at x = 0, where D is the operator (1+x)*d/dx. Cf. A003128, A052852 and A009737. - Peter Bala, Nov 25 2011
From Sergei N. Gladkovskii, Oct 11 2012 to Jan 26 2014: (Start)
Continued fractions:
G.f.: 1/U(0) where U(k) = 1 - x*(k+3) - x^2*(k+1)/U(k+1).
G.f.: 1/(U(0)-x) where U(k) = 1 - x - x*(k+1)/(1 - x/U(k+1)).
G.f.: G(0)/(1-x) where G(k) = 1 - 2*x*(k+1)/((2*k+1)*(2*x*k+2*x-1) - x*(2*k+1)*(2*k+3)*(2*x*k+2*x-1)/(x*(2*k+3) - 2*(k+1)*(2*x*k+3*x-1)/G(k+1) )).
G.f.: (G(0) - 1)/(x-1) where G(k) = 1 - 1/(1-2*x-k*x)/(1-x/(x-1/G(k+1) )).
G.f.: -G(0)/x where G(k) = 1 - 1/(1-k*x-x)/(1-x/(x-1/G(k+1) )).
G.f.: 1 - 2/x + (1/x-1)*S where S = sum(k>=0, ( 1 + (1-x)/(1-x-x*k) )*(x/(1-x))^k / prod(i=0..k-1, (1-x-x*i)/(1-x) ) ).
G.f.: (1-x)/x/(G(0)-x) - 1/x where G(k) = 1 - x*(k+1)/(1 - x/G(k+1) ).
G.f.: (1/G(0) - 1)/x^3 where G(k) = 1 - x/(x - 1/(1 + 1/(x*k-1)/G(k+1) )).
G.f.: 1/Q(0), where Q(k)= 1 - 2*x - x/(1 - x*(k+1)/Q(k+1)).
G.f.: G(0)/(1-3*x), where G(k) = 1 - x^2*(k+1)/( x^2*(k+1) - (1 - x*(k+3))*(1 - x*(k+4))/G(k+1) ). (End)
a(n) ~ exp(n/LambertW(n) + 3*LambertW(n)/2 - n - 1) * n^(n + 1/2) / LambertW(n)^(n+1). - Vaclav Kotesovec, Jun 09 2020
a(0) = 1; a(n) = 2 * a(n-1) + Sum_{k=0..n-1} binomial(n-1,k) * a(k). - Ilya Gutkovskiy, Jul 02 2020
a(n) ~ n^2 * Bell(n) / LambertW(n)^2 * (1 - LambertW(n)/n). - Vaclav Kotesovec, Jul 28 2021
a(n) = Sum_{k=0..n} 3^k*A124323(n, k). - Mélika Tebni, Jun 02 2022

Extensions

Definition revised by David Callan, Oct 11 2005

A001858 Number of forests of trees on n labeled nodes.

Original entry on oeis.org

1, 1, 2, 7, 38, 291, 2932, 36961, 561948, 10026505, 205608536, 4767440679, 123373203208, 3525630110107, 110284283006640, 3748357699560961, 137557910094840848, 5421179050350334929, 228359487335194570528, 10239206473040881277575, 486909744862576654283616
Offset: 0

Views

Author

Keywords

Comments

The number of integer lattice points in the permutation polytope of {1,2,...,n}. - Max Alekseyev, Jan 26 2010
Equals the number of score sequences for a tournament on n vertices. See Prop. 7 of the article by Bartels et al., or Example 3.1 in the article by Stanley. - David Radcliffe, Aug 02 2022
Number of labeled acyclic graphs on n vertices. The unlabeled version is A005195. The covering case is A105784, connected A000272. - Gus Wiseman, Apr 29 2024

Examples

			From _Gus Wiseman_, Apr 29 2024: (Start)
Edge-sets of the a(4) = 38 forests:
  {}  {12}  {12,13}  {12,13,14}
      {13}  {12,14}  {12,13,24}
      {14}  {12,23}  {12,13,34}
      {23}  {12,24}  {12,14,23}
      {24}  {12,34}  {12,14,34}
      {34}  {13,14}  {12,23,24}
            {13,23}  {12,23,34}
            {13,24}  {12,24,34}
            {13,34}  {13,14,23}
            {14,23}  {13,14,24}
            {14,24}  {13,23,24}
            {14,34}  {13,23,34}
            {23,24}  {13,24,34}
            {23,34}  {14,23,24}
            {24,34}  {14,23,34}
                     {14,24,34}
(End)
		

References

  • B. Bollobas, Modern Graph Theory, Springer, 1998, p. 290.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

The connected case is A000272, rooted A000169.
The unlabeled version is A005195, connected A000055.
The covering case is A105784, unlabeled A144958.
Row sums of A138464.
For triangles instead of cycles we have A213434, covering A372168.
For a unique cycle we have A372193, covering A372195.
A002807 counts cycles in a complete graph.
A006125 counts simple graphs, unlabeled A000088.
A006129 counts covering graphs, unlabeled A002494.

Programs

  • Maple
    exp(x+x^2+add(n^(n-2)*x^n/n!, n=3..50));
    # second Maple program:
    a:= proc(n) option remember; `if`(n=0, 1, add(
          binomial(n-1, j-1)*j^(j-2)*a(n-j), j=1..n))
        end:
    seq(a(n), n=0..20);  # Alois P. Heinz, Sep 15 2008
    # third Maple program:
    F:= exp(-LambertW(-x)*(1+LambertW(-x)/2)):
    S:= series(F,x,51):
    seq(coeff(S,x,j)*j!, j=0..50); # Robert Israel, May 21 2015
  • Mathematica
    nn=20;t=Sum[n^(n-1)x^n/n!,{n,1,nn}];Range[0,nn]!CoefficientList[ Series[Exp[t-t^2/2],{x,0,nn}],x] (* Geoffrey Critzer, Sep 05 2012 *)
    nmax = 20; CoefficientList[Series[-LambertW[-x]/(x*E^(LambertW[-x]^2/2)), {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Jul 19 2019 *)
  • PARI
    a(n)=if(n<0,0,sum(m=0,n,sum(j=0,m,binomial(m,j)*binomial(n-1,n-m-j)*n^(n-m-j)*(m+j)!/(-2)^j)/m!)) /* Michael Somos, Aug 22 2002 */

Formula

E.g.f.: exp( Sum_{n>=1} n^(n-2)*x^n/n! ). This implies (by a theorem of Wright) that a(n) ~ exp(1/2)*n^(n-2). - N. J. A. Sloane, May 12 2008 [Corrected by Philippe Flajolet, Aug 17 2008]
E.g.f.: exp(T - T^2/2), where T = T(x) = Sum_{n>=1} n^(n-1)*x^n/n! is Euler's tree function (see A000169). - Len Smiley, Dec 12 2001
Shifts 1 place left under the hyperbinomial transform (cf. A088956). - Paul D. Hanna, Nov 03 2003
a(0) = 1, a(n) = Sum_{j=0..n-1} C(n-1,j) (j+1)^(j-1) a(n-1-j) if n>0. - Alois P. Heinz, Sep 15 2008

Extensions

More terms from Michael Somos, Aug 22 2002

A027710 Number of ways of placing n labeled balls into n unlabeled (but 3-colored) boxes.

Original entry on oeis.org

1, 3, 12, 57, 309, 1866, 12351, 88563, 681870, 5597643, 48718569, 447428856, 4318854429, 43666895343, 461101962108, 5072054649573, 57986312752497, 687610920335610, 8442056059773267, 107135148331162767, 1403300026585387686, 18946012544520590991
Offset: 0

Views

Author

George Yuhasz (gyuhasz(AT)vt.edu) and John W. Layman

Keywords

Comments

Binomial transform of this sequence is A078940 and a(n+1) = 3*A078940(n). - Paul D. Hanna, Dec 08 2003
First column of the cube of the matrix exp(P)/exp(1) given in A011971. - Gottfried Helms, Mar 30 2007. Base matrix in A011971, second power in A078937, third power in A078938, fourth power in A078939.
The number of ways of putting n labeled balls into a set of bags and then putting the bags into 3 labeled boxes. - Peter Bala, Mar 23 2013

Crossrefs

Programs

  • Maple
    b:= proc(n, m) option remember; `if`(n=0,
          1, m*b(n-1, m)+3*b(n-1, m+1))
        end:
    a:= n-> b(n, 0):
    seq(a(n), n=0..27);  # Alois P. Heinz, Aug 03 2021
  • Mathematica
    colors=3; Array[ bell, 25 ]; For[ x=1, x<=25, x++, bell[ x ]=0 ]; bell[ 1 ]=colors;
    Print[ "1 ", colors ]; For[ n=2, n<=25, n++, bell[ n ]=colors*bell[ n-1 ];
    For[ i=1, n-i>1, i++, bell[ n-i ]=bell[ n-i ]*(n-i)+colors*bell[ n-i-1 ] ];
    bellsum=0; For[ t=0, tVaclav Kotesovec, Mar 12 2014 *)
  • PARI
    a(n)=if(n<0,0,n!*polcoeff(exp(3*(exp(x+x*O(x^n))-1)),n))
    
  • Sage
    from sage.combinat.expnums import expnums2
    expnums(22, 3) # Zerinvary Lajos, Jun 26 2008

Formula

E.g.f.: exp {3(e^x-1)}. - Michael Somos, Oct 18 2002
a(n) = exp(-3)*Sum_{k>=0} 3^k*k^n/k!. - Benoit Cloitre, Sep 25 2003
G.f.: 3*(x/(1-x))*A(x/(1-x)) = A(x) - 1; thrice the binomial transform equals the sequence shifted one place left. - Paul D. Hanna, Dec 08 2003
a(n) = Sum_{k = 0..n} 3^k*A048993(n, k); A048993: Stirling2 numbers. - Philippe Deléham, May 09 2004
PE=exp(matpascal(5))/exp(1); A = PE^3; a(n)= A[ n,1 ] with exact integer arithmetic: PE=exp(matpascal(5)-matid(6)); A = PE^3; a(n)=A[ n,1]. - Gottfried Helms, Apr 08 2007
G.f.: (G(0) - 1)/(x-1)/3 where G(k) = 1 - 3/(1-k*x)/(1-x/(x-1/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 16 2013
G.f.: T(0)/(1-3*x), where T(k) = 1 - 3*x^2*(k+1)/( 3*x^2*(k+1) - (1-3*x-x*k)*(1-4*x-x*k)/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Oct 24 2013
a(n) ~ n^n * exp(n/LambertW(n/3)-3-n) / (sqrt(1+LambertW(n/3)) * LambertW(n/3)^n). - Vaclav Kotesovec, Mar 12 2014
G.f.: Sum_{j>=0} 3^j*x^j / Product_{k=1..j} (1 - k*x). - Ilya Gutkovskiy, Apr 07 2019

Extensions

Entry revised by N. J. A. Sloane, Apr 25 2007

A078944 First column of A078939, the fourth power of lower triangular matrix A056857.

Original entry on oeis.org

1, 4, 20, 116, 756, 5428, 42356, 355636, 3188340, 30333492, 304716148, 3218555700, 35618229364, 411717043252, 4957730174836, 62045057731892, 805323357485684, 10820999695801908, 150271018666120564, 2153476417340487476
Offset: 0

Views

Author

Paul D. Hanna, Dec 18 2002

Keywords

Comments

Also, the number of ways of placing n labeled balls into n unlabeled (but 4-colored) boxes. Binomial transform of this sequence is A078945 and a(n+1) = 4*A078945(n). - Paul D. Hanna, Dec 08 2003
First column of PE^4, where PE is given in A011971, second power in A078937, third power in A078938, fourth power in A078939. - Gottfried Helms, Apr 08 2007
The number of ways of putting n labeled balls into a set of bags and then putting the bags into 4 labeled boxes. - Peter Bala, Mar 23 2013
Exponential self-convolution of A001861. - Vladimir Reshetnikov, Oct 06 2016

Crossrefs

Programs

  • Maple
    A056857 := proc(n,c) combinat[bell](n-1-c)*binomial(n-1,c) ; end: A078937 := proc(n,c) add( A056857(n,k)*A056857(k+1,c),k=0..n) ; end: A078938 := proc(n,c) add( A078937(n,k)*A056857(k+1,c),k=0..n) ; end: A078939 := proc(n,c) add( A078938(n,k)*A056857(k+1,c),k=0..n) ; end: A078944 := proc(n) A078939(n+1,0) ; end: seq(A078944(n),n=0..25) ; # R. J. Mathar, May 30 2008
    # second Maple program:
    b:= proc(n, m) option remember; `if`(n=0, 4^m,
          add(b(n-1, max(m, j)), j=1..m+1))
        end:
    a:= n-> b(n, 0):
    seq(a(n), n=0..25);  # Alois P. Heinz, Aug 03 2021
  • Mathematica
    Table[n!, {n, 0, 20}]CoefficientList[Series[E^(4E^x-4), {x, 0, 20}], x]
    Table[BellB[n,4],{n,0,20}] (* Vaclav Kotesovec, Mar 12 2014 *)
    With[{nn=20},CoefficientList[Series[Exp[4(Exp[x]-1)],{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, May 03 2022 *)
  • Sage
    expnums(20, 4) # Zerinvary Lajos, Jun 26 2008

Formula

PE=exp(matpascal(5))/exp(1); A = PE^4; a(n)= A[ n,1 ] with exact integer arithmetic: PE=exp(matpascal(5)-matid(6)); A = PE^4; a(n)=A[ n,1]. - Gottfried Helms, Apr 08 2007
E.g.f.: exp(4*(exp(x)-1)).
a(n) = exp(-4)*Sum_{k>=0} 4^k*k^n/k!. - Benoit Cloitre, Sep 25 2003
G.f.: 4*(x/(1-x))*A(x/(1-x)) = A(x) - 1; four times the binomial transform equals this sequence shifted one place left. - Paul D. Hanna, Dec 08 2003
a(n) = Sum_{k = 0..n} 4^k*A048993(n, k); A048993: Stirling2 numbers. - Philippe Deléham, May 09 2004
G.f.: (G(0) - 1)/(x-1)/4 where G(k) = 1 - 4/(1-k*x)/(1-x/(x-1/G(k+1))); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 16 2013
G.f.: T(0)/(1-4*x), where T(k) = 1 - 4*x^2*(k+1)/(4*x^2*(k+1) - (1-(k+4)*x)*(1-(k+5)*x)/T(k+1)); (continued fraction). - Sergei N. Gladkovskii, Oct 28 2013
a(n) ~ n^n * exp(n/LambertW(n/4)-4-n) / (sqrt(1+LambertW(n/4)) * LambertW(n/4)^n). - Vaclav Kotesovec, Mar 12 2014
G.f.: Sum_{j>=0} 4^j*x^j / Product_{k=1..j} (1 - k*x). - Ilya Gutkovskiy, Apr 07 2019

Extensions

More terms from R. J. Mathar, May 30 2008
Edited by N. J. A. Sloane, Jul 02 2008 at the suggestion of R. J. Mathar
Showing 1-10 of 76 results. Next