OFFSET
1,2
EXAMPLE
MATHEMATICA
f[p_, e_] := 4*(Floor[p/4] + Floor[Mod[p, 4]/2]) + (2 - Mod[p, 4])*I; f[2, e_] := 2; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[500000], SquareFreeQ[#] && Im[s[#]] == 0 &] (* Amiram Eldar, Feb 24 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 08 2002
EXTENSIONS
More terms from Amiram Eldar, Feb 24 2024
STATUS
approved