|
%I #11 Sep 08 2022 08:45:52
%S 1,1,1,1,15,1,1,65,65,1,1,175,225,175,1,1,369,529,529,369,1,1,671,
%T 1025,1135,1025,671,1,1,1105,1761,2065,2065,1761,1105,1,1,1695,2785,
%U 3391,3585,3391,2785,1695,1,1,2465,4145,5185,5681,5681,5185,4145,2465,1
%N Triangle T(n,k) = 1 + 2*k*(n-k)*(k^2 -n*k +2*n^2) read by rows.
%C This could be written T(n,k) = 1-(n-k)^4 -k^4 +n^4, the quartic analog of A176284.
%C Row sums are {1, 2, 17, 132, 577, 1798, 4529, 9864, 19329, 34954, 59345, ...} = (n+1)*(9*n^4 -9*n^3 -n^2 +n +15)/15.
%H G. C. Greubel, <a href="/A176286/b176286.txt">Rows n = 0..100 of triangle, flattened</a>
%F T(n,k) = T(n,n-k).
%e Triangle begins as:
%e 1;
%e 1, 1;
%e 1, 15, 1;
%e 1, 65, 65, 1;
%e 1, 175, 225, 175, 1;
%e 1, 369, 529, 529, 369, 1;
%e 1, 671, 1025, 1135, 1025, 671, 1;
%e 1, 1105, 1761, 2065, 2065, 1761, 1105, 1;
%e 1, 1695, 2785, 3391, 3585, 3391, 2785, 1695, 1;
%e 1, 2465, 4145, 5185, 5681, 5681, 5185, 4145, 2465, 1;
%e 1, 3439, 5889, 7519, 8449, 8751, 8449, 7519, 5889, 3439, 1;
%p seq(seq(n^4 -(n-k)^4 -k^4 +1, k=0..n), n=0..12); # _G. C. Greubel_, Nov 25 2019
%t (* First program *)
%t f[n_, m_, q_]:= f[n, m, q] = 1 -(n-m)^q -m^q +n^q;
%t Table[Flatten[Table[Table[f[n, m, q], {m, 0, n}], {n, 0, 10}]], {q,1,10}]
%t (* Second program *)
%t Table[n^4 -(n-k)^4 -k^4 +1, {n,0,12}, {k,0,n}]//Flatten (* _G. C. Greubel_, Nov 25 2019 *)
%o (PARI) T(n,k) = n^4 -(n-k)^4 -k^4 +1; \\ _G. C. Greubel_, Nov 25 2019
%o (Magma) [n^4 -(n-k)^4 -k^4 +1: k in [0..n], n in [0..12]]; // _G. C. Greubel_, Nov 25 2019
%o (Sage) [[n^4 -(n-k)^4 -k^4 +1 for k in (0..n)] for n in (0..12)] # _G. C. Greubel_, Nov 25 2019
%o (GAP) Flat(List([0..12], n-> List([0..n], k-> n^4 -(n-k)^4 -k^4 +1 ))); # _G. C. Greubel_, Nov 25 2019
%K nonn,tabl,easy
%O 0,5
%A _Roger L. Bagula_, Apr 14 2010
%E Edited by _R. J. Mathar_, May 03 2013
|
