|
%I #24 Oct 27 2023 20:50:12
%S 0,0,1,1,2,0,3,4,8,2,4,1,5,8,24,24,6,2,7,15,107,46,8,4,135,101,347,83,
%T 9,0,10,460,1019,431,1308,13,11,842,2760,214,12,2,13,1418,5124,2977,
%U 14,42,2021,720,17355,4997,15,70,21108,3674,40702,16907,16,1,17
%N The number of phi-partitions of n.
%C The number of partitions n = a1+a2+...+ak which have at least two parts and obey phi(n) = phi(a1)+phi(a2)+...+phi(ak). phi(.) = A000010(.) is Euler's totient. The trivial result with one part, n=a1, is not counted; that would induce another sequence with terms a(n)+1.
%H Alois P. Heinz, <a href="/A283528/b283528.txt">Table of n, a(n) for n = 1..1000</a> (first 100 terms from Giovanni Resta)
%H P. Jones, <a href="http://www.fq.math.ca/Scanned/29-4/jones.pdf">Phi-partitions</a>, Fib. Quart. 29 (4) (1991) 347-350.
%H C. Powell, <a href="http://www.fq.math.ca/Scanned/34-3/powell.pdf">On the uniqueness of reduced phi-partitions</a>, Fib. Quart. 34 (3) (1996) 194-199.
%H J. Wang, <a href="http://www.fq.math.ca/Scanned/31-4/wang.pdf">Reduced phi-partitions of positive integers</a>, Fib. Quart. 31 (4) (1993) 365-369.
%e a(7) = 3 counts the partitions 1+1+1+1+1+2 = 1+1+1+1+3 = 1+1+5.
%e a(8) = 4 counts the partitions 2+2+2+2 = 2+2+4 = 4+4 = 1+1+6.
%p A283528 := proc(n)
%p local a,k,phip ;
%p a := 0 ;
%p for k in combinat[partition](n) do
%p if nops(k) > 1 then
%p phip := add( numtheory[phi](p),p =k) ;
%p if phip = numtheory[phi](n) then
%p a := a+1 ;
%p end if;
%p end if;
%p end do:
%p a ;
%p end proc:
%p # second Maple program:
%p with(numtheory):
%p b:= proc(n, m, i) option remember; `if`(n=0,
%p `if`(m=0, 1, 0), `if`(i<1 or m<0, 0, b(n, m, i-1)+
%p `if`(i>n, 0, b(n-i, m-phi(i), i))))
%p end:
%p a:= n-> b(n, phi(n), n)-1:
%p seq(a(n), n=1..70); # _Alois P. Heinz_, Mar 10 2017
%t Table[ Length@ IntegerPartitions[n 10^7 + EulerPhi[n], {2, Infinity},
%t EulerPhi@ Range[n-1] + 10^7 Range[n-1]], {n, 60}] (* _Giovanni Resta_, Mar 10 2017 *)
%t b[n_, m_, i_] := b[n, m, i] = If[n == 0,
%t If[m == 0, 1, 0], If[i < 1 || m < 0, 0, b[n, m, i - 1] +
%t If[i > n, 0, b[n - i, m - EulerPhi[i], i]]]];
%t a[n_] := b[n, EulerPhi[n], n]-1;
%t Array[a, 70] (* _Jean-François Alcover_, Jun 01 2021, after _Alois P. Heinz_ *)
%Y Cf. A000010, A271384, A283530.
%K nonn
%O 1,5
%A _R. J. Mathar_, Mar 10 2017
%E a(56)-a(61) from _Giovanni Resta_, Mar 10 2017
|
