(m, n, p)
(1,2,3)
(2,1,3)
-------
p3s = [i**3 for i in range(1000) if isprime(i)]
big = 1000
print('(m, n, p)')
for m in range(1, big+1):
for n in range(-big, big+1):
lhs = (m**3 + n)*(n**3 + m)
if lhs not in p3s:
continue
p = round( lhs ** (1/3) )
print('({},{},{})'.format(m,n,p))
|
Posted by Larry
on 2025-04-09 15:00:58 |