OFFSET
0,2
FORMULA
T(n, k) = [x^k] Sum_{k=0..n} 4^(n-k)*binomial(n, k)^2*binomial(2*k, k)*k!*x^k, see A393915.
EXAMPLE
Triangle starts:
[0] 1;
[1] 4, 2;
[2] 16, 32, 12;
[3] 64, 288, 432, 120;
[4] 256, 2048, 6912, 7680, 1680;
[5] 1024, 12800, 76800, 192000, 168000, 30240;
[6] 4096, 73728, 691200, 3072000, 6048000, 4354560, 665280;
MAPLE
p := (n, x) -> 4^n * hypergeom([-n, -n, 1/2], [1], x):
T := (n, k) -> coeff(simplify(p(n, x)), x, k):
for n from 0 to 6 do seq(T(n, k), k = 0..n) od;
PROG
(SageMath)
def P(n):
x = SR.var('x')
poly = c = SR(1)
for k in range(n):
c = (c * (k - n)^2 * (2*k + 1) / 2) / (k + 1)^2
poly += c * x^(k + 1)
return 4^n * poly
for n in range(7): print(n, P(n).list())
(PARI) T(n, k) = 4^n * polcoef(hypergeom([-n, -n, 1/2], [1], x), k); \\ Michel Marcus, May 12 2026
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Peter Luschny, May 12 2026
STATUS
approved
