A071379 -id:A071379 - cp's OEIS frontend

A090210 Triangle of certain generalized Bell numbers.

1, 1, 1, 2, 1, 1, 5, 7, 1, 1, 15, 87, 34, 1, 1, 52, 1657, 2971, 209, 1, 1, 203, 43833, 513559, 163121, 1546, 1, 1, 877, 1515903, 149670844, 326922081, 12962661, 13327, 1, 1, 4140, 65766991, 66653198353, 1346634725665, 363303011071, 1395857215, 130922, 1, 1

Offset: 1

Wolfdieter Lang, Dec 01 2003

Let B_{n}(x) = sum_{j>=0}(exp(j!/(j-n)!*x-1)/j!) and
S(n,k) = k! [x^k] taylor(B_{n}(x)), where [x^k] denotes the
coefficient of x^k in the Taylor series for B_{n}(x).
Then S(n,k) (n>0, k>=0) is the square array representation of the triangle.
To illustrate the cross-references of T(n,k) when written as a square array.
n\k         A000012,A000012,A002720,A069948,A182925,A182924,...
0: A000012: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
1: A000110: 1, 1, 2, 5, 15, 52, 203, 877, 4140, ...
2: A020556: 1, 1, 7, 87, 1657, 43833, 1515903, ...
3: A069223: 1, 1, 34, 2971, 513559, 149670844, ...
4: A071379: 1, 1, 209, 163121, 326922081, ...
5: A090209: 1, 1, 1546, 12962661, 363303011071,...
6:  ...     1, 1, 13327, 1395857215, 637056434385865,...
Note that the sequence T(0,k) is not included in the data.
- Peter Luschny, Mar 27 2011

			Triangle begins:
1;
1, 1;
2, 1, 1;
5, 7, 1, 1;
15, 87, 34, 1, 1;
52, 1657, 2971, 209, 1, 1;
203, 43833, 513559, 163121, 1546, 1, 1;

P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.

P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem.
W. Lang, First 8 rows.
K. A. Penson, P. Blasiak, A. Horzela, A. I. Solomon and G. H. E. Duchamp,Laguerre-type derivatives: Dobinski relations and combinatorial identities, J. Math. Phys. 50, 083512 (2009).

Cf. A000110, A020556, A069223, A071379, A090209, A002720, A069948, A182924, A182925, A182933.

T(n,n) gives A070227.

A090210_AsSquareArray := proc(n,k) local r,s,i;
if k=0 then 1 else r := [seq(n+1,i=1..k-1)]; s := [seq(1,i=1..k-1)];
exp(-x)*n!^(k-1)*hypergeom(r,s,x); round(evalf(subs(x=1,%),99)) fi end:
seq(lprint(seq(A090210_AsSquareArray(n,k),k=0..6)),n=0..6);
# Peter Luschny, Mar 30 2011

t[n_, k_] := t[n, k] = Sum[(n+j)!^(k-1)/(j!^k*E), {j, 0, Infinity}]; t[_, 0] = 1;
Flatten[ Table[ t[n-k+1, k], {n, 0, 8}, {k, n, 0, -1}]][[1 ;; 43]] (* Jean-François Alcover, Jun 17 2011 *)

a(n, m) = Bell(m;n-(m-1)), n>= m-1 >=0, with Bell(m;k) := Sum_{p=m..m*k} S2(m;k, p), where S2(m;k, p) := (((-1)^p)/p!) * Sum_{r=m..p} ((-1)^r)*binomial(p, r)*fallfac(r, m)^k; with fallfac(n, m) := A008279(n, m) (falling factorials) and m<=p<=k*m, k>=1, m=1, 2, ..., else 0. From eqs.(6) with r=s->m and eq.(19) with S_{r, r}(n, k)-> S2(r;n, k) of the Blasiak et al. reference. [Corrected by Sean A. Irvine, Jun 03 2024]
a(n, m) = (Sum_{k>=m} fallfac(k, m)^(n-(m-1)))/exp(1), n>=m-1>=0, else 0. From eq.(26) with r->m of the Schork reference which is rewritten eq.(11) of the original Blasiak et al. reference.
E.g.f. m-th column (no leading zeros): (Sum_{k>=m} exp(fallfac(k, m)*x)/k!) + A000522(m)/m!)/exp(1). Rewritten from the top of p. 4656 of the Schork reference.

A090214 Generalized Stirling2 array S_{4,4}(n,k).

Original entry on oeis.org

1, 24, 96, 72, 16, 1, 576, 13824, 50688, 59904, 30024, 7200, 856, 48, 1, 13824, 1714176, 21606912, 76317696, 110160576, 78451200, 30645504, 6976512, 953424, 78400, 3760, 96, 1, 331776, 207028224, 8190885888, 74684104704, 253100173824

Offset: 1

Wolfdieter Lang, Dec 01 2003

The row length sequence for this array is [1,5,9,13,17,...] = A016813(n-1), n >= 1.
The g.f. for the k-th column, (with leading zeros and k >= 4) is G(k,x) = x^ceiling(k/4)*P(k,x)/Product_{p = 4..k} (1 - fallfac(p,4)*x), with fallfac(n,m) := A008279(n,m) (falling factorials) and P(k,x) := Sum_{m = 0..kmax(k)} A090221(k,m)*x^m, k >= 4, with kmax(k) := A057353(k-4)= floor(3*(k-4)/4). For the recurrence of the G(k,x) see A090221.
Codara et al., show that T(n,k) gives the number of k-colorings of the graph nK_4 (the disjoint union of n copies of the complete graph K_4). - Peter Bala, Aug 15 2013

			Table begins
n\k|   4      5      6      7      8     9   10   11   12
= = = = = = = = = = = = = = = = = = = = = = = = = = = = =
1  |   1
2  |  24     96     72     16      1
3  | 576  13824  50688  59904  30024  7200  856   48    1
...

Robert Israel, Table of n, a(n) for n = 1..10011 (rows 1 to 71, flattened)
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
P. Codara, O. M. D’Antona, and P. Hell, A simple combinatorial interpretation of certain generalized Bell and Stirling numbers, arXiv:1308.1700v1 [cs.DM], 2013.
A. Dzhumadildaev and D. Yeliussizov, Path decompositions of digraphs and their applications to Weyl algebra, arXiv preprint arXiv:1408.6764v1 [math.CO], 2014. [Version 1 contained many references to the OEIS, which were removed in Version 2. - _N. J. A. Sloane_, Mar 28 2015]
Askar Dzhumadil’daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, 22(4) (2015), #P4.10.
Wolfdieter Lang, First 4 rows.
M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.

Cf. A090215, A071379 (row sums), A090213 (alternating row sums).

S_{1, 1} = A008277, S_{2, 1} = A008297 (ignoring signs), S_{3, 1} = A035342, S_{2, 2} = A078739, S_{3, 2} = A078740, S_{3, 3} = A078741.

T:= (n,k) -> (-1)^k/k!*add((-1)^p*binomial(k,p)*(p*(p-1)*(p-2)*(p-3))^n,p=4..k):
seq(seq(T(n,k),k=4..4*n),n=1..10); # Robert Israel, Jan 28 2016

a[n_, k_] := (((-1)^k)/k!)*Sum[((-1)^p)*Binomial[k, p]*FactorialPower[p, 4]^n, {p, 4, k}]; Table[a[n, k], {n, 1, 5}, {k, 4, 4*n}] // Flatten (* Jean-François Alcover, Sep 05 2012, updated Jan 28 2016 *)

a(n, k) = (-1)^k/k! * Sum_{p = 4..k} (-1)^p * binomial(k, p) * fallfac(p, 4)^n,  with fallfac(p, 4) := A008279(p, 4) = p*(p - 1)*(p - 2)*(p - 3); 4 <= k <= 4*n, n >= 1, else 0. From eq.(19) with r = 4 of the Blasiak et al. reference.
E^n = Sum_{k = 4..4*n} a(n,k)*x^k*D^k where D is the operator d/dx, and E the operator (x^4)*d^4/dx^4.
The row polynomials R(n,x) are given by the Dobinski-type formula R(n,x) = exp(-x)*Sum_{k >= 0} (k*(k - 1)*(k - 2)*(k - 3))^n*x^k/k!. - Peter Bala, Aug 15 2013

A090209 Generalized Bell numbers (from (5,5)-Stirling2 array A090216).

Original entry on oeis.org

1, 1, 1546, 12962661, 363303011071, 25571928251231076, 3789505947767235111051, 1049433111253356296672432821, 498382374325731085522315594481036, 380385281554629647028734545622539438171, 443499171330317702437047276255605780991365151

Offset: 0

Wolfdieter Lang, Dec 01 2003

Contribution from Peter Luschny, Mar 27 2011: (Start) Let B_{m}(x) = sum_{j>=0}(exp(j!/(j-m)!*x-1)/j!) then a(n) = n! [x^n] taylor(B_{5}(x)), where [x^n] denotes the coefficient of x^k in the Taylor series for B_{5}(x).

a(n) is row 5 of the square array representation of A090210. (End)

M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.

G. C. Greubel, Table of n, a(n) for n = 0..115
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem., arXiv:quant-ph/0402027, Phys. Lett. A 309 (3-4) (2003) 198-205
K. A. Penson, P. Blasiak, A. Horzela, A. I. Solomon and G. H. E. Duchamp, Laguerre-type derivatives: Dobinski relations and combinatorial identities, J. Math. Phys. 50, 083512 (2009).

Cf. (Generalized) Bell numbers from (m,m)-Stirling2 array: A000110 (m=1), A020556 (m=2), A069223 (m=3), A071379 (m=4). Triangle A090210.

A071379 := proc(n) local r,s,i;
if n=0 then 1 else r := [seq(6,i=1..n-1)]; s := [seq(1,i=1..n-1)];
exp(-x)*5!^(n-1)*hypergeom(r,s,x); round(evalf(subs(x=1,%),99)) fi end:
seq(A071379(n),n=0..8); # Peter Luschny, Mar 30 2011

fallfac[n_, k_] := Pochhammer[n-k+1, k]; a[n_, k_] := (((-1)^k)/k!)*Sum[((-1)^p)*Binomial[k, p]*fallfac[p, 5]^n, {p, 5, k}]; a[0] = 1; a[n_] := Sum[a[n, k], {k, 5, 5*n}];  Table[a[n], {n, 0, 8}] (* Jean-François Alcover, Mar 05 2014 *)

a(n) = Sum_{k=5..5*n} A090216(n, k), n>=1. a(0) := 1.
a(n) = Sum_{k >=5} (fallfac(k, 5)^n)/k!/exp(1), n>=1, a(0) := 1. From eq.(26) with r=5 of the Schork reference.
E.g.f. with a(0) := 1: (sum((exp(fallfac(k, 5)*x))/k!, k=5..infinity)+ A000522(4)/4!)/exp(1). From the top of p. 4656 with r=5 of the Schork reference.

If it is proved that A283154 and A090209 are the same, then the entries should be merged and A283154 recycled. - N. J. A. Sloane, Mar 06 2017

A182933 Generalized Bell numbers based on the rising factorial powers; square array read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 5, 27, 13, 1, 1, 15, 409, 778, 73, 1, 1, 52, 9089, 104149, 37553, 501, 1, 1, 203, 272947, 25053583, 57184313, 2688546, 4051, 1, 1, 877, 10515147, 9566642254, 192052025697, 56410245661, 265141267, 37633, 1

Offset: 0

Peter Luschny, Mar 29 2011

These numbers are related to the generalized Bell numbers based on the falling factorial powers (A090210).
The square array starts for n>=0, k>=0:
n\k=0,1,..  A000012,A000262,A182934,...
A000012: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
A000110: 1, 1, 2, 5, 15, 52, 203, 877, 4140, ...
A094577: 1, 3, 27, 409, 9089, 272947, 10515147, ...
A182932: 1, 13, 778, 104149, 25053583, 9566642254, ...
       : 1, 73, 37553, 57184313, 192052025697, ...
       : 1, 501, 2688546, 56410245661, ...
 ....  : 1, 4051, 265141267, 89501806774945, ...

Cf. A000110, A020556, A069223, A071379, A090209, A002720, A069948, A182925, A182924, A182933.

A182933_AsSquareArray := proc(n,k) local r,s,i;
r := [seq(n+1,i=1..k)]; s := [seq(1,i=1..k-1),2];
exp(-x)*n!^k*hypergeom(r,s,x); round(evalf(subs(x=1,%),99)) end:
seq(lprint(seq(A182933_AsSquareArray(n,k),k=0..6)),n=0..6);

a[n_, k_] := Exp[-1]*n!^k*HypergeometricPFQ[ Table[n+1, {k}], Append[ Table[1, {k-1}], 2], 1.]; Table[ a[n-k, k] // Round , {n, 0, 8}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Jul 29 2013 *)

Let r_k = [n+1,...,n+1] (k occurrences of n+1), s_k = [1,...,1,2] (k-1 occurrences of 1) and F_k the generalized hypergeometric function of type k_F_k, then a(n,k) = exp(-1)*n!^k*F_k(r_k, s_k | 1).

Let B_{n}(x) = sum_{j>=0}(exp((j+n-1)!/(j-1)!*x-1)/j!) then a(n,k) = k! [x^k] series(B_{n}(x)), where [x^k] denotes the coefficient of x^k in the Taylor series for B_{n}(x).

A283153 Number of set partitions of unique elements from an n X 4 matrix where elements from the same row may not be in the same partition.

Original entry on oeis.org

1, 209, 163121, 326922081, 1346634725665, 9939316337679281, 119802044788535500753, 2205421644124274191535553, 58945667435045762187763602753, 2198513228897522394476415669503377, 110833342180980170285766876408530089329, 7356710448423295420590529054176924329802529, 628972339934967292421997567343442748145219556449

Offset: 1

Marko Riedel, Mar 01 2017

nonn

Apparently a duplicate of A071379? - R. J. Mathar, Mar 06 2017

Cf. A283154, A283155, A071379.

Table[(4!^n) * Sum[Binomial[p,4]^n/p! * Sum[(-1)^k/k!,{k,0,4n-p}],{p,1,4n}],{n,1,50}] (* Indranil Ghosh, Mar 04 2017 *)

a(n) = (4!^n) * sum(p=1, 4*n, binomial(p,4)^n/p! * sum(k=0, 4*n-p, (-1)^k/k!)); \\ Indranil Ghosh, Mar 04 2017

a(n) = m!^n * Sum_{p=1..n*m} (binomial(p,m)^n/p!) * Sum_{k=0..n*m-p} (-1)^k/k! with m=4.

If it is proved that A283153 and A071379 are the same, then the entries should be merged and A283153 recycled. - N. J. A. Sloane, Mar 06 2017

A157280 a(n) arises in the normal ordering of n-th power of the operator (d/dx)(x(d/dx))^4.

Original entry on oeis.org

1, 52, 43833, 149670844, 1346634725665, 25571928251231076, 893591647147188285577, 52327970757667659912764908, 4796836032234830356783078467969

Offset: 1

Karol A. Penson, Feb 26 2009

nonn

this sequence generates the fifth terms of the following sequences:
a(2)=52=A000110(5), a(3)=43833=A020556(5), a(4)=149670844=A069223(5),
a(5)=1346634725665=A071379(5),a(6)=25571928251231076=A070227(5)

nMax = 8; kMax = 50; seq0 = {}; seq = {1}; While[seq != seq0, seq0 = seq; seq = (1/E Sum[HypergeometricPFQ[{k+1, k+1, k+1, k+1}, {1, 1, 1, 1}, x]/k!, {k, 0, kMax}] + O[x]^(nMax+1) // CoefficientList[#, x]&) Range[0, nMax]!^5 // Round; kMax += 10; Print[kMax]]; A157280 = seq (* Jean-François Alcover, Nov 07 2016 *)

Sequence defined through the following hypergeometric-type generating function, in Maple notation:
exp(-1)*sum(hypergeom([k+1,k+1,k+1,k+1],[1,1,1,1],x)/k!,k=0..infinity)=sum(a(n)*x^n/(n!)^5,n=0..infinity),
which is itself an infinite sum of hypergeometric functions.

A090213 Alternating row sums of array A090214 ((4,4)-Stirling2).

Original entry on oeis.org

1, -15, 1169, -154079, -148969375, 778633335441, -4003896394897551, 27901641934428560705, -268555885416357907647039, 3460225909437698652973995569, -56404253763542830420650221273263, 1050004356721541004548911018674177377

Offset: 1

Wolfdieter Lang, Dec 01 2003

P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.

P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem.

Cf. A000587, A090211-2. A071379 (non-alternating sum, generalized Bell-numbers).

a[n_] := Sum[(-1)^k FactorialPower[k, 4]^n/k!, {k, 2, Infinity}]*E; Array[a, 12] (* Jean-François Alcover, Sep 01 2016 *)

a(n) := sum( A090214(n, k)*(-1)^k, k=4..4*n), n>=1. a(0) := 1 may be added.
a(n) = sum(((-1)^k)*(fallfac(k, 4)^n)/k!, k=4..infinity)*exp(1), with fallfac(k, 4)=A008279(k, 4)=k*(k-1)*(k-2)*(k-3) and n>=1. This produces also a(0)=1.
E.g.f. if a(0)=1 is added: exp(1)*(sum(((-1)^k)*exp(fallfac(k, 4)*x)/k!, k=4..infinity) + A000166(3)/3!). with the subfactorials A000166. A000166(3)/3!=1/3. Similar to derivation on top of p. 4656 of the Schork reference.

cp's OEIS Frontend

A090210 Triangle of certain generalized Bell numbers.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A090214 Generalized Stirling2 array S_{4,4}(n,k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A090209 Generalized Bell numbers (from (5,5)-Stirling2 array A090216).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A182933 Generalized Bell numbers based on the rising factorial powers; square array read by antidiagonals.

Views

Author

Keywords

Comments

Crossrefs

Programs

Maple

Mathematica

Formula

A283153 Number of set partitions of unique elements from an n X 4 matrix where elements from the same row may not be in the same partition.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A157280 a(n) arises in the normal ordering of n-th power of the operator (d/dx)(x(d/dx))^4.

Views

Author

Keywords

Comments

Programs

Mathematica

Formula

A090213 Alternating row sums of array A090214 ((4,4)-Stirling2).

Views

Author

Keywords

References

Links

Crossrefs

Programs

Mathematica

Formula