Why Fast Fourier Transform works?

I had this question a while ago, and I can’t seem to quite understand it in any kind of easy way. Let’s first start with Discrete Fourier Transform. It has something to do with Vandermonde matrix but I don’t quite get it. Why DFT is found that way? Namely: \[\hat{f}(x) = \]