,
Sharath Raghvendra,
Keegan Yao
Creative Commons Attribution 4.0 International license
We give the first (relative) (1+ε)-approximation algorithm for the continuous optimal transport (OT) problem between two axis-aligned histograms in the plane, each of which is represented as a piecewise-constant function over a rectangular subdivision. Our algorithm runs in nearly quadratic time in the complexity of the histograms in the worst case. In contrast to the discrete and semi-discrete optimal transport problems, which always admit linear-size OT plans, we demonstrate that there is a quadratic lower bound on the complexity of even a (1+ε)-approximate OT plan in the continuous setting in the worst case. This suggests that the runtime of our approximation algorithm nearly matches the worst-case lower-bound complexity of an explicit (1+ε)-approximate OT plan between two histograms. We additionally provide near-linear time algorithms with either weaker approximation guarantees or restrictions on the input histograms.
@InProceedings{agarwal_et_al:LIPIcs.ESA.2026.123,
author = {Agarwal, Pankaj K. and Raghvendra, Sharath and Yao, Keegan},
title = {{Fast Algorithms for Continuous Optimal Transport Between Histograms}},
booktitle = {34th Annual European Symposium on Algorithms (ESA 2026)},
pages = {123:1--123:23},
series = {Leibniz International Proceedings in Informatics (LIPIcs)},
ISBN = {978-3-95977-445-1},
ISSN = {1868-8969},
year = {2026},
volume = {388},
editor = {Bille, Philip and Pettie, Seth and Storandt, Sabine},
publisher = {Schloss Dagstuhl -- Leibniz-Zentrum f{\"u}r Informatik},
address = {Dagstuhl, Germany},
URL = {https://drops.dagstuhl.de/entities/document/10.4230/LIPIcs.ESA.2026.123},
URN = {urn:nbn:de:0030-drops-272595},
doi = {10.4230/LIPIcs.ESA.2026.123},
annote = {Keywords: Optimal transport, min-cost flow, minimum-weight matching, quadtrees, approximation algorithms, dynamic optimal transport algorithm}
}