A262964 - OEIS

%I #21 Feb 21 2016 20:16:20

%S 1,1,0,1,0,1,1,0,3,0,1,0,6,1,2,1,0,10,4,9,0,1,0,15,10,28,7,5,1,0,21,

%T 20,69,36,30,1

%N Lower triangular matrix. See comments for definition.

%C Handmade table created by looking at the "columns" that are separated by zeros in the Mathematica program output. First column is A000012. Second column A000004. Third column A000217. Fourth column A000292. At least a few columns onwards are recognized by the OEIS automated sequence recognition as polynomials.

%C Row sums are A239605. Alternating row sums are A011782.

%e Triangle starts and row sums are A239605:

%e    1                                    =    1

%e    1    0                               =    1

%e    1    0    1                          =    2

%e    1    0    3    0                     =    4

%e    1    0    6    1    2                =   10

%e    1    0   10    4    9    0           =   24

%e    1    0   15   10   28    7    5      =   66

%e    1    0   21   20   69   36   30    1 =  178

%e Alternating row sums are A011782:

%e    1                                    =    1

%e    1   -0                               =    1

%e    1   -0    1                          =    2

%e    1   -0    3   -0                     =    4

%e    1   -0    6   -1    2                =    8

%e    1   -0   10   -4    9   -0           =   16

%e    1   -0   15  -10   28   -7    5      =   32

%e    1   -0   21  -20   69  -36   30   -1 =   64

%t (* To get this number triangle, look at the numbers separated by zeros in the columns *)

%t (* coefficients (coeff) in power series can be changed *)Clear[t, n, k, i, nn, x];

%t Clear[x]

%t coeff = {1, 1000000000000000000, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,

%t    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1};

%t mp[m_, e_] :=

%t If[e == 0, IdentityMatrix@Length@m, MatrixPower[m, e]]; nn =

%t Length[coeff]; cc = Range[nn]*0 + 1; Monitor[

%t Do[Clear[t]; t[n_, 1] := t[n, 1] = cc[[n]];

%t   t[n_, k_] :=

%t    t[n, k] =

%t     If[n >= k,

%t      Sum[t[n - i, k - 1], {i, 1, k - 1}] -

%t       Sum[t[n - i, k], {i, 1, k - 1}], 0];

%t   A4 = Table[Table[t[n, k], {k, 1, nn}], {n, 1, nn}];

%t   A5 = A4[[1 ;; nn - 1]]; A5 = Prepend[A5, ConstantArray[0, nn]];

%t   cc = Total[

%t     Table[coeff[[n]]*mp[A5, n - 1][[All, 1]], {n, 1, nn}]];, {i, 1,

%t    nn}], i];

%t cc;

%t TableForm[A4[[All,1]]]

%Y Cf. Row sums A239605, alternating row sums A011782, columns: A000012, A000004, A000217, A000292.

%K tabl,nonn

%O 1,9

%A _Mats Granvik_, Oct 05 2015