OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Rows n = 0..100 of triangle, flattened
FORMULA
Triangle read by rows, T(n,k) = Sum_{j=k..n} j!, 0 <= k <= n.
(The above is to be read as matrix product A*D*A where A is lower triangular filled with 1's and D = diag(n!, n >= 0). - M. F. Hasler, Aug 26 2020)
T(n, k) = t(k) - t(n + 1), where t(n) = (-1)^(n + 1)*Gamma(n + 1)*Subfactorial(-(n + 1)). - Peter Luschny, Jul 11 2024
EXAMPLE
First few rows of the triangle are:
1;
2, 1;
4, 3, 2;
10, 9, 8, 6;
34, 33, 32, 30, 24;
...
T(4,2) = 32 = 4! + 3! + 2! = (24 + 6 + 2).
MAPLE
a:=proc(n, k) local j; add(factorial(j), j=k..n) end: seq(seq(a(n, k), k=0..n), n=0..8); # Muniru A Asiru, Oct 16 2018
MATHEMATICA
Table[Sum[j!, {j, k, n}], {n, 0, 15}, {k, 0, n}]//Flatten (* G. C. Greubel, Oct 15 2018 *)
PROG
(PARI) for(n=0, 15, for(k=0, n, print1(sum(j=k, n, j!), ", "))) \\ G. C. Greubel, Oct 15 2018
(PARI) (A143122(n, k)=sum(j=k, n, j!)); T(n)=matrix(n++, n, i, j, i>=j)*matrix(n, n, i, j, (i>=j)*i--!) \\ M. F. Hasler, Aug 26 2020
(Magma) [[(&+[Factorial(j): j in [k..n]]): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Oct 15 2018
(GAP) Flat(List([0..8], n->List([0..n], k->Sum([k..n], j->Factorial(j))))); # Muniru A Asiru, Oct 16 2018
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Gary W. Adamson & Roger L. Bagula, Jul 26 2008
STATUS
approved