%I #14 May 02 2022 02:58:42
%S 1,0,1,4,22,241,2696,35218,529888,8998419,169486964,3496417024,
%T 78344008779,1891733424205,48923563968087,1347813311456319,
%U 39371345548420060,1214570579814316742,39430967625404799740,1343040950675651131103,47862610677098010505554
%N Number of compositions of n^2 into n prime parts.
%H Alois P. Heinz, <a href="/A341459/b341459.txt">Table of n, a(n) for n = 0..100</a>
%F a(n) = A121303(n^2,n).
%e a(3) = 4: 333, 225, 252, 522.
%p b:= proc(n, t) option remember;
%p `if`(n=0, `if`(t=0, 1, 0), `if`(t<1, 0, add(
%p `if`(isprime(j), b(n-j, t-1), 0), j=1..n)))
%p end:
%p a:= n-> b(n^2, n):
%p seq(a(n), n=0..22);
%t b[n_, t_] := b[n, t] =
%t If[n == 0, If[t == 0, 1, 0], If[t < 1, 0, Sum[
%t If[PrimeQ[j], b[n-j, t-1], 0], {j, 1, n}]]];
%t a[n_] := b[n^2, n];
%t Table[a[n], {n, 0, 22}] (* _Jean-François Alcover_, May 02 2022, after _Alois P. Heinz_ *)
%Y Cf. A023360, A035026, A121303, A299168.
%K nonn
%O 0,4
%A _Alois P. Heinz_, Feb 12 2021
|