拆系数FFT(任意模数FFT)

拆系数FFT

对于任意模数 $mod$
设$m=\sqrt{mod}$
把多项式$A(x)$和$B(x)$的系数都拆成$a\times m+b$的形式，时$a,b$都小于$m$
提出，那么一个多项式就可以拆成两个多项式的加法
一个是$a\ast m$的，一个是$b$的
直接乘法分配律，$aa$一遍，$ab$一遍，$ba$，$bb$一遍，四遍$FFT$
乘出来不会超过取模范围
然后合并直接
$$(a\times m+b)(c\times m+d)=a\times c\times m^2+(a\times c+b\times d)m+b\times d$$
这样子要进行 $7$ 次 $DFT$

如果研究一下 $myy$ $2016$ 年的集训队论文就会发现有 $2$ 次 或者 $1.5$ 次 $DFT$ 的 $FFT$ 算法
2次的够了吧
$myy$ 巧妙的运用了复数的虚部，优化了算法
具体来说

设
$$C(x)=A(x)+iB(x)$$
$$D(x)=A(x)-iB(x)$$

假设
$c(w_n^k)$ 表示将 $C(x)$ $DFT$ 后的点值
$d(w_n^k)$ 表示将 $D(x)$ $DFT$ 后的点值
$w$ 为 $n$ 次单位复数根
设 $conj(x)$ 表示 $x$ 的共轭复数
那么
c ( w 2 n k ) = A ( w 2 n k ) + i B ( w 2 n k ) = ∑ j = 0 2 n − 1 A j w 2 n j k