OFFSET
1,1
EXAMPLE
sigma(16380) = 156*sigma(156) = 182*sigma(182).
sigma(16655184) = 4368*sigma(4368) = 4836*sigma(4836) = 5642*sigma(5642).
MATHEMATICA
Select[Range[1.5*10^6], Count[(d = Divisors[#]) * DivisorSigma[1, d], DivisorSigma[1, #]] > 1 &] (* Amiram Eldar, May 12 2021 *)
PROG
(Magma) s:=func<n|{d:d in Divisors(n)|DivisorSigma(1, n) eq DivisorSigma(1, d)*d}>; [n: n in [1..1000000]|#s(n) ge 2];
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, May 10 2021
STATUS
approved