牛顿迭代

给定多项式$G(x)$，已知$A(x)$满足
$$G(A(x))\equiv 0mod{x}^n$$
求模$x^n$意义下的$A(x)$

倍增法

首先解出
$$\begin{array}{r}[x^0]G(A(x))\equiv 0modx\end{array}$$
因为是在模$x$意义下解方程，所以$A(x)$求出来也是在模$x$意义下，所以$A(x)$是常数。
假设求出了模$x^{\lceil \frac{n}{2}\rceil }$意义下的解$A^{{}^{\prime }}(x)$，要求模$x^n$意义下的解$A(x)$
因为$G(A(x))\equiv 0mod{x}^n$和$x^{\lceil \frac{n}{2}\rceil }|x^n$，则
$$\{\begin{array}{ll}G(A(x))& \equiv 0mod{x}^{\lceil \frac{n}{2}\rceil }\\ G(A^{{}^{\prime }}(x))& \equiv 0mod{x}^{\lceil \frac{n}{2}\rceil }\end{array}\to A(x)\equiv A^{{}^{\prime }}(x)mod{x}^{\lceil \frac{n}{2}\rceil }$$
把$G(A(x))$在$A^{{}^{\prime }}(x)$处泰勒展开，得
$$G(A(x))=\sum _{i=0}^{\infty }\frac{G^{(i)}(A^{{}^{\prime }}(x))}{i!}[A(x)-A^{{}^{\prime }}(x){]}^i$$
因为$A(x)-A^{{}^{\prime }}(x)\equiv 0mod{x}^{\lceil \frac{n}{2}\rceil }$，所以$[A(x)-A^{{}^{\prime }}(x){]}^i\equiv 0mod{x}^{i\lceil \frac{n}{2}\rceil }$
当$i\ge 2$时，$i\lceil \frac{n}{2}\rceil \ge n$，所以$x^n|x^{i\lceil \frac{n}{2}\rceil }$，所以$[A(x)-A^{{}^{\prime }}(x){]}^i\equiv 0mod{x}^n$。所以泰勒展开在模$x^n$得意义下只保留前两项。
$$\begin{array}{rl}G(A(x))& \equiv G(A^{{}^{\prime }}(x))+G^{(1)}(A^{{}^{\prime }}(x))[A(x)-A^{{}^{\prime }}(x)]mod{x}^n\\ 0& \equiv G(A^{{}^{\prime }}(x))+G^{(1)}(A^{{}^{\prime }}(x))[A(x)-A^{{}^{\prime }}(x)]mod{x}^n\\ A(x)& \equiv A^{{}^{\prime }}(x)-\frac{G(A^{{}^{\prime }}(x))}{G^{(1)}(A^{{}^{\prime }}(x))}mod{x}^n\end{array}$$

应用：

• 多项式求逆元
  给定多项式$A(x)$，求$B(x)$满足$A(x)B(x)\equiv 1mod{x}^n$。
  $$\begin{array}{rl}A(x)B(x)& \equiv 1mod{x}^n\\ A(x)-\frac{1}{B(x)}& \equiv 0mod{x}^n\end{array}$$
  首先把$A(x)-\frac{1}{B(x)}$看成关于$B(x)$的函数，则设$G(B(x))=A(x)-\frac{1}{B(x)}$，所以$G(B(x))\equiv 0mod{x}^n$。和上面的形式一样，所以就有。
  $$\begin{array}{rl}B(x)& \equiv B^{{}^{\prime }}(x)-\frac{G(B^{{}^{\prime }}(x))}{G^{(1)}(B^{{}^{\prime }}(x))}mod{x}^n\\ & \equiv B^{{}^{\prime }}(x)-\frac{A(x)-\frac{1}{B^{{}^{\prime }}(x)}}{\frac{1}{B^{{}^{\prime }}(x{)}^2}}mod{x}^n\\ & \equiv B^{{}^{\prime }}(x)-A(x)B^{{}^{\prime }}(x{)}^2+B^{{}^{\prime }}(x)mod{x}^n\\ & \equiv 2B^{{}^{\prime }}(x)-A(x)B^{{}^{\prime }}(x{)}^{2}mod{x}^n\end{array}$$
  这样就可以用倍增计算
• 多项式开方
  给定多项式$A(x)$，求$B(x)$满足$A(x)\equiv B(x{)}^{2}mod{x}^n$。
  设$G(B(x))=A(x)-B(x{)}^2$，则
  $$\begin{array}{rl}B(x)& \equiv B^{{}^{\prime }}(x)-\frac{G(B^{{}^{\prime }}(x))}{G^{(1)}(B^{{}^{\prime }}(x))}mod{x}^n\\ & \equiv B^{{}^{\prime }}(x)-\frac{A(x)-B^{{}^{\prime }}(x{)}^2}{-2B^{{}^{\prime }}(x)}mod{x}^n\\ & \equiv \frac{B^{{}^{\prime }}(x)}{2}+\frac{A(x)}{2B^{{}^{\prime }}(x)}mod{x}^n\end{array}$$
• 多项式EXP
  给定多项式$A(x)$，求$B(x)$满足$e^{A(x)}\equiv B(x)mod{x}^n$。
  注意$e^{A(x)}$存在必须满足$[x^0]A(x)=0$
  两边取对数，设$G(B(x))=A(x)-\ln B(x)$，则
  $$\begin{array}{rl}B(x)& \equiv B^{{}^{\prime }}(x)-\frac{G(B^{{}^{\prime }}(x))}{G^{(1)}(B^{{}^{\prime }}(x))}mod{x}^n\\ & \equiv B^{{}^{\prime }}(x)-\frac{A(x)-\ln B^{{}^{\prime }}(x)}{-\frac{1}{B^{{}^{\prime }}(x)}}mod{x}^n\\ & \equiv B^{{}^{\prime }}(x)[1-\ln B^{{}^{\prime }}(x)+A(x)]mod{x}^n\end{array}$$
• 多项式的幂
  计算$f^k(x)mod{x}^n$

1. 当$[x^0]f(x)=1$时，$f^k(x)\equiv \exp (k\ln f(x))mod{x}^n$
2. 当$[x^0]f(x)\ne 1$时，设$f(x)$的最低次项为$a_i{x}^i$，则
   $f^k(x)\equiv a_i^k{x}^{ki}\exp (k\ln \frac{f(x)}{a_i{x}^i})mod{x}^n$