%I #22 Sep 11 2019 06:25:48
%S 1,1,2,7,10,25,68,111,208,435,1218,1773,3586,6077,12156,31961,47624,
%T 86825,151962,265525,469610,1242607,1750108,3217663,5263928,9205197,
%U 14713474,26440503,63610938,90877893,159360628,258871127,431309688,687140639,1134231986
%N "AFK" (ordered, size, unlabeled) transform of 1,2,3,4,...
%C Sum of products of parts in all compositions of n into distinct parts. - _Vladeta Jovovic_, Feb 21 2005
%C Number of compositions of n into distinct parts if there are i kinds of part i. a(3) = 7: 3, 3', 3'', 21, 2'1, 12, 12'.
%H Alois P. Heinz, <a href="/A032007/b032007.txt">Table of n, a(n) for n = 0..5000</a>
%H C. G. Bower, <a href="/transforms2.html">Transforms (2)</a>
%p b:= proc(n, i, p) option remember;
%p `if`(n=0, p!, `if`(i<1, 0, b(n, i-1, p)+
%p `if`(i>n, 0, i*b(n-i, i-1, p+1))))
%p end:
%p a:= n-> b(n$2, 0):
%p seq(a(n), n=0..40); # _Alois P. Heinz_, Sep 05 2015
%t b[n_, i_, p_] := b[n, i, p] = If[n == 0, p!, If[i < 1, 0, b[n, i - 1, p] + If[i > n, 0, i*b[n - i, i - 1, p + 1]]]];
%t a[n_] := b[n, n, 0];
%t a /@ Range[0, 40] (* _Jean-François Alcover_, Sep 11 2019, after _Alois P. Heinz_ *)
%o (PARI) seq(n)={apply(p->subst(serlaplace(p), y, 1), Vec(prod(k=1, n, 1 + k*x^k*y + O(x*x^n))))} \\ _Andrew Howroyd_, Jun 21 2018
%Y Cf. A022629.
%Y Cf. A032005, A032006, A032008.
%K nonn
%O 0,3
%A _Christian G. Bower_
%E a(0)=1 prepended by _Alois P. Heinz_, Sep 05 2015