Revision history [back]

Got it! There was indeed a missing asterisk, after all!

The Pari-GP script should have been coded as such:

for(x=1, 100000, for(y=1, 100000, if(2*x*y-sigma(x*y) >= (2*x-sigma(x))*(sigma(y)-y),print(x,"   ",y))))

I am able to run this version of the script now.