快速傅里叶变换（FFT）

贴过来做个模板。

C为复数类，newArr和at为内存池。

 

void FFT(int n, C *x)

{

    C W = C(cos(2 * M_PI / n), sin(2 * M_PI / n)), w = C(1), 

        *u = newArr(n >> 1), *v = newArr(n >> 1);

    for (int i = 0; i < n >> 1; ++i)

        u[i] = x[i << 1], v[i] = x[i << 1 | 1];

    if (n > 1) FFT(n >> 1, u), FFT(n >> 1, v);

    for (int i = 0; i < n >> 1; ++i, w *= W)

        x[i] = u[i] + w * v[i], x[i + (n >> 1)] = u[i] - w * v[i];

    at -= n + 5;

}