We present a novel analogue for finite exchangeable sequences of the de Finetti, Hewitt and Savage theorem and investigate its implications for multi-marginal optimal transport (MMOT) and Bayesian statistics. If $(Z_1,...,Z_N)$ is a finitely exchangeable sequence of $N$ random variables taking values in some Polish space $X$, we show that the law $μ_k$ of the first $k$ components has a representation of the form $μ_k=\int_{{\mathcal P}_{\frac{1}{N}}(X)} F_{N,k}(λ) \, \mbox{d} α(λ)$ for some probability measure $α$ on the set of $1/N$-quantized probability measures on $X$ and certain universal polynomials $F_{N,k}$. The latter consist of a leading term $N^{k-1}\! /{\small \prod_{j=1}^{k-1}(N\! -\! j)\, λ^{\otimes k}}$ and a finite, exponentially decaying series of correlated corrections of order $N^{-j}$ ($j=1,...,k$). The $F_{N,k}(λ)$ are precisely the extremal such laws, expressed via an explicit polynomial formula in terms of their one-point marginals $λ$. Applications include novel approximations of MMOT via polynomial convexification and the identification of the remainder which is estimated in the celebrated error bound of Diaconis-Freedman between finite and infinite exchangeable laws.