OFFSET
1,1
MATHEMATICA
s[n_] := DivisorSum[n, # &, #^2 < n && CoprimeQ[#, n/#] &]; seq={}; s1 = s[1]; s2 = s[2]; s3 = s[3]; Do[s4 = s[n]; If[s1 == s2 && s2 == s3 && s3 == s4, AppendTo[seq, n - 3]]; s1 = s2; s2 = s3; s3 = s4, {n, 4, 10^9}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 12 2020
STATUS
approved