|
%I #21 Apr 04 2022 11:41:51
%S 0,0,2,5,11,23,119,237,776,6665,16518,207953,892680,1824445,8374988,
%T 96208461,978217302,2059770725,18616884428,78013141907,158103168924,
%U 1386674113487,6734724875544,82189835767618,2013603833805429,9101106147506177,19147196940580651
%N Number of compositions of the n-th prime into a prime number of prime parts.
%H Alois P. Heinz, <a href="/A344782/b344782.txt">Table of n, a(n) for n = 1..300</a>
%e a(3) = 2: [2,3], [3,2].
%e a(4) = 5: [5,2], [2,5], [3,2,2], [2,3,2], [2,2,3].
%p b:= proc(n, c) option remember; `if`(n=0, `if`(isprime(c),
%p 1, 0), add(`if`(isprime(j), b(n-j, c+1), 0), j=2..n))
%p end:
%p a:= n-> b(ithprime(n), 0):
%p seq(a(n), n=1..31);
%t b[n_, c_] := b[n, c] = If[n == 0, If[PrimeQ[c], 1, 0],
%t Sum[If[PrimeQ[j], b[n - j, c + 1], 0], {j, 2, n}]];
%t a[n_] := b[Prime[n], 0];
%t Table[a[n], {n, 1, 31}] (* _Jean-François Alcover_, Apr 04 2022, after _Alois P. Heinz_ *)
%Y Cf. A000040, A316154, A344790.
%K nonn
%O 1,3
%A _Alois P. Heinz_, May 28 2021
|
