OFFSET
0,5
LINKS
G. C. Greubel, Rows n = 0..50 of the triangle, flattened
FORMULA
T(n, k) = c(n, q)/c(k, q) if k <= floor(n/2), otherwise c(n, q)/c(n-k, q), where c(n, q) = Product_{j=1..n} (1 - q^j) and q = 2.
T(n, n-k) = T(n, k).
EXAMPLE
Triangle begins as:
1;
1, 1;
1, -3, 1;
1, -7, -7, 1;
1, -15, 105, -15, 1;
1, -31, 465, 465, -31, 1;
1, -63, 1953, -29295, 1953, -63, 1;
1, -127, 8001, -248031, -248031, 8001, -127, 1;
1, -255, 32385, -2040255, 63247905, -2040255, 32385, -255, 1;
1, -511, 130305, -16548735, 1042570305, 1042570305, -16548735, 130305, -511, 1;
MATHEMATICA
c[n_, q_]= Product[1-q^i, {i, n}];
T[n_, k_, q_]= If[Floor[n/2]>=k, c[n, q]/c[n-k, q], c[n, q]/c[k, q]];
Table[T[n, k, 2], {n, 0, 12}, {k, 0, n}]//Flatten
(* Second program *)
T[n_, k_, q_]= With[{QP=QPochhammer}, If[k<=Floor[n/2], QP[q, q, n]/QP[q, q, n-k], QP[q, q, n]/QP[q, q, k]]];
Table[T[n, k, 2], {n, 0, 15}, {k, 0, n}]//Flatten (* G. C. Greubel, Dec 04 2022 *)
PROG
(Magma)
QPochhammer:= func< n, a, q | n eq 0 select 1 else (&*[1-a*q^j: j in [0..n-1]]) >;
T:= func< n, k, q | k le Floor(n/2) select QPochhammer(n, q, q)/QPochhammer(n-k, q, q) else QPochhammer(n, q, q)/QPochhammer(k, q, q) >;
A174387:= func< n, k | T(n, k, 2) >;
[A174387(n, k): k in [0..n], n in [0..15]]; // G. C. Greubel, Dec 04 2022
(SageMath)
from sage.combinat.q_analogues import q_pochhammer
def T(n, k, q):
if ((n//2)>k-1): return q_pochhammer(n, q, q)/q_pochhammer(n-k, q, q)
else: return q_pochhammer(n, q, q)/q_pochhammer(k, q, q)
def A174387(n, k): return T(n, k, 2)
flatten([[A174387(n, k) for k in range(n+1)] for n in range(16)]) # G. C. Greubel, Dec 04 2022
CROSSREFS
KEYWORD
sign,tabl
AUTHOR
Roger L. Bagula, Mar 18 2010
EXTENSIONS
Edited by G. C. Greubel, Dec 04 2022
STATUS
approved