login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A157272
Triangle T(n, k, m) = (m*(n-k) + 1)*T(n-1, k-1, m) + (m*k + 1)*T(n-1, k, m) + m*f(n,k)*T(n-2, k-1, m) with T(n, 0, m) = T(n, n, m) = 1, f(n, k) = 2*k+1 if k <= floor(n/2) otherwise 2*(n-k)+1, and m = 1, read by rows.
23
1, 1, 1, 1, 7, 1, 1, 20, 20, 1, 1, 47, 155, 47, 1, 1, 102, 753, 753, 102, 1, 1, 213, 3004, 7109, 3004, 213, 1, 1, 436, 10800, 48727, 48727, 10800, 436, 1, 1, 883, 36517, 280736, 551251, 280736, 36517, 883, 1, 1, 1778, 118795, 1454163, 4879214, 4879214, 1454163, 118795, 1778, 1
OFFSET
0,5
FORMULA
T(n, k, m) = (m*(n-k) + 1)*T(n-1, k-1, m) + (m*k + 1)*T(n-1, k, m) + m*f(n,k)*T(n-2, k-1, m) with T(n, 0, m) = T(n, n, m) = 1, f(n, k) = 2*k+1 if k <= floor(n/2) otherwise 2*(n-k)+1, and m = 1.
T(n, n-k, m) = T(n, k, m).
EXAMPLE
Triangle begins as:
1;
1, 1;
1, 7, 1;
1, 20, 20, 1;
1, 47, 155, 47, 1;
1, 102, 753, 753, 102, 1;
1, 213, 3004, 7109, 3004, 213, 1;
1, 436, 10800, 48727, 48727, 10800, 436, 1;
1, 883, 36517, 280736, 551251, 280736, 36517, 883, 1;
1, 1778, 118795, 1454163, 4879214, 4879214, 1454163, 118795, 1778, 1;
MATHEMATICA
f[n_, k_]:= If[k<=Floor[n/2], 2*k+1, 2*(n-k)+1];
T[n_, k_, m_]:= T[n, k, m]= If[k==0 || k==n, 1, (m*(n-k)+1)*T[n-1, k-1, m] + (m*k+1)*T[n-1, k, m] + m*f[n, k]*T[n-2, k-1, m]];
Table[T[n, k, 1], {n, 0, 12}, {k, 0, n}]//Flatten (* modified by G. C. Greubel, Feb 04 2022 *)
PROG
(Sage)
def f(n, k): return 2*k+1 if (k <= n//2) else 2*(n-k)+1
@CachedFunction
def T(n, k, m): # A157207
if (k==0 or k==n): return 1
else: return (m*(n-k) +1)*T(n-1, k-1, m) + (m*k+1)*T(n-1, k, m) + m*f(n, k)*T(n-2, k-1, m)
flatten([[T(n, k, 1) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Feb 04 2022
KEYWORD
nonn,tabl
AUTHOR
Roger L. Bagula, Feb 26 2009
EXTENSIONS
Edited by G. C. Greubel, Feb 04 2022
STATUS
approved