sympy - 多项式求逆的两种实现

Question

我不是很懂数学，但在dup_revert找到的用于计算多项式模乘逆的算法似乎有点类似于_series_inversion1。

请问有谁知道有什么区别吗？

标签： sympy

Answer 1

模块的想法rs_series是使用多项式工具来有效地处理幂级数，因此它采用类似的算法也就不足为奇了。它不简单地借用所有多项式方法的一个原因是该rs_series模块也适用于 Puiseux 系列（其中指数是有理数，不一定是整数）。例如：

from sympy.polys.domains import QQ
from sympy.polys.rings import ring
from sympy.polys.ring_series import _series_inversion1
R, x = ring('x', QQ)
p = x**(S(2)/3) + 1
_series_inversion1(p, x, 4)  

返回-x**(10/3) + x**(8/3) - x**2 + x**(4/3) - x**(2/3) + 1，而不是

R.dup_revert(p, 4)  

导致“TypeError：'Rational'对象不能被解释为整数”