Adding the three equations for gives the necessary solvability condition . Conversely, subtraction of the first two and use of the third force
If , these values also satisfy the first two equations: for example,
The analogous calculation gives . Thus existence and uniqueness hold exactly under the stated compatibility condition:
For the three-variable system the matrix and data vector are
Perform Gaussian elimination on the augmented matrix. The operations and give
Each operation is an invertible row scaling followed by a row addition, so the solution set is preserved. Next gives
Back substitution now yields
and then
All three pivots are nonzero. Thus the solution exists uniquely for every , the matrix rank is , and the kernel of a linear map represented by is . The general rank-nullity theorem for a matrix with columns is
here .
Finally, reverse the viewpoint and regard as prescribed. The necessary and sufficient condition already proved applies with , , . Their sum is . Consequently can be recovered exactly when
When this holds, their unique values are and . The allowable triples form an affine plane, not all of .
The dual vector space is , the vector space of all linear functionals on . For a vector subspace , its annihilator of a vector subspace is
Given a basis , define by and extend linearly. If , evaluation at gives , so the are linearly independent. Every satisfies
so they also span and therefore form its dual basis.
Write . The condition is
Gaussian elimination gives and . Hence
Gaussian elimination without pivoting gives
Thus and the solution is unique exactly when .
For , forward substitution in a triangular system in gives . Consistency of requires . backward substitution in a triangular system then gives
Gaussian elimination without row exchanges gives the LU decomposition
Because has unit diagonal and is upper triangular,
Solving the two triangular systems gives
and
Gaussian elimination gives the LU decomposition
The prescribed diagonal entries of are visible, and direct multiplication gives .
Apply symmetric Gaussian elimination without row exchanges. At step , let be the leading diagonal entry of the remaining symmetric Schur complement. If , stop and report that is not positive definite. If , use it to eliminate the rest of its row and column. If all steps succeed, this constructs an LDL decomposition
where is unit lower triangular and has positive diagonal.
The test is correct from first principles. If all , then for every nonzero ,
because is invertible. Conversely, if is positive definite, its first pivot is , and completing the square gives
Choosing shows that the Schur complement is positive definite. Induction forces every pivot to be positive. This is also Sylvester's criterion.
At step , updating the remaining matrix costs arithmetic operations. Hence the total is
which proves the existence of the required algorithm.
Prepare , query the oracle, and measure or discard the output register. For , the input register becomes the coset state
Apply , the quantum Fourier transform over . The two amplitudes interfere destructively unless the binary inner product satisfies
and every vector in this orthogonal subspace is sampled uniformly. Repeat until independent equations have been collected, then use Gaussian elimination over to find their one-dimensional null space; its nonzero vector is . This is Simon's algorithm. It uses oracle queries with high probability and polynomial classical work. When , the function is injective and the samples eventually span all of , which distinguishes that case with arbitrarily high probability.
Choose finite generating sets and , using the Hilbert basis theorem. The assumed inclusion says that each vanishes on . By the Strong Hilbert Nullstellensatz, for every there is an exponent such that
The coefficients in an expression solve a finite system of linear equations with rational coefficients. Since it has a complex solution, Gaussian elimination gives a rational solution. Clearing the finitely many denominators produces a nonzero integer such that
for every .
For any prime , reduce these identities modulo . At a common zero of in the algebraic closure , they give , hence , for all . Therefore
for every prime except the finitely many divisors of . This is the spreading out of an affine zero-set inclusion.
Simon's algorithm 2026-09-28
Simon's algorithm prepares hidden-subgroup coset states, applies the quantum Fourier transform over , and samples vectors satisfying over . Repeated samples and Gaussian elimination recover with oracle queries and polynomial computation.