std::vector<int> convolution(std::vector<int> a, std::vector<int> b)
{
    a.insert(a.end(), a.size(), 0);
    b.insert(b.end(), b.size(), 0);

    a = forward_dft(a);
    b = forward_dft(b);