|
%I #40 Feb 26 2024 10:34:47
%S 1,-1,1,3,-4,1,-21,31,-11,1,315,-486,196,-26,1,-9765,15381,-6562,1002,
%T -57,1,615195,-978768,428787,-69688,4593,-120,1,-78129765,124918731,
%U -55434717,9279163,-652999,19833,-247,1
%N Triangle read by rows, matrix inverse of A139382.
%C This triangle appears to be the q-analog of A008275 (Stirling numbers of the first kind) for q=2. However, A333142 has a similar definition.
%C Row sums of unsigned triangle are A006125 with offset 1.
%H John Keith, <a href="/A342186/b342186.txt">Rows n = 1..20 of triangle, flattened</a>
%F T(n,k) = T(n-1,k-1) - (2^(n-1)-1) * T(n-1,k), n, k >= 1, T(1, 1) = 1, T(n, 0) = 0.
%F For unsigned triangle, T(n, 1) = A005329(n-1); T(n, 2) = A203011(n-1); T(n, n-1) = A000295(n+1); T(n, n-2) = A203242(n-1).
%F T(n,k) = Sum_{j=k..n} (-1)^(n-j)*2^binomial(n-j,2)*qBinomial(n,j,2)*binomial(j,k), where qBinomial(n,k,2) is A022166(n,k). - _Fabian Pereyra_, Feb 08 2024
%e The triangle begins:
%e 1;
%e -1, 1;
%e 3, -4, 1;
%e -21, 31, -11, 1;
%e 315, -486, 196, -26, 1;
%e -9765, 15381, -6562, 1002, -57, 1;
%e 615195, -978768, 428787, -69688, 4593, -120, 1;
%e -78129765, 124918731, -55434717, 9279163, -652999, 19833, -247, 1;
%e ...
%p A342186 := proc(n, k) if n = 1 and k = 1 then 1 elif k > n or k < 1 then 0 else
%p A342186(n-1, k-1) - (2^(n-1) - 1) * A342186(n-1, k) fi end:
%p for n from 1 to 8 do seq(A342186(n, k), k = 1..n) od; # _Peter Luschny_, Jun 28 2022
%t T[1, 1] := 1; T[n_, k_] := T[n, k] = If[k > n || k < 1, 0, T[n - 1, k - 1] - (2^(n - 1) - 1)*T[n - 1, k]]; Table[T[n, k], {n, 1, 8}, {k, 1, n}] (* after _G. C. Greubel_'s program for A139382 *)
%o (PARI) mat(nn) = my(m = matrix(nn, nn)); for (n=1, nn, for(k=1, nn, m[n,k] = if (n==1, if (k==1, 1, 0), if (k==1, 1, (2^k-1)*m[n-1, k] + m[n-1, k-1])););); m; \\ A139382
%o tabl(nn) = 1/mat(nn); \\ _Michel Marcus_, Mar 18 2021
%Y Cf. A008275, A139382, A333142, A333143, A006125 (row sums).
%Y Columns of unsigned triangle: A005329, A203011, A000295, A203242.
%K sign,tabl
%O 1,4
%A _John Keith_, Mar 04 2021
|
