Fit a degree-n polynomial through n+1 points and you are solving a linear system whose matrix is the Vandermonde matrix — rows of powers [1, x, x², …]. Its determinant is a thing of beauty: ∏i<j(xj − xi) — the product of every pairwise gap between the nodes. That single number is load-bearing: it is nonzero if and only if the nodes are distinct, and that is exactly the condition under which the interpolating polynomial exists and is unique. Reconstruction is not merely correlated with distinct nodes; it is caused by the determinant being nonzero. Let two nodes collide and the determinant vanishes and the reconstruction is impossible — there is no polynomial to recover.
The Vandermonde determinant of a set of nodes. The demo computes it for distinct nodes (nonzero — reconstruction possible) and for a repeated node (zero — impossible): live demo
“Any set of points can be fit by a polynomial.” — only if the nodes are distinct, and the Vandermonde determinant is the number that decides it. Reconstruction has a precondition, and this determinant is it. cited
The interpolant exists and is unique because this determinant is nonzero — and it is nonzero because the nodes are distinct. Collapse a node, and there is nothing to recover. Vandermonde (misnomer)
On the canonical compiler, nodes {1,2,4} give det = (2−1)(4−1)(4−2) = 6 — reconstruction possible; nodes {1,2,2} give det = 0 — impossible: