A101845 - OEIS

%I #15 Nov 19 2024 11:47:19

%S 1,1,3,1,7,16,1,15,61,115,1,31,206,626,1056,1,63,659,2989,7554,11774,

%T 1,127,2052,13308,47349,105099,154624,1,255,6297,56935,274677,824331,

%U 1660957,2337507,1,511,19162,237862,1518478,5960818,15747154,29428654

%N Triangle formed by left half of A101842, read by rows.

%e Triangle begins:

%e   1,

%e   1,  3,

%e   1,  7,  16,

%e   1, 15,  61,  115,

%e   1, 31, 206,  626, 1056,

%e   1, 63, 659, 2989, 7554, 11774,

%e   ...

%p A101842 := proc(n,k) option remember ; if k < -n or k >= n then 0 ; elif n = 1 then 1; else (n-k)*A101842(n-1,k-1)+A101842(n-1,k)+(n+k+1)*A101842(n-1,k+1) ; fi ; end: A101845 := proc(n,k) A101842(n,-n+k-1) ; end: for n from 1 to 10 do for k from 1 to n do printf("%d, ",A101845(n,k)) ; od: od: # _R. J. Mathar_, Aug 07 2007

%t (* T is A101842 *)

%t T[n_, k_] /; -n <= k <= n-1 := T[n, k] = (n-k)*T[n-1, k-1]+T[n-1, k]+(n+k+1)* T[n-1, k+1];

%t T[1, -1] = T[1, 0] = 1; T[_, _] = 0;

%t A101845[n_, k_] := T[n, k-n-1];

%t Table[A101845[n, k], {n, 1, 10}, {k, 1, n}] // Flatten (* _Jean-François Alcover_, Nov 17 2024 *)

%Y Cf. A101842.

%K nonn,tabl

%O 1,3

%A _David Applegate_, Jun 19 2007

%E More terms from _R. J. Mathar_, Aug 07 2007