OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..100
MATHEMATICA
f[p_, e_] := If[Abs[p] == 1, 1, (p^e + 1)]; normUsigma[n_] := Abs[Times @@ f @@@ FactorInteger[n, GaussianIntegers -> True]]^2; seq = {}; u1 = normUsigma[1]; Do[u2 = normUsigma[n]; If[u1 == u2, AppendTo[seq, n - 1]]; u1 = u2, {n, 2, 10^6}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 13 2020
STATUS
approved