CountSketch Input-Sparsity Subspace Embedding and Regression
Algorithms AIMech:CountSketchInputSparsitySubspaceEmbeddingAndRegression DRAFT
Aliases: Clarkson-Woodruff Input-Sparsity Transform (CW Sketch), Sparse CountSketch Oblivious Subspace Embedding, Input-Sparsity Time Low-Rank Approximation and Regression
CountSketch Input-Sparsity Subspace Embedding and Regression is an oblivious sparse random projection and numerical linear algebra framework introduced by Clarkson and Woodruff (2013, 2017) that constructs an O(n)-column, t-row sketching matrix S = Phi D in R^{t x n} with exactly ONE nonzero entry per column (S_{h(i), i} = sigma_i in {-1, +1}, where h: [n] -> [t] and sigma: [n] -> {-1, +1} are random hash and sign maps) so that S A can be computed in true input-sparsity time O(nnz(A)) for any matrix A in R^{n x d}, where nnz(A) is the number of nonzero entries of A. Theorem 2.3 and Lemma 2.5 prove that for t = O(d^2 / epsilon^2) (or t = O(k^2 / epsilon^2) for a k-dimensional subspace), S is a (1 +/- epsilon)-l_2 subspace embedding satisfying ||S A x||_2 = (1 +/- epsilon)||A x||_2 simultaneously for all x in R^d by decomposing the basis U of range(A) into heavy rows (||u_i||_2^2 >= tau, which are isolated into distinct hash buckets without collision and preserved isometrically) and light rows (bounded via high-moment trace inequalities), yielding O(nnz(A) + d^3 epsilon^{-2} log(d/epsilon)) algorithms for (1+epsilon)-approximate overconstrained least-squares regression (Theorem 2.6) and O(nnz(A) + n poly(k/epsilon)) low-rank matrix approximation.
Identity Notes & Concept Scope
Distinct from Subsampled Randomized Hadamard Transform Subspace Embedding (`AIMech:SubsampledRandomizedHadamardTransformSubspaceEmbedding`), which requires O(n d log t) dense Walsh-Hadamard butterfly operations regardless of matrix sparsity, and from Frequent Directions Deterministic Streaming Matrix Sketching (`AIMech:FrequentDirectionsDeterministicStreamingMatrixSketching`), which uses O(n d l) deterministic SVD shrinkage. Specifically captures the Clarkson-Woodruff single-nonzero-per-column CountSketch matrix S = Phi D, the heavy/light leverage row decomposition (Lemma 2.5), and the O(nnz(A)) input-sparsity time algorithms for subspace embedding, least-squares regression, and rank-k approximation.
Claims & Grounded Evidence
3 Claims
Mechanism
#claim-01
For the sparse embedding matrix S = Phi D in R^{t x n} with a single uniformly hashed bucket h(i) in [t] and random Rademacher sign sigma_i in {-1, +1} per column i, the matrix product S A is computable in O(nnz(A)) arithmetic operations, satisfies the exact unbiased second-moment identity E[(S^T S)_{ij}] = delta_{ij} (so E[||S x||_2^2] = ||x||_2^2 and E[||S A||_F^2] = ||A||_F^2), and preserves the Euclidean norm of any subset of heavy rows isolated into distinct hash buckets with exact isometry.
Scope
#claim-02
Theorem 2.3 and Theorem 2.6 prove that for t = O(d^2 / epsilon^2) rows, S = Phi D is a (1 +/- epsilon)-l_2 subspace embedding for any d-dimensional subspace V in R^n with constant probability, and combining CountSketch preconditioning with iterative refinement solves min_x ||A x - b||_2 to (1 + epsilon) relative residual error in O(nnz(A) log(1/epsilon) + d^3 log^2 d + d^2 log(1/epsilon)) time.
Limitation
#claim-03
Because S = Phi D hashes each row of A to a single bucket without prior Hadamard spreading, avoiding collisions among up to d high-leverage (spiky canonical basis) rows by the birthday paradox requires target dimension t = Omega(d^2 / epsilon^2) rather than O(d / epsilon^2), and hash collisions between two heavy basis rows with opposite effective signs cause destructive cancellation that breaks the subspace embedding property.
⚠️ Curated Limitations & Negative Results
- Single-pass sparse CountSketch subspace embedding requires O(d^2 / epsilon^2) target rows (quadratic in subspace dimension d due to birthday-paradox collision avoidance on high-leverage rows) unless composed in a two-stage sketch with SRHT.
Artificial Intelligence Ontology (AIO) Grounding
Mapped AIO Classes
| AIO Target Class | Mapping Relation | Mapping Rationale |
|---|---|---|
AIO:DimensionalityReduction
Dimensionality Reduction
|
narrow_match | Sparse O(nnz(A)) oblivious subspace embedding S = Phi D reducing n rows to t = O(d^2/epsilon^2) rows. |
AIO:LinearRegression
Linear Regression
|
related_match | Yields input-sparsity time (1+epsilon)-approximate overconstrained least-squares regression and preconditioned iterative solvers. |
AIO:PrincipalComponentAnalysis
Principal Component Analysis
|
related_match | Provides input-sparsity time O(nnz(A) + n poly(k/epsilon)) low-rank matrix approximation. |
Evidence Sources & Primary Literature
3 Sources| Source ID | Publication / Document | Type / Access | Reference & Checksum |
|---|---|---|---|
r331-primary-pdf
|
Kenneth L. Clarkson, David P. Woodruff
|
primary paper
full_text
|
ID: https://doi.org/10.1145/2488608.2488620
sha256:8005ec41c281…
|
r331-arxiv-abs
|
Kenneth L. Clarkson, David P. Woodruff
|
primary paper
abstract_only
|
ID: https://arxiv.org/abs/1207.6365
sha256:5b39d1b4f68d…
|
r331-metadata-cr
|
Kenneth L. Clarkson, David P. Woodruff
|
primary paper
abstract_only
|
ID: https://doi.org/10.1145/2488608.2488620
sha256:2ca807efd904…
|
Ontology & Grounding Topology
Open in interactive network graph →Curation Provenance
Curator: Antigravity
Method: ai_assisted
Model: Gemini
Created: 2026-10-19T06:54:00Z
Updated: 2026-10-19T06:54:00Z
Source File:
data/curated/CountSketchInputSparsitySubspaceEmbeddingAndRegression.yaml
· Schema:
curated_concept.yaml
Explore across X-Mech
Related computational & biological Mechs:
Embedding Map
Network Graph
ProteinTraitsMech (ESM-2 pLM)
DUFMech (Domain pLM)
TraitMech
CultureMech
All 14 X-Mech Repositories