Revision history [back]
 | 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).
