我想在 Python 中计算多元多项式的余数，我发现 sympy 中的 div() 应该可以解决问题（我还需要 sympy 来进行 Gröbner 计算）。但我一直发现的问题是 div() 似乎只检查除法的前导项，因为
q, r = div(x**2 + y, 2*x)
给出r=y，而
q, r = div(x**2 + y, 2*y)
给r=x**2+y.
我想做诸如理想成员资格之类的事情，因此在除以 G = {g_1，...，g_s} 时找到一些多项式 f 的余数，其中我现在不能依赖 div()。
在使用 Sage 时，我没有遇到这个问题（使用 (x^2+y)%y 给出 x^2），但我更熟悉 Python，更喜欢通过 Python 来完成。
如果我做错了什么，有人可以告诉我吗？或者有人知道用于余数的更好功能吗？