The Discrete Fourier Transform, which sends a vector $x=\left(x_j\right)_{j=0}^{N-1}$ to $y=\mathrm{DFT}(x)$ such that $$y_k=\frac{1}{\sqrt{N}}\sum_{j=0}^{N-1}e^{2\pi i \times jk/N}x_j$$ has a matrix representation $$\mathrm{DFT}_{jk}=e^{2\pi i \times jk/N}=\left(e^{2\pi i /N}\right)^{j\times k},$$ which is in fact a doubly Vandermonde matrix: both it and its transpose are Vandermonde matrices. With this you can use the Vandermonde determinant to prove that $\mathrm{DFT}$ is nonsingular, and if you prove using other means that it is unitary (rather easy) then you will get, I think, a nontrivial expression for 1 as a product of differences of roots of unity.