快速数论变换（NTT）

开头

NTT就是对FFT进行了改进，上一篇中讲了FFT，可以看到FFT所用的单位根 $w_n^k$ 的实部和虚部都是通过正弦和余弦函数计算而来的，所以不可避免地会有很多浮点数运算

所以NTT就是在整数范围内寻找和单位根有相同性质的那些数，可以提升计算的精度

这种整数被找到了，就是原根
        

欧拉函数

对于一个正整数 n，小于等于 n 且与 n 互质的正整数的个数，记作 $\phi (n)$，有
$$\phi (x)=x(1-\frac{1}{P(1)})(1-\frac{1}{P(2)})......(1-\frac{1}{P(n)})$$其中 $P(k)$ 为 $x$ 的质因数，且相同的 $P(k)$ 只会出现一次

如：
$$\begin{array}{rl}\phi (99)& =99\times (1-\frac{1}{3})\times (1-\frac{1}{11})\\ & =60\\ \phi (15)& =99\times (1-\frac{1}{3})\times (1-\frac{1}{5})\\ & =8\end{array}$$

阶

定义：假设 $m>1$，且 $gcd(a,m)=1$ （$a$、$m$ 互质），那么使得 $a^r\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}m$ 成立的最小的整数 $r$，称为 $a$ 对模 $m$ 的阶，记作
$${\delta }_m(a)$$

原根

如果有
$${\delta }_m(a)=\phi (m)$$则称 $a$ 为模 $m$ 的一个原根

---

举个例子：

设 $m=7$，则 $\phi (m)=7\times (1-\frac{1}{7})=6$

一个一个取 $a$

$$\begin{array}{rl}a=2，& 2^3=8\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}7\\ & 得到{\delta }_7(2)=3\\ & \because 3<6，\therefore {\delta }_7(2)̸=\phi (7)\\ & 2不是模7的原根\\ a=3，& 3^k̸=1\mathrm{m}\mathrm{o}\mathrm{d}7,k\in [1,5]\\ & 3^6=729\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}7\\ & 得到{\delta }_7(3)=6\\ & \because 6=6，\therefore {\delta }_7(3)=\phi (7)\\ & 3是模7的原根\end{array}$$
原根就具有单位根的性质，使得它可以和单位根一样进行运算

原根存在的条件是：
$m=2,4,p^n,2p^n$，这里 $p$ 为奇素数， $n$ 为正整数
        

NNT

这里我们取
$$g^{\frac{P-1}{n}}\mathrm{m}\mathrm{o}\mathrm{d}P=w_n$$其中 $g$ 是 $P$ 的原根，这样就可以用左式代替右式进行类FFT计算，即用整数代替了浮点数

这里 $P$ 通常取998244353、1004535809，它们的原根 $g$ 均为 3

那NNT就是将FFT中的单位根换成了原根

如果说FFT的核心代码如下的话：

x = a[j + k];
y = w * a[j + k + mid];
a[j + k] = x + y;
a[j + k + mid] = x - y;

那NNT的核心代码就是这样：

x = a[j + k];
y = w * a[j + k + mid];
a[j + k] = (x + y) % P;
a[j + k + mid] = (x - y + P) % P;