我想在 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 来完成。
如果我做错了什么,有人可以告诉我吗?或者有人知道用于余数的更好功能吗?