|
%I #10 Apr 11 2022 08:36:10
%S 1,1,4,40,355,4425,69521,1162951,22259414,478096938,11614132907,
%T 299700810545,8456607358157,255883964141333,8275199908539114,
%U 287869753459458468,10564476589147507523,409845503129745719513,16777378294629533764699,720626728499888159724831
%N Number of colored compositions of n using all colors of an initial interval of the color palette such that all parts have different color patterns and the patterns for parts i are sorted and have i colors (in arbitrary order).
%H Alois P. Heinz, <a href="/A327675/b327675.txt">Table of n, a(n) for n = 0..200</a>
%p b:= proc(n, i, k, p) option remember;
%p `if`(n=0, p!, `if`(i<1, 0, add(binomial(k^i, j)*
%p b(n-i*j, min(n-i*j, i-1), k, p+j)/j!, j=0..n/i)))
%p end:
%p a:= n-> add(add(b(n$2, i, 0)*(-1)^(k-i)*
%p binomial(k, i), i=0..k), k=0..n):
%p seq(a(n), n=0..23);
%t b[n_, i_, k_, p_] := b[n, i, k, p] =
%t If[n == 0, p!, If[i < 1, 0, Sum[ Binomial[k^i, j]*
%t b[n - i*j, Min[n - i*j, i - 1], k, p + j]/j!, {j, 0, n/i}]]];
%t a[n_] := Sum[Sum[b[n, n, i, 0]*(-1)^(k - i)*
%t Binomial[k, i], {i, 0, k}], {k, 0, n}];
%t Table[a[n], {n, 0, 23}] (* _Jean-François Alcover_, Apr 11 2022, after _Alois P. Heinz_ *)
%Y Row sums of A327673.
%K nonn
%O 0,3
%A _Alois P. Heinz_, Sep 21 2019
|
