reviewed
approved
reviewed
approved
editing
proposed
editing
proposed
# Assuming offset 2:
# Alternatively, as coefficients using the generating function of the row polynomials:
mpa rgf := (n, x) -> ((sqrt(x) - 1)^(2*n)*(2*n*sqrt(x) + x + 1) - (sqrt(x) + 1)^(2*n)*(-2*n*sqrt(x) + x + 1))/(16*sqrt(x)):
T := (n, k) -> coeff(expand(mpargf(n, x)), x, k):
proposed
editing
editing
proposed
editing
proposed
proposed
reviewed
proposed
reviewed
T(n, k) = binomial(n, k)*hypergeom([3, k - n], [k + 1], -1). - Peter Luschny, Sep 23 2024
[0] 1
[1] 4, 1
[2] 13, 5, 1
{1}; {4,1}; {13,5,1}; {[3] 38, 18, 6, 1};...
[4] 104, 56, 24, 7, 1
[5] 272, 160, 80, 31, 8, 1
[6] 688, 432, 240, 111, 39, 9, 1
[7] 1696, 1120, 672, 351, 150, 48, 10, 1
Fourth row polynomial (n = 3): p(3, x) = 38 + 18*x + 6*x^2 + x^3.
T := (n, k) -> binomial(n, k)*hypergeom([3, k - n], [k + 1], -1):
for n from 0 to 7 do seq(simplify(T(n, k)), k = 0..n) od; # Peter Luschny, Sep 23 2024
approved
editing