login
Triangle read by rows: T(n, k) = t(n-k, k), where t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0, else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), and f(n) = 5*n + 2.
10

%I #14 Mar 01 2022 05:32:37

%S 1,2,2,4,28,4,8,244,244,8,16,1844,5856,1844,16,32,13260,101620,101620,

%T 13260,32,64,93684,1511160,3455080,1511160,93684,64,128,657836,

%U 20663388,91981880,91981880,20663388,657836,128,256,4609588,269011408,2121603436,4047202720,2121603436,269011408,4609588,256

%N Triangle read by rows: T(n, k) = t(n-k, k), where t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0, else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), and f(n) = 5*n + 2.

%H G. C. Greubel, <a href="/A257614/b257614.txt">Rows n = 0..50 of the triangle, flattened</a>

%F T(n, k) = t(n-k, k), where t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0, else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), and f(n) = 5*n + 2.

%F Sum_{k=0..n} T(n, k) = A008546(n).

%F From _G. C. Greubel_, Mar 01 2022: (Start)

%F t(k, n) = t(n, k).

%F T(n, n-k) = T(n, k).

%F t(0, n) = T(n, 0) = A000079(n). (End)

%e Array t(n,k) begins as:

%e 1, 2, 4, 8, 16, ... A000079;

%e 2, 28, 244, 1844, 13260, ...;

%e 4, 244, 5856, 101620, 1511160, ...;

%e 8, 1844, 101620, 3455080, 91981880, ...;

%e 16, 13260, 1511160, 91981880, 4047202720, ...;

%e 32, 93684, 20663388, 2121603436, 146321752612, ...;

%e 64, 657836, 269011408, 44675623468, 4648698508440, ...;

%e Triangle T(n,k) begins as:

%e 1;

%e 2, 2;

%e 4, 28, 4;

%e 8, 244, 244, 8;

%e 16, 1844, 5856, 1844, 16;

%e 32, 13260, 101620, 101620, 13260, 32;

%e 64, 93684, 1511160, 3455080, 1511160, 93684, 64;

%e 128, 657836, 20663388, 91981880, 91981880, 20663388, 657836, 128;

%t t[n_, k_, p_, q_]:= t[n, k, p, q] = If[n<0 || k<0, 0, If[n==0 && k==0, 1, (p*k+q)*t[n-1,k,p,q] + (p*n+q)*t[n,k-1,p,q]]];

%t T[n_, k_, p_, q_]= t[n-k, k, p, q];

%t Table[T[n,k,5,2], {n,0,12}, {k,0,n}]//Flatten (* _G. C. Greubel_, Mar 01 2022 *)

%o (Sage)

%o @CachedFunction

%o def t(n,k,p,q):

%o if (n<0 or k<0): return 0

%o elif (n==0 and k==0): return 1

%o else: return (p*k+q)*t(n-1,k,p,q) + (p*n+q)*t(n,k-1,p,q)

%o def A257614(n,k): return t(n-k,k,5,2)

%o flatten([[A257614(n,k) for k in (0..n)] for n in (0..12)]) # _G. C. Greubel_, Mar 01 2022

%Y Cf. A000079, A008546 (row sums), A142460, A257623.

%Y Cf. A038208, A256890, A257609, A257610, A257612, A257616, A257617, A257618, A257619.

%Y Similar sequences listed in A256890.

%K nonn,tabl

%O 0,2

%A _Dale Gerdemann_, May 09 2015