1 | initial version |
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))))
2 | No.2 Revision |
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))))
(the original had 2x
instead of 2*x
).