editing
approved
editing
approved
with(numtheory); P:=proc(q) local a, n;
for n from 1 to q do a:=sigma(n)/n; if type(a, integer) and a>2 then print(n);
fi; od; end: P(10^6); # Paolo P. Lava, Oct 24 2013
approved
editing
proposed
approved
editing
proposed
Achim Flammenkamp, <a href="http://wwwwwwhomes.uni-bielefeld.de/~achim/mpn.html">The Multiply Perfect Numbers Page</a>
approved
editing
proposed
approved
editing
proposed
Achim Flammenkamp, <a href="http://www.uni-bielefeld.de/~achim/mpn.html">The Multiply Perfect Numbers Page</a>
(PARI) isok(n) = sn = sigma(n)/n ; (type(sn) == "t_INT") && (sn > 2); \\ Michel Marcus, Oct 24 2013
proposed
editing
editing
proposed
with(numtheory); P:=proc(q) local a, n;
for n from 1 to q do a:=sigma(n)/n; if type(a, integer) and a>2 then print(n);
fi; od; end: P(10^6); # Paolo P. Lava, Oct 24 2013
approved
editing