Faster Algorithms for Multimarginal Optimal Transport
2026-08-10 • Data Structures and Algorithms
Data Structures and Algorithms
AI summaryⓘ
The authors study how to approximate the multimarginal optimal transport (MOT) distance, which measures how to optimally move mass between multiple distributions. They provide a classical algorithm that improves previous methods by balancing accuracy and problem size better than before. They also propose two quantum algorithms that run faster for large problem dimensions but have worse accuracy trade-offs. Additionally, they establish lower bounds showing the minimum number of queries any algorithm must make to achieve a certain precision. Their work advances understanding of both classical and quantum approaches to solving MOT problems.
multimarginal optimal transportoptimal transport distanceclassical algorithmsquantum algorithmslinear programming dualSinkhorn algorithmentropy regularizationadditive approximationquery complexitysubgradient method
Authors
Brandon Augustino, Yue Sun, Atithi Acharya, Shouvanik Chakrabarti, Junhyung Lyle Kim, Shree Hari Sureshbabu, Charlie Che
Abstract
We study algorithms for approximating the multimarginal optimal transport (MOT) distance, a generalization of the classic optimal transport distance, between $m$ discrete probability distributions each supported on at most $n$ points. We give a classical algorithm that computes a coupling between these marginals whose expected transportation cost is within an additive $\varepsilon > 0$ of the MOT distance in time $O(m^2 n^m \varepsilon^{-1}\mathrm{polylog}(m,n,\varepsilon^{-1}))$. This is, to our knowledge, the first bound for general MOT problems with simultaneous linear dependence on the dimension $n^m$ and on the accuracy parameter $\varepsilon^{-1}$, improving the prior state of the art. On the quantum side, we give two algorithms that achieve speedups in dimension, though with worse accuracy dependence than classical approaches. First, we construct a quantum projected subgradient method for estimating the MOT distance within an additive $\varepsilon >0$ with runtime $O( m^3 n^{\frac{m}{2}+1} \varepsilon^{-2} \mathrm{polylog}(m,n,\varepsilon^{-1}))$. This algorithm works with the linear programming dual of the MOT problem, and does not return a coupling. We also give a quantum multimarginal Sinkhorn algorithm for entropy-regularized MOT. This algorithm returns an implicit description of an approximately optimal coupling with runtime $O(m^8n^{\frac{m+1}{2}} \varepsilon^{-5} \mathrm{polylog}(m,n,\varepsilon^{-1})))$ after the usual reduction from entropic MOT to unregularized MOT. We also record query lower bounds: for any precision $\varepsilon<1/2$, randomized classical algorithms require $Ω(n^m/(1+\varepsilon n))$ queries and quantum algorithms require $Ω(\sqrt{n^m/(1+\varepsilon n)})$ queries.