OFFSET
1,2
COMMENTS
The unitary version of A260653.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Nicholas Schwab and Lola Thompson, A generalization of the practical numbers, arXiv:1701.08504 [math.NT], 2017.
EXAMPLE
The unitary divisors of 12 are 1, 3, 4 and 12, and the set of their uphi values is {1, 2, 3, 6}. Each number below 12 is the sum of a subset, e.g., 11 = 2 + 3 + 6, 10 = 1 + 3 + 6, etc.
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[10000], uPhiPracticalQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 15 2017
STATUS
approved