math.CONov 4, 2026

How many sorted projections determine n points in k dimensions?

Authors: Christopher Gorham Lester, Raymond Hau Sørensen

Organizations: Cavendish Laboratory, Department of Physics, University of Cambridge, UK

Abstract

Let d1,…,dmd_1,\dots,d_m be directions in Rk\mathbb{R}^k, where k≥2k\ge 2. A multiset XX of nn points of Rk\mathbb{R}^k projects to the mm multisets of numbers {{x⋅dj:x∈X}}\{\{x\cdot d_j : x\in X\}\}, or equivalently to the mm sorted lists of these numbers. How many directions are needed for these projections to determine every multiset of nn points? Two simple bounds are known in every dimension: more than (k−1)log⁡2n(k-1)\log_2 n directions are always needed, and (k−1)n+1(k-1)n+1 directions in general position always suffice. In the plane much more is known: Matousek, Privetivy and Skovron showed that almost every choice of O(log⁡nlog⁡log⁡n)O(\log n \log\log n) directions suffices. Whether anything similar holds in higher dimensions was left open by Dym, Wellershoff, Tsoukanis, Levy and Balan. We observe that two lemmas of Hair and Sahai, proved for list recovery of random linear codes, answer that question, and by tightening their argument we show that in every dimension almost every choice of (4k−2)ln⁡(2n)ln⁡(klog⁡22n)+k+1(4k-2)\ln(2n)\ln(k\log_2 2n)+k+1 directions suffices. For k≤10k\le 10 this beats the linear bound once nn is larger than about 100: for example, 144 directions suffice for multisets of 1000 points in the plane, and an exact computation lowers this to 133. Letting some of the random blocks in their argument fail removes a factor log⁡k\log k: almost every choice of O(klog⁡nlog⁡log⁡n)O(k\log n\log\log n) directions suffices, uniformly in kk and nn. We also show that 2k−12k-1 directions never determine three points, so that the logarithmic lower bound is not tight; that every mm directions in the plane fail on some pair of multisets of size O(1.78m)O(1.78^m), improving the bound O(1.80m)O(1.80^m) of Matousek, Privetivy and Skovron; and that for each kk, mm and nn, either almost every choice of mm directions works or almost none does.

Figures & tables

Appendix figures & tables5 assets

Supplementary material from the paper’s appendix.

Appendix

Explore similar work

Aug 21, 2026

Sorting from Counterexamples

Consider the following problem of learning an unknown linear order on nn items. In each round, the learner guesses a complete ordering of the items and receives either confirmation that the guess is correct or a counterexample: a pair of items in the wrong order. The goal is to identify the unknown order using as few queries as possible. We study this problem when up to kk of the returned counterexamples may be untruthful, where kk is not known in advance. We determine the optimal query complexity up to constant factors: Θ(nlog⁡n+nk).Θ(n\log n + nk). Thus, while the noiseless complexity matches the classical complexity of sorting, each untruthful counterexample incurs an additional cost of order nn. The upper bound is based on a geometric representation of permutations and Grünbaum's theorem, while the lower bound combines sorting arguments with a Condorcet-type construction. We also study the case where the target ranking has a low-dimensional geometric representation: each item is represented by a point in Rd\mathbb{R}^d, and the ranking is obtained by projecting the points onto an unknown direction. For these classes we give an upper bound of O(d2log⁡n+dk)O(d^2\log n+dk) and a lower bound of Ω(dlog⁡n+dk)Ω(d\log n+dk), leaving a factor of dd gap in the noiseless term.
Jul 3, 2026

Dimension Reduction for Curves: Simplified and Generalized

We revisit random projections for reducing the dimension of high-dimensional polygonal curves. Drawing from the toolbox of randomized linear algebra, we give a considerably simplified proof of the known O(ε−2log⁡(nm))O(\varepsilon^{-2}\log(nm)) bound on the target dimension of a random projection that preserves the continuous Fréchet distance of polygonal curves up to a factor (1±ε)(1\pm\varepsilon). Our proof is based on the concept of sparse oblivious subspace embeddings. While previous techniques were limited to the case of the Fréchet distance, our techniques are fairly general and extend to all possible distance measures that involve the maximum, a sum or an integral over Euclidean distances between pairs of points on both input curves. We define a generalized dissimilarity measure for curves that includes several popular measures such as Fréchet, qq-DTW, Hausdorff, etc. as special cases and show that the same dimension reduction technique works for this generalized dissimilarity measure. Finally, we apply the same framework for dimension reduction to piecewise linear surfaces, after extending the distance measure suitably to such surfaces.
Jul 9, 2026

High-Dimensional Procrustes Matching via Tree Counts

Suppose we observe two sets of nn Gaussian vectors in Rd\mathbb{R}^d, with the promise that, after applying a permutation of [n][n] and a rotation of Rd\mathbb{R}^d, the two sets are ρρ-correlated. The Procrustes matching problem asks us to recover the unknown permutation of [n][n] that aligns the two sets. The problem is well-studied in the low-dimensional regime d=O(log⁡n)d=O(\log n), but the high-dimensional regime d≫log⁡nd\gg \log n has remained largely uncharted: prior matching guarantees require nearly perfect correlation ρ=1−o(1)ρ=1-o(1), even for information-theoretic recovery. Our main result is a polynomial-time algorithm for exact recovery at constant correlation. The algorithm works by computing and comparing weighted counts of a specially chosen family of ``wide'' trees. So long as d≥polylog(n)d\ge \mathrm{polylog}(n), the algorithm succeeds with high probability for any ρ2>αρ^2>\sqrtα, where α≈0.338α\approx 0.338 is Otter's tree-counting constant. We complement this algorithmic result with an improved information-theoretic guarantee, showing that exact recovery is possible when ρ2≳max⁡{log⁡n/d,log⁡n/n}ρ^2 \gtrsim \max\{\log n/d,\sqrt{\log n/n}\}. We also carry out a low-degree advantage calculation, which suggests that the condition ρ2>αρ^2 > \sqrtα is necessary for any tree-counting algorithm.