[学习笔记][省选算法]多项式乘法之 FFT & NTT 及卷积

一、开头

（Place ：山东省神犇协会第 998244353 会议厅）
SD 神犇 998244353 号（会长）：所有神犇，洛谷黑题都切完了吗？
全体神犇：切完了！
神犇 1004535809 号： 998244353 ，我昨天看到一个弱鸡，洛谷名叫 xyz32768 ，他一道黑题都没做，您认为应该怎样呢？
神犇 998244353 号：这个问题您们考虑一下，应该用怎样的方式把 xyz32768 D 一遍。时刻记住，我们协会的口号：见一弱鸡， D 一弱鸡！
神犇 167772161 号：要不，考他 FFT 吧。这是我们协会里每个人都会的算法！
神犇 998244353 号：Good idea ！
（Place ： FFT 星球 DFT 国 IDFT 省 NTT 市）
神犇 167772161 号：哈哈，好久没见，我们搞 OI 多年，能回忆起的比赛，至少也从三年前的 SDOI 2015 开始吧，那次比赛有一道十分模板的题。题目名称：序列统计。题面你自己去 BZOJ 或洛谷上去查，如果你 1h 内写不出来我会对整个 FFT 星球说你菜！
xyz32768 ：BZOJ 3992 / 洛谷 P3321 ？
神犇 1004535809 号：对。这是一道紫题，但此大佬已经将之视为宇宙之水题了。
（蒟蒻 xyz32768 当然不会，看了看算法标签）
xyz32768 ： 什么？？？？？？ 快速傅里叶变换,DFT,FFT ？？？？？？
神犇 998244353 号：哈哈哈哈哈，你真是菜啊， FFT 这么水的算法都不会，再见蒟蒻！
神犇 998244353 号 & 神犇 1004535809 号 & 神犇 167772161 号（大声喊）： xyz32768 菜， xyz32768 好菜， xyz32768 最菜！

二、预备知识

1、复数

下面的 $i$ ，除非作为 $\sum$ 求和的变量，其余都表示虚数单位 $\sqrt {-1}$ 。
一个复数可以表示成 $a+bi$ 的形式，也可以表示成 $r(\cos\theta+i\sin\theta)$ 的形式。
上面 $a$ 为实部， $b$ 为虚部， $r=\sqrt{a^2+b^2}$ 为模长， $\theta=\arctan\frac b a$ 为辐角。
C++ 中，

#include <bits/stdc++.h>
using namespace std;

复数类为（ T 为实部和虚部的数据类型）：

complex<T>

定义一个复数（ $a$ 为实部， $b$ 为虚部）

complex <T> x(a, b)

实部

x.real()

一个复数 $a+bi$ 对应复平面上一个点 $(a,b)$ ，模长 $r$ 为点到原点的距离，辐角 $\theta$ 为终边为 $(0,0)-(a,b)$ 的角（即该边与 x 轴正半轴形成的角）。

2、单位根

$$\omega_n=\cos\frac{2\pi}n+i\sin\frac{2\pi}n$$
那么
$$\omega_n^0,\omega_n^1,\omega_n^2,...,\omega_n^{n-1}$$ 称为 $n$ 次单位根。
如：
$$\omega_4^0=1,\omega_4^1=i,\omega_4^2=-1,\omega_4^3=-i$$
$n$ 次单位根的性质 ：
（1） $n$ 次单位根所表示的点分布在复平面的单位圆上，并且把圆周等分成 $n$ 段。
（2） $n$ 次单位根的 $n$ 次方等于 $1$ 。
（3） $$\omega_{2n}^{2i}=\omega_n^i,\omega_{2n}^n=-1,\omega_{2n}^{n+i}=-\omega_{2n}^i$$

3、原根

在模 $p$ （为质数）意义下，如果存在数 $g$ ，满足

$$g^0,g^1,g^2,...,g^{p-2}$$
互不相同，那么 $g$ 为模 $p$ 的原根。
原根都比较小，可以暴力枚举。将 $p-1$ 分解质因数，对于一个 $g$ ，如果对于 $p-1$ 的任何一个质因子 $x$ ，都有
$$g^{\frac{p-1}x}\ne1$$
那么 $g$ 是模 $p$ 的一个原根。

三、多项式乘法

步入正题。一个 $n$ 次多项式 $A$ 和一个 $n$ 次多项式 $B$ 相乘，设 $A$ 的系数为 $a_{0...n}$ （ $x^i$ 的项系数为 $a_i$ ）， $B$ 的系数为 $b_{0...n}$ ，那么

$$ab_k=\sum_{i+j=k}a_ib_j$$
复杂度显然是 $O(n^2)$ 的。快速傅里叶变换（ FFT ）则能做到 $O(n\log n)$ 的复杂度。

四、多项式的点值表示及 DFT

多项式的点值表示：选取 $n+1$ 个互不相同的值 $x_{0...n}$ 代入 $A$ 对多项式进行求值，得到 $n+1$ 个多项式值，那么这 $n+1$ 个自变量和这 n+1