For the first output qubit,
A Clifford circuit maps the Pauli observable under conjugation to a tensor-product Pauli operator , found by propagating it backward through the circuit in polynomial time. Since the input is a product state,
Each factor is efficiently computable, so both one-bit output probabilities are strongly simulatable. This is Heisenberg propagation of a Pauli observable through a Clifford circuit.
Propagate each output observable backwards through the Clifford circuit and write
The are mutually commuting Pauli operators, and measuring them on has exactly the required joint output distribution. Initially the first qubits are constrained by the stabilizer generators . Process the in order. If anticommutes with a current generator, part b gives a uniform outcome and a Clifford operation that replaces that generator by ; this step needs no measurement on . If commutes with every current generator, multiply it by known generators to remove its action on the first register. What remains is a Pauli observable on the resource qubits and is measured there. The nontrivial are independent and mutually commuting. An independent commuting family of Pauli operators on qubits has at most members, so . All effective observables are fixed by the original commuting family; the sampled outcomes merely update the classical Clifford frame. The resulting nonadaptive Pauli-based computation, followed by the stated outputs, is therefore a weak classical simulation of a quantum circuit with the same joint distribution.
The Gottesman--Knill theorem states that stabilizer-state preparation, Clifford circuits, and adaptive measurements of Pauli observables can be simulated in classical polynomial time. For , choose generators , , and . With each row written as , its sign-free stabilizer tableau is
Conjugating these generators successively by , , and gives , , and . Hence the output tableau, again ignoring signs, is
A strong classical simulation of a quantum circuit computes any requested output probability
in polynomial time, to the prescribed inverse-polynomial accuracy. A weak classical simulation of a quantum circuit instead produces classical samples from the circuit's output distribution, with exact or suitably small total-variation error.
The Extended Gottesman--Knill theorem states that a unitary Clifford circuit with an arbitrary product state input and final computational-basis measurements is weakly classically simulable. It is strongly simulable when only output qubits are measured. Indeed, each joint output projector expands into Pauli operators, and Clifford conjugation maps every such operator to another Pauli operator whose expectation factors over the input qubits. The factor is polynomial precisely for .