|
%I #27 Sep 08 2022 08:45:30
%S 1,6,1,6,12,1,6,18,18,1,6,24,36,24,1,6,30,60,60,30,1,6,36,90,120,90,
%T 36,1,6,42,126,210,210,126,42,1,6,48,168,336,420,336,168,48,1,6,54,
%U 216,504,756,756,504,216,54,1
%N T(n,k) = 6*binomial(n,k) - 5*I(n,k), where I is the identity matrix; triangle T read by rows (n >= 0 and 0 <= k <= n).
%C Row sums give A048488.
%H G. C. Greubel, <a href="/A131114/b131114.txt">Rows n = 0..100 of triangle, flattened</a>
%F T(n,k) = 6*A007318(n,k) - 5*I(n,k), where A007318 = Pascal's triangle and I = Identity matrix.
%F Bivariate o.g.f.: Sum_{n,k>=0} T(n,k)*x^n*y^k = (1 + 5*x - x*y)/((1 - x*y)*(1 - x - x*y)).
%e Triangle T(n,k) (with rows n >= 0 and columns k = 0..n) begins:
%e 1;
%e 6, 1;
%e 6, 12, 1;
%e 6, 18, 18, 1;
%e 6, 24, 36, 24, 1;
%e 6, 30, 60, 60, 30, 1;
%e 6, 36, 90, 120, 90, 36, 1;
%e ...
%p seq(seq(`if`(k=n, 1, 6*binomial(n,k)), k=0..n), n=0..10); # _G. C. Greubel_, Nov 18 2019
%t Table[If[k==n, 1, 6*Binomial[n, k]], {n,0,10}, {k,0,n}]//Flatten (* _G. C. Greubel_, Nov 18 2019 *)
%o (PARI) T(n,k) = if(k==n, 1, 6*binomial(n,k)); \\ _G. C. Greubel_, Nov 18 2019
%o (Magma) [k eq n select 1 else 6*Binomial(n,k): k in [0..n], n in [0..10]]; // _G. C. Greubel_, Nov 18 2019
%o (Sage)
%o def T(n, k):
%o if (k==n): return 1
%o else: return 6*binomial(n, k)
%o [[T(n, k) for k in (0..n)] for n in (0..10)] # _G. C. Greubel_, Nov 18 2019
%o (GAP)
%o T:= function(n,k)
%o if k=n then return 1;
%o else return 6*Binomial(n,k);
%o fi; end;
%o Flat(List([0..10], n-> List([0..n], k-> T(n,k) ))); # _G. C. Greubel_, Nov 18 2019
%Y Cf. A007318, A048488, A131110, A131111, A131112, A131113, A131115.
%K nonn,tabl,easy,less
%O 0,2
%A _Gary W. Adamson_, Jun 15 2007
|
