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.
Algorithm Composition, Precision Budgets, and Dequantization
1 · Prerequisites
- Absolute and Conditional Convergence; Rearrangement; Products
- Algebraic Closure, Embeddings, and Separability
- Algebraic Extensions, Extension Degree, and Finite Fields
- Approximate Universality and Solovay–Kitaev
- Binary Operations, Monoids, Groups and Subgroups
- Bounded Linear Operators and Quotient Spaces
- Channels, Kraus, Choi, and Stinespring Representations
- Circuit Syntax, Denotational Semantics, and Measurement
- Compactness in Metric Spaces
- Completeness, Completion, and Uniform Continuity
- Complex Coordinate and Path Transport
- Complex Differentiability and the Cauchy–Riemann Equations
- Complex Norms and Realification
- Complex Power Series and Analytic Functions
- Complex Subspace Realification
- Complexification, Realification and Real Structures
- Composition Series, the Jordan–Hölder Theorem and Solvable Groups
- Congruences, the Integers Modulo n and the Chinese Remainder Theorem
- Connectedness
- Construction of the Natural Numbers
- Construction of the Real Numbers via Cauchy Sequences
- Construction of the Real Numbers via Dedekind Cuts
- Continuity, IVT, EVT, and Uniform Continuity
- Controlled Gates, Phase, and Entanglement
- Cosets, Index and Lagrange's Theorem
- Countability and Uncountability
- Cyclic Groups and Direct Products
- Determinants of Matrices over a Commutative Ring
- Diagonalisation and the Minimal Polynomial
- Divisibility, Euclidean Domains, Principal Ideal Domains and Unique Factorisation
- Divisibility, Greatest Common Divisors and Bézout's Identity
- Dual Spaces, Bilinear and Quadratic Forms, and Sylvester's Law of Inertia
- Eigenvalues, Eigenvectors and the Characteristic Polynomial
- Filters and Ultrafilters
- Finite Counting, Factorials and Binomial Coefficients
- Finite Dimensional Normed Spaces and Riesz Lemma
- Finite Probability and the Probabilistic Method
- Finite Probability Spaces and Random Variables
- Finite-Dimensional Complex Inner-Product Spaces and Matrix Interfaces
- Formal Languages, Encodings, and Decision Problems
- Foundations of the Real Numbers for Analysis
- Free Modules, Exact Sequences, Projective and Injective Modules
- Fundamental Trigonometric Identities
- Group Actions, Orbits, Stabilisers and Cayley's Theorem
- Group Homomorphisms and the Isomorphism Theorems
- Hilbert Space Geometry and Riesz Representation
- Ideals, Quotient Rings and the Isomorphism Theorems for Rings
- Independence Borel Cantelli and Zero One Laws
- Infinite Product Measures and Kolmogorov Extension
- Inner Product Spaces, Gram-Schmidt, Projections and Adjoints
- Limits of Real Functions
- limsup, liminf, and Subsequential Limits
- Linear Independence, Bases and Dimension
- Linear Recurrences and Rational Generating Functions
- Linear Transformations, Rank-Nullity and Quotient Spaces
- Matrices, the Matrix of a Linear Map, and Change of Basis
- Matrix Norms, Condition Numbers and Numerical Stability
- Measurable Functions and Simple Approximation
- Measurements, POVMs, and Instruments
- Measures and Their Basic Properties
- Metric Spaces
- Modules, Submodules, Quotient Modules and the Isomorphism Theorems
- Monotone Functions, Discontinuities, and Continuity Sets
- Monotone Sequences, Bolzano-Weierstrass, and Cauchy Completeness
- Normal Subgroups and Quotient Groups
- Normed and Banach Spaces
- Order, Zorn's Lemma, and the Axiom of Choice
- Outer Measure and the Caratheodory Extension Theorem
- Partitions of Unity and Paracompactness
- Polynomial Rings, the Division Algorithm and Roots
- Power Series and Real-Analytic Functions
- Primes, Euclid's Lemma and the Fundamental Theorem of Arithmetic
- Product Measures and the Fubini Tonelli Theorems
- Properties of the Integral and the Working FTC
- Pure and Mixed States, Ensembles, and Schmidt Decomposition
- Qubits, Bloch Sphere, Pauli Algebra, and One-Qubit Gates
- Registers, Tensor Products, and Partial Trace
- Relations, Functions, and Quotients
- Resource Bounds and Machine Invariance
- Reversible Computation, Ancillas, and Coherent Oracles
- Rings, Subrings, Integral Domains and Fields
- Rⁿ as a Normed Space; Vector-Valued Functions
- Roots, Rational Powers, and Classical Inequalities
- Separation Axioms: the Hierarchy
- Sequences and Limits
- Sequences and Series of Functions; Uniform Convergence
- Series: Convergence and the Nonnegative Tests
- Sigma Algebras and Borel Sets
- Simple Field Extensions and the Construction of the Complex Numbers
- Sine, Cosine, and the Definition of Pi
- Splitting Fields
- State Preparation, Isometry Synthesis, and Limits
- Strongly Continuous Semigroups and Hille Yosida
- Subspaces, Products, and Quotients
- Suprema and Infima
- Sylow's Theorems, p-Groups and Nilpotent Groups
- Symmetric Groups, Cycle Decomposition and the Sign Homomorphism
- Tensor Products of Modules
- The Cantor Set, Baire Category, and Measure Zero in ℝ
- The Complex Exponential and Euler's Formula
- The Derivative and the Mean Value Theorems
- The Determinant of a Linear Operator, Cofactors and Cramer's Rule
- The Exponential Function
- The Fundamental Theorem of Algebra
- The Fundamental Theorem of Finite Abelian Groups
- The Galois Correspondence
- The Lebesgue Integral and the Convergence Theorems
- The Logarithm and General Powers
- The Polynomial Hierarchy and Relativization
- The Riemann Integral: Definition and Integrability
- The Spectral Theorem, Positive Operators and Singular Value Decomposition
- The Topology of Euclidean Space
- The Total Derivative in ℝᵐ → ℝⁿ
- The ZFC Axioms and the Basic Set Constructions
- Topological Spaces and Continuity
- Topology of ℝ
- Trace Distance, Fidelity, and Information Limits
- Triangularisation, Generalised Eigenspaces and Jordan Canonical Form
- Turing Machines, Configurations, and Computation
- Uniform Circuit Families, Promise Problems, and Complexity
- Vector Spaces, Linear Subspaces, Span and Direct Sums
2 · Summary
This page defines end-to-end channel and output contracts and proves product-state and tensor-channel error bounds, including final-measurement variation for arbitrary reference inputs. The finite hybrid argument proves the composed-channel bound using the current channel-norm interface. The page allocates precision and failure budgets, composes declared cost coordinates, and separates classical output, quantum-state output, input-access costs, and heralded success. Its vector sampling interface uses squared moduli for complex vectors; the ModFKV and Tang matrix algorithms specialize to real matrices. The ModFKV approximation is proved from exponential Gram concentration and a deterministic threshold-band argument. Tang's recommendation sampler uses a computable adaptive row-refinement rule, an explicit finite accepted branch, and a countable probability model for its almost-sure termination claim.
3 · Logical flowchart
4 · Definitions, theorems and proofs
Threshold-band singular-value approximation
Definition
Let have a singular value decomposition as in Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition, let , and let . Write for the orthogonal projection onto the span of right singular vectors whose singular values are at least , and for the span of right singular vectors whose singular values lie in . For any orthogonal projection with range contained in , a threshold-band target is
The choice inside the threshold band is allowed to vary; all singular directions above the band are retained and all directions below it are discarded. The definition is independent of the choice of orthonormal bases within singular subspaces because it is stated using their subspaces. If a singular value equals the upper endpoint it belongs to ; equality at the lower endpoint belongs to . A range contained in contributes zero to . However, an allowed can mix nullspace directions with nonzero singular directions when both lie in , so changing that projection can change the target.
Truncated singular values give the best Frobenius rank-k approximation
Statement
Let have singular values , padded by zeros to . For an integer let be the truncation to the first singular terms of an SVD of . Then , and for every matrix with ,
Facts & Assumptions
Given: The finite real matrix and integer in the Statement.
Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition: the full SVD supplies orthonormal bases of and of , with and for zero singular values; in particular, the right singular vectors for positive singular values are orthogonal to and can be completed within .
The Frobenius norm on real or complex matrices: the Frobenius norm is the Euclidean norm of the full entry list.
The orthogonal projection is the -component in : every finite-dimensional subspace has an orthogonal projection onto it.
Proof
By [F1], write for its positive singular values, and set . Its image lies in the span of at most left singular vectors, so its rank is at most . The residual is the sum of the positive singular terms with index greater than ; padding the singular-value list by zeros, [F2]–[F3] give .
Let be the orthogonal projection onto the row space of , which exists by [F4]. Its rank is at most , every row of lies in , and hence . Rowwise, , where the first summand has each row in and the second has each row in . These subspaces are orthogonal, so summing the rowwise Pythagorean identities gives .
Let be the right singular vectors for the positive singular values, and complete them to an orthonormal basis of using vectors in ; set for . The SVD gives , so orthonormality of and [F2] give , where . If is an orthonormal basis of , then Parseval in the basis gives . For , every . For , put and use for and for to obtain Thus in either case.
Substitute the bound from step 1.3 into step 1.2 and use from [F3]. This gives , proving the claim.
Tensor-product state errors add under trace distance
Statement
For a positive integer and density operators on finite-dimensional , , Consequently every finite POVM on the product output changes its outcome law by at most this sum.
Facts & Assumptions
Given: The finite density-operator pairs in the Statement.
Density operator: each state is nonnegative of trace one.
Measurement Variation At Most Trace Distance: a finite POVM changes total variation by at most the input trace distance.
A non-negative operator is equivalently self-adjoint with nonnegative eigenvalues, a positive semidefinite matrix in an orthonormal basis, or an operator of the form S^*S: every finite-dimensional nonnegative self-adjoint operator has an orthonormal eigenbasis with nonnegative eigenvalues.
Tensor products of finitely many linear maps: tensor maps act factorwise on product vectors; Ordered tensor registers, their inner product, and wire permutations supplies the finite product-basis convention.
Trace Norm Duality And Norm Properties: the trace norm satisfies the triangle inequality.
Trace Norm: the trace norm is the sum of singular values.
The basis-independent trace of an endomorphism of a finite-dimensional vector space: the trace is the sum of diagonal entries in any orthonormal basis.
Proof
By [F5], write each in an orthonormal eigenbasis with nonnegative eigenvalues. By [F9] their sum is . The tensor-product basis from [F6] diagonalizes with eigenvalues all products of those nonnegative eigenvalues. Their sum is , so the product is a density operator; the same argument applies to .
Add and subtract the products with their first factors changed, for , to obtain This is a finite telescoping identity.
Apply [F7] to the trace-norm triangle inequality for the finite sum in step 1.2. By [F1], [F5], [F8], and [F9], each unchanged state factor has trace norm one; by [F3], the th summand has trace norm exactly . Divide by and use [F2] to obtain the asserted sum.
Apply [F4] to the product output states verified in step 1.1 and the declared finite POVM. The resulting outcome total variation is at most their trace distance, hence at most the sum from step 2.1.
End-to-end channel and output error contract
Definition
A finite -stage composition contract consists of nonzero finite-dimensional complex Hilbert spaces and, for each , an ideal channel and an implemented channel with the same domain and codomain. Channels are completely positive and trace-preserving on the full operator space, as in Finite-dimensional quantum channel. Products act right to left, so applies stage first, following Denotational semantics of a finite layered circuit.
For channels on the same nonzero finite-dimensional complex Hilbert input and output spaces define the normalized diamond distance where is the unhalved supremum of induced trace norms over all finite reference spaces, as in Diamond Norm. For density operators define , as in Trace Distance. For probability distributions on a common finite output alphabet define .
The local channel error at stage is . The final channel error is . A final output-distribution error is total variation distance between the finite laws obtained by applying a declared final POVM to the two output states. A failure probability is the probability of a declared failure event in the finite execution history. These four quantities have different definitions and must be reported separately.
When an input is entangled with a finite reference register, include that reference in the density operator of Density operator and retain it through the channel call; the diamond distance already takes the supremum over such references. Channel, state, and classical output contracts are distinct. Cost coordinates for input description and preprocessing, state preparation, gates, queries including controlled and inverse access, workspace, measurement, postprocessing, and classical output are specified separately by Classical output, quantum-state output, and access costs and Circuit-family and computation cost model.
Hybrid bound for a composed channel circuit
Statement
Let , let be nonzero finite-dimensional complex Hilbert spaces, and let be channels for . Then where is the normalized distance in End-to-end channel and output error contract. For any finite reference space , input density operator on , and finite final POVM on , the total variation distance between the two final outcome laws is at most the same sum. The claim is uniform over the reference dimension; no product-input assumption is made.
Facts & Assumptions
Given: The finite channel sequence and the finite input/reference/measurement specified in the Statement.
End-to-end channel and output error contract: , where the diamond norm takes the supremum over every finite reference space.
Channel Composition Hybrid Bound: for two matching pairs of channels between nonzero finite-dimensional complex Hilbert spaces, replacing both stages changes the raw diamond norm by at most the sum of the two local raw diamond errors.
Finite-dimensional channels compose and tensor: the composition of finite-dimensional channels is a channel.
Measurement Variation At Most Trace Distance: the output total variation of a finite POVM is at most the trace distance between its input density operators.
Trace Norm Duality And Norm Properties: the trace norm satisfies the triangle inequality.
Density operator: a density operator is positive semidefinite and has trace one.
Trace Norm: the trace norm is the sum of singular values.
A non-negative operator is equivalently self-adjoint with nonnegative eigenvalues, a positive semidefinite matrix in an orthonormal basis, or an operator of the form S^*S: positive semidefinite operators have orthonormal eigenbases with nonnegative eigenvalues; those eigenvalues are their singular values.
The basis-independent trace of an endomorphism of a finite-dimensional vector space: the trace is the sum of diagonal entries in any orthonormal basis.
For every nonzero finite-dimensional space , is a channel: each block amplification fixes every operator, so it is completely positive, and it preserves trace because it is the identity (Completely positive maps, Finite-dimensional quantum channel).
Proof
For any two maps with matching domains, the triangle inequality for the induced trace norm holds at each finite reference space by [F5]. Taking the supremum over reference spaces in [F1] gives .
For , let and . If either product is empty, it is the identity on the corresponding nonzero finite-dimensional space and is a channel by [F10]; otherwise its factors compose to a channel by [F3]. Apply [F2] first to the identical pair followed by , and then to the pair followed by the identical pair . After dividing by , the two applications give The finite telescoping identity follows by changing stages one at a time.
Apply the diamond-norm triangle inequality from step 1.1 to the finite sum in step 1.2. Each summand has normalized distance at most its local value, so the final normalized diamond distance is at most .
The identity is a channel by [F10], so tensor-channel closure [F3] makes and channels. Their complete positivity at ancilla dimension one sends the nonnegative to nonnegative operators, and trace preservation keeps trace one, so the outputs are density operators by [F6] and the channel definition. By the supremum in [F1], [F7]–[F9] give because .
Apply [F4] to and the declared final POVM. Its outcome-law total variation is at most , hence at most the sum from step 2.1.
Tensor-product channel error bound
Statement
For a positive integer and a finite family of matching channel pairs , , on nonzero finite-dimensional complex Hilbert spaces, For every nonzero finite-dimensional reference space , every density operator on , and every finite POVM on , the total variation distance between the outcome laws after applying and to is at most the same sum. No product-input assumption is made; may be entangled across the component inputs and reference.
Facts & Assumptions
Given: The finite channel pairs, product spaces, nonzero finite reference space , input density operator , and final finite POVM specified in the Statement.
End-to-end channel and output error contract: and each include a factor ; the raw diamond norm does not.
Channel Tensor Error Bound: for two matching channel pairs, the raw diamond distance between tensor products is at most the sum of the raw local diamond errors; the diamond norm is multiplicative under tensor products.
Finite-dimensional channels compose and tensor: tensor products and compositions of finite-dimensional channels are channels.
Hybrid bound for a composed channel circuit: for matching finite-dimensional channels, every finite reference input and final POVM has outcome total variation at most the normalized diamond distance between the channels.
Proof
For the channel inequality is equality by [F1].
Assume the bound for factors. By [F3], the first channel factors in each tensor product form channels. Apply [F2] to this prefix channel pair and the final channel pair. Dividing the raw diamond inequality by gives the normalized distance of the full tensor products at most the prefix normalized distance plus . The induction hypothesis bounds the prefix term by the first summands.
Steps 1.1 and 1.2 prove the channel inequality for every finite by induction.
By [F3], the product maps and are channels. Apply [F4] to the one-stage pair , with the declared reference , input , and final POVM. Its outcome total variation is at most , which step 2.1 bounds by .
Allocate local error budgets through a pipeline
Definition
Inputs: a finite list of stage compilers , each accepting a rational local precision and, when it uses randomized verification, a rational local failure budget; rational targets ; and declared stage input/output interfaces , including encodings or register spaces. Stage must accept and return at every admissible budget, so adjacent interfaces match. For a channel pipeline these are the matching spaces and channel interfaces of End-to-end channel and output error contract. Each compiler is assumed to terminate with an implemented stage, its output type, and its actual cost record under Classical output, quantum-state output, and access costs. The sequential pipeline is run under the same declared cost-accounting model.
Output: a pipeline from to , together with the budget lists and cost records below. If , a single interface is supplied; return its identity procedure, empty budget lists, no stage records, and an auxiliary record charging any identity-interface and allocator work. If , set and for each . For in order, invoke once when stage uses randomized verification, and otherwise invoke once and mark its failure coordinate not applicable. Retain each returned implementation, output type, and actual stage cost record, then compose the stages chronologically. Record shared input/preprocessing, allocation arithmetic, and any interstage wiring, control, postprocessing, or final-output costs in a separate auxiliary record, counting each operation once. Return the composed implementation, local budget lists, per-stage records, and auxiliary record. The correctness theorem gives an additive resource bound from these records under its explicit accounting assumptions; a requested precision alone does not determine an implementation or its cost.
The correctness claim is separate: Allocated local errors meet the global budget proves the global channel and failure bounds and composes the recorded cost coordinates under explicit stage and pipeline accounting assumptions. This algorithm does not create implementations or amplify a compiler's error; each compiler's correctness method and cost remain part of its own guarantee.
Classical output, quantum-state output, and access costs
Definition
An algorithm result contract declares exactly one output type:
- Classical string: a finite bit string in a stated encoding. Its cost includes every output bit written and any required formatting.
- Classical sample: an element of a stated finite set, with the output law and the cost of writing the sample encoding.
- Quantum state: a density operator, as in Density operator, on an explicitly named finite register, produced by a declared state-preparation procedure. This does not mean that its amplitude list or a classical description of the state is output.
The accompanying cost record keeps separate: input encoding length and preprocessing; each data-structure or oracle operation; state preparation and its inverse or controlled forms; ordinary gates, controlled gates, and inverse gates; register width and clean/dirty/other workspace; measurement and classical postprocessing; finite-gate approximation precision; and classical output bits. A coordinate that the model does not use is marked not applicable. No state-preparation oracle, QRAM, controlled access, inverse access, or full classical state description is assigned zero cost without an explicit interface.
For the exact continuous-rotation table model, QC24 State Preparation Cost Model and QC24 State Preparation specify the table input and circuit procedure; QC24 Explicit State-Preparation Cost Guarantee gives the corresponding gate/ancilla upper bound. That bound uses exact continuous rotations and gives no finite-gate or uniform bit-time claim. For a fixed finite gate set, Solovay–Kitaev Unitary Compiler for SU(2) and Solovay–Kitaev Error and Resource Bound supply a separate approximation and resource contract. The circuit-family coordinates are those of Circuit-family and computation cost model.
Heralded success and conditional verification contract
Definition
A heralded procedure has an explicit finite outcome set containing success and failure, a success probability , and a declared output type. For an input density operator as in Density operator, its success branch is the subnormalized state of the instrument Finite instruments and conditional states, with . If , its conditional output state is ; if , no normalized conditional state is assigned and the branch is denoted . The outcome law is specified by the instrument and the POVM conventions of Finite POVMs and their outcome law.
The contract states separately (i) the success probability, (ii) the guarantee conditional on success, (iii) any unconditional guarantee including failure, and (iv) the cost of verification and of outputting or recording the success flag. An acceptance probability from Quantum acceptance probability for a finite promise computation is not itself a conditional-state guarantee. A claim about a quantum state output uses the type and cost distinctions in Classical output, quantum-state output, and access costs.
ℓ2 sample/query access to vectors and matrices
Definition
For a finite complex vector , put , where is the complex modulus of Real and imaginary parts, complex conjugation, and modulus. Its squared-entry distribution is on when . A vector has ℓ2 sample/query access when the interface supplies: (a) an entry query returning ; (b) its norm ; and (c) a sample distributed according to . For any finite adaptive sequence of sample calls, conditioned on the preceding history and the current requested sampler, each returned index has the declared law and uses fresh randomness independent of that history; in particular, repeated calls to a fixed sampler are mutually independent in the sense of Pairwise and mutual independence of finite-valued random variables. Every draw is charged. The interface states the cost and precision of every query, norm computation, and sample, and any preprocessing or storage cost. If , is undefined and the sampler returns the explicit status ZERO; it does not return a fictitious normalized sample.
For a finite real matrix , put , the Frobenius norm of The Frobenius norm on real or complex matrices. Row and column ℓ2 access supplies entry queries, row and column norms, and samplers for each nonzero row and column according to its squared-entry distribution. These real rows and columns are instances of the vector interface above. It also declares preprocessing, storage, and whether their costs are amortized. The common row sampler has law for , and the column sampler has law ; zero rows and columns have probability zero. An all-zero matrix returns ZERO for a request for a norm-normalized sample. All row-vector, column-vector, matrix-level row, and matrix-level column samplers obey the same conditional fresh-randomness rule; every draw is charged. Access to entries alone does not imply either sampler. The access and output costs are recorded separately under Classical output, quantum-state output, and access costs, and the distributions are finite probability laws as in Finite probability spaces, outcome weights, events, and event probabilities.
This is an input-data structure assumption. It is not a quantum circuit, state-preparation oracle, QRAM, or an assertion that the preprocessing needed to implement it is sublinear. The Tang matrix algorithms on this page specialize to real matrices; no complex-matrix extension of those algorithms is asserted.
Writing an N-bit output costs at least N bit operations
Statement
In an output-interface model with designated output positions, where one bit operation can write at most one position and every output position must be written, emitting an -bit classical string costs at least bit operations. This lower bound concerns the output interface; it is not a lower bound on preparing a quantum state or on producing a succinct circuit description.
Facts & Assumptions
Given: An integer and the output-interface convention in the Statement.
Classical output, quantum-state output, and access costs: a classical string output is a finite bit string in a declared encoding, and its cost includes every output bit written.
Proof
If , the empty string requires no output positions and the lower bound is . If , the output has distinct designated positions by [F1] and the stated interface.
Each bit operation writes at most one position, while all positions must be written. Hence any execution that emits the string uses at least bit operations.
ModFKV low-rank description procedure
Definition
Inputs: a nonzero real matrix with the exact row/column ℓ2 sample/query interface of ℓ2 sample/query access to vectors and matrices, a threshold , tolerances and , and an integer sample count . The procedure returns a compact description; it does not write all entries of the output matrix.
- Sample independently from the row distribution . Form the conceptual matrix whose th row is . The entries of can be queried through the stored row indices; no zero-probability row is selected.
- Define a column distribution by choosing uniformly in and then sampling a column from . Draw independently from . Form with . A sampled column has ; otherwise return INVALID.
- Compute a singular value decomposition of . Retain the left singular vectors whose singular values are strictly greater than , and set and . If , return the zero matrix description. Otherwise return the row indices and the matrices needed to evaluate
The sample count , matrix-entry and norm queries, sample operations, arithmetic, SVD, storage, and any finite-precision implementation are charged in the declared access model. Correctness, success probability, approximation error, and the parameter-dependent choice of are separate claims of ModFKV low-rank approximation supplier.
Finite rejection sampler for a linear combination
Definition
Inputs: an integer , vectors with ℓ2 sample/query access, real coefficients , and a positive integer cap . Set and . Output: an index or FAIL.
If , return FAIL. Otherwise, for each of at most trials:
- Using fresh randomness, choose with probability , then sample from . A zero-weight term is never chosen, so every requested vector distribution is defined.
- Query all and set By Cauchy--Schwarz, , so this is a valid acceptance probability.
- With a fresh independent acceptance coin, accept index with probability and return upon acceptance.
If all trials reject, return FAIL. When but , every defined acceptance probability is zero and the cap still forces a finite FAIL output. No normalized target distribution is assigned when . All samples, entry queries, random bits, cap, and output writing are charged. Correctness and cost of the linear-combination sampler ↗ proves the target law on the accepted branch and the failure and query bounds under its separate hypotheses on , the cancellation ratio, and .
Two-stage ℓ2 sampling controls Gram matrices
Statement
Let be nonzero, let be a probability distribution satisfying for some , and let be a positive integer. Draw independent row indices from this distribution and form with row equal to ; a zero row with is never sampled. Then and for every , For column sampling, assume the analogous probability lower bound on columns and apply the result to .
In particular, let be nonzero with , and use samples at each of two stages. First draw independent with probabilities and set . Conditional on , draw independent with probabilities and set . These are the exact normalized two-stage sketches; zero-probability indices are never selected. Then , both Gram estimators are unbiased (the second conditional on ), and For every , the probability that either Gram error is at least is at most ; the two errors need not be independent.
Facts & Assumptions
Given: The finite matrices, sampling probabilities, positive sample count, and independent samples within each stage in the Statement; second-stage sampling is conditional on the realized first stage.
Finite probability spaces, outcome weights, events, and event probabilities: a finite probability law assigns nonnegative weights summing to one; here the assumed row law assigns weight to index .
Pairwise and mutual independence of finite-valued random variables: independent sampled indices make their entrywise matrix-valued contributions independent.
The Frobenius norm on real or complex matrices: is the finite sum of squared entry moduli.
Markov's inequality on a finite probability space: for nonnegative and , .
Proof
The assumed lower bound on forces whenever . For one sampled row , define the random matrix when , and on zero-probability zero rows. For each matrix entry , Thus is unbiased.
The squared Frobenius norm of one contribution is . Therefore using the lower bound on each .
Expand entrywise. For distinct , independence [F2] and zero means make each cross expectation zero. Hence its expectation is by step 1.2.
Apply [F4] to the nonnegative random variable with threshold . Step 2.1 yields the probability bound. Transposing interchanges rows and columns, proving the column assertion.
In the exact two-stage construction every sampled row of has squared norm , so and . Every sampled column of likewise has squared norm , giving . Apply steps 1.1–2.1 with , , and row law for the first Gram estimator. For each fixed , apply the transposed result with , , and row law for the second estimator. This proves both unbiasedness assertions and both second-moment bounds, with the second conditional on .
Step 3.1 bounds the first-stage tail by and bounds the second-stage tail by the same number conditional on every . Averaging the latter bound over the finite first-stage histories gives the same unconditional bound. Pointwise the indicator of either failure is at most the sum of the two failure indicators, so summing their finite probability weights gives . No independence between the stages' errors is used.
Squared-vector distributions are close in total variation
Statement
Let be nonzero, and define and . Then The bound is valid even when its right side exceeds .
Facts & Assumptions
Given: Nonzero real vectors .
ℓ2 sample/query access to vectors and matrices: for nonzero , and likewise for .
End-to-end channel and output error contract: for distributions on a common finite alphabet.
Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation: Cauchy–Schwarz and the triangle inequality hold for the Euclidean norm.
Proof
Put , , , and . Then by their definitions.
Using and [F3],
Since , [F3] gives where is the reverse triangle inequality. Combining this with step 2.1 proves the result.
Estimate an inner product from ℓ2-sampling access
Statement
Let be nonzero with exact norm, ℓ2 sample access, and entry-query access, and let have entry-query access and known norm. For , a positive integer , and a positive integer with , there is a finite estimator using the lower median of block means, with samples from and at most entry queries to each of and , such that If , return without sampling. All sample and entry-query costs are charged.
Facts & Assumptions
Given: The access interfaces, and vectors in the Statement.
ℓ2 sample/query access to vectors and matrices: sampling index from nonzero gives probability .
Chebyshev's inequality on a finite probability space: for a finite real variable and , .
Pairwise and mutual independence of finite-valued random variables: independent samples give independent copies of any deterministic function of the sampled index.
The finite union bound: the probability of a finite union is at most the sum of the event probabilities.
Proof
If , the prescribed output equals and has zero error. Suppose henceforth . For one sample , query and and define when and when . Since a zero coordinate has sampling probability zero, [F1] gives .
Directly from the sampling law [F1], . Hence .
Partition independent samples into disjoint blocks of size . Let be the mean of in block and define to be the -th order statistic of (the lower median). Independence gives each block mean variance at most . Apply [F2] with the positive threshold ; each block mean is outside the target interval with probability at most .
If the median lies outside the target interval, at least block means lie outside it. For each fixed subset of that size the probability all its blocks fail is at most by [F3]; there are at most such subsets. By [F4], the failure probability is at most .
Adaptive conditional failure budgets compose by a union bound
Statement
Consider a procedure with finitely many stages and a finite probability space of execution histories. Let for every stage, and let be the event that stage is the first failing stage. If, after every successful history of positive probability that reaches stage , the conditional probability of failure at stage is at most , then No independence hypothesis is required. If fresh trials are mutually independent and each succeeds with probability at least , their miss probability is at most . For a common success probability , the expected number of trials capped at is , which tends to as .
Facts & Assumptions
Given: The finite staged procedure and its conditional failure bounds.
Finite probability spaces, outcome weights, events, and event probabilities: event probabilities are sums of weights in a finite probability space.
Normalization, nonnegativity, monotonicity, complements, and differences in a finite probability space: probabilities are nonnegative, and probabilities of disjoint events add.
Pairwise and mutual independence of finite-valued random variables: events determined by distinct members of a mutually independent family factor.
for every real , hence : for and natural , .
Proof
The stage-failure events are part of the given execution model. The events are pairwise disjoint: one execution cannot have two different first failing stages. Their union is exactly the event that some stage fails, because every failing finite history has a least index at which failure first occurs.
For mutually independent trials with failure events , [F3] gives . Apply [F4] to obtain the stated exponential bound, including and .
For fixed , partition by the finite successful histories that reach stage . For each such history, the hypothesis gives ; histories of probability zero contribute zero. Summing the conditional probabilities weighted by gives . By [F2] and step 1.1, .
Let be the number of trials used by the capped repeat-until-success procedure, with if all trials fail. For , is the event that the first trials fail, so . For a finite nonnegative integer-valued variable, by summing its finite outcome weights; hence for the finite geometric sum is . By step 1.2, as , so this capped expectation tends to .
Allocated local errors meet the global budget
Statement
Suppose stage compilers called by Allocate local error budgets through a pipeline return implemented channels with local normalized diamond errors at most and, for every successful history reaching stage , conditional failure probability at most . Then the final channel error is at most and the probability of any stage failure is at most . Use the finite execution-history model of End-to-end channel and output error contract. No independence between stage failures is assumed.
For resource composition, let be each stage's actual cost record and let record, exactly once, shared input/preprocessing, allocator arithmetic, and every interstage wiring/control, postprocessing, and final-output operation not already included in a . For each additive count, length, or sequential-time coordinate of Classical output, quantum-state output, and access costs and Circuit-family and computation cost model—including input/circuit-description length, preprocessing and sequential compiler time, each named data-structure or oracle operation, state-preparation and its controlled/inverse calls, gate occurrences by type, measurements, classical predicate/postprocessing, and output-bit writes—assume the pipeline performs precisely the operations in these records and has no unrecorded work. Then its total is If any records are upper bounds rather than exact counts, replace the equality by the corresponding upper bound. For serial circuit layers, write the actual gate word as , with the complete stage- word and its chronological auxiliary blocks, allowing empty blocks. Put and , retaining each auxiliary block at its position rather than minimizing the depth after deleting stages. Then . For each wire/workspace category (total wires, clean ancillas, dirty ancillas, or other workspace), let be the stage allocation and the shared/interface allocation. Identify these allocations with physical wire sets covering all wires of the total circuit, with each category inherited from the total circuit's declared classification. For clean or dirty ancilla counts, require the declared data/ancilla partitions in the stage and auxiliary records to be restrictions of the corresponding total partition under these wire identifications, and require all corresponding exact return contracts; apply this condition separately to each ancilla coordinate. A coordinate without these compatible supplied contracts is not covered by the bound. If stage workspaces are allocated disjointly, then . The local precision entries and error/failure probabilities remain separate fields; they are not resource counts to add.
Facts & Assumptions
Given: The compilers and local guarantees in the Statement, with rational .
Allocate local error budgets through a pipeline: for it assigns and to every stage and calls each compiler once.
Hybrid bound for a composed channel circuit: the final normalized channel error is at most the sum of local normalized channel errors.
Adaptive conditional failure budgets compose by a union bound: conditional failure bounds after every successful history imply total failure at most , without independence.
Classical output, quantum-state output, and access costs and Circuit-family and computation cost model: cost records separate input/preprocessing, named data access, state preparation, gate types, width/workspace, measurement, classical processing, output writing, and finite-gate precision; unprovided access or connecting operations are not free.
Proof
By [F1], each stage error is at most , so the sum of the local errors is at most .
By [F1], each stage has conditional failure bound . Apply [F3] to the finite execution history to obtain total failure at most .
Apply [F2] to those stage channels and use step 1.1; the final normalized diamond error is at most .
For each count-valued additive coordinate , the stated pipeline accounting partitions its operations into the auxiliary operations in and the stage operations in , each counted once. Counting this disjoint union gives ; replacing exact stage counts by upper bounds gives the corresponding inequality. Choose a minimum layering separately for each finite word , using no layers for an empty word. Concatenate these layerings in that order. Within each block every shared-wire precedence is respected, and between blocks every earlier gate is scheduled first. This gives an admissible layering of the actual full word with layers, proving the minimum-depth upper bound. For a wire/workspace category , let be the total set counted and let be the sets counted in the auxiliary and stage records. Coverage and the compatible category restrictions give . Thus finite cardinality gives , as asserted. In particular a total ancilla cannot be recorded only as stage data; the clean and dirty bounds each use their own compatible partitions and supplied return contracts. Precision, error, and failure fields remain separate by [F4].
Correctness and cost of the linear-combination sampler
Statement
Run Finite rejection sampler for a linear combination with , , and positive integer cap . Put , and suppose a proof-side bound and satisfy . The routine itself receives but does not test or require a certificate for ; these are hypotheses for this theorem's guarantee. A proposal is accepted with probability , and conditional on acceptance its index has law . The probability of returning FAIL is at most . The worst-case entry-query count is at most , with at most mixture draws and vector-sample calls in addition. If an entry query, a mixture draw, and a vector-sample call have uniform cost upper bounds , , and in the declared model, their total access cost is at most , where charges the initial norm queries and mixture setup. Arithmetic, acceptance randomness, and output writing are charged separately; no unit-cost sampling or bit-time bound is assumed. As , the capped expected number of trials tends to .
Facts & Assumptions
Given: The access and coefficient hypotheses of the Statement, with .
Finite rejection sampler for a linear combination: with , the proposal chooses term proportionally to and then samples from ; its acceptance probability is the displayed ratio. The routine's finite cap does not depend on detecting the truth of .
ℓ2 sample/query access to vectors and matrices: each nonzero term sampler has law with fresh conditional randomness; entry queries, norm queries, and sample calls have separately declared costs.
for every real , hence : for , .
Adaptive conditional failure budgets compose by a union bound: for independent trials of success probability , the expected capped number of trials tends to as the cap tends to infinity.
Proof
Put . The proposal probability of coordinate is , because choosing has weight and the conditional coordinate law is . Its target probability is .
When , [F3] shows the acceptance ratio lies in . If the denominator is zero, then every , so and the prescribed acceptance is .
Multiplying proposal and acceptance probabilities gives . Summing over shows total acceptance probability is , and conditioning on acceptance yields . Also by [F3] applied to each coordinate and summed.
Independent trials miss with probability . Since , this is at most by [F4] and the cap hypothesis. Each executed trial queries the entries at its sampled coordinate, makes one mixture draw, and calls one vector sampler. At most trials therefore use at most entry queries and calls of each other type. Initial norms and mixture setup contribute ; multiplying the call counts by the declared uniform upper bounds gives . This count excludes arithmetic, acceptance randomness, and output writing, which retain their separate charges. With per-trial success probability , [F5] shows the capped expected trial count tends to as .
Exponential Frobenius concentration for normalized Gram sampling
Statement
Let , , and let be its exact normalized two-stage sketches in Two-stage ℓ2 sampling controls Gram matrices, using samples at each stage. For and , if then, with probability at least , both Deterministically . The confidence bound has no matrix-dimension factor.
Facts & Assumptions
Given: The finite sampling laws, exact sketch normalization, and parameters in the Statement.
Two-stage ℓ2 sampling controls Gram matrices proves the normalized Gram estimator is unbiased and has squared Frobenius mean error at most , at both sampling stages, with the second statement conditional on .
Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation gives on a finite probability space. Markov's inequality on a finite probability space bounds a nonnegative variable's upper tail by its mean divided by the threshold.
Proof
We first prove a finite bounded-differences estimate. If a centered finite real variable lies in an interval of length , let . Differentiating finite sums shows is the variance under the weights proportional to . Any variable in has variance at most . Thus ; since , integration gives for all real .
Suppose a real function of independent finite samples changes by at most if sample is replaced. Reveal samples in order and let be the finite conditional mean of . For each fixed earlier history, the possible values of lie in an interval of length at most , since averaging the future preserves the coordinate oscillation bound. Therefore is conditionally centered in an interval of that length. Applying step 1.1 conditionally and summing over successive histories gives . Markov applied to this exponential, followed by , gives when and the denominator is positive. If all , is constant on the product support and its upper tail is zero.
For a finite random vector with identically, set and . Independence cancels the mixed centered terms, and , so . By [F2], . Replacing one changes by at most , by the triangle inequality. Step 2.1 now gives .
The row sampling is step 3.1 with and . Every row of has squared norm , so . Its column mixture probability is ; hence every normalized sampled column of has squared norm , and . Conditional on every realized , step 3.1 applies to the column sampling with the same and failure bound . Averaging that bound over the finite row histories and adding the first-stage failure bound proves the assertion. Zero-probability rows and columns are never selected and require no division by zero.
A deterministic threshold-band bound for two Gram sketches
Statement
Let be nonzero, , , , and . Suppose and satisfy and Let be the orthonormal left singular vectors of with singular values strictly greater than , and let contain those singular values. Write for their number, , and , using empty matrices when . Then , the eigenvalues of lie in , and some threshold-band target of Threshold-band singular-value approximation obeys
Facts & Assumptions
Given: The matrices and bounds in the Statement. Matrices are real and all projections below are orthogonal.
The Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition identifies the nonzero eigenvalues of and . The Real spectral theorem: a self-adjoint endomorphism of a finite-dimensional real inner product space has an orthonormal eigenbasis diagonalizes every real symmetric matrix used below.
Cauchy-Schwarz with its equality case, the triangle inequality for , the parallelogram law and polarisation applies to matrix entries, giving and .
A Threshold-band singular-value approximation is , where projects onto singular values at least and is any projection supported in the band .
Proof
Set , , , , and . Since and each retained eigenvalue exceeds , , so . For , . Therefore is the projection onto the range of . Also has Frobenius norm at most . By the SVD, the nonzero eigenvalues of are those of , giving . For , use and the same bounds with empty matrices.
We record the spectral comparisons explicitly. For a symmetric with decreasing eigenvalues, : the span of its first eigenvectors attains the value, and every -dimensional intersects the span of the eigenvectors numbered onward. Thus . Moreover over projections : in an eigenbasis the diagonal entries of are in and sum to , so moving their weight to the first indices maximizes the sum. By [F2] the top- sums for differ by at most . Pad matrices by zero blocks when dimensions differ. Applying these comparisons to and , and using [F1] between , gives and a top- sum difference at most .
If , is positive semidefinite. Hence . Also . Combining these inequalities with the top- comparisons in step 1.2 gives . The variational formula gives . If , .
Let count eigenvalues of strictly greater than . Between indices and , every crossing eigenvalue lies within of , by step 1.2 and the strict threshold rule for . There are at most crossing indices: if , then since each of its first eigenvalues exceeds and ; if , then . Thus the difference between the sums of the largest values of and all positive values of is at most . Put in an eigenbasis of . Since , the nonnegative threshold deficit is For the same crossing argument applies, and zero eigenvalues can be padded without changing the sums.
Let project onto eigenvalues respectively in , , and . Write for the part of on the first and third sets, and . Both outside gaps from are at least , so and . This uses no strict separation of eigenvalues within the band.
On the band space, is a positive contraction. Let be its spectral projection onto eigenvalues at least . In an eigenbasis of , . For the last inequality, the high cross block satisfies and the low block satisfies . These squares sum to at most .
Evaluate by summing its rows in an eigenbasis of . Its high contribution is , because there. Its low contribution is at most , and its band contribution is at most . Thus the sum is at most . By [F3], is a permitted target. Finally ; the triangle inequality gives the result. The same block formulas cover empty high, band or low spaces, , repeated eigenvalues, , and both endpoint conventions.
ModFKV low-rank approximation supplier
Statement
Let be nonzero, , and suppose it has the exact row/column ℓ2 sample/query access of ℓ2 sample/query access to vectors and matrices. Let , , , and . Set and . For a sufficiently large universal constant in the sample count the procedure ModFKV low-rank description procedure outputs a description of a rank-at-most- matrix such that, with probability at least , there is a threshold-band target from Threshold-band singular-value approximation satisfying The ModFKV description itself queries matrix entries and computes an SVD of a matrix at arithmetic cost, before any later sampling from 's rows. These operation counts use the stated exact-access arithmetic model and do not include unprovided preprocessing for the input access structure.
Facts & Assumptions
Given: The matrix, exact access, and parameters of the Statement.
Exponential Frobenius concentration for normalized Gram sampling controls both Gram errors with probability at least , and proves exactly.
A deterministic threshold-band bound for two Gram sketches gives the retained rank, near-isometry, and the threshold-band error from these two Gram errors.
Proof
Set . Choose the integer . This has the asserted order, uniformly for , and [F1] gives both Gram errors at most except on an event of probability at most . The algorithm's mixture probability equals , so a zero-probability column is never drawn and its INVALID branch has probability zero.
Since , and , . On the event in step 1.1, [F2] applies to the actual strictly retained singular values and to . It gives rank at most and a threshold-band target with error at most . Empty retention returns zero and is included in [F2]. No gap between actual singular values is assumed.
The row samples, mixture samples and the probabilities of the sampled columns can be obtained using at most input queries: each mixture probability sums known row probabilities and queried squared entries. Forming uses entries; its exact-arithmetic SVD and the stored matrices require arithmetic. These costs give the stated counts for the description and make no claim to materialize , or in full.
Tang's conditional recommendation sampler
Statement
Let have the row/column entry-query and ℓ2-sampling access of ℓ2 sample/query access to vectors and matrices, let , and fix a row with . Here is the ModFKV matrix description from ModFKV low-rank description procedure and its success event is specified by ModFKV low-rank approximation supplier. Assume Assume countable choice and dependent choice, and model the finite sketch history together with the countably many potential fresh row-sample calls in the adaptive refinement by the joint product probability law supplied by Assuming countable and dependent choice, countable products of arbitrary probability spaces. After a finite row-refinement stop, the capped rejection sampler uses its finite conditional law; the joint execution law is the countable mixture over finite stopping histories together with any residual nontermination event, on which no rejection call occurs. There is a Tang recommendation procedure whose sketch succeeds unconditionally with probability at least . For every fixed successful description having , its conditional probability of a successful row estimate and an accepted sample is at least . Thus the unconditional probability of a nonzero realized target row together with failure of these guarantees is at most . If the nonzero-row promise holds almost surely, this implies overall success probability at least . For every fixed successful sketch and successful row estimate, the law of the accepted classical sample is within total variation (where as in End-to-end channel and output error contract) of the normalized squared-entry distribution of . A capped sampling failure returns FAIL; no conditional distribution is assigned to a zero target row. The nonzero-row hypothesis is a promise on the realized description, not an assumed ability to compute its full row norm. The failure bound is unconditional; it is not a claim that conditioning on a nonzero realized row preserves ModFKV's confidence. At constant precision its query and arithmetic-time bounds are respectively For the final total-variation guarantee, the finite row refinement below adds a query and arithmetic-time factor . These bounds hold on the successful event; refinement terminates almost surely for every realized description satisfying the nonzero-row promise. Hidden logarithms, access-structure preprocessing, arithmetic precision, and the cost of writing the returned sample remain charged. The result is conditional on this exact input model and does not output as a full matrix or a full classical amplitude list.
Facts & Assumptions
Given: The exact input access and parameters, and a realized nonzero target row as promised in the Statement.
ModFKV low-rank approximation supplier and its proof give a successful sketch with failure probability at most its assigned budget; A deterministic threshold-band bound for two Gram sketches gives and eigenvalues of in for , where and . The retained singular values in lie in .
With exact norm, fresh ℓ2 samples and entry queries to , and entry queries and a known norm for , Estimate an inner product from ℓ2-sampling access estimates to error with failure at most using samples and queries.
Apply Squared-vector distributions are close in total variation to the ordered pair and use symmetry of total variation: for nonzero vectors, .
Correctness and cost of the linear-combination sampler proves exact sampling on acceptance and a miss probability at most the assigned budget under a proof-side cancellation bound and a finite cap.
Under countable choice and dependent choice, the finite sketch history and the countable stream of repeated fresh samples from the fixed row distribution admit the joint product probability law stated in the Statement (Assuming countable and dependent choice, countable products of arbitrary probability spaces). Every finite stopping history has a finite conditional sampler law by [F4], so the stopped branches form a countable mixture; paths that never stop remain the residual nontermination event in the product row-sample space.
A countable union of events has probability at most the sum of their probabilities (Finite and countable subadditivity of measures).
If a sequence of events has summable probabilities, only finitely many occur almost surely (First Borel-Cantelli lemma for events).
By the singular-value decomposition and the operator-norm characterization, for the real matrix (The operator norm is zero on the zero domain and otherwise is max_{||v||=1} ||Tv||, Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition, The operator norm is 0 on the zero domain and otherwise equals the largest singular value, attained at a right-singular vector).
Proof
Put . If , return FAIL: then and the promised input is absent. Otherwise is available from exact access. Run ModFKV with failure budget , retaining its rank , sample count, and of [F1]. If the retained rank is , return FAIL because again . The sketch event of [F1] has probability at least ; all subsequent estimates use fresh randomness conditional on the actual sketch. Write and , so .
At stages , set , , and . Put and choose the least positive integer with . For each of the entries of , use [F2]'s lower-median estimator with samples per block, blocks, tolerance , and failure budget ; use mutually independent fresh sample batches for these estimates. Each row of has norm . Thus on that stage's simultaneous success, , and the explicit vector obeys . Stop at the first stage with . All quantities in this rule are computable from stored coefficients and exact input row norms; no query of the unknown full output-row norm occurs.
Fix a realized sketch and use [F5] to predefine the disjoint finite sample batches for every potential stage . Let be the cylinder event that at least one of that stage's inner-product estimates fails its tolerance. Fresh samples, [F2], and [F6] give . Since , [F6] gives , and [F7] gives . Thus with probability at least every stage estimate succeeds, and almost surely all sufficiently late stage estimates succeed. On those late stages , so . The nonzero-row promise implies ; hence the stopping threshold tends to zero while , and the procedure terminates almost surely, including for a sketch outside the successful sketch event whenever its realized target row is nonzero. On the event that every stage estimate succeeds, it stops no later than a stage with , because . Hence the stopped stage satisfies .
On the successful sketch and estimation event, at the stopping stage . With and , [F1] gives . In particular , and [F3] yields . This holds for every fixed successful sketch and estimate, so any mixture of their accepted output laws has the same bound when compared with the corresponding mixture of target laws.
Form the coefficient row , so . Its numerator in the cancellation ratio is . On the successful sketch, , so the ratio is at most . Take as a proof-side bound and integer cap , then invoke the finite sampler [F4] with terms and cap . The bound is established only on the successful sketch event; the procedure does not test or require a certificate for it. On every other history the sampler is still executable: if its component-square sum is zero it returns FAIL immediately, and otherwise its finite proposal/acceptance loop returns an accepted index or FAIL by the cap. Conditional on any fixed good coefficients, [F4] gives miss probability at most and accepted law exactly .
Sketch failure has probability at most . Conditional on every fixed sketch, any row-estimation failure has probability at most by step 3.1. Conditional on a successful sketch and row estimate with nonzero target row, capped rejection failure has probability at most by [F4]. By [F5], the finite conditional rejection laws are attached to their countably many finite stopping histories, so the failure event is measurable in the joint execution model; [F6] bounds the union of the sketch, estimation, and applicable rejection failures by . Conditioning on rejection success does not manufacture confidence: for each fixed successful estimate its accepted law is already by [F4]. Step 4.1 supplies the asserted total-variation conclusion on that branch. Outside it the routine may return FAIL or an uncontrolled sample; no assertion treats a failed event as successful.
Stage uses queries and arithmetic for inner products, plus coefficient arithmetic. Summing geometrically to the stopping stage gives on simultaneous success. Since , this is . Rejection uses trials with queries per trial; forming its weights costs . Even the looser bound suffices.
Substituting , the base query bound is , and the SVD time is . Since , the refined estimate cost in step 6.2 is at most these base bounds times . On the successful sketch, and [F1], [F8] give . Hence , so this last factor is . Replacing by gives the stated powers, with all access preprocessing, precision and writing costs charged separately.
Dequantization is relative to the input access and output task
Remarks
Under the nonzero-row promise and countable-product sampling model stated in the theorem, Tang's real-matrix result matches a particular input structure: entry queries together with efficient row/column ℓ2-norm sampling, including the preprocessing and storage that provide those samplers. Its task is to sample from a specified recommendation row derived from a thresholded low-rank approximation, and its output is a sample rather than a full matrix or a classical list of amplitudes. It therefore establishes a classical counterpart for that access/task/output contract. It does not assert efficient classical simulation of arbitrary quantum circuits, arbitrary quantum-state output, or algorithms whose input model lacks the required ℓ2 samplers. Compare algorithms only after matching input access, output type, approximation metric, and charged preprocessing.
5 · Examples, counterexamples and false statements
None yet.
Sources
- Ewin Tang, A quantum-inspired classical algorithm for recommendation systems, version 3
- Alan Frieze, Ravi Kannan, and Santosh Vempala, Fast Monte-Carlo Algorithms for Finding Low-Rank Approximations, author manuscript (22 October 1998), section 4.1, PDF page 6
- John Watrous, The Theory of Quantum Information
- Andrew M. Childs, Lecture Notes on Quantum Algorithms
- Ronald de Wolf, Quantum Computing: Lecture Notes
- Alan Frieze, Ravi Kannan, and Santosh Vempala, Fast Monte-Carlo Algorithms for Finding Low-Rank Approximations, full author manuscript
- Alan Frieze, Ravi Kannan, and Santosh Vempala, Fast Monte-Carlo Algorithms for Finding Low-Rank Approximations
- J. Matoušek and J. Vondrák, The Probabilistic Method