OFFSET
1,2
COMMENTS
LINKS
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
W. Lang, First 6 rows.
M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.
FORMULA
Recursion: a(n, k) = sum(binomial(2, p)*fallfac(2*(n-1)-p+k, 2-p)*a(n-1, k-p), p=0..2), n>=2, 2<=k<=2*n, a(1, 2)=1, else 0. Rewritten from eq.(19) of the Schork reference with r=4, s=2. fallfac(n, m) := A008279(n, m) (falling factorials triangle).
a(n, k) = (((-1)^k)/k!)*sum(((-1)^p)*binomial(k, p)*product(fallfac(p+2*(j-1), 2), j=1..n), p=2..k), n>=1, 2<=k<=2*n, else 0. From eq. (12) of the Blasiak et al. reference with r=4, s=2.
a(n, k) = ((2*n)!/k!)*binomial(2*(n-1), k-2), n>=1, 2<=k<=2*n.
E.g.f. column k>=2 (with leading zeros): (((-1)^k)/k!)*(sum(((-1)^p)*binomial(k, p)*hypergeom([(p-1)/2, p/2], [], 4*x), p=2..k)-(k-1)).
Coefficient triangle of the polynomials (2*n+2)!*hypergeom([-2*n],[3],-x)/2. - Peter Luschny, Apr 08 2015
Coefficient triangle of Laguerre polynomials (2*n)!*L(2*n,2,-x)). - Peter Luschny, Apr 08 2015
MAPLE
with(PolynomialTools):
p := n -> (2*n+2)!*hypergeom([-2*n], [3], -x)/2:
seq(CoefficientList(simplify(p(n)), x), n=0..5); # Peter Luschny, Apr 08 2015
MATHEMATICA
a[n_, k_] := (-1)^k/k!*Sum[(-1)^p*Binomial[k, p]*Product[FactorialPower[p + 2*(j-1), 2], {j, 1, n}], {p, 2, k}]; Table[a[n, k], {n, 1, 8}, {k, 2, 2 n}] // Flatten (* Jean-François Alcover, Sep 01 2016 *)
CROSSREFS
KEYWORD
nonn,easy,tabf
AUTHOR
Wolfdieter Lang, Dec 23 2003
STATUS
approved