login
Triangle T(n,k) read by rows with q-e.g.f.: 1/Product_{k>0} (1-x^k/faq(k,q)).
4

%I #28 Mar 11 2022 12:31:41

%S 1,1,2,1,3,3,3,1,5,7,11,11,8,4,1,7,13,25,36,44,42,36,24,13,5,1,11,24,

%T 54,93,142,184,215,222,208,172,126,81,44,19,6,1,15,39,98,195,344,532,

%U 753,964,1150,1264,1294,1226,1082,880,661,451,278,151,70,26,7,1

%N Triangle T(n,k) read by rows with q-e.g.f.: 1/Product_{k>0} (1-x^k/faq(k,q)).

%H Alois P. Heinz, <a href="/A152534/b152534.txt">Rows n = 0..50, flattened</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/q-ExponentialFunction.html">q-Exponential Function</a>.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/q-Factorial.html">q-Factorial</a>.

%F Sum_{k=0..binomial(n,2)} T(n,k)*q^k = Sum_{pi} faq(n,q)/Product_{i=1..n} faq(i,q)^e(i), where pi runs over all nonnegative integer solutions to e(1) + 2*e(2) + ... + n*e(n) = n and faq(i,q) = Product_{j=1..i} (q^j-1)/(q-1), i = 1..n. Sum_{k=0..binomial(n,2)} T(n,k)*exp(2*Pi*I*k/n)) = 1.

%F Sum_{k=0..binomial(n,2)} (-1)^k*T(n,k) = A152536(n). - _Alois P. Heinz_, Aug 09 2021

%e Triangle begins:

%e 1;

%e 1;

%e 2, 1;

%e 3, 3, 3, 1;

%e 5, 7, 11, 11, 8, 4, 1;

%e 7, 13, 25, 36, 44, 42, 36, 24, 13, 5, 1;

%e ...

%p multinomial2q := proc(n::integer,k::integer,nparts::integer)

%p local lpar ,res, constrp;

%p res := [] ;

%p if n< 0 or nparts <= 0 then

%p ;

%p elif nparts = 1 then

%p if n = k then

%p return [[n]] ;

%p end if;

%p else

%p for lpar from 0 do

%p if lpar*nparts > n or lpar > k then

%p break;

%p end if;

%p for constrp in procname(n-nparts*lpar,k-lpar,nparts-1) do

%p if nops(constrp) > 0 then

%p res := [op(res),[op(constrp),lpar]] ;

%p end if;

%p end do:

%p end do:

%p end if ;

%p return res ;

%p end proc:

%p multinomial2 := proc(n::integer,k::integer)

%p local res,constrp ;

%p res := [] ;

%p for constrp in multinomial2q(n,k,n) do

%p if nops(constrp) > 0 then

%p res := [op(res),constrp] ;

%p end if ;

%p end do:

%p res ;

%p end proc:

%p faq := proc(i,q)

%p mul((q^j-1)/(q-1),j=1..i) ;

%p end proc;

%p A152534 := proc(n,k)

%p pi := [] ;

%p for sp from 0 to n do

%p pi := [op(pi),op(multinomial2(n,sp))] ;

%p end do;

%p tqk := 0 ;

%p for p in pi do

%p faqe :=1 ;

%p for i from 1 to nops(p) do

%p faqe := faqe* faq(i,q)^op(i,p) ;

%p end do:

%p tqk := tqk+faq(n,q)/faqe ;

%p end do;

%p tqk ;

%p coeftayl(tqk,q=0,k) ;

%p end proc:

%p for n from 1 to 8 do

%p for k from 0 to binomial(n,2) do

%p printf("%d,",A152534(n,k)) ;

%p end do;

%p printf("\n") ;

%p end do: # _R. J. Mathar_, Sep 27 2011

%p # second Maple program:

%p f:= proc(n) option remember; `if`(n<2, 1, f(n-1)*(q^n-1)/(q-1)) end:

%p b:= proc(n, i) option remember; simplify(`if`(n=0 or i=1, 1,

%p add(b(n-i*j, i-1)/f(i)^j, j=0..n/i)))

%p end:

%p T:= n-> (p-> seq(coeff(p, q, i), i=0..degree(p)))(simplify(f(n)*b(n$2))):

%p seq(T(n), n=0..10); # _Alois P. Heinz_, Aug 09 2021

%t f[n_] := f[n] = If[n < 2, 1, f[n - 1]*(q^n - 1)/(q - 1)];

%t b[n_, i_] := b[n, i] = Simplify[If[n == 0 || i == 1, 1,

%t Sum[b[n - i*j, i - 1]/f[i]^j, {j, 0, n/i}]]];

%t T[n_] := CoefficientList[Simplify[f[n]*b[n, n]], q];

%t Table[T[n], {n, 0, 10}] // Flatten (* _Jean-François Alcover_, Mar 11 2022, after _Alois P. Heinz_ *)

%Y Cf. A005651 (row sums), A000041 (first column), A076276 (second column), A152474, A152536.

%Y T(n,n) gives A346980.

%K nonn,tabf

%O 0,3

%A _Vladeta Jovovic_, Dec 06 2008

%E T(0,0)=1 prepended by _Alois P. Heinz_, Aug 09 2021