1 Introduction
A linear projection can concentrate a measure even when the original measure has no atoms. The amount of concentration depends on how much mass can lie near the affine subspaces that a projection collapses. This paper studies that dependence for random exact projections and uses it to analyze what a finite memory state can retain from noiseless Gaussian equations.
We consider three questions. The first asks for an integrated moment of a projected density. Let \nu be a finite positive measure in \R^d whose mass on every ball of radius r is at most Kr^\beta, and let P be a random matrix with k orthonormal rows. When P_\#\nu has density g_P, an integrated estimate controls the average over P of \int g_P^q. Such an estimate measures concentration throughout the label space. Theorem 3.3 proves a bound with explicit dependence on the dimension, support radius, and local-mass exponent.
The second question asks for an actual positive dominating measure after a finite routing decision. Let \sigma be uniform probability on S^{d-1}, and let \mu_B be \sigma conditioned on a spherical cap B. A routing rule observes a Gaussian matrix X and the exact label Xs, and selects an output w with probability K_w(X,Xs). Its average likelihood f_w(s)=\E_XK_w(X,Xs) weights the parent measure to f_w\mu_B. We seek a positive sum of normalized cap measures that dominates this measure, with a controlled cost for small child radii. The positive-domination theorem in Section 4 constructs a countable such sum and bounds its cost by a power of a reference routing probability.
The third question concerns evaluation at a specified label. For a measurable function h:S^{d-1}\to[0,1], the ambient Riesz functional
is the supremum of the raw Riesz potentials of h\sigma over ambient centers, and it controls mass near every center. If h\sigma is restricted to a ball about z, Section 5 estimates a Gaussian projected density at Xz+b for every fixed offset b. It uses one explicitly chosen density version for all offsets. This qualification matters: an integrated density norm is unchanged by modifications on a matrix-dependent graph, so it cannot by itself assign or bound the values on that graph.
All three estimates give a precision bound for finite-state regression. The signal is uniform on S^{d-1}, and each fresh observation is a standard Gaussian feature together with its exact inner product with the signal. Write T for the learner’s deterministic sample horizon. In the bit-memory model of Section 2, a learner with M(d)=o(d^2) retained bits that achieves angular error at most \epsilon(d), where 0<\epsilon(d)\le1/10, with average probability at least 2/3 requires
for an absolute c>0 and all sufficiently large d, as stated in Theorem 2.4. The dimension threshold may depend on the memory sequence. Under the same accuracy and success requirements, the integrated projection argument also gives the explicit tradeoff
for every nonnegative integer M and all d above an absolute threshold (Corollary 2.5). In particular, the same precision bound holds with absolute constants throughout M\le d^2. The underlying success-probability bounds (17), (30), and Proposition 5.4 retain explicit dependence on arbitrary M.
1.1 Projection geometry and earlier methods
The local growth condition belongs to classical potential theory. Mattila’s work on Hausdorff dimension and orthogonal projections uses inverse-distance energy averaging and compares positive-dimensional Hausdorff measure with the existence of a measure satisfying a ball-growth bound [8, Lemma 5.1, p. 236, and Theorem 5.4, p. 237]. Mattila describes the potential-theoretic projection method as following Kaufman’s work [8, Introduction, p. 227]. Here the growth measure is supplied by a hypothesis or constructed from a cap test. The quantitative estimates require successive distances to affine spans and moments whose order grows with d.
Affine simplex volumes also have an established role in integral geometry. Drury’s affine-plane change of variables contains a simplex-volume Jacobian [3, Lemma 1, pp. 497–498]. In our frame calculation, the product of successive affine heights is the volume of the associated parallelepiped, and the k projected rows raise its inverse to the kth power. We derive that row-Jacobian and its probability normalization directly. The exact density of a rectangular block of a Haar orthogonal matrix is classical; see Jiang [7, Lemma 2.5]. Our derivation keeps its normalization explicit before the affine-height change of variables. Classical beta descriptions of a uniform direction projected onto a subspace are discussed by Frankl and Maehara [4]; Lemma 2.1 gives the vector density and conditional direction law used here.
The regression application sits within a broader study of learning with limited retained information. Shamir’s finite-message protocols include a bounded-memory online subclass [11, Section 2]. Steinhardt, Valiant, and Wager related memory, communication, and statistical queries [14, Section 1.1]; Raz subsequently proved a qualitative quadratic-memory versus exponential-sample separation for parity learning [9, Theorem 1]. For regression, Steinhardt and Duchi established memory-dependent minimax risk bounds in a sparse noisy model [13, Section 1.1].
Sharan, Sidford, and Valiant developed a nonuniform finite-width branching-program analysis for continuous regression. Their principal theorem concerns independent isotropic Gaussian features with independent uniform additive noise of half-width 2^{-d/5}. With at most d^2/4 bits of memory, success probability at least 2/3 for a uniformly distributed unit signal, and Euclidean accuracy d^{-r} in its stated range r\le O(d/\log d), it gives an \Omega(d\log r) sample bound [12, Theorem 1]. Their projection argument expands high moments into independent copies and uses successive orthogonalization [12, Section 7, Lemma 12]. That is a close methodological predecessor of our tuple calculations. Their lemma treats one noisy interval under a global L^2 density condition. The estimates here treat exact multirow labels under local ball or Riesz-potential conditions. Dagan, Kur, and Shamir obtained related quadratic-space obstructions for approximate null vectors of Gaussian rows and for empirical residual error in a different equation model [2, Theorems 3 and 10]. These comparisons concern distinct statistical targets and observation laws.
1.2 How the three arguments work
For the integrated estimate, expand the qth density moment using q points of the measure. The collision of their labels is controlled by the projection of the q-1 differences from the first point. Successive coordinates of an orthonormal frame have ellipsoidal densities; their determinant powers telescope, leaving one nonnegative determinant power and an explicit \sqrt d scale. Gram–Schmidt converts the remaining change of variables into inverse affine heights. A tube about each preceding affine span is covered by balls, so the local mass hypothesis integrates these heights in reverse order. This proves Theorem 3.3.
The first application separates a Gaussian matrix into its positive radial factor and an independent orthonormal frame. The routing rule may use the whole radial factor; averaging that dependence produces weights on frame and projected label. The projected support occupies a k-ball. Its volume cancels the dimensional scale in the density norm and returns the same radius exponent as the incoming local-mass bound. The backward argument therefore preserves the ambient dimension through every block. Section 3 completes this route with local stopping-index counting and the common full-data bound of Section 2.
For positive cap domination, the order is different. Bounded nonnegative tests on a countable cap dictionary imply a local ball bound for a rescaled test measure. The normalized instance of the integrated theorem bounds the test of f_w\mu_B. To turn all these inequalities into a positive dominating measure, the proof uses a weighted sum of cap densities. Finitely many caps occur at each scale and the weighted cost makes the fine-scale L^1 tail vanish. This gives norm compactness. Strict separation of the closed downward convex set then produces exactly the bounded tests already controlled. The resulting positive sums propagate the unnormalized joint laws of signal and state, with a separate reference probability for each parent cap.
The specified-label argument uses q vectors anchored at z, rather than q-1 differences anchored at a sampled point. Their images under the Gaussian rows have covariance equal to the anchored Gram matrix. The resulting ball-average estimate is uniform in every fixed offset. Taking a single lower-limit version and applying Fatou’s lemma on each offset graph preserves that estimate for exact densities. Operator-norm bins and radial shells then control the Riesz potential of earlier continuations. This route retains a factor of order 1+\log(1/\epsilon) from its shell count and a global terminal clock; its own short-run argument controls both.
The three analytic conclusions and the costs of their regression applications are compared below. Here T is the deterministic sample horizon, and a stopping index is charged only where its value enters the corresponding argument.
| Estimate | Analytic guarantee | Regression argument |
| Integrated projection | A density moment integrated over the frame and label space | Backward local-mass propagation; stopping indices counted within each block |
| Positive cap domination | An actual dominating sum of cap measures, with weighted radius cost | Forward propagation of joint signal–state laws; one final sum over stopping indices |
| Fixed-offset density | One density version with a moment bound at every fixed affine offset | Backward Riesz-potential propagation; global terminal tags and a logarithmic shell count |
Section 2 supplies the probability normalizations, finite-state model, Borel-kernel reduction, and common full-data endpoint. Section 3 proves the integrated estimates and their backward application; these two sections suffice for the regression consequence and the explicit memory tradeoff. Section 4 uses the normalized projection moment to construct positive dominating measures and propagate them forward. Section 5 first proves its specified-offset estimate, then develops the Riesz recurrence and its parameter argument. The cap dual constraints test averages on spherical caps; they are distinct from integration on affine subspheres.
2 Spherical probabilities and exact observations
We collect the geometric facts used by all three estimates and then specify the finite-state experiment to which they will be applied. The geometric lemmas keep their normalizations visible: all measures on spheres below are probability measures, whereas densities on Euclidean label spaces are relative to Lebesgue measure unless another reference measure is stated.
2.1 Coordinates and caps
Write \sigma_d for uniform probability on the unit sphere S^{d-1}\subset\R^d, and write B(z,r) for the closed Euclidean ball with center z and radius r. When the ambient dimension is fixed we abbreviate \sigma_d to \sigma. The norm is Euclidean, and \norm A_{\op} denotes the operator norm of a matrix. Constants denoted by C are positive and absolute and may increase from one occurrence to the next. A subscript records any permitted additional dependence. All logarithms are natural unless a base is displayed.
Lemma 2.1 (Coordinates of a uniform direction). Let D\ge2 and 1\le k<D be integers. If U is uniform on S^{D-1} and P is any k-by-D matrix satisfying PP^{\mathsf T}=I_k, then PU has Lebesgue density
The direction of the component of U in \ker P is uniform on its unit sphere and is independent of the pair consisting of PU and the length of that complementary component. In particular, its conditional law given PU is uniform.
Proof. Rotation invariance reduces the projection to the first k coordinates. Let G\in\R^k and H\in\R^{D-k} be independent standard Gaussian vectors. Then (G,H)/\sqrt{\norm G^2+\norm H^2} is uniform on the sphere. The random variables \norm G^2/2 and \norm H^2/2 are independent unit-scale gamma variables of shapes k/2 and (D-k)/2. Changing variables from these two positive numbers to their sum and their first fraction shows that the fraction has density
The direction of G is uniform and independent of these lengths. Polar coordinates in \R^k, with \vol(S^{k-1})=2\pi^{k/2}/\Gamma(k/2), now give (1). The surface-area identity itself follows by integrating e^{-\norm x^2} in polar coordinates. The direction of H is uniform and independent of both lengths and of G, which proves the conditional assertion. The event H=0 is null.
The beta distribution in this lemma is a classical description of random orthogonal projections; see Frankl and Maehara [4]. We use the displayed derivation, including its conditional direction law.
Write v_k for the Lebesgue volume of the unit ball in \R^k. The following estimate will cancel the dimensional density scale in both projection arguments.
Lemma 2.2 (Unit-ball volume). For every integer k\ge1,
Proof. On the ball of radius \sqrt k, the function e^{-\norm z^2/2} is at least e^{-k/2}. Comparing its integral there with the full Gaussian integral gives
which is the claimed bound.
Lemma 2.3 (Spherical cap mass). Put n=d-1 for d\ge2. There is an absolute a>0 such that, for every c\in S^{d-1} and 0<r\le2,
For every ambient center z\in\R^d and radius r>0,
An output direction independent of the signal has probability at most \epsilon^n of angular error at most \epsilon, averaged over \sigma.
Proof. For r\le1, rotate c to the north pole. A sphere point in B(c,r) has last coordinate at least 1-r^2/2\ge1/2. Its surface element as a graph over c^\perp is at most twice Lebesgue measure, and its projection lies in the n-ball of radius r. The full sphere has area at least twice the volume of the unit n-ball, by projecting both hemispheres. This proves the upper bound for r\le1; total mass one proves it for r>1.
For the lower bound when r\le1, the angular radius 2\arcsin(r/2) is at least r. In polar angle the cap mass is the integral of \sin^{n-1}\theta over this interval divided by its integral over [0,\pi]. Use \sin\theta\ge2\theta/\pi for 0\le\theta\le r and bound the denominator by \pi. Since n\le2^{n-1}, the result is at least (r/\pi)^n. For 1\le r\le2, monotonicity therefore gives the lower bound with a=1/(2\pi).
If an ambient ball meets the sphere, choose c in the intersection. Its sphere portion is contained in B(c,2r), proving (4); an empty intersection is immediate. Finally, angular error at most \epsilon implies chordal error at most \epsilon, because 2\sin(\theta/2)\le\theta. Apply the sphere-centered upper bound to each fixed output and then integrate over any independent output randomization.
2.2 The finite-state regression experiment
An unknown direction s\in S^{d-1} is observed through independent standard Gaussian features and their exact labels:
At each sample index the learner retains a state from a set of at most 2^M values, where M is a nonnegative integer. A transition may use the current state, the entire current pair (X_t,Y_t), the sample index, and fresh independent randomness. Computation within the transition is unrestricted. Every piece of information depending on earlier samples that remains available later must be encoded by the finite state. The learner neither chooses a feature nor revisits an earlier one.
The learner stops after at most a prescribed integer T\ge0 samples. Its stopping decisions obey the same information restriction. For fixed rules, its output is a unit vector whose law depends only on the terminal state and stopping index, with fresh randomness allowed. In particular, the last observed pair influences the output only through the terminal state. The sample index is available without being stored in the state; the proofs below account for the information in the stopping index.
A shared random seed may select the rules before sampling. The pair consisting of this seed and the initial finite state is jointly independent of the signal and the entire sequence of feature rows. The seed and initial state may depend on each other. Fresh transition and output randomness is independent of these variables and of initialization. Nonuniform rules may also depend on d,\epsilon, and the sample index. This is a bit-memory model with a deterministic finite horizon.
We impose the following measurability convention. The original experiment is jointly measurable in the signal, pre-generated feature rows, shared seed, initial state, and fresh randomness, including the resulting state path, stopping index, and output. For almost every fixed shared seed, the finite vector of transition and stopping probabilities at each fixed state and index is measurable in the completion of \gamma_d(dx)\,dy, where \gamma_d=N(0,I_d). A rule given instead as a completed-measurable function of the current pair and an independent fresh seed is also allowed; Lemma 2.6 justifies this presentation. Each terminal output law is a probability measure on the Borel sphere.
For the lower bounds, let S\sim\sigma be independent of the feature rows. We call the probability of \arccos\inner{\widehat S}{S}\le\epsilon, including all learner randomness, the uniform-sphere average success probability.
Theorem 2.4 (Finite-memory precision bound). Let M(d)=o(d^2) be a sequence of nonnegative integers and let 0<\epsilon(d)\le1/10. There is an absolute c>0 such that, for all sufficiently large d, every learner in the experiment above whose uniform-sphere average success satisfies
must use
The dimension threshold may depend on the memory sequence. In particular, the conclusion holds when the same success probability is required separately for every s\in S^{d-1}.
An every-signal guarantee implies the average guarantee by integration against \sigma. All three proofs use the average guarantee itself. The first proof also yields a uniform dependence on the memory budget.
Corollary 2.5 (Finite-dimensional memory dependence). There are absolute constants c>0 and d_0 such that, for every integer d\ge d_0, every nonnegative integer M, and every 0<\epsilon\le1/10, a learner in the experiment above with uniform-sphere average success at least 2/3 and deterministic finite horizon T satisfies
The proof appears at the end of Section 3. In particular, the same precision bound holds with absolute constants throughout M\le d^2 once d exceeds an absolute threshold.
2.3 Measurable continuations and shared randomness
The block arguments will start fresh executions from specified states, including states the learner never reaches. The next lemma supplies Borel rules for all such continuations while preserving the actual uniform-prior experiment. Its proof uses the marginal law of each observation pair, so it does not require a density conditional on a selected state.
Lemma 2.6 (Borel versions and counterfactual continuations). Fix d\ge2, finite T, and finite state sets. Suppose each transition and stopping probability, for each fixed state and index, is measurable in the completion of the Borel sigma field of \R^d\times\R under \lambda(dx,dy)=\gamma_d(dx)\,dy, where \gamma_d is standard Gaussian probability on \R^d. The terminal output law depends only on the terminal state and index. There are everywhere Borel probability kernels with the same state sets and output laws that induce the same uniform-prior joint law of signal, state path, stopping index, and output. For the Borel learner, starting from any state and index and supplying fresh future features defines a measurable continuation success function, whether or not that state is reachable in the original execution. The same conclusion holds when a rule is given before seed integration by a finite probability vector measurable in the completion of \lambda\otimes\kappa, where \kappa is the law of an independent fresh seed.
Proof. First suppose the rule includes a fresh seed of law \kappa. Each completed-product measurable coordinate has a product-measurable version: approximate it by simple functions and replace the countably many level sets by product-measurable sets modulo null sets. Replace the resulting finite vector by a fixed point mass wherever it fails to be a probability vector. This repair occurs on a product-measurable null set. Integrating against \kappa now gives a Borel probability vector in the current pair, and Fubini shows agreement with the original integrated rule for \lambda-almost every pair. It therefore suffices to consider the coordinate probabilities in the first statement.
For each fixed state and index, combine the possible transitions and stops into a finite set of destinations, tagged as continuing or stopping. Choose Borel versions of their coordinate probabilities by the same simple-function construction. On the Borel null set where the versions fail to be nonnegative or to sum to one, replace the vector by a fixed point mass. Finitely many states and indices are involved, so all changes lie in one Borel \lambda-null set Z.
For every nonzero feature x, Lemma 2.1 gives a Lebesgue density for \inner{x}{S}. Since a Gaussian feature is nonzero almost surely, every pre-generated pair (X_t,\inner{X_t}{S}) has law absolutely continuous with respect to \lambda. Thus each pair avoids Z almost surely. Even if the current state depends on previous observations, the probability of being in a particular state and seeing a pair in Z is bounded by the marginal probability of seeing that pair in Z, which is zero.
Almost surely none of the finitely many pre-generated pairs lies in Z. Couple the original and modified finite destination kernels by the same fresh uniform random numbers. Induction gives identical state paths and stopping decisions; couple their unchanged terminal output laws identically as well. This proves equality of the specified joint laws. Finally, finite composition and integration of the Borel kernels against independent future features and coins defines a measurable continuation from every specified state and index, including unreachable states.
Convention for the three applications. We first condition on a shared seed for which the model assumptions hold. The joint independence of seed and initial state from signal and feature rows ensures that the conditional initial-state law still contains no signal or sample information. We then use the Borel learner supplied by Lemma 2.6. The lemma need only preserve the uniform-prior law; preservation at every fixed signal is unnecessary. All bounds below are uniform in the fixed rules and are averaged over the shared seed at the end. The original conditional success probabilities are measurable in that seed because the original experiment is jointly measurable. Neither the replacement kernels nor any subsequent proof decompositions need to be selected jointly across seeds.
2.4 A common full-data endpoint
The three moment arguments produce bounds involving T+d. The next elementary obstruction supplies the linear lower bound on T needed to remove that additive d, without using any memory restriction.
Lemma 2.7 (Antipodal full-data bound). In the experiment above, if \epsilon\le1/10 and T\le d/8, even an estimator given all T features and labels has uniform-prior success probability at most 5/8.
Proof. Let F be the row span of the pre-generated feature matrix and put u=P_FS. Its dimension is T almost surely, with the evident interpretation for T=0. The labels determine u. For T=0, the residual is the original uniform signal. For T>0, Lemma 2.1, applied after fixing the feature matrix, shows that conditional on all the data the residual is \sqrt{1-\norm u^2}\,\omega, with \omega uniform on the unit sphere of F^\perp. Before conditioning on the labels, isotropy gives \E\norm u^2=T/d. Consequently \Prob\{\norm u^2\ge1/2\}\le2T/d.
When \norm u^2<1/2, opposite residual directions give two signals at distance 2\sqrt{1-\norm u^2}>\sqrt2>2\epsilon. They have the same conditional law, and at most one can be within chordal distance \epsilon of a fixed output. Angular success implies this chordal condition, so conditional success is at most 1/2, also after averaging independent estimator randomness. Total success is at most
Giving the estimator all pre-generated data also covers a learner that stops at a prefix, with its independent shared and fresh randomness.
3 Integrated moments of orthonormal projections
This section answers the first question: how does a bound on the mass of every Euclidean ball control the density after a random orthonormal projection, when the label is integrated over its whole space? The proof has two geometric parts. A frame calculation expresses the cost of projecting several points through their affine Gram determinant. A covering argument then integrates that determinant using the ball-mass bound. We retain the dimensional factor in the density norm because it cancels against the projected support volume in the regression application.
For integers 1\le k<d, let \mathcal V_{d,k} be the space of k-by-d matrices P satisfying PP^{\mathsf T}=I_k. Write \vartheta_{d,k} for the law of the first k rows of a uniform orthogonal matrix. One direct construction fixes the probability normalization. For a square standard Gaussian matrix G, which is invertible almost surely, put O=G(G^{\mathsf T}G)^{-1/2}. Then O is orthogonal. For fixed orthogonal U,V, the same formula applied to UGV gives UOV; the Gaussian law is unchanged by both multiplications. Thus the law of O is a probability invariant on both sides, the Haar probability used here. Its columns can equivalently be sampled successively uniformly in the orthogonal complements of the preceding columns. To verify the equivalence, average either invariant frame law under an independent O; transitivity makes the averaged law independent of the initial frame. As in Section 2, v_k is the unit-ball volume.
3.1 The projected frame and its determinant
We first record the classical truncated Haar density, together with the affine-height form used below. Jiang [7, Lemma 2.5] records this formula and refers to earlier formulations by Diaconis, Eaton, and Lauritzen and by Eaton. We include its derivation to retain the dimensional normalization.
Lemma 3.1 (Projected frame density). Let k,m\ge1 and k+m+1\le d. Fix a d-by-m matrix E with orthonormal columns, and draw P\sim\vartheta_{d,k}. The k-by-m matrix PE has Lebesgue density
Here >0 means positive definite. In particular, this density is bounded by (C\sqrt d)^{km}.
If x_1,\ldots,x_{m+1}\in\R^d are affinely independent and
then the joint Lebesgue density of the columns P(x_i-x_1), 2\le i\le m+1, is bounded by
Proof. Right rotational invariance allows us to replace the fixed columns of E by an independent uniform orthonormal frame. Conditional on P, its projection has the same law as the first k coordinates of that frame. Sample the full columns successively and let Z_j be the matrix of the first k coordinates of its first j columns. Set
Given the preceding full columns, the next column is uniform on the unit sphere in their d-j-dimensional orthogonal complement. The coordinate map from that complement to the first k coordinates has row Gram matrix S_j. When S_j>0, multiplication by S_j^{-1/2} makes its rows orthonormal. Lemma 2.1 and a linear change of variables therefore give the next projected column z conditional density
on z^{\mathsf T}S_j^{-1}z<1. This expression depends only on Z_j, so it is also the conditional density given those projected columns. The strict inequality implies S_{j+1}=S_j-zz^{\mathsf T}>0 almost surely. The determinant identity
shows that multiplication of the m conditional densities cancels every intermediate determinant. The remaining factor is \det S_m^{(d-k-m-1)/2}, proving (5). Its exponent is nonnegative and 0<\det S_m\le1.
The normalizing constants have the stated scale. For a unit-scale gamma variable G of shape a=(d-j-k)/2>0, Cauchy–Schwarz and the gamma recurrence give
Together with the factor \pi^{-k/2}, this bounds each c(d-j,k) by (C\sqrt d)^k.
For the last assertion, Gram–Schmidt writes [x_2-x_1,\ldots,x_{m+1}-x_1]=ER, where E has orthonormal columns and R is triangular with \abs{\det R}=\prod_{i=2}^{m+1}\delta_i. Right multiplication by R on k-by-m matrices has Jacobian \abs{\det R}^k, one factor for each row. Applying this change of variables to the preceding density proves (6).
The product of the affine heights in (6) is m! times the m-dimensional volume of the simplex spanned by the points. Related simplex-volume factors occur in Drury’s affine-plane Jacobian [3, Lemma 1, pp. 497–498]. Here the exponent and normalization come from the displayed k-row change of variables.
3.2 From local mass to an exact density norm
The local growth hypothesis below is of Frostman type, and inverse-distance averaging is a classical projection method; compare Mattila [8, Lemma 5.1, p. 236, and Theorem 5.4, p. 237]. We prove the quantitative affine-tube and growing-order moment bounds needed here.
Lemma 3.2 (Affine inverse distances). Let 0<\beta\le d, K\ge0, and R>0. Let \nu be a finite positive Borel measure supported in B(z,R)\subset\R^d such that
If F is an affine j-plane, k\ge1 is an integer, and \beta-j-k\ge1, then \nu(F)=0 and
provided j,k\le d.
Proof. For 0<t\le R, the projection of B(z,R) onto F is contained in a j-ball of radius R. A maximal t-separated subset is a t-net with at most (3R/t)^j points, by packing disjoint j-balls of radius t/2. For j=0, one point suffices. The corresponding ambient balls of radius 2t cover the points of B(z,R) within distance t of F. Thus
Because \beta-j>0, continuity from above as t\downarrow0 gives \nu(F)=0. Distances greater than R contribute at most R^{-k}\nu(B(z,R))\le KR^{\beta-k}. On the shell 2^{-h-1}R<\dist(x,F)\le2^{-h}R, (9) gives a contribution at most
The sum over h\ge0 is at most twice its first term. Since \beta,j,k\le d, the resulting factor is bounded by C^d.
Theorem 3.3 (Integrated orthonormal-projection moment). Let d,k,q be integers and \beta a real number with k\ge1, q\ge2, and
Let \nu be a finite positive Borel measure supported in B(z,R), R>0, satisfying (7) for some K\ge0. Assume that P_\#\nu is absolutely continuous with respect to Lebesgue measure for \vartheta_{d,k}-almost every P. There is a jointly measurable version g(P,u) of its density such that
For every P, this version vanishes when \norm{u-Pz}>R.
Proof. When K=0, the measure is zero and the assertion is immediate. For K>0, let v_k be the unit-ball volume and define
These are jointly measurable nonnegative functions. Expanding their qth powers and using Tonelli’s theorem expresses the integrated moment through q points x_1,\ldots,x_q drawn from \nu. The intersection of the q label balls of radius \delta has volume at most v_k\delta^k; it is empty unless \norm{P(x_i-x_1)}\le2\delta for i\ge2.
Put \delta_i=\dist(x_i,\aff(x_1,\ldots,x_{i-1})). Lemma 3.2 shows that every affine plane encountered here has zero \nu-mass: its dimension is at most q-2, so the last condition in (10) applies. Consequently the tuple is affinely independent for \nu^{\otimes q}-almost every choice. For such a tuple, Lemma 3.1 applies with m=q-1; the condition k+q\le d is exactly its determinant condition. The probability of all q-1 inequalities is bounded by
The label-ball volumes cancel those in the definition of g_\delta. Absorbing 2^{k(q-1)} in C, we obtain
The degenerate tuples form a null set and can be assigned an infinite upper bound in this nonnegative integral.
Integrate first in x_q, then in x_{q-1},\ldots,x_2. At each step Lemma 3.2 bounds the current factor by C^d K R^{\beta-k}. The remaining first-point mass is at most KR^\beta. Hence the last integral in (12) is at most
This proves the qth power of (11) uniformly in \delta.
Define g(P,u)=\liminf_{j\to\infty}g_{1/j}(P,u). The assumed absolute continuity and Lebesgue differentiation [15, Theorem 1.6.19, pp. 146–147] imply that, for almost every P, this is a density of P_\#\nu for almost every u. Fatou’s lemma on the product space [15, Corollary 1.4.47, p. 110] gives (11). If \norm{u-Pz}>R, sufficiently small label balls miss P(B(z,R)), since P is a contraction; the displayed version is then zero. This limiting calculation concerns the exact projected measure. It introduces no noise into an observation.
We now spell out the two parameter choices that will be used. Their normalizations differ, even though their frame calculation is shared.
Corollary 3.4 (All-radii projection of spherical local mass). For sufficiently large d, put
Let 0\le f\le1 be Borel on S^{d-1}, and suppose for K\ge0 that
For z\in\R^d, R>0, and E=B(z,R), the projection of \1_E f\sigma by P\sim\vartheta_{d,k} has a jointly measurable Lebesgue density g(P,u) satisfying
Proof. Restriction preserves the stated all-ball bound. For every P, the projection of \sigma is absolutely continuous by Lemma 2.1, so the same holds for the restricted weighted measure. The exact margins are
Thus Theorem 3.3 applies. This also records the nonnegative frame-determinant exponent (d-k-q)/2=(d-2\lfloor d/4\rfloor)/2.
Theorem 3.5 (Normalized moment for unit-ball measures). For sufficiently large d, put
Let X be a k-by-d standard Gaussian matrix and define
on its probability-one full-row-rank event, and set P to a fixed frame on the null complement. Let m_k be uniform probability on the unit ball of \R^k. Suppose \nu is a finite positive Borel measure fixed independently of X, of mass at most one, supported in the unit ball of \R^d, such that for a fixed D>0,
Assume P_\#\nu is absolutely continuous for almost every X. There is a jointly measurable nonnegative g_X(z) such that, for almost every X, P_\#\nu(dz)=g_X(z)\,dz/v_k. This version vanishes when \norm z>1 and satisfies
Proof. Right rotational invariance of X makes P have law \vartheta_{d,k}; the independence assertion about P,H used later is proved below. The exact determinant and inverse-distance margins are
Apply Theorem 3.3 with q=p, \beta=\lambda, K=D^d, and R=1. If \ell(P,z) is its Lebesgue density, then g_X(z)=v_k\ell(P,z) is the required normalized density. Moreover
Use the unit-ball volume estimate (2). Since k\ge d/32 for sufficiently large d,
Theorem 3.3 now bounds the normalized L^p norm by \exp(C_Dd), proving (14). The density version and its support are inherited from that theorem. The mass-at-most-one and absolute-continuity hypotheses have been retained explicitly.
3.3 Radial averaging and dimension-preserving propagation
We next apply the all-radii estimate to the actual Gaussian data. A learner may use both the row directions and their lengths and correlations. The following factorization retains all of that radial information before averaging it into the routing weights.
Lemma 3.6 (Gaussian radial factor and frame). Let X be a k-by-d standard Gaussian matrix with k<d. Almost surely,
The frame P has law \vartheta_{d,k} and is independent of H.
Proof. Full row rank holds almost surely because the rank-deficient matrices are the zero set of nonzero minors and the Gaussian law has a density. Right multiplication of X by a fixed orthogonal matrix leaves its law unchanged, fixes H, and right-rotates P. Let O be an independent uniform orthogonal matrix. For bounded measurable \phi,\psi, invariance and averaging give
where \mu is uniform orthogonal probability. The inner integral is independent of P: the right action on \mathcal V_{d,k} is transitive, and \mu is invariant. Its value is \int\psi\,d\vartheta_{d,k}. This factorizes the expectation for all bounded measurable tests, which proves the stated law and independence.
Proposition 3.7 (One-block all-radii propagation). For sufficiently large d, use k=q=\lfloor d/4\rfloor and \beta=d-1-q. Let f_1,\ldots,f_N:S^{d-1}\to[0,1] be Borel functions satisfying
Let Q_w(X,y) be nonnegative Borel functions on \R^{k\times d}\times\R^k with \sum_wQ_w\le1, and define
Then
Proof. Write X=HP as in Lemma 3.6 and average the full radial dependence of the routing rule:
These are measurable nonnegative weights with \sum_w\tau_w\le1. Independence gives the exact identity
Fix E=B(z,R), and let g_w(P,u) be the exact density in Corollary 3.4 for \1_Ef_w\sigma. After integrating (16) over E, pushforward integration yields \int\!\int\sum_w\tau_w g_w\,du\,d\vartheta_{d,k}. Each g_w is supported in \norm{u-Pz}\le R. On this domain the measure with density \tau_w, relative to frame probability, label Lebesgue measure, and counting measure in w, has total mass at most v_kR^k. Hölder’s inequality therefore gives
Substitute (13). The radius powers add exactly to \beta. Also k\ge d/8, so (2) gives
This proves (15). The original Q_w could use all entries of H; only after that dependence was retained did we average it into \tau_w.
3.4 Backward propagation of continuation success
We now apply Proposition 3.7 to the success probability of a fresh continuation from each state. Its local-mass exponent survives every block, and the resulting finite-dimensional estimate retains the full dependence on M.
Proposition 3.8 (Success bound from orthonormal projections). There are absolute constants C\ge1 and d_0 such that the following holds for every integer d\ge d_0, every nonnegative integer M, every 0<\epsilon\le1/10, and every learner of deterministic finite horizon T in Section 2. Put
Then
The constant and dimension threshold are independent of M,\epsilon,T and the learner’s rules.
Proof. Use the fixed-seed and Borel-kernel convention of Section 2, and set n=d-1 and \beta=n-q. For each terminal state and stopping index, let h(s) be the success probability of its signal-independent output law. The cap estimates of Lemma 2.3 give
The last inequality follows from R\le\epsilon or R\ge\epsilon, respectively. Thus the all-radii coefficient for every terminal function is at most K_0=C^d\epsilon^q.
Divide the horizon into b blocks of k rows, padding the last block with independent rows ignored by its rules. Starting at a specified state and boundary index, define its continuation by a fresh future execution. This definition does not condition on that state being reached. The continuation functions are Borel and take values in [0,1] by Lemma 2.6.
Within a block, a destination is either a continuing state at the next boundary or a terminal state tagged by its local stopping index. The possible local stopping indices are 0,1,\ldots,k, including an immediate stop. There are at most (k+1)2^M terminal destinations and 2^M continuing destinations, hence at most N_* in total. The global stopping index is determined by the fixed block start and this local index, so the terminal output retains its full allowed dependence on time.
Finite composition of the transition kernels produces Borel routing probabilities Q_w(X,y) for these destinations, with \sum_wQ_w\le1. Their continuation functions f_w use only future independent rows and coins; terminal destinations use the functions h above. The incoming continuation therefore has precisely the form
from Proposition 3.7. The Gaussian block is fresh: conditional on the fixed shared seed, its rows are jointly independent of the signal and entering state, which depends only on initialization, past pairs, and past coins.
Backward induction through the blocks multiplies the all-radii coefficient by at most C^dN_*^{1/q} per block. This factor is at least one, so the coefficient also bounds any earlier terminal option. Integrate the initial continuation over B(0,1) and average over the conditional initial-state law, which is independent of the signal and rows. Conditional success is at most the right side of (17). The bound is uniform in the fixed rules, so the shared-seed convention gives the same bound for the original learner. When T=0, the terminal bound supplies the conclusion directly.
Proof of Corollary 2.5 and first proof of Theorem 2.4. Put L=\log(1/\epsilon). Choose an absolute dimension threshold so that Proposition 3.8 applies and k=q\ge d/8. If success is at least 2/3, (17) implies
Since b\le8T/d+1 and \log(k+2)/q is bounded by an absolute constant, this gives an absolute C_1 such that
No bound on M was used. Lemma 2.7 forces T>d/8, and hence T+d<9T. Substitution into (18) proves Corollary 2.5 with an absolute constant.
Finally, if M(d)=o(d^2), then 1+M(d)/d^2\le2 for all sufficiently large d. The corollary proves Theorem 2.4, with an absolute lower-bound constant and a dimension threshold that may depend on the memory sequence. This proof is uniform over the accuracy sequence and uses only within-block stopping indices.
4 Positive domination by cap measures
This section proves the sample bound by propagating the joint law of the signal and each memory state forward through the stream. We bound each such law by a positive sum of uniform measures on spherical caps. A term is charged for its mass and for a power of the reciprocal of its radius. The main step shows that a block of observations increases the total charge by a controlled factor. The passage from bounds on test functions to an actual dominating measure requires compactness in the norm of L^1.
We prove Theorem 2.4 by a forward argument in the finite-state experiment of Section 2. Use that section’s fixed-seed convention and Borel kernels from Lemma 2.6. All bounds below are uniform in the fixed rules; their scalar probability conclusions are averaged over the shared seed at the end.
4.1 The spherical experiment and the projection input
Write \sigma for uniform probability on S^{d-1}, and set
We take d sufficiently large that p\ge2 and k,p\ge d/32. For c\in S^{d-1} and 0<r\le2, define
These are spherical caps and their normalized probability measures. For an ambient ball in \R^d we use the distinct notation D(b,t)=\{u\in\R^d:\norm{u-b}\le t\}.
Lemma 2.3 supplies the two-sided estimate
In particular every \mu_{c,r} is well defined.
Let X be a k\times d matrix with independent standard Gaussian entries and put
on the probability-one event that X has full row rank. Define P arbitrarily on the remaining null set. We may make this definition measurable. On the full-rank event PP^{\mathsf T}=I_k and X=HP; the law of P is invariant under right orthogonal multiplication. Let v_k be the Lebesgue volume of the unit ball of \R^k and let m_k(dz)=v_k^{-1}\1_{\{\norm z\le1\}}\,dz be uniform probability on that ball.
The use of independent replicas for high projection moments is related to the argument of Sharan, Sidford, and Valiant [12, Section 7, Lemma 12]. We use the normalized moment theorem (Theorem 3.5) with the parameters in (19). Its inverse-distance margin satisfies \lambda-(p-2)-k\ge3d/8. Below we verify its unit support, mass, local-growth, independence, and absolute-continuity hypotheses for each finite positive test measure before invoking its normalized density bound.
4.2 One block: from bounded tests to a dominating measure
Fix a parent cap B=B(c,r). Let \mathcal W be a finite set and let K_w(X,Y)\in[0,1], for w\in\mathcal W, be Borel functions of a k\times d matrix and a vector in \R^k with \sum_wK_w\le1. Define
The functions f_w are measurable. The variables z in the second line are independent of X and are used only to define a reference law. The kernel still receives the full matrix X and the radial factor H through the displayed label. In particular, \sum_wq_w\le1.
Proposition 4.1 (Positive domination after one block). For each w\in\mathcal W, there are countably many caps B(c_i,\rho_i) with 0<\rho_i\le r and coefficients a_i\ge0 such that
The first inequality is an inequality of finite positive measures on S^{d-1}: it holds after evaluating both sides on every measurable set. The constant C is independent of the parent cap and of the kernels.
Proof. We prove a bound for nonnegative tests and then convert it into a countable positive domination.
Cap tests and their local mass. Put \rho_j=r2^{-j} and choose a finite \rho_j/2-net N_j of the sphere for every j\ge0. Take as a dictionary the parent cap, indexed by *, together with all B(c',\rho_j), c'\in N_j. Repetitions may have separate indices. For this countable index set I, write B_i=B(c_i,\rho_i) and
Thus \phi_*\,d\sigma=\mu_{c,r}, every \phi_i is a probability density, b_*=1, and the weights at level j are 2^{\alpha j}. If h\in L^\infty(\sigma) is nonnegative and satisfies
then its integral against any positive dictionary sum is at most that sum’s weighted cost. We will bound \int h f_w\,d\mu_{c,r} for every such test.
Let \nu be the image of h\mu_{c,r} under s\mapsto(s-c)/r. It has unit-ball support and mass at most one by the parent constraint. We claim that, for an absolute D\ge1,
For 0<t\le1/4, suppose the ball meets the rescaled cap, and choose a point s_0\in B whose image lies in it. The corresponding sphere points are all within 2rt of s_0. Set
A center in N_j within \rho/2 of s_0 gives a dictionary cap containing those points, since 2rt+\rho/2\le\rho. Its constraint and the two-sided cap estimate imply
An empty intersection has zero mass. For t>1/4, use the mass bound and 1\le4^\lambda t^\lambda. These estimates prove (24).
The quantitative test estimate. All hypotheses of Theorem 3.5 now hold for \nu. Indeed, the test, parent cap, and measure are fixed before the fresh matrix X is drawn. For every orthonormal k-frame P, the coordinate density in Lemma 2.1 gives P_\#\sigma\ll dz. Because h\mu_{c,r}\ll\sigma, its projection is also absolutely continuous, and the affine rescaling
gives the same conclusion for P_\#\nu. The theorem therefore supplies a jointly measurable density satisfying, for almost every X,
The density is relative to dm_k=dz/v_k on the unit ball; D, and hence the new constant, is absolute. Since s=c+ru gives Xs=Xc+rHPu, Tonelli’s theorem and Hölder’s inequality on the product law of X and z yield
The last line uses 0\le K_w\le1. Thus every test satisfying (23) obeys (25).
Representation by positive cap sums. The geometric estimate is complete; the remaining step uses only the dictionary and its test bound. If q_w=0, the feasible test h=1 shows that f_w\mu_{c,r}=0, and the empty sum suffices. Otherwise set K=\exp(Cd)q_w^\theta>0, with the constant from (25), and define in real L^1(\sigma)
Each series converges in L^1, since b_i\ge1 and \norm{\phi_i}_{L^1}=1. This convex set is norm compact. To see this, take a sequence and choose a coefficient array for each term. A diagonal subsequence converges coefficientwise, since 0\le a_i\le K/b_i. The limiting array has cost at most K by Fatou’s lemma. There are finitely many indices through every level J, and uniformly over all admissible arrays,
The same estimate holds for the limiting array. Finite-level convergence followed by (26) gives convergence in L^1. This proves sequential compactness, and hence compactness in the metric space L^1(\sigma) [10, Section 1.3, Theorem 1.37, p. 29].
Its signed downward closure
is convex, contains every nonpositive L^1 function, and is norm closed. In fact, if u_m\to u and u_m\le v_m\in\mathcal C_K, compactness provides a subsequence v_m\to v\in\mathcal C_K; then
Put F=f_w\phi_*. If F\notin\mathcal D_K, strict separation of the closed convex set \mathcal D_K and the compact singleton \{F\} [10, Section 5.2, Theorem 5.20(b), pp. 98–99] gives a continuous real functional \ell with
A probability measure is localizable [5, Theorem 211L], so real L^1–L^\infty duality [6, Theorem 243G(b)] represents \ell(u)=\int hu\,d\sigma for a real h\in L^\infty(\sigma). The finite supremum forces h\ge0: if h\le-\delta on a set E of positive measure, the functions -t\1_E\in\mathcal D_K would make it infinite. Define
Then
The cost bound gives the upper inequality, and the single terms (K/b_i)\phi_i give the reverse inequality by taking the supremum over i. This does not require the supremum to be attained. If S_h=0, the parent term gives \int h\phi_*\,d\sigma=0, and 0\le f_w\le1 implies \ell(F)=0, contrary to strict separation. If S_h>0, the function h/S_h satisfies (23), so (25) again contradicts separation. Therefore F\in\mathcal D_K. An admissible series dominates its density almost everywhere; integrating gives the asserted measure inequality and cost bound. Its total mass is at most K.
4.3 Propagation, bounded stopping, and the accuracy endpoint
The preceding proposition applies to each cap in a joint-state bound. The positivity of all measures makes the iteration linear before any estimate is taken.
Lemma 4.2 (Cost of a forward block). Let \mathcal V and \mathcal W be finite state sets with \abs{\mathcal W}\le2^M. For each v\in\mathcal V, let K_{v,w} be block kernels of the type in (21), with likelihoods f_{v,w}. Let \eta_v be finite positive measures on the sphere and suppose
where each sum is countable, a_{v,i}\ge0, and 0<r_{v,i}\le2. Define the outgoing measures by
Then the \eta'_w admit countable positive cap dominations whose total cost \mathcal P' satisfies
Proof. Multiplication by a nonnegative measurable function preserves a measure inequality, as follows first for simple functions and then by monotone convergence. Summing over incoming states and cap terms therefore gives
Apply Proposition 4.1 separately to every term. For a fixed incoming cap, write q_{v,i,w} for its reference probabilities. They depend on that cap as well as on the starting state, and satisfy \sum_wq_{v,i,w}\le1. Hölder’s inequality for the finite sum gives
If the outgoing terms for this parent have coefficients \beta_j and radii \rho_j, their contribution to the global cost is bounded by
Sum first over w and then over the countably many incoming terms to obtain (28). Tonelli’s theorem applies throughout because every summand is nonnegative.
At each of finitely many block boundaries there are only countably many state–cap terms. Choosing a dominating sum for each produces a deterministic proof representation, with another countable sum at the next boundary. There is no choice indexed by the random signal or by the fresh matrix. The state measures here remain unnormalized joint measures; the argument never divides them by their total masses.
Lemma 4.3 (Terminal success under a cap). Let Q be any probability law on unit-vector outputs, independent of the signal, and put h_Q(s)=Q\{z:\arccos\inner{z}{s}\le\epsilon\} for 0<\epsilon\le1/10. For every cap B(c,r),
Proof. For a fixed unit output z, write \vartheta=\arccos\inner{z}{s}. Angular success implies \norm{s-z}=2\sin(\vartheta/2)\le\epsilon. If \epsilon<r, Lemma 2.3 bounds its probability under \mu_{c,r} by
If \epsilon\ge r, use probability at most one and (\epsilon/r)^\alpha\ge1. Integration over Q proves the result.
Proposition 4.4 (A fixed stopping index). For all sufficiently large d, let a learner obey the observation and information restrictions in Section 2, with 0<\epsilon\le1/10, any nonnegative integer M bits of state, and horizon T, and let \mathsf T be its stopping index. For every fixed integer 0\le t\le T,
In particular, for a sequence M(d)=o(d^2) the right side is eventually at most
with an absolute C and a dimension threshold that may depend on the memory sequence.
Proof. Fix t>0 throughout the argument and split its first t indices into b=\lceil t/k\rceil consecutive blocks of length at most k. At a prefinal boundary retain only paths that have not stopped and are eligible to continue. In the final block retain exactly the paths stopping at t, with their terminal state. For a fixed starting state and a supplied block of pairs, integrating the internal randomness gives the probability of each retained ending state. These are Borel kernels with sum at most one, since all omitted paths are discarded. A shorter block is padded by independent unused Gaussian rows, which the kernels ignore. Its likelihood for a signal s is therefore f_{v,w}(s) of (21). Conditional on the fixed seed, the fresh Gaussian feature rows in the block are jointly independent of the signal and entering state, indeed of the signal, initialization, past rows, and past coins. Given the signal, the past affects this likelihood only through the starting state and the fixed index.
At any counted boundary, let \eta_v(E) be the probability that S\in E and the path is counted in state v. These are unnormalized joint measures and their outgoing laws have exactly the formula in Lemma 4.2. Initially, B(c,2)=S^{d-1} and \mu_{c,2}=\sigma. Write V_0 for the initial finite state. For the fixed shared seed \xi, joint initialization independence gives
Thus the initial total cost is one before any discarded paths are removed, and at most one afterward. After b applications of that lemma, the counted terminal measures admit positive cap dominations of total cost at most [\exp(Cd)2^{M/p}]^b.
At each fixed terminal state and index t, the output is drawn from a fixed law using fresh randomness. Apply Lemma 4.3 to each term of its dominating cap sum. A term of coefficient a and radius r contributes at most \exp(Cd)(\epsilon/2)^\alpha a(2/r)^\alpha. Summing proves the conditional version of (30). At t=0, use the initial domination and the same terminal estimate directly. The index has been fixed throughout, so only the at most 2^M states at each ending boundary were counted. No additional state is charged for the value of t.
The countable dominating sums are chosen only after fixing the shared seed. The bound is uniform in that seed, so the common convention of Section 2 allows us to average the original conditional success probabilities. No measurable choice of cap decompositions across shared seeds is required. Finally, for M(d)=o(d^2) we eventually have M(d)\le d^2. Since p\ge d/32, the logarithm of the extra factor 2^{M/p} is then at most 32(\log2)d. Absorbing it into the absolute exponential constant proves (31).
We now sum the fixed-index estimate over all possible stopping times and use the common full-data bound to absorb the resulting lower-order terms.
Proof of Theorem 2.4 by positive cap domination. Put L=\log(1/\epsilon).
For the given sequence M(d)=o(d^2), use (31) in all sufficiently large dimensions. The stopping events are disjoint. Since k\ge d/32, \alpha\ge d/3, and \lceil t/k\rceil\le t/k+1, summing over 0\le t\le T yields
Here (\epsilon/2)^\alpha\le\exp(-dL/3), and C_1,C_2 are absolute. If success is at least 2/3, taking logarithms gives
since \log(T+1)\le T for every T\ge0. Lemma 2.7 also gives T>d/8, and therefore T\ge1. Consequently
This proves the claimed bound with an absolute constant, uniformly over the full accuracy range. Only the dimension at which the memory factor is absorbed may depend on the sequence M(d).
5 Riesz estimates at fixed affine offsets
An integrated density norm does not specify the value of a density on a matrix-dependent graph of labels. We now construct a single density version whose moment is controlled at every fixed affine offset. The bound uses the supremum of ambient Riesz potentials of the original measure. We first prove this analytic estimate, then use it on distance shells to propagate finite-state continuation probabilities backward. The shell argument retains an explicit logarithmic loss.
Throughout this section, fix an integer d\ge16 and put
These choices give
The matrix X will always denote a k by d standard Gaussian matrix. Constants denoted by C\ge1 are absolute and may increase between uses.
5.1 Ambient potentials and affine distances
For a measurable h:S^{d-1}\to[0,1], define the supremum of its ambient Riesz potentials by
For a measure \mu, the raw Riesz potential is U_\alpha^\mu(z)=\int\norm{s-z}^{-\alpha}\,d\mu(s) [1, Section 1, p. 2]. Thus V(h) is its supremum over ambient centers for \mu=h\sigma; it is not the Riesz energy. The kernel is understood as a nonnegative extended-valued function at s=z. Write B_r(z)=\{x\in\R^d:\norm{x-z}\le r\}, and intersect such balls with the sphere when integrating against \sigma.
Lemma 5.1 (Ambient potential bounds). For every z\in\R^d and 0<r\le1/10,
Every measurable h:S^{d-1}\to[0,1] has finite potential. Moreover, for 0<\epsilon\le1/10,
Proof. The ambient estimate (35) follows from Lemma 2.3, with 2^{d-1} absorbed into C^d. For dyadic outer radii r_m=\epsilon2^{-m}, the shell r_m/2<\norm{s-z}\le r_m contributes at most (2/r_m)^\alpha C^d r_m^{d-1}. The point s=z has zero \sigma-measure, and \sum_{m\ge0}2^{-m(d-1-\alpha)}\le2. Summing proves (36). The same argument at radius 1/10, together with the kernel bound 10^\alpha outside that ball, proves V(h)<\infty uniformly over 0\le h\le1.
The potential bounds the measure of a tube around any affine subspace through an ambient center. We use that fact to integrate successive inverse distances when computing the Gaussian projection moment.
Lemma 5.2 (Inverse distance to an anchored affine space). Let h:S^{d-1}\to[0,1] be measurable, z\in\R^d, \rho>0, and let F be an affine subspace through z of dimension \ell\le q-1. Then F\cap B_\rho(z) has zero h\,d\sigma mass, and
Proof. For every ambient center c and radius r>0, definition (34) gives
Fix 0<r\le\rho. A maximal r-separated subset of F\cap B_\rho(z) has at most (1+2\rho/r)^\ell \le(3\rho/r)^\ell points, by comparing disjoint radius-r/2 balls in F with the radius-(\rho+r/2) ball. It is an r-net. The orthogonal projection onto F of any point in B_\rho(z) lies in F\cap B_\rho(z), since z\in F. Consequently the part of B_\rho(z) within distance r of F is covered by the corresponding ambient radius-2r balls. Equation (38) yields
This includes \ell=0, with a one-point net. Since \alpha-\ell>0, letting r decrease to zero shows that F\cap B_\rho(z) has zero mass.
Every point of B_\rho(z) is within distance \rho of F. Apply (39) on the shells with outer distance r_m=\rho2^{-m} and bound the inverse power there by (2/r_m)^k. Their sum is at most
Here \alpha-\ell-k\ge\alpha-(q-1)-k\ge2, so the series is at most 2. Since \ell\le d/8, \alpha=d/2, and k\le d/8, the remaining constant is at most C^d. This proves (37).
5.2 A density at fixed affine offsets
The next estimate uses q points anchored at the chosen center z. Their images under each Gaussian row have an anchored Gram matrix as covariance. The inverse-distance bound just proved will integrate the negative determinant power in the joint Gaussian density.
Lemma 5.3 (Projection density at deterministic offsets). Let z\in\R^d, \rho>0, let H\subset S^{d-1}\cap B_\rho(z) be measurable, and let h:S^{d-1}\to[0,1] be measurable. Fix these objects before drawing a k by d standard Gaussian matrix X. The pushforward of h\mathbf 1_H\,d\sigma by s\mapsto Xs admits a single jointly Borel function
such that R(X,\cdot) is its Lebesgue density for every full-row-rank X, and, for every fixed b\in\R^k,
The same function R is used for all b. Moreover,
Proof. Let \mu be the finite Borel measure given by \mu(E)=\int_E h\mathbf1_H\,d\sigma on Borel subsets of the sphere. For \eta>0, let B_\eta^k=\{u\in\R^k:\norm u\le\eta\} and let \operatorname{vol}_k denote k-dimensional Lebesgue measure. Define
The integrand is Borel in (X,y,t), so the parameter integral is Borel in (X,y). We first prove (40) with R_\eta in place of R, uniformly in \eta.
Expand its qth moment at Xz+b by Tonelli’s theorem, using t_1,\ldots,t_q from the q copies of \mu. For a fixed tuple, the columns X(t_j-z) form an array whose k rows are independent centered Gaussian vectors with covariance
If G is nonsingular, their joint density is bounded everywhere by (2\pi)^{-kq/2}(\det G)^{-k/2}. Therefore the probability that all columns belong to the radius-\eta ball about b is at most this bound times \operatorname{vol}_k(B_\eta^k)^q.
The singular tuples have zero \mu^q mass. Indeed, at the first linearly dependent vector t_j-z, the point t_j lies in z+\operatorname{span}(t_1-z,\ldots,t_{j-1}-z), an affine space through z of dimension at most q-1. Lemma 5.2 assigns zero \mu mass to it. The degeneracy set is measurable, and Fubini’s theorem applied at each j proves the claim. For a nonsingular tuple, Gram–Schmidt triangularizes the matrix with columns t_j-z with diagonal entries equal to the successive orthogonal distances. Taking the determinant of its Gram matrix gives
Integrate this product first in t_q, then in t_{q-1}, and so on. For every nonsingular prefix, Lemma 5.2 bounds the next integral by C^dV(h)\rho^{\alpha-k}. The exceptional prefixes have zero product measure by the same degeneracy argument. It follows that the integral of (43) is at most (C^dV(h)\rho^{\alpha-k})^q. The ball volumes cancel in the moment expansion, and the Gaussian prefactor is at most one. Thus
for every \eta>0 and every fixed b.
Now choose the function once and for all:
It is jointly Borel. To see that it is a density, fix a full-row-rank X. Set D=(XX^{\mathsf T})^{1/2} and P=D^{-1}X, so D is invertible, P has orthonormal rows, and X=DP. Lemma 2.1 makes P_\#\sigma absolutely continuous, and the invertible map D preserves absolute continuity. Since \mu\le\sigma, the pushforward X_\#\mu is also absolutely continuous. The Lebesgue differentiation theorem [15, Theorem 1.6.19], applied to its L^1 density, shows that (44) equals that density almost everywhere in y. A Gaussian X has full row rank almost surely, since each successive Gaussian row avoids the span of the preceding rows.
For each fixed b, the identity R(X,Xz+b)=\liminf_m R_{1/m}(X,Xz+b) holds by the definition of the single function R. Fatou’s lemma [15, Corollary 1.4.47] over X now gives (40). These are separate moment statements for each deterministic b; no simultaneous almost-sure bound in X over all offsets is asserted. Finally, if \norm{y-Xz}>\norm X_{\mathrm{op}}\rho, then every Xt with t\in H is at distance at least \norm{y-Xz}-\norm X_{\mathrm{op}}\rho>0 from y. The averages in (42) are eventually zero, proving (41).
The explicit choice (44) is the feature needed for the next step. A density chosen only up to a Lebesgue-null set in y for each X would not justify the moment estimate at the matrix-dependent point Xz+b. An integrated L^q(dX\,dy) estimate controls a density only up to product-null sets, while the graph y=Xz+b has zero product measure. Positive cap domination instead controls the measure of every Borel set by a positive sum of cap measures. Neither conclusion supplies the specified-offset moment in Lemma 5.3; its anchored determinant estimate and explicit limit are therefore kept as a separate argument.
5.3 Continuation probabilities and the terminal potential
We apply the specified-offset estimate to the finite-state experiment of Section 2. Fix nonnegative integers M,T and an accuracy 0<\epsilon\le1/10, and write
The success probability averages over the uniform signal, Gaussian features, and learner randomness. The factor W will count the retained states together with their possible terminal stopping indices.
Proposition 5.4 (Riesz success estimate). There is an absolute constant C_0\ge1 such that every learner just specified satisfies
This includes T=0. No asymptotic assumption on M is needed for this estimate.
We first bound the potential of each terminal output law. To control its increase across a block, we apply Lemma 5.3 at Xz+b and then integrate over b. The factor 1+L in (45) counts the distance shells and remains in the final parameter argument.
Use the fixed-seed and Borel-kernel convention of Section 2: first fix an admissible shared seed, under which the initial state is independent of the signal and rows, and then apply Lemma 2.6 to choose Borel rules. These rules define continuations from every state. All bounds below are uniform in the fixed rules; after proving the conditional bound, we average the original conditional success probabilities over the seed.
Complete unused rules and terminal laws arbitrarily, and force a halt by T. After a halt, retain a tag consisting of the stopping index and terminal state and ignore further samples. At any time there are at most 2^M live states and at most (T+1)2^M such tags, including a possible immediate halt. Hence the stated bound W applies. Add independent unused Gaussian rows after T to complete the last block of length k. Let \mathcal V_i be the state set at the boundary after block i, 0\le i\le N; thus \lvert\mathcal V_i\rvert\le W.
For every v\in\mathcal V_i and s\in S^{d-1}, let h_i(v,s) be the success probability of the remaining execution started at that boundary state. This definition also applies to states not reached by the original run. At i=N it is the success probability of the fixed output law associated with the terminal state. These functions are Borel and take values in [0,1]: the terminal assertion follows by integrating the Borel success indicator against the output law, and the assertion at earlier boundaries follows inductively from the formula below. If X is a fresh k by d standard Gaussian matrix, then for 1\le i\le N
Here K^{(i)}_{v,w}(X,y) is the probability of the block transition from v to w on the exact row pairs (X,y). These kernels are Borel, nonnegative, and sum to one in w for every (X,y). They incorporate internal randomization and the ignored rows after a halt. The functions h_i(w,\cdot) are fixed independently of X. The fresh Gaussian block is independent of the signal, initialization, past rows, and past coins. Starting the continuation at the specified state therefore gives exactly (46); the state and boundary index contain all retained sample-dependent information.
Lemma 5.5 (Terminal potential). Let Q be a signal-independent probability law on unit-vector outputs, and let h(s)=Q\{u:\arccos\inner{u}{s}\le\epsilon\}. Then
Proof. Lemma 2.3 gives \int h\,d\sigma\le\epsilon^{d-1}, also for a randomized output law. For any ambient center z, use h\le1 and (36) inside B_\epsilon(z). Outside that ball, bound the kernel by \epsilon^{-\alpha} and use the total mass bound. The sum is at most C^d\epsilon^{d-1-\alpha} uniformly in z, proving (47).
Set
The maximum is finite by Lemma 5.1. Lemma 5.5 gives V_N\le C^d\epsilon^{d-1-\alpha}. We next estimate the contribution of a single distance shell in (46).
5.4 One block and the distance shells
Lemma 5.6 (One shell across a block). For 1\le i\le N, v\in\mathcal V_{i-1}, z\in\R^d, and \rho>0, let H_\rho(z)=\{s\in S^{d-1}:\rho/2<\norm{s-z}\le\rho\}. Then
Proof. For each successor w, use the particular density R_w(X,y) of Lemma 5.3 for the measure h_i(w,\cdot)\mathbf1_{H_\rho(z)}\,d\sigma. Tonelli’s theorem, (46), and pushforward integration give
The pushforward identity used in the first inequality holds for every full-row-rank X, hence almost surely. The second inequality uses the probability-vector property of K^{(i)} and (41).
We need a deterministic integration region after averaging over X. Choose a sufficiently large absolute D_0 and partition the matrices into
These bins obey
For completeness, take 1/4-nets on the unit spheres of \R^d and \R^k. A maximal separated set has at most 9^d and 9^k points, respectively, by the disjoint-ball volume argument. Approximating each unit vector by its net point incurs an error at most \norm X_{\mathrm{op}}/2 in a bilinear form. Hence \norm X_{\mathrm{op}} is at most twice the largest absolute bilinear form on the two nets. Each such form is standard Gaussian, so its elementary tail bound and a union bound give
Putting t=D_0(j-1)\sqrt d proves (50) after increasing D_0, since k\le d/8.
For every deterministic b\in\R^k, Hölder’s inequality and (40) give
No conditional version of the density lemma is being used: Hölder’s inequality pays for the bin event. For a fixed matrix in E_j, translate y=Xz+b in (49); the remaining region is \norm b\le D_0j\sqrt d\,\rho. Joint measurability permits Tonelli’s theorem to integrate (52) over this deterministic ball.
The volume of that ball is at most C^d j^k\rho^k. Indeed, the unit-ball bound (2) gives v_k\le(2\pi e/k)^{k/2}. Since d/k\le16, v_k(D_0j\sqrt d\,\rho)^k\le C^d j^k\rho^k. Also, by q\ge2 and (50),
For the last bound use k\log j\le(d/8)(j-1)^2 for j\ge2 and sum the resulting geometric upper bound. Combining the volume and bin estimates with (49) proves (48). The powers of the shell radius cancel exactly:
Proof of Proposition 5.4. Outside B_\epsilon(z) use the shells
For every ambient z, at most C(1+L) of these meet the sphere. If \norm z\le3, all distances are at most 4, giving at most \lceil\log_2(4/\epsilon)\rceil nonempty shells. If \norm z>3, the distances lie between \norm z-1 and \norm z+1, whose ratio is less than 2; only a bounded number of dyadic shells can then meet that interval.
Inside B_\epsilon(z) use h_{i-1}\le1 and (36). On each remaining shell use Lemma 5.6. Taking the supremum over z and the maximum over v gives
Choose C_0 large enough that the first term here and the terminal value V_N are at most B in (45), and that the second coefficient is at most J-1. With U_i=\max\{V_i,B\}, equation (53) gives U_{i-1}\le J U_i, while U_N=B. Hence V_0\le BJ^N, also when N=0.
At z=0 the kernel in (34) equals one on the unit sphere, so \int h_0(v,s)\,d\sigma(s)\le V_0 for each initial state. For almost every fixed shared seed, averaging over its conditional initial-state law preserves the bound, because that law is independent of the signal and rows. Averaging the original conditional success probabilities over the seed then preserves it again. This proves (45).
5.5 The logarithmic lower bound
We now prove Theorem 2.4 from the Riesz success estimate. The parameter argument retains the stopping-index factor in W and the shell-count factor 1+L. Terminal-tag counting controls both in the only case still requiring proof, namely T<dL.
Proof of Theorem 2.4 by the Riesz method. There is a useful bound that counts possible terminal tags. For fixed rules as above, write p_{\mathrm{succ},\mathrm{rule}} for the conditional success probability. For a tag a, let f_a(s) be the probability of reaching it given s, and let g_a(s) be its output law’s success probability. This output law depends only on the tag, and 0\le f_a\le1. Lemma 2.3 gives \int g_a\,d\sigma\le C^d\epsilon^{d-1}, including a randomized output law. Sum over at most W tags to obtain
The right side is independent of the chosen rules, so averaging also gives p_{\mathrm{succ}}\le W C^d\epsilon^{d-1}.
Write L=\log(1/\epsilon). If T\ge dL, the conclusion already holds. In the remaining case T<dL, assume p_{\mathrm{succ}}\ge2/3. Taking logarithms in the averaged form of (54) gives
For the given memory sequence, eventually M\le d^2. Since L\ge\log10 and T<dL, we also have
for large d. It follows that (d-2)L\le d^2\log2+d\log C+\log(2d)+\log(3/2). Consequently, only in this remaining case,
Thus the multiplier from Proposition 5.4 obeys
using q\ge d/16. In particular the shell-count term has been bounded rather than removed.
Now (45) and p_{\mathrm{succ}}\ge2/3 imply
because N=\lceil T/k\rceil\le16T/d+1. The common full-data bound, Lemma 2.7, gives T>d/8; hence the additive d on the right of (56) is at most 8T. For d\ge4 the left coefficient is at least d/4, so T\ge c\,dL for an absolute c>0. This combines with the case T\ge dL already separated.
References
- [1]Matthew T. Calef and Douglas P. Hardin. Riesz s-Equilibrium Measures on d-Rectifiable Sets as s Approaches d. 2008. arXiv:0808.3802v1, 28 August 2008; the potential definition is in Section 1, p. 2. https://arxiv.org/abs/0808.3802v1.
- [2]Yuval Dagan, Gil Kur, and Ohad Shamir. Space Lower Bounds for Linear Prediction in the Streaming Model. In Proceedings of the Thirty-Second Conference on Learning Theory, pages 929–954, 2019. Theorem locators refer to the published PMLR version. https://proceedings.mlr.press/v99/dagan19b.html.
- [3]S. W. Drury. Generalizations of Riesz Potentials and L^p Estimates for Certain k-Plane Transforms. Illinois Journal of Mathematics, 28(3):495–512, 1984.
- [4]Peter Frankl and Hiroshi Maehara. Some Geometric Applications of the Beta Distribution. Annals of the Institute of Statistical Mathematics, 42(3):463–474, 1990. DOI: https://doi.org/10.1007/BF00049302.
- [5]D. H. Fremlin. Measure Theory, Chapter 21: Taxonomy of Measure Spaces. University of Essex, author-hosted online development extract. Section 211 version 20 November 2003; accessed 27 September 2026. https://www1.essex.ac.uk/maths/people/fremlin/chap21.pdf.
- [6]D. H. Fremlin. Measure Theory, Chapter 24: Function Spaces. University of Essex, author-hosted online development extract. Chapter version 15 November 2013; Section 243 version 30 April 2004; accessed 27 September 2026. https://www1.essex.ac.uk/maths/people/fremlin/chap24.pdf.
- [7]Tiefeng Jiang. How Many Entries of a Typical Orthogonal Matrix Can Be Approximated by Independent Normals?. The Annals of Probability, 34(4):1497–1529, 2006. Lemma locators refer to the published-article reprint, arXiv:math/0601457v2, 25 September 2006. https://arxiv.org/abs/math/0601457v2. DOI: https://doi.org/10.1214/009117906000000205.
- [8]Pertti Mattila. Hausdorff Dimension, Orthogonal Projections and Intersections with Planes. Annales Academiae Scientiarum Fennicae, Series A I, Mathematica, 1(2):227–244, 1975. DOI: https://doi.org/10.5186/aasfm.1975.0110.
- [9]Ran Raz. Fast Learning Requires Good Memory: A Time-Space Lower Bound for Parity Learning. 2016. arXiv:1602.05161v1, 16 February 2016. https://arxiv.org/abs/1602.05161v1.
- [10]Roland Schnaubelt. Functional Analysis: Lecture Notes of Winter Semester 2017/18. Karlsruhe Institute of Technology, author-hosted lecture notes, 2025. Title-page version 21 July 2025; accessed 27 September 2026. https://iana.math.kit.edu/downloads/iana3/schnaubelt/Skripten/fa17-skript.pdf.
- [11]Ohad Shamir. Fundamental Limits of Online and Distributed Algorithms for Statistical Learning and Estimation. In Advances in Neural Information Processing Systems, pages 163–171, 2014. https://proceedings.neurips.cc/paper_files/paper/2014/hash/cc427d934a7f6c0663e5923f49eba531-Abstract.html.
- [12]Vatsal Sharan, Aaron Sidford, and Gregory Valiant. Memory-Sample Tradeoffs for Linear Regression with Small Error. In Proceedings of the 51st Annual ACM SIGACT Symposium on Theory of Computing, pages 890–901, 2019. Full version arXiv:1904.08544v2, 10 October 2020; theorem and section locators refer to that version. https://arxiv.org/abs/1904.08544v2. DOI: https://doi.org/10.1145/3313276.3316403.
- [13]Jacob Steinhardt and John Duchi. Minimax Rates for Memory-Bounded Sparse Linear Regression. In Proceedings of the 28th Conference on Learning Theory, pages 1564–1587, 2015. https://proceedings.mlr.press/v40/Steinhardt15.html.
- [14]Jacob Steinhardt, Gregory Valiant, and Stefan Wager. Memory, Communication, and Statistical Queries. In 29th Annual Conference on Learning Theory, pages 1490–1516, 2016. https://proceedings.mlr.press/v49/steinhardt16.html.
- [15]Terence Tao. An Introduction to Measure Theory. Graduate Studies in Mathematics 126, American Mathematical Society, 2011. Theorem and page locators refer to the author-hosted preliminary version; accessed 27 September 2026. https://terrytao.wordpress.com/wp-content/uploads/2012/12/gsm-126-tao5-measure-book.pdf.