OFFSET
1,2
COMMENTS
EXAMPLE
10 is a term since its sum of unitary divisors in Gaussian integers is -12 + 16*i, whose norm (-12)^2 + 16^2 = 400 is divisible by 10^2 = 100.
MATHEMATICA
f[p_, e_] := If[Abs[p] == 1, 1, (p^e + 1)]; Select[Range[21000], Divisible[Abs[ Times @@ f @@@ FactorInteger[#, GaussianIntegers -> True]]^2, #^2] &]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Feb 13 2020
STATUS
approved