fix(math/poly/fft): 修改递归版 FFT 代码实现

docs/math/poly/fft.md

@@ -247,16 +247,13 @@ $$
       Comp *g = f, *h = f + n / 2;
       // 递归 DFT
       DFT(g, n / 2, rev), DFT(h, n / 2, rev);
-      // cur 是当前单位复根，对于 k = 0 而言，它对应的单位复根 omega^0_n = 1。
-      // step 是两个单位复根的差，即满足 omega^k_n = step*omega^{k-1}*n，
-      // 定义等价于 exp(I*(2*M_PI/n*rev))
-      Comp cur(1, 0), step(cos(2 * M_PI / n), sin(2 * M_PI * rev / n));
       for (int k = 0; k < n / 2;
            ++k) {  // F(omega^k_n) = G(omega^k*{n/2}) + omega^k*n\*H(omega^k*{n/2})
+        // cur 是当前单位复根，定义等价于 exp(I*(2*M_PI/n*rev))
+        Comp cur(cos(2 * k * M_PI / n), sin(2 * k * M_PI * rev / n));
         tmp[k] = g[k] + cur * h[k];
         // F(omega^{k+n/2}*n) = G(omega^k*{n/2}) - omega^k_n*H(omega^k\_{n/2})
         tmp[k + n / 2] = g[k] - cur * h[k];
-        cur *= step;
       }
       for (int i = 0; i < n; ++i) f[i] = tmp[i];
     }