Необработанное исключение по адресу 0x00007FFC1E08AFEC (ucrtbased.dll) в lab work 4.exe: Недопустимый параметр был

#include <algorithm>
#include <cstdio>
#include <ctime>
#include <vector>
#include <complex>
#include <iostream>
#define M_PI 3.14159265
 
using namespace std;
 
typedef complex<double> cd;
typedef vector<cd> vcd;
 
vcd traditional(const vcd& as)
{
    int n = as.size();
    if (n == 1) return vcd(1, as[0]);
    vcd w(n);
    for (int i = 0; i < n; i++) 
    {
        double alpha = 2 * M_PI * i / n;
        w[i] = cd(cos(alpha), sin(alpha));
    }
    vcd A(n / 2), B(n / 2);
    for (int i = 0; i < n / 2; i++) 
    {
        A[i] = as[i * 2];
        B[i] = as[i * 2 + 1];
    }
    vcd Av = traditional(A);
    vcd Bv = traditional(B);
    vcd res(n);
    for (int i = 0; i < n; i++)
        res[i] = Av[i % (n / 2)] + w[i] * Bv[i % (n / 2)];
    return res;
}
 
vcd FFT(const vcd& as)
{
    int n = as.size();
    int k = 0; 
    while ((1 << k) < n) k++;
    vector<int> rev(n);
    rev[0] = 0; int high1 = -1;
    for (int i = 1; i < n; i++) {
        if ((i & (i - 1)) == 0) 
            high1++;
        rev[i] = rev[i ^ (1 << high1)]; 
        rev[i] |= (1 << (k - high1 - 1)); 
    }
    vcd roots(n);
    for (int i = 0; i < n; i++) {
        double alpha = 2 * M_PI * i / n;
        roots[i] = cd(cos(alpha), sin(alpha));
    }
    vcd cur(n);
    for (int i = 0; i < n; i++)
        cur[i] = as[rev[i]];
    for (int len = 1; len < n; len <<= 1) {
        vcd ncur(n);
        int rstep = roots.size() / (len * 2);
 
        for (int pdest = 0; pdest < n;) {
            int p1 = pdest;
            for (int i = 0; i < len; i++) {
                cd val = roots[i * rstep] * cur[p1 + len];
                ncur[pdest] = cur[p1] + val;
                ncur[pdest + len] = cur[p1] - val;
                pdest++, p1++;
            }
            pdest += len;
        }
        cur.swap(ncur);
    }
    return cur;
}
 
void combo(int n, int n0, double kof, vcd as, clock_t stime) 
{
    cout << "\n Combined method :\n";
    stime = clock();
    for (int i = 0; i <= 100; i++)
    {
        stime = clock();
        if (n >= n0)
        {
            vcd res = FFT(as);
        }
        else if (n < n0)
        {
            vcd res2 = traditional(as);
        }
        kof += (clock() - stime);
    }
    kof = kof / 100;
    cout << "Time (in ms.) :" << kof << endl;
    kof = 0;
    for (int i = 0; i < n; i++)
    {
        if (i >= n0)
        {
            vcd res = FFT(as);
            printf("%.4lf %.4lf\n", res[i].real(), res[i].imag());
        }
        else if (i <= n0)
        {
            vcd res2 = traditional(as);
            printf("%.4lf %.4lf\n", res2[i].real(), res2[i].imag());
        }
    }
}
 
int main()
{
    int n;
    int n0 = 0;
    double kof = 0;
    cout << "Enter dimension: ";
    cin >> n;
    vcd as(n);
    for (int i = 0; i < n; i++) 
    {
        as[i] = rand() % 10;
        cout << "X" << i + 1 << " = " << as[i] << endl;
    }
    cout << "\n FFT algorithm: \n";
    clock_t stime = clock();
    for (int i = 0; i <= 100; i++)
    {
        stime = clock();
        vcd res = FFT(as);
        kof += (clock() - stime);
    }
    kof = kof / 100;
    cout << "Time (in ms.) :" << kof << endl;
    kof = 0;
    vcd res = FFT(as);
    for (int i = 0; i < n; i++)
        printf("%.4lf %.4lf\n", res[i].real(), res[i].imag());
    cout << "\n Traditional algorithm: \n";
    for (int i = 0; i <= 100; i++)
    {
        stime = clock();
        vcd res2 = traditional(as);
        kof += (clock() - stime);
    }
    kof = kof / 100;
    cout << "Time (in ms.) :" << kof << endl;
    kof = 0;
    vcd res2 = traditional(as);
    for (int i = 0; i < n; i++)
        printf("%.4lf %.4lf\n", res2[i].real(), res2[i].imag());
    combo(n, n0, kof, as, stime);
    return 0;
}

 for (int i = 0; i < n; i++)
    cur[i] = as[rev[i]];