OFFSET
0,4
COMMENTS
Define an integer-valued polynomial p_n(x) = Product_{i=1..n} ((n+1)*x+i) * (2*x+1) / n!. Then the n-th row describes the numerator of the generating function Sum_{k >= 0} p_n(k) t^k as a rational function.
EXAMPLE
The first few rows are:
1, 1;
1, 6, 1;
1, 26, 26, 1;
1, 100, 310, 100, 1;
1, 372, 2752, 2752, 372, 1;
1, 1379, 21259, 48034, 21259, 1379, 1;
PROG
(SageMath)
x = polygen(QQ, 'x')
t = x.parent()[['t']].gen()
def row(n):
pn = prod((n+1)*x+i for i in range(1, n+1)) * (2*x+1) / factorial(n)
poly = (1 - t)**(n + 2) * sum(pn(k) * t**k for k in range(n + 2)).O(n + 2)
return poly.coefficients()
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
F. Chapoton, Apr 26 2026
STATUS
approved
