多项式求逆 与 多项式开方

介绍与基本概念

多项式求逆元是多项式除法、多项式开根和多项式取模的必要过程，用快速傅里叶变换及倍增算法可以做到用$O$($n$ $log$ $n$)的时间复杂度求出一个多项式的逆元。

对于一个多项式 $A{(x)}$，称其最高项的次数为这个多项式的度（$degree$），记作 $deg_A$
对于多项式 $A{(x)}$,$B{(x)}$，存在唯一的$Q{(x)}$,$R{(x)}$满足 $A{(x)}$=$Q{(x)}$*$B{(x)}$+$R{(x)}$，其中$deg_R$<$deg_B$，我们称$Q{(x)}$为$B{(x)}$除$A{(x)}$的商，$R{(x)}$为$B{(x)}$除$A{(x)}$的余数。

多项式求逆

对于一个多项式$A{(x)}$，存在$B{(x)}$满足$deg_B$≤$deg_A$并且

$A{(x)}$$B{(x)}$≡$1$ ($mod$ $x^n$)
那么称$B{(x)}$为$A{(x)}$在($mod$ $x^n$)意义下的逆元，记作${A^{−1}}{(x)}$。

接下来讲一下如何求多项式的逆元。
当$n$=$1$时，直接求常数项的逆元即可。
当$n$>$1$时，考虑如何从($mod$ $x^{n}$)意义下的逆元推到($mod$ $x^{2n}$)意义下的逆元。

假设$A{(x)}在$($mod$ $x^{n}$)意义下的逆元为$H{(x)}$，$A{(x)}$在($mod$ $x^{2n}$)意义下的逆元为$G{(x)}$，那我们就是已知$H{(x)}$求$G{(x)}$。

由定义可知，

$$A{(x)}H{(x)}≡1 (mod_{ }\ x^n)$$
移项得，
$$A{(x)}H{(x)}-1≡0 (mod_{ }\ x^n)$$
考虑对此式子的两边平方，这样就能变成在( $mod$ $x^{2n}$)意义下了，便有
$$A^2(x)H^2(x)-2A(x)*H(x)+1≡0 (mod_\ x^{2n})$$
移项便得
$$2A(x)*H(x)-A^2(x)H^2(x)≡1 (mod_\ x^{2n})$$
$$A(x)[2H(x)-A(x)H^2(x)]≡1 (mod_\ x^{2n})$$
所以
$$G{(x)}≡2H(x)-A(x)H^2(x)$$
这样就求出 $G{(x)}$了，时间复杂度为 $T(n)$= $O$($n$ $log$ $n$)+$T$( $n \over 2$)= $O$($n$ $log$ $n$)。

多项式开根

给定$A{(x)}$，求 $B{(x)}$满足 $B^2{(x)}$≡ $1$ ($mod$ $x^n$)
当 $n$=$1$时，这个比较麻烦，待会再说。
当 $n$>$1$时，假设已求出了 $B^2{(x)}$≡ $1$ ($mod$ $x^n$)，现在要求出 $G^2{(x)}$≡ $1$ ($mod$ $x^{2n}$)
同样两式相减可得

$$G^2{(x)}-B^2{(x)}≡0 (mod \ x^{n})$$
一样的套路，为了使其变成在( $mod$ $x^{2n}$)意义下，对式子两边平方，得
$$G^4{(x)}-2G^2{(x)}B^2{(x)}+B^4{(x)}≡0 (mod \ x^{2n})$$
$$G^4{(x)}+2G^2{(x)}B^2{(x)}+B^4{(x)}≡4G^2{(x)}B^2{(x)} (mod \ x^{2n})$$
此时两边都能化成平方的形式，妙哉
$$(G^2{(x)}+B^2{(x)})^2≡(2G{(x)}B{(x)})^2 (mod \ x^{2n})$$
$$G^2{(x)}+B^2{(x)}≡2G{(x)}B{(x)} (mod \ x^{2n})$$
由 $G^2{(x)}$≡ $1$ ($mod$ $x^{2n}$)可得
$$B^2{(x)}+A{(x)}≡2G{(x)}B{(x)} (mod \ x^{2n})$$
两边同时除掉 $2G{(x)}$，可得
$${B{(x)}\over 2}+{A{(x)}\over 2B{(x)}}≡G{(x)} (mod \ x^{2n})$$

那只用求$B^{-1}{(x)}$即可，求多项式逆元的复杂度是$O$($n$ $log$ $n$)的，因而总的时间复杂度为$T(n)$=$O$($n$ $log$ $n$)+$T$($n \over 2$)=$O$($n$ $log$ $n$)。

现在回到问题当$n$=$1$时该怎么求，其实问题等同于已知$b$,$c$，求$a$满足$a^2$≡$b$ ($mod$ $c$)，这里推介一位大佬的$blog$二次剩余Cipolla算法，这里就不细讲了。

若有地方写错了，还望各大佬指出。