A072678 Generalized Bell numbers B_{4,2}.

1, 21, 1045, 93289, 12975561, 2581284541, 693347907421, 241253367679185, 105394372192969489, 56410454014314490981, 36271084122927079387941, 27567930377271475039277881, 24435533594428382909107147225

Offset

1,2

Links

Robert Israel, Table of n, a(n) for n = 1..220

P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem., arXiv:quant-phys/0402027, 2004.

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.

Formula

a(n) = (2*n)!*hypergeom([2*n+1], [3], 1)/(2*exp(1)), n=1, 2, ... Special values of the confluent hypergeometric function 1F1.

a(n) = sum(A090438(n, k), k=2..2*n)= sum((1/k!)*product(fallfac(k+(j-1)*(4-2), 2), j=1..n), k=1..infinity)/exp(1), n>=1. From eq.(9) of the Blasiak et al. reference with r=4, s=2. fallfac(n, m) := A008279(n, m) (falling factorials triangle). a(0) := 1 may be added.

8*n*(2*n-1)*(2*n+1)*(n+1)^2*(n+3)*(n+2)*a(n)+(2*(n+1))*(8*n^3+32*n^2+42*n+13)*a(n+1)*(n+3)*(n+2)-(8*n^2+38*n+51)*(n+3)*(n+2)*a(n+2)+(n+3)*(n+2)*a(n+3) = 0. - Robert Israel, May 23 2016

a(n) = A052852(2*n-1). - Mark van Hoeij, Sep 05 2022

Maple

f:= n -> simplify((2*n)!*hypergeom([2*n+1], [3], 1)/(2*exp(1))):
map(f, [$1..30]); # Robert Israel, May 23 2016

Mathematica

a[n_] := n*(2n-1)!*Hypergeometric1F1[2-2n, 3, -1]; Array[a, 30] (* Jean-François Alcover, Sep 01 2016 *)

Crossrefs

Cf. A090439 (alternating row sums of A090438).

Sequence in context: A186392 A360505 A263994 * A218821 A012153 A231521

Adjacent sequences: A072675 A072676 A072677 * A072679 A072680 A072681

Keyword

nonn

Author

Karol A. Penson, Jul 01 2002

Extensions

Edited by Wolfdieter Lang, Dec 23 2003

Status

approved