【learning】多项式开根详解+模板

概述

多项式开跟是一个非常重要的知识点,许多多项式题目都要用到这一算法。

用快速数论变换,多项式求逆元和倍增法可以在$O(n log n)$的时间复杂度下求出一个$n$次多项式的开根。

 

前置技能

快速数论变换(NTT),多项式求逆元,二次剩余。

 

多项式的开根

给定一个多项式$A(x)$,其次数为$deg_A$,若存在一个多项式$B(x)$,使其满足$deg_B≤deg_A$,且$ B^2(x) \equiv A(x) (mod\ x^n)$,则$B(x)$即为$A(x)$在模$x^n$意义下的的开根。

 

求多项式的开根

我们不妨假设,$n=2^k,k∈N$。

若$n=1$,则$b_0^2 \equiv a_0(mod\ x)$。其中$a_0$,$b_0$表示多项式$A$和多项式$B$的常数项。

若需要求出$b_0$,直接用二次剩余求出即可。(不过由于大部分题目公式特别优美,常数项通常为0)

当$n>1$时:

我们假设在模$x^{\frac{n}{2}}$的意义下$A(x)$的开根$B'(x)$我们已经求得。

依据定义,则有

$B^2(x)\equiv A(x) (mod\ x^{n})$          $(1)$

因为多项式$B'(x)$满足

$B'^2(x)\equiv A(x) (mod\ x^{\frac{n}{2}})$          $(2)$

将$(1)-(2)$,得

$B^2(x)-B'^2(x)\equiv 0(mod\ x^{\frac{n}{2}})$          $(3)$

我们用平方差公式将$(3)$展开

$(B(x)+B'(x))(B(x)-B'(x))\equiv 0(mod\ x^{\frac{n}{2}})$          $(4)$

然后我们会发现多项式$B(x)$有两个解,下面考虑其中一种情况

$B(x)-B'(x)\equiv 0(mod\ x^{\frac{n}{2}})$          $(5)$

我们将$(5)$式同余符号两边平方,得

$B^2(x)-2B(x)B'(x)+B'^2(x)\equiv 0(mod\ x^n)$          $(6)$

考虑到式$(1)$,我们将式$(1)$代入式$(6)$,得

$A(x)-2B(x)B'(x)+B'^2(x)\equiv 0(mod\ x^n)$          $(7)$

最后对式 $(7)$进行移项,得到

$B(x)\equiv \dfrac{A(x)+B'^2(x)}{2B'(x)} (mod\ x^n)$          $(8)$

 

显然,我们可以用上述式子,通过倍增,求出n次项的$B(x)$。

这一步的计算我们可以使用多项式求逆元和$NTT$,时间复杂度为$O(n log n)$。

我们可以通过递归的方法,求解出最终的$N$次项的$B(x)$

时间复杂度$T(n)=T(\dfrac{n}{2})+O(n log n)=O(n log n)$。

 

模板如下:(暂时还没有)

 

posted @ 2018-07-26 17:41  AlphaInf  阅读(2220)  评论(4编辑  收藏  举报