%I #43 May 24 2020 14:58:20
%S 1,1,2,1,2,3,1,3,4,5,1,3,5,6,7,1,4,7,9,10,11,1,4,8,10,12,13,14,1,5,9,
%T 14,16,18,19,20,1,5,12,17,21,23,25,26,27,1,6,13,21,26,30,32,34,35,36,
%U 1,6,16,25,33,37,41,43,45,46,47,1,7,19,32,42,50,54,58,60,62,63,64
%N Triangle T(n,k) of number of distinct values of multinomial coefficients corresponding to sequence A026820 (n >= 1, 1 <= k <= n).
%C Differs from A026820 after position 24.
%C Includes sequence A070289 when k = n.
%H Alois P. Heinz, <a href="/A213008/b213008.txt">Rows n = 1..45, flattened</a>
%H Katsuhisa Yamanaka, Shin-ichiro Kawano, Yosuke Kikuchi, and Shin-ichi Nakano, <a href="https://doi.org/10.1093/ietfec/e90-a.5.888">Constant Time Generation of Integer Partitions</a>, IEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences, vol. E90-A, no.5, pp. 888-895, (May-2007).
%H Sergei Viznyuk, <a href="http://phystech.com/ftp/A213008.c">C-Program for this sequence</a>, 2012.
%H Sergei Viznyuk, <a href="/A213008/a213008.c.txt">C-Program for sequences A026820, A070289, and A213008 (local copy)</a>, 2012.
%e Triangle T(n,k) begins:
%e 1;
%e 1, 2;
%e 1, 2, 3;
%e 1, 3, 4, 5;
%e 1, 3, 5, 6, 7;
%e 1, 4, 7, 9, 10, 11;
%e 1, 4, 8, 10, 12, 13, 14;
%e ...
%e Thus, for n = 7 and k = 6 there are 13 distinct values of multinomial coefficients corresponding to partitions of n = 7 into at most k = 6 parts. The corresponding number of partitions from sequence A026820 is 14. That is because partitions 7 = 4 + 1 + 1 + 1 and 7 = 3 + 2 + 2 produce the same value of multinomial coefficient 7!/(4!*1!*1!*1!) = 7!/(3!*2!*2!).
%p b:= proc(n, i, k) option remember; if n=0 then {1} elif i<1
%p then {} else {b(n, i-1, k)[], seq(map(x-> x*i!^j,
%p b(n-i*j, i-1, k-j))[], j=1..min(n/i, k))} fi
%p end:
%p T:= (n, k)-> nops(b(n, n, k)):
%p seq(seq(T(n,k), k=1..n), n=1..14); # _Alois P. Heinz_, Aug 14 2012
%t b[n_, i_, k_] := b[n, i, k] = If[n == 0, {1}, If[i<1, {}, Join[b[n, i-1, k], Table[ Function[#*i!^j] /@ b[n-i*j, i-1, k-j], {j, 1, Min[n/i, k]}] // Flatten] // Union] ]; T[n_, k_] := Length[b[n, n, k]]; Table[Table[T[n, k], {k, 1, n}], {n, 1, 10}] // Flatten (* _Jean-François Alcover_, Mar 12 2015, after _Alois P. Heinz_ *)
%Y Cf. A026820, A070289.
%K nonn,tabl
%O 1,3
%A _Sergei Viznyuk_, Jun 01 2012