python - この Wolfram Solve の例を sympy で複製する方法

Question

Wolfram言語で得られた結果をsympyで再現したい.

Wolframcloud を使用すると、この式

Solve[m^2+m*n==500 && m>n,{m,n},PositiveIntegers]

私が探している結果が得られます：

{{m->20,n->5}}

sympy を使用してこれを再現するにはどうすればよいですか?

私が試してみました

import sympy as sp
m,n = sp.symbols('m n',integer=True)
sp.solve(m**2 + m*n - 500, m,n)

を与える

[(m, -m + 500/m)]

これは正しいですが、特に役に立ちません。

この質問は Project Euler Problem 9 に触発されていることに注意してください。