OFFSET
1,2
COMMENTS
The unitary arithmetic numbers are in A103826.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
a(8) = 5 because the eighth unitary arithmetic number is A103826(8) = 12, the unitary divisors of 12 are 1, 3, 4 and 12 and (1 + 3 + 4 + 12)/4 = 5.
MAPLE
with(numtheory):unitdiv:=proc(n) local A, k: A:={}: for k from 1 to tau(n) do if gcd(divisors(n)[k], n/divisors(n)[k])=1 then A:=A union {divisors(n)[k]} else A:=A fi od end:utau:=n->nops(unitdiv(n)):usigma:=n->add(unitdiv(n)[j], j=1..nops(unitdiv(n))): p:=proc(n) if type(usigma(n)/utau(n), integer)=true then usigma(n)/utau(n) else fi end:seq(p(n), n=1..109);
MATHEMATICA
Select[Table[Mean[Select[Divisors[n], GCD[#, n/#]==1&]], {n, 150}], IntegerQ] (* Harvey P. Dale, May 20 2012 *)
Select[Times @@ ((1 + Power @@@ FactorInteger[#])/2) & /@ Range[100], IntegerQ] (* Amiram Eldar, Jun 14 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Emeric Deutsch, Feb 17 2005
STATUS
approved