%I #12 May 30 2022 08:14:58
%S 0,0,1,3,6,15,33,74,160,344,731,1544,3237,6753,14022,29009,59819,
%T 123010,252341,516560,1055476,2153115,4385889,8922556,18131000,
%U 36805009,74643126,151255021,306267833,619719217,1253191291,2532750315,5116124712,10329574480
%N Number of prime parts, counted without multiplicity, in all compositions of n.
%H Alois P. Heinz, <a href="/A336632/b336632.txt">Table of n, a(n) for n = 0..1000</a>
%e a(4) = 0 + 1 + 1 + 1 + 1 + 1 + 1 + 0 = 6: 1111, 11(2), 1(2)1, (2)11, (2)2, 1(3), (3)1, 4.
%p b:= proc(n, i, p) option remember; `if`(n=0, [p!, 0],
%p `if`(i<1, 0, add((h-> [0, `if`(j>0 and isprime(i),
%p h[1], 0)]+h)(b(n-i*j, i-1, p+j)/j!), j=0..n/i)))
%p end:
%p a:= n-> b(n$2, 0)[2]:
%p seq(a(n), n=0..38);
%t b[n_, i_, p_] := b[n, i, p] = If[n == 0, {p!, 0},
%t If[i < 1, 0, Sum[Function[h, {0, If[j > 0 && PrimeQ[i],
%t h[[1]], 0]} + h][b[n - i*j, i - 1, p + j]/j!], {j, 0, n/i}]]];
%t a[n_] := b[n, n, 0][[2]];
%t Table[a[n], {n, 0, 38}] (* _Jean-François Alcover_, May 30 2022, after _Alois P. Heinz_ *)
%Y Cf. A000040, A102291, A336579.
%K nonn
%O 0,4
%A _Alois P. Heinz_, Jul 28 2020