A182930 -id:A182930 - cp's OEIS frontend

A020556 Number of oriented multigraphs on n labeled arcs (without loops).

1, 1, 7, 87, 1657, 43833, 1515903, 65766991, 3473600465, 218310229201, 16035686850327, 1356791248984295, 130660110400259849, 14177605780945123273, 1718558016836289502159, 230999008481288064430879, 34208659263890939390952225, 5549763869122023099520756513

Offset: 0

Gilbert Labelle (gilbert(AT)lacim.uqam.ca) and Simon Plouffe

nonn

Generalized Bell numbers: a(n) = Sum_{k=2..2*n} A078739(n,k), n >= 1.
Let B_{m}(x) = Sum_{j>=0} exp(j!/(j-m)!*x-1)/j! then
a(n) = n! [x^n] taylor(B_{2}(x)), where [x^n] denotes the coefficient of x^n in the Taylor series for B_{2}(x). a(n) is row 2 of the square array representation of A090210. - Peter Luschny, Mar 27 2011
Also the number of set partitions of {1,2,...,2n+1} such that the block |n+1| is a part but no block |m| with m < n+1. - Peter Luschny, Apr 03 2011

			Example: For n = 2 the a(2) = 7 are the number of set partitions of 5 such that the block |3| is a part but no block |m| with m < 3: 3|1245, 3|4|125, 3|5|124, 3|12|45, 3|14|25, 3|15|24, 3|4|5|12. - _Peter Luschny_, Apr 05 2011

G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004.

Alois P. Heinz, Table of n, a(n) for n = 0..288
P. Blasiak, K. A. Penson and A. I. Solomon, The Boson Normal Ordering Problem and Generalized Bell Numbers, arXiv:quant-ph/0212072, 2002.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
P. Codara, O. M. D'Antona, P. Hell, A simple combinatorial interpretation of certain generalized Bell and Stirling numbers, arXiv preprint arXiv:1308.1700 [cs.DM], 2013.
G. Labelle, Counting enriched multigraphs according to the number of their edges (or arcs), Discrete Math., 217 (2000), 237-248.
Peter Luschny, Set partitions
G. Paquin, Dénombrement de multigraphes enrichis, Mémoire, Math. Dept., Univ. Québec à Montréal, 2004. [Cached copy, with permission]
M. Riedel, Set partitions of unique elements from an n-by-m matrix where elements from the same row may not be in the same partition, Mathematics Stack Exchange.
M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.

Cf. A020554, A014500, A020558, A090210, A095149, A106436, A182930.

A020556 := proc(n) local k;
add((-1)^(n+k)*binomial(n,k)*combinat[bell](n+k),k=0..n) end:
seq(A020556(n),n=0..17); # Peter Luschny, Mar 27 2011
# Uses floating point arithmetic, increase working precision for large n.
A020556 := proc(n) local r,s,i;
if n=0 then 1 else r := [seq(3,i=1..n-1)]; s := [seq(1,i=1..n-1)];
exp(-x)*2^(n-1)*hypergeom(r,s,x); round(evalf(subs(x=1,%),99)) fi end:
seq(A020556(n),n=0..15); # Peter Luschny, Mar 30 2011
T := proc(n, k) option remember;
  if n = 1 then 1
elif n = k then T(n-1,1) - T(n-1,n-1)
else T(n-1,k) + T(n, k+1) fi end:
A020556 := n -> T(2*n+1,n+1);
seq(A020556(n), n = 0..99); # Peter Luschny, Apr 03 2011

f[n_] := f[n] = Sum[(k + 2)!^n/((k + 2)!*(k!^n)*E), {k, 0, Infinity}]; Table[ f[n], {n, 1, 16}]
(* Second program: *)
a[n_] := Sum[(-1)^k*Binomial[n, k]*BellB[2n-k], {k, 0, n}]; Table[a[n], {n, 0, 17}] (* Jean-François Alcover, Jul 11 2017, after Vladeta Jovovic *)

a(n)={my(bell=serlaplace(exp(exp(x + O(x^(2*n+1)))-1))); sum(k=0, n, (-1)^k*binomial(n,k)*polcoef(bell, 2*n-k))} \\ Andrew Howroyd, Jan 13 2020

a(n) = e*Sum_{k>=0} ((k+2)!^n/(k+2)!)*(k!^n), n>=1.
a(n) = (1/e)*Sum_{k>=2} (k*(k-1))^n/k!, n >= 1. a(0) := 1. (From eq.(26) with r=2 of the Schork reference.)
E.g.f.: (1/e)*(2 + Sum_{k>=2} ((exp(k*(k-1)*x))/k!)) (from top of p. 4656 of the Schork reference).
a(n) = Sum_{k=0..n} (-1)^k*binomial(n, k)*Bell(2*n-k). - Vladeta Jovovic, May 02 2004
a(n) = A095149(2n,n). - Alois P. Heinz, Dec 20 2018
a(n) = A106436(2n,n) = A182930(2n+1,n+1). - Alois P. Heinz, Jan 29 2019

Edited by Robert G. Wilson v, Apr 30 2002

A106436 Difference array of Bell numbers A000110 read by antidiagonals.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 1, 2, 3, 5, 4, 5, 7, 10, 15, 11, 15, 20, 27, 37, 52, 41, 52, 67, 87, 114, 151, 203, 162, 203, 255, 322, 409, 523, 674, 877, 715, 877, 1080, 1335, 1657, 2066, 2589, 3263, 4140, 3425, 4140, 5017, 6097, 7432, 9089, 11155, 13744, 17007, 21147

Offset: 0

Philippe Deléham, May 29 2005

Essentially Aitken's array A011971 with first column A000296.

Mirror image of A182930. - Alois P. Heinz, Jan 29 2019

			   1;
   0,  1;
   1,  1,  2;
   1,  2,  3,  5;
   4,  5,  7, 10, 15;
  11, 15, 20, 27, 37, 52;
  ...

Alois P. Heinz, Rows n = 0..140, flattened
Jocelyn Quaintance and Harris Kwong, A combinatorial interpretation of the Catalan and Bell number difference tables, Integers, 13 (2013), #A29.

Cf. A000110, A182930.
Diagonals give A005493, A011965-A011967, A191099, A000298, A011968-A011970.
T(2n,n) gives A020556.

b:= proc(n) option remember; `if`(n=0, 1, add(
      b(n-j)*binomial(n-1, j-1), j=1..n))
    end:
T:= proc(n, k) option remember; `if`(k=0, b(n),
      T(n+1, k-1)-T(n, k-1))
    end:
seq(seq(T(n, d-n), n=0..d), d=0..12);  # Alois P. Heinz, Jan 29 2019

bb = Array[BellB, m = 12, 0];
dd[n_] := Differences[bb, n];
A = Array[dd, m, 0];
Table[A[[n-k+1, k+1]], {n, 0, m-1}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 26 2019 *)
a[0,0]:=1; a[n_,0]:=a[n-1,n-1]-a[n-1,0]; a[n_,k_]/;0Oliver Seipel, Nov 23 2024 *)

Double-exponential generating function: sum_{n, k} a(n-k, k) x^n/n! y^k/k! = exp(exp{x+y}-1-x). a(n,k) = Sum_{i=k..n} (-1)^(n-i)*binomial(n-k,i-k)*Bell(i). - Vladeta Jovovic, Oct 14 2006

cp's OEIS Frontend

A020556 Number of oriented multigraphs on n labeled arcs (without loops).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions