[拉格朗日反演][FFT][NTT][多项式大全]详解

1、目录

2、多项式的两种表示法

1.系数表示法

我们最常用的多项式表示法就是系数表示法，一个次数界为$n$的多项式$S(x)$可以用一个向量$s=(s_0,s_1,s_2,\cdots,s_n-1)$系数表示如下：

$$S(x)=\sum_{k=0}^{n-1}s_kx^k$$

系数表示法很适合做加法，可以在$O(n)$的时间复杂度内完成，表达式为：

$$S(x)=A(x)+B(x)=\sum_{k=0}^{n-1}(a_k+b_k)x^k$$

当中

$$s_k=a_k+b_k$$

但是，系数表示法不适合做乘法，时间复杂度为$O(n^2)$，表达式为：

$$S(x)=A(x)B(x)=\sum_{k=0}^{n-1}\left(\sum_{j=0}^{n-1}a_j b_{k-j}\right)x^k$$

当中

$$s_k=\sum_{j=0}^k a_j b_{k-j}$$

这就是卷积的一般形式，记$s=a\otimes b$，我们要想办法加速这个过程。

2.点值表示法

顾名思义，点值就是多项式在一个点处的值。多项式$A(x)$的点值表达是一个集合：

$$\{(x_0,y_0),(x_1,y_1),(x_2,y_2),\cdots,(x_{n-1},y_{n-1})\}$$

使得对于$k=0,1,2,\cdots,n-1$有$x_k$两两不相同且$y_k=A(x_k)$。

$n$个点可以确定唯一一个$n$次多项式。

点值表达有很多优良的性质，加法和乘法都可以在$O(n)$的时间复杂度内完成。

现有$A(x)$的点值表达

$$\{(x_0,y_0),(x_1,y_1),(x_2,y_2),\cdots,(x_{n-1},y_{n-1})\}$$ 和 $B(x)$的点值表达 $$\{(x_0,y'_0),(x_1,y'_1),(x_2,y'_2),\cdots,(x_{n-1},y'_{n-1})\}$$

则$C(x)=A(x)+B(x)$的点值表达为：

$$\{(x_0,y_0+y'_0),(x_1,y_1+y'_1),(x_2,y_2+y'_2),\cdots,(x_{n-1},y_{n-1}+y'_{n-1})\}$$

$C(x)=A(x)B(x)$的点值表达为：

$$\{(x_0,y_0 y'_0),(x_1,y_1 y'_1),(x_2,y_2 y'_2),\cdots,(x_{n-1},y_{n-1} y'_{n-1})\}$$

可见，点值表示可以帮助我们更快地进行卷积，可是如何在系数表示法和点值表示法之间相互转化呢？

3、复数

当$x$为实数时，无法很好地对转换方法进行优化。为了优化计算$x^n$所浪费的时间，我们需要$x$有循环的性质。但点值表示法需要$n$个两两不同的值，而在实数域中只有$1$和$-1$，因此，我们需要复数的帮助。

1.复数、复平面的定义

我们把形如$a+bi$的数称为复数$z$，其中$a$为实部(Real)，记为$\Re z$；$b$为虚部(Imaginary)，记为$\Im z$。

每一点都对应唯一复数的平面叫复平面，相当于一个把$\Re z$作为横坐标，把$\Im z$作为纵坐标的笛卡尔坐标系。如图：

模长：复平面上原点到复数$z$的距离，记为$|z|$。根据勾股定理，$|z|=|a+bi|=\sqrt{a^2+b^2}$

辐角：复平面上$x$轴与复数$z$所对应向量之间的夹角，在$(-\frac{\pi}{2},\frac{\pi}{2})$之间的记为辐角主值$\arg z$。

2.欧拉公式

大名鼎鼎的欧拉公式：

$$e^{i t}=\cos t+i\sin t$$

根据三角函数在单位圆上的几何意义，公式是容易理解的。

几何意义：

当中$\theta$为角度，$t$为弧长。

则根据欧拉公式，可将一个复数表示为一个二元组$(a,\theta)$，即模长和辐角（相当于复平面上极坐标系的表示方法）。值为:$a(\cos\theta+i\sin\theta)$

特殊情况：欧拉恒等式$e^{i\pi}+1=0$

3.复数的运算

（1）复数加法

运算规则：实部、虚部分别相加

$$(a+bi)+(c+di)=a+c+bi+di=(a+c)+(b+d)i$$

几何意义：如图

结果相当于两个向量所构成的平行四边形的对角线。如果把一个复数所对应的向量视为一个移动的变换，那么向量加法就是连续运用这两个变换相当于的新变换。

（2）复数乘法

运算规则：展开

$$(a+bi)(c+di)=ac+adi+bci+bdi^2=(ac-bd)+(ad+bc)i$$

几何意义：如图

如图，$\arg a+\arg b=\arg a\times b$,$|a|\times|b|=|a\times b|$

总结就是：模长相乘，辐角相加。

因此，如果模长为$1$，那么它的$n$次方一定还在单位圆上。

证明：

根据欧拉公式，已知$x=(a_1,\theta_1)=a_1(\cos\theta_1+i\sin\theta_1),y=(a_2,\theta_2)=a_2(\cos\theta_2+i\sin\theta_2)$

则

$$\begin{align*} x\times y&=a_1 a_2(\cos\theta_1+i\sin\theta_1)(\cos\theta_2+i\sin\theta_2)\\ &=a_1 a_2\left[(\cos\theta_1\cos\theta_2-\sin\theta_1\sin\theta_2)+i(\cos\theta_1\sin\theta_2+\sin\theta_1\cos\theta_2)\right]\\ &=a_1 a_2\left[\left(\frac{\cos(\theta_1+\theta_2)+\cos(\theta_1-\theta_2)}{2}+\frac{\cos(\theta_1+\theta_2)-\cos(\theta_1-\theta_2)}{2}\right)\right.\\ &\left.+i\left(\frac{\sin(\theta_1+\theta_2)-\sin(\theta_1-\theta_2)}{2}+\frac{\sin(\theta_1+\theta_2)+\sin(\theta_1-\theta_2)}{2}\right)\right]\tag{积化和差公式}\\ &=a_1 a_2\left[\cos(\theta_1+\theta_2)+i\sin(\theta_1+\theta_2)\right] \end{align*}$$

$\therefore |x\times y|=|x|\times |y|,\arg(x\times y)=\arg x+\arg y$

证毕。

4.单位复数根

（1）基本性质

单位复数根是方程$\omega^n=1$的解，第$k$个解记为$\omega_n^k$(这里的$k$事实上是乘方的含义)

$n=16$的解在复平面上的位置如下：

可以看到，$n$个解把单位圆分成了$n$等弧，交点即为根。而且，$\omega_n^k$实际上是$\omega_n$的$n$次方，模长仍为$1$,辐角翻倍！

为什么呢？

$\because |x^n|=|x|^n,\arg x^n=n\arg x$

$\therefore |\omega|^n=|\omega^n|,\arg\omega^n=n\arg\omega$

$\therefore |\omega|^n=1(|\omega|\in\mathbb{R}^+),\arg\omega=\frac{360^\circ}{n}$

$\therefore |\omega|=1,\arg\omega=\frac{360^\circ}{n}$

这就很明显了。

所以，$\omega_n^k$事实上表示的是$\omega_n$的$k$次幂。为什么选择单位复数根呢？因为它有循环的优良性质，即$\omega_n^n=1$。由于其他的都可以由$\omega_n^1$得到，因此称为主$n$次单位根，又记为$\omega_n$。

根据单位复数根的平分圆的意义和欧拉公式，$\omega_n^k=e^\frac{2\pi i k}{n}=\cos\frac{2\pi k}{n}+i\sin\frac{2\pi k}{n}$。

（2）计算引理

显然，由于单位复数根循环($\omega_n^{zn}=e^{2\pi iz}=\left[\left(e^{\pi i}\right)^2\right]^z=1^z=1$)，有变换恒等式：

$$\omega_n^k=\omega_n^{k+wn}(w\in\mathbb{Z})$$

每一份再分成$k$份，编号也变成$k$倍，位置自然不变($\omega_{d n}^{d k}=e^\frac{2\pi i d k}{dn}=e^\frac{2\pi i k}{n}=\omega_n^k$)，所以有消去引理：

ωd