快速数论变换（NTT）

前置知识：

阶：对于互质整数 $a,n$ ，满足 $a^r\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}n$ 的最小的 $r$ 就是 $a$ 模 $n$ 的阶。
原根：对于正整数 $n$ ，若整数 $a$ 模 $n$ 的阶等于 $\phi (n)$ ，则称 $a$ 为模 $n$ 的一个原根。

快速数论变换

$NTT$，实际上就是利用一些特殊的质数（例如 $998244353=119\times 2^{23}+1$，这些质数的特征都是可以可以表示 $q\times 2^k+1$ 的形式）进行膜意义下的运算代替浮点复数运算来保证精度，原理是质数原根和复数单位根在DFT运算中具有相同的性质。
下列出我们在 $FFT$ 中利用的单位复根的性质：

1. ${\omega }_n^0={\omega }_n^n=1$
2. 若 $i\ne j\mathrm{m}\mathrm{o}\mathrm{d}n$ ，有 ${\omega }_n^i\ne {\omega }_n^j$
3. ${\omega }_{dn}^{dk}={\omega }_n^k$
4. ${\omega }_n^{k+\frac{n}{2}}=-{\omega }_n^k$
5. ${\sum }_{i=1}^{n-1}({\omega }_n^{j-k}{)}^i=[k==j]n$

对于质数 $p$ ，设 $g$ 是 $p$ 的一个原根，则有 $g^0\equiv g^{\phi (p)}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$ ，故对于 $i\not\equiv j\mathrm{m}\mathrm{o}\mathrm{d}\phi (p)$ ，有 $g^i\not\equiv g^j\mathrm{m}\mathrm{o}\mathrm{d}p$
对于质数 $p=k\times 2^N+1$ ，设 $g_n\equiv g^{\frac{p-1}{n}}\mathrm{m}\mathrm{o}\mathrm{d}p$ ，这个 $g_n$ 同样拥有着上述单位复根的 $5$ 条重要性质。
我们之所以取模数 $p=k\times 2^N+1$ ，是因为我们要保证 $\frac{p-1}{n}$ 这个指数是一个整数，由于 $FFT$ 的过程是一个不断把区间大小缩小一半来处理的过程，所以我们可以保证 $n$ 总是 $2$ 的幂次，所以模数必须足够多的 $2$ ，所以 $p$ 需要取成 $k\times 2^N+1$ 而且 $N$ 需要尽可能的大一些
所以在写 $NTT$ 时，只需要把 $FFT$ 中的单位复根换成原根即可
具体见代码：

int n,m,len = 2,L,rev[maxn<<1],a[maxn<<1],b[maxn<<1],c[maxn<<1];
const int mod = 998244353;//mod原根为3 
const int g = 3;//原根 
ll qpow(ll a,ll b)
{
	ll res = 1;
	while(b)
	{
		if(b&1) res = res*a%mod;
		a = a*a%mod;
		b >>= 1;
	}
	return res;
}
void get_rev()
{
	for(int i = 0;i < len;i++)
		rev[i] = (rev[i>>1]>>1)|(len>>1)*(i&1);
//	for(int i = 0;i < len;i++) rev[i] = (rev[i>>1]>>1)|((i&1)<<(L-1));
}

void ntt(int *a,int dft)
{
	for(int i = 0;i < len;i++)
		if(i < rev[i]) swap(a[i],a[rev[i]]);//不加这条if会交换两次（就是没交换）
	for(int mid = 1;mid < len;mid <<= 1)//mid是准备合并序列的长度的二分之一
	{
		ll W = qpow(g,(mod-1)/(mid<<1));
		if(dft == -1) W = qpow(W,mod-2);
		for(int i = 0;i < len;i += mid<<1)//mid*2是准备合并序列的长度，i是合并到了哪一位
		{
			ll w = 1;
			for(int j = i;j < mid+i;j++,w = w*W%mod)//只扫左半部分，得到右半部分的答案
			{
				int x = a[j];
				int y = w*a[j+mid]%mod;
				a[j] = (x+y)%mod;
				a[j+mid] = ((x-y)%mod+mod)%mod;
			}
		}
	}
	if(dft == -1) 
	{
		int tmp = qpow(len,mod-2);
		for(int i = 0;i < len;i++) a[i] = (ll)a[i]*tmp%mod;
	}
}

任意模数的快速数论变换 $MTT$ （待更）