多项式——多项式牛顿迭代

多项式——多项式牛顿迭代

多项式的牛顿迭代法是多项式函数的基础算法。考虑方程 $F(P(x))\equiv 0\mathrm{mod}{x}^n$ 其中 $F(x)$ 是一个已知的多项式插值函数（其实是可以泰勒展开的任意函数）定义为：

$$F(x)=\sum _{i=0}^{\infty }{\alpha }_i(x-\beta {)}^i$$

${\alpha }_i$ 称为系数， $\beta$ 称为展开点，都是定值。

$P(x)$ 是一个多形式，现在请求出 $P(x)$ 。

将 $F(x)$ 在 $y$ 点泰勒展开得到：

$$F(x)=F(y)+(x-y)F^{\prime }(y)+(x-y{)}^{2}G(x,y)$$

其中 $F^{\prime }(x)$ 是 $F(x)$ 的导数，定义为：

$$F^{\prime }(x)=\sum _{i=0}^{\infty }(i+1){\alpha }_{i+1}(x-\beta {)}^i$$

而 $G(x,y)$ 是关于 $x$ 和 $y$ 的形式幂级数。

令 $F(Q_k)\equiv 0\mathrm{mod}{x}^a$ 我们想找到 $F(Q_{k+1})\equiv 0\mathrm{mod}{x}^{2a}$ 满足 $Q_{k+1}\equiv Q_k+x^{a}C\mathrm{mod}{x}^{2a}$ 。

带入 $x=Q_{k+1}$ 和 $y=Q_k$ 我们得到：

$$F(Q_{k+1})\equiv F(Q_k)+(Q_{k+1}-Q_k)F^{\prime }(Q_k)+(Q_{k+1}-Q_k{)}^{2}G(x,y)\mathrm{mod}{x}^{2a}.$$

因为 $Q_{k+1}-Q_k\equiv 0\mathrm{mod}{x}^a$ 那么 $(Q_{k+1}-Q_k{)}^2\equiv 0\mathrm{mod}{x}^{2a}$ 。因此：

$$0\equiv F(Q_{k+1})\equiv F(Q_k)+(Q_{k+1}-Q_k)F^{\prime }(Q_k)\mathrm{mod}{x}^{2a}$$

等式的第三部分告诉了我们 $Q_k$ 和 $Q_{k+1}$ 的关系，即为：

$$Q_{k+1}\equiv Q_k-\frac{F(Q_k)}{F^{\prime }(Q_k)}\mathrm{mod}{x}^{2a}$$

你已经学会如何对多项式取逆和如何计算 $F(Q_k)$ 了，那么我们就能在 $O(n\log n)$ 的情况下计算出多项式 $P(x)$ 。

上述过程就称为多项式的牛顿迭代法，可用于计算多项式函数。