a(m,n)=forprime(p=2,5000,if(ispseudoprime((m^p-n^p)/(m-n)),print1(p, ", ")))
for(m=2,128,for(n=-m,m,if(gcd(m,n)==1,print1(m, ", ", n, ": ");a(m,n);print())))