OFFSET
1,2
COMMENTS
A unitary phi-practical number k is a number k such that each number in the range 1..k is a subsum of a the multiset {uphi(d) : d | k, gcd(d, k/d) = 1}. This sequence is restricted to cases in which all the values in this multiset are distinct.
Are all the terms above 3 divisible by 5?
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
uphi[n_] := If[n == 1, 1, (Times @@ (Table[#[[1]]^#[[2]] - 1, {1}] & /@ FactorInteger[n]))[[1]]];
uDivisors[n_] := Select[Divisors[n], GCD[#, n/#] == 1 &]; uPhiPracticalQ[n_] := If[n < 1, False, If[n == 1, True, (lst = Sort @ Map[uphi, uDivisors[n]]; ok = True; Do[If[lst[[m]] > Sum[lst[[l]], {l, 1, m - 1}] + 1, (ok = False; Break[])], {m, 1, Length[lst]}]; ok)]];
Select[Range[9000], UnsameQ @@ uphi /@ Divisors[#] && uPhiPracticalQ[#] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Dec 31 2022
STATUS
approved