OFFSET
1,1
COMMENTS
a(7) > 10^12. - Giovanni Resta, Mar 11 2014
EXAMPLE
If k = 17296 then k' = sigma(k) = 35712. If k = 24016 then k' = sigma(k) = 49600.
MAPLE
with(numtheory); P:= proc(q) local a, n, p;
for n from 1 to q do a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
if sigma(n)=a then print(n); fi; od; end: P(10^9);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Oct 14 2013
EXTENSIONS
a(3)-a(5) from Giovanni Resta, Oct 14 2013
a(6) from Giovanni Resta, Mar 11 2014
STATUS
approved