OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..75
EXAMPLE
MATHEMATICA
sigma[p_, e_] := If[Abs[p] == 1, 1, (p^e + 1)]; tau[p_, e_] := If[Abs[p] == 1, 1, 2]; unitaryHarmoincQ[n_] := Divisible[n * Times @@ tau @@@ (f = FactorInteger[n, GaussianIntegers -> True]), Times @@ sigma @@@ f]; Select[Range[10^6], unitaryHarmoincQ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 13 2020
STATUS
approved