login
Number of prime parts, counted without multiplicity, in all compositions of n.
3

%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