Write r:R→R, r(x)=−x, so that r∘r is the identity and −A=r−1(A) for every A⊆R.
Claim 1. Let A⊆R and let ((am,bm))m∈N be a sequence of open intervals covering A as in Lebesgue Outer Measure on the Real Line. Then the intervals (−bm,−am) cover −A: if x∈−A then −x∈A, so am<−x<bm for some m, i.e., −bm<x<−am; and each has the same length (−am)−(−bm)=bm−am. Thus every covering of A induces a covering of −A with the same total length, so λ∗(−A)≤λ∗(A); applying this inequality to −A in place of A and using −(−A)=A gives the reverse inequality, hence λ∗(−A)=λ∗(A).
Next, the family D={B⊆R: −B∈B(R)} is a σ-algebra: −R=R, −(R∖B)=R∖(−B), and −⋃mBm=⋃m(−Bm). It contains every open subset U of R: if x∈−U then −x∈U, so there is ε>0 with (−x−ε,−x+ε)⊆U, and then (x−ε,x+ε)=−(−x−ε,−x+ε)⊆−U, so −U is open, hence Borel. Since B(R) is the σ-algebra generated by the open sets, B(R)⊆D; that is, −B is Borel whenever B is. Since λ is the restriction of λ∗ to B(R) by Existence of Lebesgue Measure on the Real Line, λ(−B)=λ∗(−B)=λ∗(B)=λ(B).
Claim 2. Let g(x)=exp(−x2/2) as in Standard Normal Distribution, and note g(−x)=g(x) for every x, since (−x)2=x2. We first show that for every measurable f:R→R (with respect to B(R)) with f≥0,
∫Rf∘rdλ=∫Rfdλ.
First, f∘r is measurable: (f∘r)−1(B)=−f−1(B) is Borel by Claim 1. For a nonnegative simple function s=∑ici1Ai with Borel Ai and ci≥0, the composition s∘r=∑ici1−Ai is again a nonnegative simple function, and its integral is ∑iciλ(−Ai)=∑iciλ(Ai) by Claim 1; that is, simple functions and their reflections have equal integrals. Moreover s↦s∘r is a bijection from the set of simple functions s with 0≤s≤f∘r pointwise onto the set of simple functions s′ with 0≤s′≤f pointwise: it is its own inverse because r∘r is the identity, and 0≤s≤f∘r holds pointwise if and only if 0≤s∘r≤f does. Hence the two suprema defining the integrals of f∘r and of f range over equal sets of real numbers, and the displayed identity follows.
Now fix a Borel set B and apply the identity to f=1Bg, which is measurable (a product of measurable functions, by the preliminaries of Square-Integrable Random Variables and the Mean-Square Inner Product applied on the measurable space (R,B(R))) and nonnegative. Since 1B(−x)=1−B(x) and g(−x)=g(x), we get f∘r=1−Bg, so
∫R1−Bgdλ=∫R1Bgdλ,
and dividing by the normalizing constant c of Standard Normal Distribution gives N(−B)=N(B).
Claim 3. The function −Z=(−1)Z is a random variable by the closure preliminaries of Square-Integrable Random Variables and the Mean-Square Inner Product. For every Borel set B, the events {−Z∈B} and {Z∈−B} are equal, since −Z(ω)∈B holds exactly when Z(ω)∈−B. Because the distribution of Z is N and −B is Borel by Claim 1,
P(−Z∈B)=P(Z∈−B)=N(−B)=N(B)
by Claim 2. Hence the distribution of −Z is N, i.e., −Z is standard normal by Standard Normal Distribution. ■