cp's OEIS Frontend

A141760 Triangle T, read by rows, where the g.f. of column k in matrix power T^m is given by: 1/(1-x)^m = Sum_{n>=k} [T^m](n,k) * x^(n-k)/(1+x)^{n(n-1)/2 - k(k-1)/2} for k>=0.

%I A141760 #14 Jun 14 2017 01:08:08
%S A141760 1,1,1,1,1,1,2,2,1,1,6,6,3,1,1,26,26,13,4,1,1,154,154,77,23,5,1,1,
%T A141760 1188,1188,594,175,36,6,1,1,11474,11474,5737,1678,336,52,7,1,1,134432,
%U A141760 134432,67216,19579,3863,576,71,8,1,1,1863168,1863168,931584,270683,52944,7731
%N A141760 Triangle T, read by rows, where the g.f. of column k in matrix power T^m is given by: 1/(1-x)^m = Sum_{n>=k} [T^m](n,k) * x^(n-k)/(1+x)^{n(n-1)/2 - k(k-1)/2} for k>=0.
%F A141760 Matrix powers satisfy: T^m = P(i)^-1 * P(m+i) for all m and i, where P(m) is given by:
%F A141760 [P(m)](n,k) = [x^(n-k)] 1/(1-x)^m * (1+x)^{n(n-1)/2 - k(k-1)/2} for n>=k>=0.
%F A141760 Let U = unsigned matrix inverse (T^-1) with leftmost column dropped, then U = A107876 where [U^k](n,k) = U(n,k-1) for n>=k>0.
%F A141760 G.f. for column k of T: 1/(1-x) = Sum_{n>=0} T(n,k)*x^(n-k)/(1+x)^{n(n-1)/2 - k(k-1)/2}.
%F A141760 T(n,k) = 1 - Sum_{j=k..n-1} T(j,k)*(-1)^(n-j)*C(j(j-1)/2 - k(k-1)/2 + n-j-1, n-j) for n>k with T(k,k)=1 for k>=0.
%F A141760 G.f. for column k of matrix power T^m:
%F A141760 1/(1-x)^m = Sum_{n>=0} [T^m](n,k)*x^(n-k)/(1+x)^{n*(n-1)/2 - k*(k-1)/2}.
%F A141760 [T^m](n,k) = C(m+n-1,n) - Sum_{j=k..n-1} [T^m](j,k)*(-1)^(n-j)*C(j(j-1)/2 - k(k-1)/2 + n-j-1,n-j) for n>k with [T^m](k,k)=1 for k>=0.
%e A141760 Triangle T begins:
%e A141760 1;
%e A141760 1, 1;
%e A141760 1, 1, 1;
%e A141760 2, 2, 1, 1;
%e A141760 6, 6, 3, 1, 1;
%e A141760 26, 26, 13, 4, 1, 1;
%e A141760 154, 154, 77, 23, 5, 1, 1;
%e A141760 1188, 1188, 594, 175, 36, 6, 1, 1;
%e A141760 11474, 11474, 5737, 1678, 336, 52, 7, 1, 1;
%e A141760 134432, 134432, 67216, 19579, 3863, 576, 71, 8, 1, 1;
%e A141760 1863168, 1863168, 931584, 270683, 52944, 7731, 911, 93, 9, 1, 1; ...
%e A141760 Matrix square, T^2, begins:
%e A141760 1;
%e A141760 2, 1;
%e A141760 3, 2, 1;
%e A141760 7, 5, 2, 1;
%e A141760 23, 17, 7, 2, 1;
%e A141760 105, 79, 33, 9, 2, 1;
%e A141760 641, 487, 205, 55, 11, 2, 1;
%e A141760 5034, 3846, 1626, 433, 83, 13, 2, 1; ...
%e A141760 where g.f. for column k of matrix square T^2 is:
%e A141760 1/(1-x)^2 = Sum_{n>=0} [T^2](n,k)*x^(n-k)/(1+x)^{n(n-1)/2 - k(k-1)/2}.
%e A141760 Matrix inverse, T^-1, begins:
%e A141760 1;
%e A141760 -1, 1;
%e A141760 0, -1, 1;
%e A141760 0, -1, -1, 1;
%e A141760 0, -2, -2, -1, 1;
%e A141760 0, -7, -7, -3, -1, 1;
%e A141760 0, -37, -37, -15, -4, -1, 1;
%e A141760 0, -268, -268, -106, -26, -5, -1, 1; ...
%e A141760 Let U = unsigned T^-1 with leftmost column dropped,
%e A141760 then U = A107876 where [U^k](n,k) = U(n,k-1) for n>=k>0.
%e A141760 The g.f. for column k of matrix inverse T^-1 is:
%e A141760 1-x = Sum_{n>=0} [T^-1](n,k) * x^(n-k)/(1+x)^{n(n-1)/2 - k(k-1)/2}.
%e A141760 MATRIX PRODUCTS:
%e A141760 T = P(1)^-1 * P(2) = P(2)^-1 * P(3) = P(m)^-1 * P(m+1);
%e A141760 P(1) begins:
%e A141760 1;
%e A141760 1, 1;
%e A141760 2, 2, 1;
%e A141760 8, 7, 3, 1;
%e A141760 57, 42, 16, 4, 1;
%e A141760 638, 386, 130, 29, 5, 1;
%e A141760 9949, 4944, 1471, 299, 46, 6, 1; ...
%e A141760 where [P(1)](n,k) = [x^(n-k)] 1/(1-x)*(1+x)^{n(n-1)/2-k(k-1)/2};
%e A141760 P(2) begins:
%e A141760 1;
%e A141760 2, 1;
%e A141760 5, 3, 1;
%e A141760 20, 12, 4, 1;
%e A141760 129, 72, 23, 5, 1;
%e A141760 1268, 630, 187, 38, 6, 1;
%e A141760 17548, 7599, 2063, 392, 57, 7, 1; ...
%e A141760 where [P(2)](n,k) = [x^(n-k)] 1/(1-x)^2*(1+x)^{n(n-1)/2-k(k-1)/2};
%e A141760 P(3) begins:
%e A141760 1;
%e A141760 3, 1;
%e A141760 9, 4, 1;
%e A141760 38, 18, 5, 1;
%e A141760 240, 111, 31, 6, 1;
%e A141760 2223, 955, 256, 48, 7, 1;
%e A141760 28672, 11124, 2794, 500, 69, 8, 1; ...
%e A141760 where [P(3)](n,k) = [x^(n-k)] 1/(1-x)^3*(1+x)^{n(n-1)/2-k(k-1)/2}.
%t A141760 T[n_, k_, m_] := T[n, k, m] = If[n<k || k<0, 0, If[n == k, 1, Binomial[m+n- 1, n] - Sum[T[j, k, m]*(-1)^(n-j)*Binomial[j*(j-1)/2 - k*(k-1)/2 + n-j-1, n-j], {j, k, n-1}]]]; Table[T[n, k, 1], {n, 0, 10}, {k, 0, n}] // Flatten (* _Jean-François Alcover_, Sep 19 2016, adapted from PARI *)
%o A141760 (PARI) T(n,k,m=1)=if(n<k || k<0,0,if(n==k,1, binomial(m+n-1,n) - sum(j=k,n-1,T(j,k,m)*(-1)^(n-j)*binomial(j*(j-1)/2-k*(k-1)/2+n-j-1,n-j))))
%Y A141760 Cf. columns: A141761, A141762, A141763; A107876 (unsigned inverse).
%K A141760 nonn,tabl
%O A141760 0,7
%A A141760 _Paul D. Hanna_, Jul 18 2008