Section 13.3 A Lemma About Square Roots Modulo n
ΒΆWe'll continue our formal investigation of what numbers are sums of two squares by taking a look at a seemingly unrelated lemma about square roots. In Section 14.1 we'll see that square roots of negative one (thinking of β1βZ, not Zn) are connected to sums of squares as well, so it is not completely implausible to connect roots and these sums. Before we do this, let's codify something we already have discussed since Question 7.1.1 at various times, e.g. in Fact 7.3.1 or Section 7.6.Definition 13.3.1.
We say that a number a has a square root modulo n if there is some number x with
Fact 13.3.2.
For an odd prime p, the only way there is a square root of β1 modulo p is if pβ‘1 (mod 4).
Proof.
We will use group theory to prove this.
Assume there is a square root \(f\text{,}\) so that
Then the order of \(f\) in \(U_p\) is four, since
We know that the order
but then Lagrange's (group theory) Theorem 8.3.12 says that four divides \(p-1\text{.}\)
Given that, the only possible kind of prime \(p\) solving this is the form \(4k+1\text{.}\)
Lemma 13.3.3.
For an odd prime pβ‘1 (mod 4), there actually does exist a square root of β1 modulo p. That is, there is an f such that
Proof.
Before we start the proof, recall Wilson's Theorem, which states that
Do you remember our proof? We paired up all the numbers from \(2\) to \(p-2\) in pairs of multiplicative inverses (mod \(p\)), thus:
Our strategy for this proof will be similar.
This time, pair up the numbers from \(1\) to \(p-1\) in a different way, in pairs of additive inverses (mod \(p\)):
This makes sense because \((p-1)/2\) is an integer halfway between \(1\) and \(p\text{,}\) as \(p\) is odd.
If we rethink things (mod \(p\)), we can rewrite this in a more suggestive way. Let \(\left(1\cdot 2\cdot 3\cdots \frac{p-1}{2}\right)\) be called \(f\text{.}\) This is also \(\left(\frac{p-1}{2}\right)!\text{,}\) of course. Then, keeping in mind that \(\frac{p+1}{2}=p-\frac{p-1}{2}\text{,}\)
Remember that our hypothesis is \(p\equiv 1\) (mod \(4\)). Then \(p=4k+1\) for integer \(k\text{,}\) so \(\frac{p-1}{2}=2k\) is even and by Wilson's Theorem
xxxxxxxxxx
def _(p=(13,[q for q in prime_range(200) if q%4==1])):
f=mod(factorial((p-1)/2),p)
pretty_print(html(r"The potential square roots of $-1$ are $\pm \left(\frac{%s-1}{2}\right)!=%s,%s\text{ (mod }%s)$"%(p,f,-f,p)))
pretty_print(html(r"And we can compute that ${0}^2\equiv{1}$ and ${2}^2\equiv {3}$ modulo ${4}$".format(f,f^2,-f,(-f)^2,p)))