Alphabeta Math
TheoremStatement: AI-adaptedProof: AI-adaptedPipeline-generatedprecheck passjudge pass (gpt-6-luna)
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Sparse constrained least-squares risk bound

Statement

For every admissible input to Constrained support-search least squares, its finite scan terminates and returns a support minimizing ∥y−PVSy∥22 over all supports of size at most s. Among tied minimizers, it returns the earliest support in the prescribed order.

Moreover, let n,d be positive integers, let X∈Rn×d be fixed, and let θ∗∈Rd have at most s nonzero coordinates, where 1≤s≤d/2. Suppose Y=Xθ∗+ε and ε is centered sub-Gaussian with proxy σ2≥0 in the linear-functional sense of Sub-Gaussian variables and vectors, and sub-exponential variables. Let μ^s be the fitted vector returned by Constrained support-search least squares. For 0<δ<1, with C:=32(4+2log⁡6),

E[1n∥μ^s−Xθ∗∥22]≤Cσ2slog⁡(ed/s)n,

and

P ⁣(1n∥μ^s−Xθ∗∥22≤Cσ2[slog⁡(ed/s)+log⁡(1/δ)]n)≥1−δ. No design-rank, incoherence, or column-normalization condition is assumed. The claim concerns fitted values; it gives no coefficient-error bound.

Facts & Assumptions

Given: For the scan assertion, any admissible input listed in the Definition of Constrained support-search least squares. For the risk bounds, the fixed-design model and noise assumptions in the statement, with 1≤s≤d/2 and 0<δ<1.

[F1]

The procedure enumerates the finite nonempty list of supports of size at most s in its prescribed order and replaces its incumbent only when it finds a strictly smaller residual (Constrained support-search least squares, ∣P(A)∣=2∣A∣ for finite A).

[F2]

For a fixed subspace of rank r and centered sub-Gaussian noise with proxy σ2>0, P(∥PVε∥22>t)≤2⋅6re−t/(8σ2) for t>0; proxy zero gives zero projection error (Sub-Gaussian projection risk bounds).

[F3]

For natural numbers 1≤m≤d, (dm)≤(ed/m)m (A binomial support-count bound).

[F5]

The Euclidean inner product satisfies Cauchy--Schwarz (The Cauchy-Schwarz inequality for finite sums).

[F6]

Orthogonal projection has its component in the subspace and its residual orthogonal to that subspace (The orthogonal projection PWv is the W-component in V=W⊕W⊥).

[F7]

Probability is countably subadditive (Basic identities for a probability measure).

[F8]

Risk is the expected loss of a decision rule under the sampling law (Loss functions and risk).

[F9]

The selected support index and fitted-vector maps are Borel measurable for this finite score family (Finite tie-broken projection searches are Borel measurable).

[F10]

Orthogonal projection minimizes distance to its subspace (The orthogonal projection is the unique nearest point in the subspace).

[F11]

The natural logarithm is the inverse of exp⁡, is strictly increasing, has the product, quotient, and reciprocal laws and log⁡1=0; repeated application gives its integer-power law, and e=exp⁡(1) (The natural logarithm as the inverse of the exponential function, Order, continuity, range, and the product, quotient, and reciprocal laws for the natural logarithm, The real exponential function and the number e by a power series).

[F12]

For T≥0, c>0, x0≥0, h=c/2, and N∈N, min⁡(T,x0+Nh)≤x0+h∑j<N1{T>x0+jh}. If P(T>x)≤Aexp⁡(−x/c) and x0=clog⁡A, monotone convergence and expectation linearity, with the Archimedean property giving x0+Nh→∞, bound ET by x0+h∑j≥0exp⁡(−j/2). The exponential series gives exp⁡(1/2)≥13/8. Set q=exp⁡(−1/2)=1/exp⁡(1/2)≤8/13<1; the addition law gives exp⁡(−j/2)=qj by induction, and the geometric-series formula gives ∑j≥0exp⁡(−j/2)=1/(1−q)≤13/5 (Monotone convergence for the integral, Every complete ordered field is Archimedean, Linearity, monotonicity, and the modulus bound for expectation, The exponential addition formula exp⁡(x+y)=exp⁡(x)exp⁡(y), The exponential is positive and satisfies exp⁡(−x)=1/exp⁡(x), The real exponential function and the number e by a power series, Integer powers am, For ∣r∣<1, ∑k≥0rk=1/(1−r), and for ∣r∣≥1 the series diverges).

[F13]

Expectation is linear for integrable random variables (Linearity, monotonicity, and the modulus bound for expectation).

Proof

1.1F1algebra

Let S1,…,SM be the finite nonempty ordered candidate list from [F1]. After the first support is processed, the incumbent has least residual and earliest index among minimizers in the first one-element list. If this is true after j<M supports, then the next residual either is strictly smaller, in which case replacing the incumbent gives the new least residual, or is at least the incumbent residual, in which case retaining it preserves the minimum and keeps the earlier index in a tie. Induction proves the invariant through M, so the scan terminates and returns a least-residual support with the stated tie rule.

1.2F1F10given

Write μ=Xθ∗, let S∗={j:θj∗≠0}, and put Δ=μ^s−μ. Since ∣S∗∣≤s, S∗ is a candidate support and μ∈VS∗. Its candidate fitted vector is PVS∗Y. By [F10], ∥Y−PVS∗Y∥22≤∥Y−μ∥22.

2.1step 1.1step 1.2algebra

By the scan-minimizer result in step 1.1, the selected residual is no larger than that of PVS∗Y, which step 1.2 bounds by the residual at μ: ∥Y−μ^s∥22≤∥Y−μ∥22=∥ε∥22. Expanding Y=μ+ε gives ∥Δ∥22≤2⟨ε,Δ⟩.

3.1step 2.1F5F6algebra

If Δ≠0, its support difference lies in the span of columns indexed by S∗∪S^, whose size is at most 2s. Since 2s≤d, that set is contained in some T⊆{1,…,d} with ∣T∣=2s. Thus Δ∈VT, and orthogonality of the projection gives ⟨ε,Δ⟩=⟨PVTε,Δ⟩ by [F6]. By [F5] and step 2.1, ∥Δ∥22≤2∥PVTε∥2∥Δ∥2, hence ∥Δ∥22≤4∥PVTε∥22. If Δ=0, the last inequality holds for every such T as well.

4.1F1F2F7F9step 1.1step 1.2algebrastep 3.1

There are (d2s) possible sets T. For each, its projection rank is at most 2s, so [F2] and the finite union bound give, for t>0 and σ2>0, P(∥Δ∥22>4t)≤2(d2s)62se−t/(8σ2)=:Ae−t/(8σ2). If σ2=0, apply [F2] to the identity subspace: for every integer j≥1, P(∥ε∥22>1/j)=0. Since {∥ε∥22>0}=⋃j≥1{∥ε∥22>1/j}, [F7] gives ε=0 almost surely. Then [F1] and the scan result in step 1.1 show that the selected residual is zero, so its fitted vector equals Y=μ; both claims hold with zero error.

5.1F3F9F11step 4.1algebra

For σ2>0, take t=8σ2[log⁡A+log⁡(1/δ)] in step 4.1. By [F11] the threshold is positive. The resulting tail is at most δ, and log⁡A=log⁡2+2slog⁡6+log⁡(d2s). By [F3] and [F11], log⁡(d2s)≤2slog⁡(ed/(2s))=2slog⁡(ed/s)−2slog⁡2. Hence log⁡A≤2slog⁡(ed/s)+2slog⁡6+(1−2s)log⁡2≤(2+2log⁡6)slog⁡(ed/s). Dividing the threshold 4t by n and using C=32(4+2log⁡6) proves the stated high-probability bound.

6.1F3F8F9F11F12F13step 4.1step 5.1algebra∎

Put Z=∥Δ∥22, c=32σ2, h=c/2, and x0=clog⁡A. Step 4.1 gives P(Z>x0+jh)≤exp⁡(−j/2) for every j≥0. Apply [F12] to Z, using [F13] for expectation linearity and monotone convergence, to get EZ≤x0+h∑j≥0exp⁡(−j/2)≤c(log⁡A+13/10). By step 5.1, log⁡A≤(2+2log⁡6)slog⁡(ed/s); since slog⁡(ed/s)≥1 and 13/10<2, this yields EZ/n≤Cσ2slog⁡(ed/s)/n. By [F8], this is the expected risk under squared fitted-value loss. The zero-proxy case was settled in step 4.1.

Depends on

Used by

Cited to discharge well-definedness by Constrained support-search least squares.

Dependency tree · two levels

120 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources