OFFSET
1,1
COMMENTS
The unitary version of A181595.
All the terms up to a(23) are divisible by 2^8 * 3 * 5. - Giovanni Resta, Apr 26 2018
EXAMPLE
295680 is in the sequence since usigma(295680) - 2*295680 = 592128 - 591360 = 768 and 768 is a unitary divisor of 295680.
MATHEMATICA
usigma[n_] := If[n == 1, 1, Times @@ (1 + Power @@@ FactorInteger[n])]; aQ[n_] :=
Module[{d}, d = usigma[n] - 2 n; If[d <= 0, False, Divisible[n, d] && GCD[d, n/d] == 1]]; n = 1; seq={}; Do[ If[aQ[n], AppendTo[seq, n]]; n++, {k, 1, 300000}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 22 2018
EXTENSIONS
a(9)-a(23) from Giovanni Resta, Apr 26 2018
STATUS
approved