A157638 Triangle of the elementwise product of binomial coefficients with q-binomial coefficients [n,k] for q = 2.
1, 1, 1, 1, 6, 1, 1, 21, 21, 1, 1, 60, 210, 60, 1, 1, 155, 1550, 1550, 155, 1, 1, 378, 9765, 27900, 9765, 378, 1, 1, 889, 56007, 413385, 413385, 56007, 889, 1, 1, 2040, 302260, 5440680, 14055090, 5440680, 302260, 2040, 1, 1, 4599, 1563660, 66194940
Offset: 0
Examples
Triangle begins: 1; 1, 1; 1, 6, 1; 1, 21, 21, 1; 1, 60, 210, 60, 1; 1, 155, 1550, 1550, 155, 1; 1, 378, 9765, 27900, 9765, 378, 1; 1, 889, 56007, 413385, 413385, 56007, 889, 1; 1, 2040, 302260, 5440680, 14055090, 5440680, 302260, 2040, 1; 1, 4599, 1563660, 66194940, 417028122, 417028122, 66194940, 1563660, 4599, 1;
Links
- Andrew Howroyd, Rows n=0..49 of triangle, flattened
Programs
-
Magma
q:=2; [[k le 0 select 1 else Binomial(n,k)*(&*[(1-q^(n-j))/(1-q^(j+1)): j in [0..(k-1)]]): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Nov 17 2018
-
Mathematica
t[n_, m_] = Product[Sum[k*(m + 1)^i, {i, 0, k - 1}], {k, 1, n}]; b[n_, k_, m_] = t[n, m]/(t[k, m]*t[n - k, m]); Flatten[Table[Table[b[n, k, 1], {k, 0, n}], {n, 0, 10}]]
-
PARI
T(n,k) = {binomial(n,k)*polcoef(x^k/prod(j=0, k, 1-2^j*x+x*O(x^n)), n)} \\ Andrew Howroyd, Nov 14 2018
-
PARI
q=2; for(n=0,10, for(k=0,n, print1(binomial(n,k)*prod(j=0,k-1, (1-q^(n-j))/(1-q^(j+1))), ", "))) \\ G. C. Greubel, Nov 17 2018
-
Sage
[[ binomial(n,k)*gaussian_binomial(n,k).subs(q=2) for k in range(n+1)] for n in range(10)] # G. C. Greubel, Nov 17 2018
Formula
T(n,k) = t(n)/(t(k)*t(n-k)) where t(n) = Product_{k=1..n} Sum_{i=0..k-1} k*2^i.
T(n,k) = binomial(n,k)*A022166(n,k) for 0 <= k <= n. - Werner Schulte, Nov 14 2018
Extensions
Edited and simpler name by Werner Schulte and Andrew Howroyd, Nov 14 2018
Comments