%I #16 Mar 10 2023 09:10:11
%S 1,3,5,6,7,9,11,12,13,14,15,17,19,21,22,23,24,25,27,28,29,30,31,33,35,
%T 37,38,39,41,42,43,44,45,46,47,48,49,51,53,54,55,56,57,59,60,61,62,63,
%U 65,66,67,69,70,71,73,75,76,77,78,79,81,83,84,85,86,87,88,89,91,92,93
%N Unitary arithmetic numbers (those for which the arithmetic mean of the unitary divisors is an integer).
%C The arithmetic means of the unitary arithmetic numbers are in A103827.
%C From _Amiram Eldar_, Mar 10 2023: (Start)
%C Union of the odd numbers (A005408) and twice the numbers that are not the sum of 2 squares (A022544).
%C The asymptotic density of this sequence is 1. (End)
%H Charles R Greathouse IV, <a href="/A103826/b103826.txt">Table of n, a(n) for n = 1..10000</a>
%e 12 is a unitary arithmetic number because the unitary divisors of 12 are 1,3,4 and 12 and (1+3+4+12)/4=5 is an integer.
%p 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 n else fi end:seq(p(n),n=1..103);
%t udiQ[n_]:=IntegerQ[Mean[Select[Divisors[n],GCD[#,n/#]==1&]]]; Select[ Range[ 100],udiQ] _Harvey P. Dale_, May 20 2012
%t Select[Range[100], IntegerQ[Times @@ ((1 + Power @@@ FactorInteger[#])/2)] &] (* _Amiram Eldar_, Jun 14 2022 *)
%o (PARI) is(n)=my(f=factor(n)); prod(i=1,#f~, f[i,1]^f[i,2]+1)%2^#f~==0 \\ _Charles R Greathouse IV_, Sep 01 2015
%Y Cf. A005408, A022544, A103827, A034444, A034448.
%K nonn
%O 1,2
%A _Emeric Deutsch_, Feb 17 2005