A108087 -id:A108087 - cp's OEIS frontend

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

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

N. J. A. Sloane, Simon Plouffe

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

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

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).

Chai Wah Wu, Table of n, a(n) for n = 0..574 (first 101 terms from T. D. Noe)
V. E. Adler, Set partitions and integrable hierarchies, arXiv:1510.02900 [nlin.SI], 2015.
M. Brookes, J. East, C. Miller, J.D. Mitchell, and N. Ruskuc, Heights of one- and two-sided congruence lattices of semigroups, arXiv:2310.08229 [math.GR], 2023.
Giulio Cerbai, Modified ascent sequences and Bell numbers, arXiv:2305.10820 [math.CO], 2023. See p. 11.
Martin Cohn, Shimon Even, Karl Menger, Jr., and Philip K. Hooper, On the Number of Partitionings of a Set of n Distinct Objects, Amer. Math. Monthly 69 (1962), no. 8, 782--785. MR1531841.
Martin Cohn, Shimon Even, Karl Menger, Jr., and Philip K. Hooper, On the Number of Partitionings of a Set of n Distinct Objects, Amer. Math. Monthly 69 (1962), no. 8, 782--785. MR1531841. [Annotated scanned copy]
R. B. Corcino and C. B. Corcino, On generalized Bell polynomials, Discrete Dyn. Nat. Soc. Article ID: 623456 (2011).
C. B. Corcino and R. B. Corcino, An asymptotic formula for r-Bell numbers with real arguments, ISRN Discrete Math, 2013 (2013), Article ID 274697.
Gábor Czédli, Lattices with lots of congruence energy, arXiv:2205.02294 [math.RA], 2022.
A. Dil and V. Kurt, Polynomials related to harmonic numbers and evaluation of harmonic number series II, Appl. Anal. Discrete Math. 5 (2011), 212-229.
Eldar Fischer, Johann A. Makowsky, and Vsevolod Rakita, MC-finiteness of restricted set partition functions, arXiv:2302.08265 [math.CO], 2023.
A. L. L. Gao, S. Kitaev, and P. B. Zhang, On pattern avoiding indecomposable permutations, arXiv:1605.05490 [math.CO], 2016.
S. Getu et al., How to guess a generating function, SIAM J. Discrete Math., 5 (1992), 497-499.
S. K. Ghosal and J. K. Mandal, Stirling Transform Based Color Image Authentication, Procedia Technology, 2013 Volume 10, 2013, Pages 95-104.
Alain Hertz, Hadrien Mélot, Sébastien Bonte, Gauvain Devillez, and Pierre Hauweele, Upper bounds on the average number of colors in the non-equivalent colorings of a graph, arXiv:2105.01120 [math.CO], 2021.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 152
R. Jakimczuk, Successive Derivatives and Integer Sequences, J. Int. Seq. 14 (2011) # 11.7.3.
S. Kitaev, Generalized pattern avoidance with additional restrictions, Sem. Lothar. Combinat. B48e (2003).
S. Kitaev and T. Mansour, Simultaneous avoidance of generalized patterns, arXiv:math/0205182 [math.CO], 2002.
J. W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Philippe Leroux, L-algebras, triplicial-algebras, within an equivalence of categories motivated by graphs (formerly "An equivalence of categories motivated by weighted directed graphs"), arXiv:0709.3453 [math-ph], 2007-2008.
Toufik Mansour and Mark Shattuck, A recurrence related to the Bell numbers, INTEGERS 11 (2011), #A67.
I. Mezo, The r-Bell numbers, J. Integer Seq. 14(1) (2011), Article 11.1.1.
I. Mezo, Periodicity of the last digits of some combinatorial sequences, arXiv preprint arXiv:1308.1637 [math.CO], 2013.
A. M. Odlyzko, Asymptotic enumeration methods, pp. 1063-1229 of R. L. Graham et al., eds., Handbook of Combinatorics, 1995; see Example 12.16 (pdf, ps)
Jocelyn Quaintance and Harris Kwong, A combinatorial interpretation of the Catalan and Bell number difference tables, Integers, 13 (2013), #A29.
J. Riordan, Letter, Oct 31 1977
Eric Weisstein's World of Mathematics, Stirling Transform.
Eric Weisstein's World of Mathematics, Bell Triangle
E. G. Whitehead, Jr., Stirling number identities from chromatic polynomials, J. Combin. Theory, A 24 (1978), 314-317.
David W. K. Yeung, Recursive sequence identifying the number of embedded coalitions, International Game Theory Review 10 (1) (2008), 129-136.

Cf. A000110, A005494, A008277, A011971, A011968, A049020, A106436, A124323, A137650, A152433, A159573.
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

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

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 *)

{a(n) = if( n<0, 0, n! * polcoeff( exp( exp( x + x * O(x^n)) + 2*x - 1), n))}; /* Michael Somos, Oct 09 2002 */

{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 */

# 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

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

Definition revised by David Callan, Oct 11 2005

A005494 3-Bell numbers: E.g.f.: exp(3*z + exp(z) - 1).

Original entry on oeis.org

1, 4, 17, 77, 372, 1915, 10481, 60814, 372939, 2409837, 16360786, 116393205, 865549453, 6713065156, 54190360453, 454442481041, 3952241526188, 35590085232519, 331362825860749, 3185554606447814, 31581598272055879, 322516283206446897, 3389017736055752178

Offset: 0

N. J. A. Sloane, Simon Plouffe

nonn

For further information, references, programs, etc. for r-Bell numbers see A005493. - N. J. A. Sloane, Nov 27 2013
From expansion of falling factorials (binomial transform of A005493).
Row sums of Sheffer triangle (exp(3*x), exp(x)-1). - Wolfdieter Lang, Sep 29 2011

			G.f. = 1 + 4*x + 17*x^2 + 77*x^3 + 372*x^4 + 1915*x^5 + 10481*x^6 + 60814*x^7 + ...

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Eldar Fischer, Johann A. Makowsky, and Vsevolod Rakita, MC-finiteness of restricted set partition functions, arXiv:2302.08265 [math.CO], 2023.
R. Jakimczuk, Successive Derivatives and Integer Sequences, J. Int. Seq. 14 (2011) # 11.7.3.
J. W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Toufik Mansour and Mark Shattuck, A recurrence related to the Bell numbers, INTEGERS 11 (2011), #A67.
I. Mezo, The r-Bell numbers, J. Int. Seq. 14 (2011) # 11.1.1.
J. Riordan, Letter, Oct 31 1977
N. J. A. Sloane, Transforms
Earl Glen Whitehead Jr., Stirling number identities from chromatic polynomials, J. Combin. Theory, A 24 (1978), 314-317.

Cf. A000110, A005493, A045379, A108087, A124323, A196834.

A row or column of the array A108087.

A005494:= func< n | (&+[Binomial(n,j)*3^(n-j)*Bell(j): j in [0..n]]) >;
[A005494(n): n in [0..30]]; // G. C. Greubel, Dec 01 2022

seq(add(3^(n-i)*combinat:-bell(i)*binomial(n,i),i=0..n), n=0..50); # Robert Israel, Dec 16 2014
# second Maple program:
b:= proc(n, m) option remember; `if`(n=0,
      m^2, m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n+1, 0)-b(n, 0):
seq(a(n), n=0..23);  # Alois P. Heinz, Aug 03 2025

Range[0, 40]! CoefficientList[Series[Exp[3 x + Exp[x] - 1], {x, 0, 40}], x] (* Vincenzo Librandi, Mar 04 2014 *)

def A005494(n): return sum( 3^(n-j)*bell_number(j)*binomial(n,j) for j in range(n+1))
[A005494(n) for n in range(31)] # G. C. Greubel, Dec 01 2022

a(n) = Sum_{i=0..n} 3^(n-i)*B(i)*binomial(n,i) where B(n) = Bell numbers A000110(n). - Fred Lunnon, Aug 04 2007
a(n) = exp(-1)*Sum_{k>=0} ((k+3)^n)/k!. - Gerald McGarvey, Jun 03 2004. May be rewritten as a(n) = Sum_{k>=3} (k^n*(k-1)*(k-2)/k!)/exp(1), which is a Dobinski-type relation for this sequence. - Karol A. Penson, Aug 18 2006
Define f_1(x), f_2(x), ... such that f_1(x) = x^2*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
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,-3). - Milan Janjic, Jul 08 2010
a(n) = Sum_{k=3..n+3} A143495(n+3,k), n >= 0. - Wolfdieter Lang, Sep 29 2011
G.f.: 1/U(0) where U(k)= 1 - x*(k+4) - x^2*(k+1)/U(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Oct 11 2012
G.f.: Sum_{k>0} x^(k-1) / ((1 - 3*x) * (1 - 4*x) * ... * (1 - (k+2)*x)). - Michael Somos, Feb 26 2014
G.f.: Sum_{k>0} k * x^(k-1) / ((1 - 2*x) * (1 - 3*x) * ... * (1 - (k+1)*x)). - Michael Somos, Feb 26 2014
a(n) ~ exp(n/LambertW(n) - n - 1) * n^(n + 3) / LambertW(n)^(n + 7/2). - Vaclav Kotesovec, Jun 10 2020
a(0) = 1; a(n) = 3 * a(n-1) + Sum_{k=0..n-1} binomial(n-1,k) * a(k). - Ilya Gutkovskiy, Jul 02 2020
a(n) = Sum_{k=0..n} 4^k*A124323(n, k). - Mélika Tebni, Jun 10 2022

A045379 Expansion of e.g.f.: exp(4*z + exp(z) - 1).

Original entry on oeis.org

1, 5, 26, 141, 799, 4736, 29371, 190497, 1291020, 9131275, 67310847, 516369838, 4116416797, 34051164985, 291871399682, 2588914083065, 23733360653955, 224592570163192, 2191466128865567, 22024934452712437, 227771488390279260

Offset: 0

N. J. A. Sloane

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..500
R. Jakimczuk, Successive Derivatives and Integer Sequences, J. Int. Seq. 14 (2011) # 11.7.3.
J. W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
I. Mezo, The r-Bell numbers, J. Int. Seq. 14 (2011) # 11.1.1.

Cf. A000110, A005493, A005494, A108087, A196834.

Equals the row sums of triangle A143496. - Wolfdieter Lang, Sep 29 2011

A045379:= func< n | (&+[Binomial(n,j)*4^(n-j)*Bell(j): j in [0..n]]) >;
[A045379(n): n in [0..30]]; // G. C. Greubel, Dec 01 2022

a[0]= 1; a[n_]:= a[n]= 4*a[n-1] +Sum[Binomial[n-1, k]*a[k], {k,0,n-1}]; Array[a, 21, 0] (* Amiram Eldar, Jul 03 2020 *)

def A045379(n): return sum( 4^(n-j)*bell_number(j)*binomial(n,j) for j in range(n+1))
[A045379(n) for n in range(31)] # G. C. Greubel, Dec 01 2022

a(n) = exp(-1)*Sum_{k>=0} ((k+4)^n)/k!. - Gerald McGarvey, Jun 03 2004
A recursive formula to compute some integer sequences (including A000110, A005493, A005494 and the present sequence). Define G(n, m), where n, m >= 0, as follows: G(0, m) = 1; G(n, m) = G(n-1, m) * (m+1) + G(n-1, m+1), where n > 0. Then G(n, 0) = A000110(n+1); G(n, 1) = A005493(n+1); G(n, 2) = A005494(n+1); G(n, 3) = A045379(n+1). - Alexey Andreev (ava12(AT)nm.ru), Jan 05 2006
Define f_1(x), f_2(x), ... such that f_1(x)=x^3*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
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,-4). - Milan Janjic, Jul 08 2010
G.f.: 1/U(0) where U(k) = 1 - x*(k+5) - x^2*(k+1)/U(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Oct 11 2012
a(n) ~ exp(n/LambertW(n) - n - 1) * n^(n + 4) / LambertW(n)^(n + 9/2). - Vaclav Kotesovec, Jun 10 2020
a(0) = 1; a(n) = 4 * a(n-1) + Sum_{k=0..n-1} binomial(n-1,k) * a(k). - Ilya Gutkovskiy, Jul 02 2020
a(n) = Sum_{j=0..n} binomial(n, j)*4^(n-j)*A000110(j). - G. C. Greubel, Dec 01 2022

A134980 a(n) = Sum_{k=0..n} binomial(n,k)n^(n-k)A000110(k).

Original entry on oeis.org

1, 2, 10, 77, 799, 10427, 163967, 3017562, 63625324, 1512354975, 40012800675, 1166271373797, 37134022033885, 1282405154139046, 47745103281852282, 1906411492286148245, 81267367663098939459, 3683790958912910588623, 176937226305157687076779

Offset: 0

Vladeta Jovovic, Feb 04 2008

Main diagonal of array "The first r-Bell numbers", p.3 of Mezo, A108087. - Jonathan Vos Post, Sep 25 2009

Number of partitions of [2n] where at least n blocks contain their own index element. a(2) = 10: 134|2, 13|24, 13|2|4, 14|23, 1|234, 1|23|4, 14|2|3, 1|24|3, 1|2|34, 1|2|3|4. - Alois P. Heinz, Jan 07 2022

Robert Israel, Table of n, a(n) for n = 0..385
R. Jakimczuk, Successive Derivatives and Integer Sequences, J. Int. Seq. 14 (2011) # 11.7.3.
Istvan Mezo, The r-Bell numbers, arXiv:0909.4417 [math.CO], 2009-2010.
I. Mezo, The r-Bell numbers, J. Int. Seq. 14 (2011) # 11.1.1.

Main diagonal of A108087.

Cf. A000110.

with(combinat): a:= n-> add(binomial(n, k)*n^(n-k)*bell(k), k=0..n):
seq(a(n), n=0..20); # Emeric Deutsch, Mar 02 2008
# Alternate:
g:= proc(n) local S;
  S:= series(exp(exp(x)+n*x-1),x,n+1);
n!*coeff(S,x,n);
end proc:
map(g, [$0..30]); # Robert Israel, Sep 29 2017
# third Maple program:
b:= proc(n, k) option remember; `if`(n=0, 1,
      k*b(n-1, k)+ b(n-1, k+1))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..20);  # Alois P. Heinz, Aug 04 2021

a[n_] := n!*SeriesCoefficient[Exp[Exp[x] + n*x - 1], {x, 0, n}]; Array[a, 19, 0] (* Jean-François Alcover, Sep 28 2017, after Ilya Gutkovskiy *)
Join[{1}, Table[Sum[Binomial[n,k]*n^(n-k)*BellB[k], {k,0,n}], {n,1,20}]] (* Vaclav Kotesovec, Jun 09 2020 *)

def A134980(n):
    return add(binomial(n, k)*n^(n-k)*bell_number(k) for k in (0..n))
[A134980(n) for n in (0..18)]  # Peter Luschny, May 05 2013

a(n) = exp(-1)*Sum_{k>=0} (n+k)^n/k!.
E.g.f.: A(x) = exp(-1)*Sum_{k>=0} (1+k*x)^k/k!.
a(n) = Sum_{k=0..n} Stirling1(n,k)*A000110(n+k). - Vladeta Jovovic, Nov 08 2009
a(n) = n! * [x^n] exp(exp(x) + n*x - 1). - Ilya Gutkovskiy, Sep 26 2017
a(n) ~ exp(exp(1) - 1) * n^n. - Vaclav Kotesovec, Jun 09 2020

More terms from Emeric Deutsch, Mar 02 2008

A259691 Triangle read by rows: T(n,k) number of arrangements of non-attacking rooks on an n X n right triangular board where the top rook is in row k (n >= 0, 1 <= k <= n+1).

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 5, 6, 3, 1, 15, 20, 12, 4, 1, 52, 74, 51, 20, 5, 1, 203, 302, 231, 104, 30, 6, 1, 877, 1348, 1116, 564, 185, 42, 7, 1, 4140, 6526, 5745, 3196, 1175, 300, 56, 8, 1, 21147, 34014, 31443, 18944, 7700, 2190, 455, 72, 9, 1

Offset: 0

N. J. A. Sloane, Jul 05 2015

Another version of A056857.
See Becker (1948/49) for precise definition.
The case of n=k+1 corresponds to the empty board where there is no top rook. - Andrew Howroyd, Jun 13 2017
T(n-1,k) is the number of partitions of [n] where exactly k blocks contain their own index element.  T(3,2) = 6: 134|2, 13|24, 13|2|4, 14|23, 1|234, 1|23|4. - Alois P. Heinz, Jan 07 2022

			Triangle begins:
    1;
    1,   1;
    2,   2,   1;
    5,   6,   3,   1;
   15,  20,  12,   4,  1;
   52,  74,  51,  20,  5, 1;
  203, 302, 231, 104, 30, 6, 1;
  ...
From _Andrew Howroyd_, Jun 13 2017: (Start)
For n=3 the 5 solutions with the top rook in row 1 are:
  x      x      x      x      x
  . .    . .    . .    . x    . x
  . . .  . . x  . x .  . . .  . . x
For n=3 the 6 solutions with the top rook in row 2 are:
  .      .      .      .      .      .
  x .    x .    x .    . x    . x    . x
  . . .  . x .  . . x  . . .  x . .  . . x
(End)

Alois P. Heinz, Rows n = 0..140, flattened (first 49 rows from Andrew Howroyd)
H. W. Becker, Rooks and rhymes, Math. Mag., 22 (1948/49), 23-26. See Table II.
H. W. Becker, Rooks and rhymes, Math. Mag., 22 (1948/49), 23-26. [Annotated scanned copy]
Fufa Beyene, Jörgen Backelin, Roberto Mantaci, and Samuel A. Fufa, Set Partitions and Other Bell Number Enumerated Objects, J. Int. Seq., Vol. 26 (2023), Article 23.1.8.
Giulio Cerbai, Modified ascent sequences and Bell numbers, arXiv:2305.10820 [math.CO], 2023. See p. 17.

First column is A000110.
Row sums are A000110(n+1).
Cf. A056857, A259697, A108087, A347420, A350589.

b:= proc(n, m) option remember; `if`(n=0, 1,
     `if`(n<0, 1/m, m*b(n-1, m)+b(n-1, m+1)))
    end:
T:= (n, k)-> k*b(n-k, k):
seq(seq(T(n, k), k=1..n+1), n=0..10);  # Alois P. Heinz, Jan 07 2022

T[n_, k_] := If[k>n, 1, k*Sum[Binomial[n-k, i]*k^i*BellB[n-k-i], {i, 0, n - k}]];
Table[T[n, k], {n, 0, 10}, {k, 1, n+1}] // Flatten (* Jean-François Alcover, Jul 03 2018, after Andrew Howroyd *)

bell(n) = sum(k=0, n, stirling(n, k, 2));
T(n,k) = if(k>n, 1, k*sum(i=0,n-k, binomial(n-k,i) * k^i * bell(n-k-i)));
for(n=0,6, for(k=1,n+1, print1(T(n,k),", ")); print) \\ Andrew Howroyd, Jun 13 2017

T(n,n+1) = 1, T(n,k) = k*Sum_{i=0..n-k} binomial(n-k,i) * k^i * Bell(n-k-i) for k<=n. - Andrew Howroyd, Jun 13 2017
From Alois P. Heinz, Jan 07 2022: (Start)
T(n,k) = k * A108087(n-k,k) for 1 <= k <= n.
Sum_{k=1..n+1} k * T(n,k) = A350589(n+1).
Sum_{k=1..n+1} (k+1) * T(n,k) = A347420(n+1). (End)

Name edited and terms a(28) and beyond from Andrew Howroyd, Jun 13 2017

A347420 Number of partitions of [n] where the first k elements are marked (0 <= k <= n) and at least k blocks contain their own index.

Original entry on oeis.org

1, 2, 5, 14, 45, 164, 667, 2986, 14551, 76498, 430747, 2582448, 16403029, 109918746, 774289169, 5715471606, 44087879137, 354521950932, 2965359744447, 25749723493074, 231719153184019, 2157494726318234, 20753996174222511, 205985762120971168, 2106795754056142537

Offset: 0

Alois P. Heinz, Jan 05 2022

nonn

			a(3) = 14 = 5 + 5 + 3 + 1: 123, 12|3, 13|2, 1|23, 1|2|3, 1'23, 1'2|3, 1'3|2, 1'|23, 1'|2|3, 1'3|2', 1'|2'3, 1'|2'|3, 1'|2'|3'.

Alois P. Heinz, Table of n, a(n) for n = 0..575
Wikipedia, Partition of a set

Antidiagonal sums of A108087.

Cf. A000110, A005490, A059841, A259691, A350589, A361380.

b:= proc(n, m) option remember;
     `if`(n=0, 1, b(n-1, m+1)+m*b(n-1, m))
    end:
a:= n-> add(b(i, n-i), i=0..n):
seq(a(n), n=0..25);

b[n_, m_] := b[n, m] = If[n == 0, 1, b[n - 1, m + 1] + m*b[n - 1, m]];
a[n_] := Sum[b[i, n - i], {i, 0, n}];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Jan 11 2022, after Alois P. Heinz *)

a(n) = Sum_{k=0..n} A108087(n-k,k).
a(n) = 1 + A005490(n).
a(n) = A000110(n) + Sum_{k=1..n} k * A259691(n-1,k).
a(n) = Sum_{k=1..n} (k+1) * A259691(n-1,k).
a(n) = A000110(n) + A350589(n).
a(n) mod 2 = A059841(n).

A005490 Number of partitions of [n] where the first k elements are marked (0 <= k <= n-1) and at least k blocks contain their own index.

Original entry on oeis.org

1, 4, 13, 44, 163, 666, 2985, 14550, 76497, 430746, 2582447, 16403028, 109918745, 774289168, 5715471605, 44087879136, 354521950931, 2965359744446, 25749723493073, 231719153184018, 2157494726318233, 20753996174222510, 205985762120971167, 2106795754056142536

Offset: 1

N. J. A. Sloane, Simon Plouffe

nonn

Old name was: From expansion of falling factorials.

			a(3) = 13 = 5 + 5 + 3: 123, 12|3, 13|2, 1|23, 1|2|3, 1'23, 1'2|3, 1'3|2, 1'|23, 1'|2|3, 1'3|2', 1'|2'3, 1'|2'|3.

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

E. G. Whitehead, Jr., Stirling number identities from chromatic polynomials, J. Combin. Theory, A 24 (1978), 314-317.

Cf. A000035, A000110, A108087, A259691, A347420.

b:= proc(n, m) option remember;
     `if`(n=0, 1, b(n-1, m+1)+m*b(n-1, m))
    end:
a:= n-> add(b(n-k, k), k=0..n-1):
seq(a(n), n=1..24);  # Alois P. Heinz, Jan 05 2022

b[n_, m_] := b[n, m] = If[n == 0, 1, b[n - 1, m + 1] + m*b[n - 1, m]];
a[n_] := Sum[b[n - k, k], {k, 0, n - 1}];
Table[a[n], {n, 1, 24}] (* Jean-François Alcover, Apr 24 2022, after Alois P. Heinz *)

a(n) = Sum_{i=1..n} b(n, i) where b(n, 1) = n and b(n+1, i+1) = (n-i) * b(n, i) + b(n+1, i) [From Whitehead]. - Sean A. Irvine, Jul 01 2016
From Alois P. Heinz, Jan 05 2022: (Start)
a(n) = Sum_{k=0..n-1} A108087(n-k,k).
a(n) = A000110(n) + Sum_{k=1..n-1} A259691(n,k)/k.
a(n) = A347420(n) - 1.
a(n) mod 2 = n mod 2 = A000035(n). (End)

More terms from Sean A. Irvine, Jul 01 2016

New name from Alois P. Heinz, Jan 07 2022

A350589 Sum over all partitions of [n] of the number of blocks containing their own index.

Original entry on oeis.org

0, 1, 3, 9, 30, 112, 464, 2109, 10411, 55351, 314772, 1903878, 12189432, 82274309, 583389847, 4332513061, 33607736990, 271657081128, 2283282938288, 19916981288017, 179994994948647, 1682624910161483, 16247280435775188, 161833756265886822, 1660836884761337248

Offset: 0

Alois P. Heinz, Jan 07 2022

nonn

Also the number of partitions of [n] where the first k elements are marked (1 <= k <= n) and at least k blocks contain their own index: a(3) = 9 = 5 + 3 + 1: 1'23, 1'2|3, 1'3|2, 1'|23, 1'|2|3, 1'3|2', 1'|2'3, 1'|2'|3, 1'|2'|3'.

			a(3) = 9 = 1 + 1 + 2 + 2 + 3: 123, 12|3, 13|2, 1|23, 1|2|3.

Alois P. Heinz, Table of n, a(n) for n = 0..575
Wikipedia, Partition of a set

Cf. A000110, A005490, A088911, A108087, A259691, A347420, A350648.

b:= proc(n, m) option remember;
     `if`(n=0, 1, b(n-1, m+1)+m*b(n-1, m))
    end:
a:= n-> add(b(n-i, i), i=1..n):
seq(a(n), n=0..25);

b[n_, m_] := b[n, m] = If[n == 0, 1, b[n - 1, m + 1] + m*b[n - 1, m]];
a[n_] := Sum[b[n - i, i], {i, 1, n}];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Jan 11 2022, after Alois P. Heinz *)

a(n) = Sum_{k=1..n} A108087(n-k,k).
a(n) = Sum_{k=1..n} k * A259691(n-1,k).
a(n) = Sum_{k=1..n} A259691(n,k)/k.
a(n) = A347420(n) - A000110(n).
a(n) = 1 + A005490(n) - A000110(n).
a(n) mod 2 = A088911(n+5).

A361781 A(n,k) is the n-th term of the k-th inverse binomial transform of the Bell numbers (A000110); square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 0, 2, 1, -1, 1, 5, 1, -2, 2, 1, 15, 1, -3, 5, -3, 4, 52, 1, -4, 10, -13, 7, 11, 203, 1, -5, 17, -35, 36, -10, 41, 877, 1, -6, 26, -75, 127, -101, 31, 162, 4140, 1, -7, 37, -139, 340, -472, 293, -21, 715, 21147, 1, -8, 50, -233, 759, -1573, 1787, -848, 204, 3425, 115975

Offset: 0

Alois P. Heinz, Mar 23 2023

			Square array A(n,k) begins:
    1,   1,   1,    1,     1,      1,       1,       1, ...
    1,   0,  -1,   -2,    -3,     -4,      -5,      -6, ...
    2,   1,   2,    5,    10,     17,      26,      37, ...
    5,   1,  -3,  -13,   -35,    -75,    -139,    -233, ...
   15,   4,   7,   36,   127,    340,     759,    1492, ...
   52,  11, -10, -101,  -472,  -1573,   -4214,   -9685, ...
  203,  41,  31,  293,  1787,   7393,   23711,   63581, ...
  877, 162, -21, -848, -6855, -35178, -134873, -421356, ...

Alois P. Heinz, Antidiagonals n = 0..150, flattened

Columns k=0-5 give: A000110, A000296, A126617, A346738, A346739, A346740.
Rows n=0-2 give: A000012, A024000, A160457.
Main diagonal gives A290219.
Antidiagonal sums give A361380.
Cf. A108087.

T:= func< n,k | (&+[(-k)^j*Binomial(n,j)*Bell(n-j): j in [0..n]]) >;
A361781:= func< n,k | T(k, n-k) >;
[A361781(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Jun 12 2024

A:= proc(n, k) option remember; uses combinat;
      add(binomial(n, j)*(-k)^j*bell(n-j), j=0..n)
    end:
seq(seq(A(n, d-n), n=0..d), d=0..10);
# second Maple program:
b:= proc(n, m) option remember;
     `if`(n=0, 1, b(n-1, m+1)+m*b(n-1, m))
    end:
A:= (n, k)-> b(n, -k):
seq(seq(A(n, d-n), n=0..d), d=0..10);

T[n_, k_]:= T[n, k]= If[k==0, BellB[n], Sum[(-k)^j*Binomial[n,j]*BellB[n-j], {j,0,n}]];
A361781[n_, k_]= T[k, n-k];
Table[A361781[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Jun 12 2024 *)

def T(n,k): return sum( (-k)^j*binomial(n,j)*bell_number(n-j) for j in range(n+1))
def A361781(n, k): return T(k, n-k)
flatten([[A361781(n, k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jun 12 2024

E.g.f. of column k: exp(exp(x) - k*x - 1).

A(n,k) = Sum_{j=0..n} (-k)^j*binomial(n,j)*Bell(n-j).

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A005494 3-Bell numbers: E.g.f.: exp(3*z + exp(z) - 1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

SageMath

Formula

A045379 Expansion of e.g.f.: exp(4*z + exp(z) - 1).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

SageMath

Formula

A134980 a(n) = Sum_{k=0..n} binomial(n,k)*n^(n-k)*A000110(k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Sage

Formula

Extensions

A259691 Triangle read by rows: T(n,k) number of arrangements of non-attacking rooks on an n X n right triangular board where the top rook is in row k (n >= 0, 1 <= k <= n+1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A347420 Number of partitions of [n] where the first k elements are marked (0 <= k <= n) and at least k blocks contain their own index.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

A134980 a(n) = Sum_{k=0..n} binomial(n,k)n^(n-k)A000110(k).