(EDIT: Also, there is a qualitative version of the above identity: just from knowing that the Vandermonde determinant is non-vanishing when the $x_i$ are distinct, one can already deduce that polynomial interpolation is theoretically possible, though to get the precise formula one still needs to go through the above identity.)