大数乘法（快速傅立叶变换）下

博客搬家：最爱午后红茶 

      大数乘法（快速傅立叶变换）上

        上篇已经已经讲了多项式乘法由系数表示法转化为点值表示法（即求值）的FFT算法的过程；接下来讲插值算法，它需不需要用新的算法写一遍呢？并不用这么麻烦！。

        我们把DFT写成，其中是由（即主单位复根，还记得吗？）的适当幂组成的一个范德蒙德矩阵：

        对 j, k = 0, 1, ... , n-1 ，  的第 k 行，第 j 列（即(k, j)，理解这个很重要！）处的元素值为；

        则求系数向量 a 的运算为 。

        其中这里有个定理：

     

        不要浮躁，仔细看，能看懂的~

        在这里再引入一定理：

        求和引理：对任何整数 n >= 1 和不能被 n 整除的非负整数 k，有：

        证明：

        等比数列求和嘛！给你公式：

      

        之后回来看看这个：

    

        如果 j = j'，则和式的值为1，由求和引理可知在其他情况下和式的值为0.注意，我们依赖-(n - 1) < j' - j < n - 1，以遍 j' - j 不能被 n 整除，这样才能引用求和引理。因此两个矩阵的乘积是单位矩阵。

        之后我们就能通过得出：；j = 0, 1, ..., n - 1

        看看这样式子跟之前求值的时候的求向量 y 的式子有何区别？

        求向量 y 的式子：；j = 0, 1, ..., n - 1