%I #41 Sep 08 2022 08:45:16
%S 1,1,1,1,3,1,1,6,5,1,1,9,15,7,1,1,12,33,28,9,1,1,15,60,81,45,11,1,1,
%T 18,96,189,161,66,13,1,1,21,141,378,459,281,91,15,1,1,24,195,675,1107,
%U 946,449,120,17,1,1,27,258,1107,2349,2673,1742,673,153,19,1
%N Triangle, read by rows, where the terms are generated by the rule: T(n,k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k-1) + T(n-3,k-1), with T(0,0)=1.
%C Row sums form A077939. This sequence was inspired by Luke Hanna.
%C Diagonal sums are A000078(n+3). - _Philippe Deléham_, Feb 16 2014
%C Riordan array (1/(1-x), x*(1+x+x^2)/(1-x)). - _Philippe Deléham_, Feb 16 2014
%H G. C. Greubel, <a href="/A102036/b102036.txt">Rows n = 0..100 of triangle, flattened</a>
%H Kuhapatanakul, Kantaphon; Anantakitpaisal, Pornpawee <a href="https://doi.org/10.1080/23311835.2017.1333293">The k-nacci triangle and applications</a>. Cogent Math. 4, Article ID 1333293, 13 p. (2017).
%H J. L. Ramírez, V. F. Sirvent, <a href="http://www.combinatorics.org/ojs/index.php/eljc/article/view/v22i1p38">A Generalization of the k-Bonacci Sequence from Riordan Arrays</a>, The Electronic Journal of Combinatorics, 22(1) (2015), #P1.38.
%F G.f.: 1/(1-y-x*(1+y+y^2)). - _Vladimir Kruchinin_, Apr 21 2015
%F T(n,k) = Sum_{m=0..(n-k)} Sum_{j=0..k}(C(j,m-j)*C(k,j))*C(n-m,k)). - _Vladimir Kruchinin_, Apr 21 2015
%F From _Werner Schulte_, Dec 07 2018: (Start)
%F G.f. of column k: Sum_{n>=0} T(n+k,k) * x^n = (1+x+x^2)^k / (1-x)^(k+1) = (1-x^3)^k / (1-x)^(2*k+1).
%F Let k >= 0 be some fixed integer and a_k(n) be multiplicative with a_k(p^e) = T(e+k,k) for prime p and e >= 0. Then we have the Dirichlet g.f.: Sum{n>0} a_k(n) / n^s = (zeta(s))^(2*k+1) / (zeta(3*s))^k. (End)
%e Generated by adding preceding terms in the triangle at positions that form the letter 'L':
%e T(n,k) =
%e T(n-3,k-1) +
%e T(n-2,k-1) +
%e T(n-1,k-1) + T(n-1,k).
%e Rows begin:
%e [1],
%e [1, 1],
%e [1, 3, 1],
%e [1, 6, 5, 1],
%e [1, 9, 15, 7, 1],
%e [1, 12, 33, 28, 9, 1],
%e [1, 15, 60, 81, 45, 11, 1],
%e [1, 18, 96, 189, 161, 66, 13, 1],
%e [1, 21, 141, 378, 459, 281, 91, 15, 1], ...
%p T:=(n,k)->add(add((binomial(j,m-j)*binomial(k,j))*binomial(n-m,k),j=0..k),m=0..n-k): seq(seq(T(n,k),k=0..n),n=0..10); # _Muniru A Asiru_, Dec 11 2018
%t T[n_, k_] := If[n < k || k < 0, 0, If[n == 0, 1, T[n - 1, k] + T[n - 1, k - 1] + T[n - 2, k - 1] + T[n - 3, k - 1]]]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* _Amiram Eldar_, Dec 07 2018 *)
%t Table[Sum[Binomial[n-m, k]*Sum[Binomial[j, m-j]*Binomial[k, j], {j, 0, k}], {m, 0, n-k}], {n, 0, 15}, {k, 0, n}]//Flatten (* _G. C. Greubel_, Dec 11 2018 *)
%o (PARI) {T(n,k)=if(n<k||k<0,0,if(n==0,1,T(n-1,k)+T(n-1,k-1)+T(n-2,k-1)+T(n-3,k-1)))}
%o (Maxima) T(n,k):=sum((sum(binomial(j,m-j)*binomial(k,j),j,0,k))*binomial(n-m,k),m,0,n-k); /* _Vladimir Kruchinin_, Apr 21 2015 */
%o (Magma) [[(&+[Binomial(n-m,k)*(&+[Binomial(j,m-j)*Binomial(k,j):j in [0..k]]): m in [0..n-k]]): k in [0..n]]: n in [0..15]]; // _G. C. Greubel_, Dec 11 2018
%o (Sage) [[sum(binomial(n-m,k)*sum(binomial(j,m-j)*binomial(k,j) for j in (0..k)) for m in (0..n-k)) for k in (0..n)] for n in range(15)] # _G. C. Greubel_, Dec 11 2018
%Y Cf. A077939, A103141.
%K nonn,tabl
%O 0,5
%A _Paul D. Hanna_, Dec 30 2004
|