First time here? Check out the FAQ!
answered 5 years ago
Blimey! I figured it out.
Pari/GP script should have been:
for(x=1, 1000000, for(y=1, 1000000, if(floor((sigma(x)/y)+(sigma(y)/x))==2*x-sigma(x),print(x," ",y))))
updated 5 years ago
(the original had 2x instead of 2*x).
2x
2*x