Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #14 Dec 22 2024 11:10:58
%S 1,0,1,2,3,3,5,5,7,8,10,12,15,15,18,21,23,26,31,33,36,41,43,48,52,58,
%T 62,72,72,82,85,95,97,112,112,125,127,142,142,161,159,181,180,200,196,
%U 222,217,243,239,269,261,291,284,316,308,341,332,370,358,394,381,427,414,456
%N Number of partitions of n into 5 distinct prime powers (including 1).
%H Alois P. Heinz, <a href="/A341142/b341142.txt">Table of n, a(n) for n = 15..5000</a>
%p q:= proc(n) option remember; nops(ifactors(n)[2])<2 end:
%p b:= proc(n, i, t) option remember; `if`(n=0,
%p `if`(t=0, 1, 0), `if`(i<1 or t<1, 0, b(n, i-1, t)+
%p `if`(q(i), b(n-i, min(n-i, i-1), t-1), 0)))
%p end:
%p a:= n-> b(n$2, 5):
%p seq(a(n), n=15..78); # _Alois P. Heinz_, Feb 05 2021
%t q[n_] := q[n] = PrimeNu[n] < 2;
%t b[n_, i_, t_] := b[n, i, t] = If[n == 0,
%t If[t == 0, 1, 0], If[i < 1 || t < 1, 0, b[n, i - 1, t] +
%t If[q[i], b[n - i, Min[n - i, i - 1], t - 1], 0]]];
%t a[n_] := b[n, n, 5];
%t Table[a[n], {n, 15, 78}] (* _Jean-François Alcover_, Jul 13 2021, after _Alois P. Heinz_ *)
%t Table[Count[IntegerPartitions[n,{5}],_?(Max[PrimeNu[#]]<2&&Length[#]==Length[Union[#]]&)],{n,15,80}] (* _Harvey P. Dale_, Dec 22 2024 *)
%Y Cf. A000961, A010055, A341123, A341132, A341134, A341140, A341141, A341143, A341144, A341145.
%K nonn
%O 15,4
%A _Ilya Gutkovskiy_, Feb 05 2021