多项式求逆的两种实现

时间:2018-06-28 11:38:41

标签: sympy

我不太懂数学,但是似乎在dup_revert处发现的用于计算多项式的模乘逆的算法似乎有点类似于_series_inversion1

请,有人知道有什么区别吗?

1 个答案:

答案 0 :(得分:0)

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'对象无法解释为整数”