|
%I #17 Mar 09 2015 06:07:15
%S 1,1,1,1,2,1,2,4,3,1,6,12,9,4,1,24,48,36,16,5,1,120,240,180,80,25,6,1,
%T 720,1440,1080,480,150,36,7,1,5040,10080,7560,3360,1050,252,49,8,1,
%U 40320,80640,60480,26880,8400,2016,392,64,9,1
%N Triangle by rows relating to the factorials, generated from A002260.
%C Row sums = A054091: (1, 2, 4, 10, 32, 130, 652,...)
%C Left border = the factorials, A000142 prefaced with a 1.
%H Alois P. Heinz, <a href="/A208058/b208058.txt">Rows n = 0..140, flattened</a>
%F Inverse of:
%F 1;
%F -1, 1;
%F 1, -2, 1;
%F -1, 2, -3, 1;
%F 1, -2, 3, -4, 1;
%F ..., where triangle A002260 = (1; 1,2; 1,2,3;...)
%e First few rows of the triangle =
%e 1;
%e 1, 1;
%e 1, 2, 1;
%e 2, 4, 3, 1;
%e 6, 12, 9, 4, 1;
%e 24, 48, 36, 16, 5, 1;
%e 120, 240, 180, 80, 25, 6, 1;
%e 720, 1440, 1080, 480, 150, 36, 7, 1;
%e 5040, 10080, 7560, 3360, 1050, 252, 49, 8, 1;
%e ...
%p T:= proc(n) option remember; local M, k;
%p M:= Matrix(n+1, (i, j)->
%p `if`(i=j, 1, `if`(i>j, j*(-1)^(i+j), 0)))^(-1);
%p seq(M[n+1, k], k=1..n+1)
%p end:
%p seq(T(n), n=0..14); # _Alois P. Heinz_, Feb 24 2012
%t T[n_] := T[n] = Module[{M}, M = Table[If[i == j, 1, If[i>j, j*(-1)^(i+j), 0]], {i, 1, n+1}, {j, 1, n+1}] // Inverse; M[[n+1]]]; Table[T[n], {n, 0, 14}] // Flatten (* _Jean-François Alcover_, Mar 09 2015, after _Alois P. Heinz_ *)
%Y Cf. A000142, A002260, A054091.
%K nonn,tabl
%O 0,5
%A _Gary W. Adamson_, Feb 22 2012
|
