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`

).

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.