STATUS
proposed
approved
proposed
approved
editing
proposed
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
approved
editing
editing
approved
A342186 := proc(n, k) if n = 1 and k = 1 then 1 elif k > n or k < 1 then 0 else
A342186(n-1, k-1) - (2^(n-1) - 1) * A342186(n-1, k) fi end:
for n from 1 to 8 do seq(A342186(n, k), k = 1..n) od; # Peter Luschny, Jun 28 2022
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
(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
approved
editing
editing
approved