cs.DSMay 24, 2026

Algorithms with Polynomially-Improved Approximation Factors for the 2 \rightarrow q Norm, and Applications

Authors: Samuel B. HopkinsStefan Tiegel

Abstract

The 2q2 \rightarrow q norm of a matrix XRn×dX \in \mathbb{R}^{n \times d} is defined as X2q=supv2=1Xvq\lVert X \rVert_{2 \rightarrow q} = \sup_{\lVert v \rVert_2 = 1} \lVert Xv \rVert_q. We give polynomial-time multiplicative approximation algorithms for this norm when q>2q > 2 (i.e. in the hypercontractive setting). This problem either directly captures or is closely related to long-standing open problems in combinatorial optimization and hardness of approximation (e.g. Small Set Expansion), quantum information (e.g. Best Separable State), and algorithmic statistics. Very little is known about what approximation factors we can achieve for this problem in polynomial time, even though such approximations have significant downstream consequences. Barak, Brandão, Harrow, Kelner, Steurer, and Zhou showed that no polynomial-time algorithm can achieve an approximation factor better than 2logn2^{\sqrt{\log n}}, assuming the Exponential Time Hypothesis (FOCS'12). On the other hand, a simple spectral algorithm gives a d1/4d^{1/4}-approximation as a baseline. We give, to the best of our knowledge, the first polynomial-time approximation algorithm beating this baseline by polynomial factors. For the important special case of q=4q = 4 it achieves a d1/8d^{1/8}-approximation. All previous algorithms required additional assumptions on XX, or only surpassed the baseline for small values of nn. Moreover, we construct sum-of-squares certificates for the 2q2 \rightarrow q norm. This directly implies improved algorithms for robust mean and covariance estimation, robust regression, and clustering, when the data only satisfies a bound on its qq-th moment.

Explore similar work

Aug 12, 2026cs.DS

Fast Length-Squared Sampling for Positive-Semidefinite Matrices

We describe a simple rejection-sampling-based algorithm to perform length-squared sampling on an n×nn \times n positive-semidefinite (psd) matrix: that is, to sample a column with probability proportional to its squared 2\ell_2-norm. The algorithm runs in just O(n)O(n) expected time, which is significantly sublinear in the input matrix size. The runtime is optimal, even when the input is assumed to be diagonal. Our result has several applications. Length-squared sampling is used by a number of sublinear time algorithms for matrix problems, like low-rank approximation and eigenvalue approximation. Often, it is assumed that the algorithm is given access to the matrix column norms, and thus can perform length-squared sampling efficiently. Our result shows that, at least for psd matrices, we can remove this assumption. We also discuss an application to an asymptotically optimal algorithm for estimating the Frobenius norm of a psd matrix to relative error. Finally, we show that our sampling algorithm yields a very simple sublinear time algorithm for the robust psd low-rank approximation problem introduced by Bakshi et al. (FOCS, 2020), which nearly matches the more complex method developed there.
Rajarshi Bhattacharjee, Ethan N. Epperly, Cameron Musco +1
May 20, 2025cs.LG

A Private Approximation of the 2nd-Moment Matrix of Any Subsamplable Input

We study the problem of differentially private second moment estimation and present a new algorithm that achieve strong privacy-utility trade-offs even for worst-case inputs under subsamplability assumptions on the data. We call an input (m,α,β)(m,α,β)-subsamplable if a random subsample of size mm (or larger) preserves w.p 1β\geq 1-β the spectral structure of the original second moment matrix up to a multiplicative factor of 1±α1\pm α. Building upon subsamplability, we give a recursive algorithmic framework similar to Kamath et al 2019, that abides zero-Concentrated Differential Privacy (zCDP) while preserving w.h.p. the accuracy of the second moment estimation upto an arbitrary factor of (1±γ)(1\pmγ). We then show how to apply our algorithm to approximate the second moment matrix of a distribution D\mathcal{D}, even when a noticeable fraction of the input are outliers.
Bar Mahpud, Or Sheffet
Sep 27, 2024math.NA

Probabilistic Analysis of Least Squares, Orthogonal Projection, and QR Factorization Algorithms Subject to Gaussian Noise

We consider the effect of Gaussian perturbations on least-squares residuals, orthogonal projections, and QR-type algorithms. The problem that motivated our investigations is as follows: suppose that a full column-rank matrix BRm×nB\in\mathbb{R}^{m\times n} has already been computed, and suppose that a new normalized column q=(x+y)/x+y2q=(x+y)/\|x+y\|_2 is to be appended to BB, where xspan(B)x\perp\operatorname{span}(B) is the ideal orthogonal component and yy represents the orthogonalization error. How large can the condition number κ([B,q])κ([B,q]) of the resulting matrix [B,q][B,q] become? While we provide a Weyl-type bound on the singular values of [B,q][B,q], in terms of the extremal singular values of BB and the quantity BTy2/x+y2\|B^T y\|_2/\|x+y\|_2, we also derive exact probability laws for norms and projection residuals under Gaussian perturbations. Finally, we use these probability laws to derive probabilistic condition-number bounds for QR-type processes with imperfect orthogonalization and exact normalization.
Ali Lotfi, Julien Langou, Mohammad Meysami