Quadratic residue

Fermat, Euler, Lagrange, Legendre, and other number theorists of the 17th and 18th centuries established theorems and formed conjectures about quadratic residues, but the first systematic treatment is § IV of Gauss's Disquisitiones Arithmeticae (1801). Article 95 introduces the terminology "quadratic residue" and "quadratic nonresidue", and says that if the context makes it clear, the adjective "quadratic" may be dropped. For a given n, a list of the quadratic residues modulo n may be obtained by simply squaring all the numbers 0, 1, ..., . Since a ≡ b (mod n) implies a2≡b2 (mod n), any other quadratic residue is congruent (mod n) to some in the obtained list. But the obtained list is not composed of mutually incongruent quadratic residues (mod n) only. Since a2≡(n−a)2 (mod n), the list obtained by squaring all numbers in the list 1, 2, ..., (or in the list 0, 1, ..., n) is symmetric (mod n) around its midpoint, hence it is actually only needed to square all the numbers in the list 0, 1, ..., \lfloor n/2 \rfloor. The list so obtained may still contain mutually congruent numbers (mod n). Thus, the number of mutually noncongruent quadratic residues modulo n cannot exceed n/2 + 1 (n even) or (n + 1)/2 (n odd). The product of two residues is always a residue. Prime modulus Modulo an odd prime number p there are (p + 1)/2 residues (including 0) and (p − 1)/2 nonresidues, by Euler's criterion. In this case, it is customary to consider 0 as a special case and work within the multiplicative group of nonzero elements of the field (\mathbb{Z}/p\mathbb{Z}). In other words, every congruence class except zero modulo p has a multiplicative inverse. This is not true for composite moduli. Following this convention, the multiplicative inverse of a residue is a residue, and the inverse of a nonresidue is a nonresidue. Following this convention, modulo an odd prime number there is an equal number of residues and nonresidues. to the law of quadratic reciprocity is that if p ≡ 1 (mod 4) then −1 is a quadratic residue modulo p, and if p ≡ 3 (mod 4) then −1 is a nonresidue modulo p. This implies the following: If p ≡ 1 (mod 4) the negative of a residue modulo p is a residue and the negative of a nonresidue is a nonresidue. If p ≡ 3 (mod 4) the negative of a residue modulo p is a nonresidue and the negative of a nonresidue is a residue. Prime power modulus All odd squares are ≡ 1 (mod 8) and thus also ≡ 1 (mod 4). If a is an odd number and m = 8, 16, or some higher power of 2, then a is a residue modulo m if and only if a ≡ 1 (mod 8). For example, mod (32) the odd squares are : 12 ≡ 152 ≡ 1 : 32 ≡ 132 ≡ 9 : 52 ≡ 112 ≡ 25 : 72 ≡ 92 ≡ 49 ≡ 17 and the even ones are : 02 ≡ 82 ≡ 162 ≡ 0 : 22 ≡ 62≡ 102 ≡ 142≡ 4 : 42 ≡ 122 ≡ 16. So a nonzero number is a residue mod 8, 16, etc., if and only if it is of the form 4k(8n + 1). A number a relatively prime to an odd prime p is a residue modulo any power of p if and only if it is a residue modulo p. If the modulus is pn, : then ''p'k'a'' :: is a residue modulo pn if k ≥ n :: is a nonresidue modulo pn if k < n is odd :: is a residue modulo pn if k < n is even and a is a residue :: is a nonresidue modulo pn if k < n is even and a is a nonresidue. Notice that the rules are different for powers of two and powers of odd primes. Modulo an odd prime power n = pk, the products of residues and nonresidues relatively prime to p obey the same rules as they do mod p; p is a nonresidue, and in general all the residues and nonresidues obey the same rules, except that the products will be zero if the power of p in the product ≥ n. Modulo 8, the product of the nonresidues 3 and 5 is the nonresidue 7, and likewise for permutations of 3, 5 and 7. In fact, the multiplicative group of the non-residues and 1 form the Klein four-group. Composite modulus not a prime power The basic fact in this case is : if a is a residue modulo n, then a is a residue modulo pk for every prime power dividing n. : if a is a nonresidue modulo n, then a is a nonresidue modulo pk for at least one prime power dividing n. Modulo a composite number, the product of two residues is a residue. The product of a residue and a nonresidue may be a residue, a nonresidue, or zero. For example, from the table for modulus 6   1, 2, 3, 4, 5 (residues in bold). The product of the residue 3 and the nonresidue 5 is the residue 3, whereas the product of the residue 4 and the nonresidue 2 is the nonresidue 2. Also, the product of two nonresidues may be either a residue, a nonresidue, or zero. For example, from the table for modulus 15   1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 (residues in bold). The product of the nonresidues 2 and 8 is the residue 1, whereas the product of the nonresidues 2 and 7 is the nonresidue 14. This phenomenon can best be described using the vocabulary of abstract algebra. The congruence classes relatively prime to the modulus are a group under multiplication, called the group of units of the ring (\mathbb{Z}/n\mathbb{Z}), and the squares are a subgroup of it. Different nonresidues may belong to different cosets, and there is no simple rule that predicts which one their product will be in. Modulo a prime, there is only the subgroup of squares and a single coset. The fact that, e.g., modulo 15 the product of the nonresidues 3 and 5, or of the nonresidue 5 and the residue 9, or the two residues 9 and 10 are all zero comes from working in the full ring (\mathbb{Z}/n\mathbb{Z}), which has zero divisors for composite n. For this reason some authors add to the definition that a quadratic residue a must not only be a square but must also be relatively prime to the modulus n. (a is coprime to n if and only if a2 is coprime to n.) Although it makes things tidier, this article does not insist that residues must be coprime to the modulus. == Notations ==

Gauss used and to denote residuosity and non-residuosity, respectively; : for example, and , or and . Although this notation is compact and convenient for some purposes, a more useful notation is the Legendre symbol, also called the quadratic character, which is defined for all integers and positive odd prime numbers as : \left(\frac{a}{p}\right) = \begin{cases}\;\;\,0&\text{ if }p \text { divides } a\\+1&\text{ if } a \operatorname{R} p \text{ and }p \text { does not divide } a\\-1&\text{ if }a \operatorname{N} p \text{ and }p \text{ does not divide } a\end{cases} There are two reasons why numbers ≡ 0 (mod ) are treated specially. As we have seen, it makes many formulas and theorems easier to state. The other (related) reason is that the quadratic character is a homomorphism from the multiplicative group of nonzero congruence classes modulo to the complex numbers under multiplication. Setting (\tfrac{np}{p}) = 0 allows its domain to be extended to the multiplicative semigroup of all the integers. One advantage of this notation over Gauss's is that the Legendre symbol is a function that can be used in formulas. It can also easily be generalized to cubic, quartic and higher power residues. There is a generalization of the Legendre symbol for composite values of , the Jacobi symbol, but its properties are not as simple: if is composite and the Jacobi symbol (\tfrac{a}{m}) = -1, then , and if then (\tfrac{a}{m}) = 1, but if (\tfrac{a}{m}) = 1 we do not know whether or . For example: (\tfrac{2}{15}) = 1 and (\tfrac{4}{15}) = 1, but and . If is prime, the Jacobi and Legendre symbols agree. ==Distribution of quadratic residues==

Although quadratic residues appear to occur in a rather random pattern modulo n, and this has been exploited in such applications as acoustics and cryptography, their distribution also exhibits some striking regularities. Using Dirichlet's theorem on primes in arithmetic progressions, the law of quadratic reciprocity, and the Chinese remainder theorem (CRT) it is easy to see that for any M > 0 there are primes p such that the numbers 1, 2, ..., M are all residues modulo p. For example, if p ≡ 1 (mod 8), (mod 12), (mod 5) and (mod 28), then by the law of quadratic reciprocity 2, 3, 5, and 7 will all be residues modulo p, and thus all numbers 1–10 will be. The CRT says that this is the same as p ≡ 1 (mod 840), and Dirichlet's theorem says there are an infinite number of primes of this form. 2521 is the smallest, and indeed 12 ≡ 1, 10462 ≡ 2, 1232 ≡ 3, 22 ≡ 4, 6432 ≡ 5, 872 ≡ 6, 6682 ≡ 7, 4292 ≡ 8, 32 ≡ 9, and 5292 ≡ 10 (mod 2521). Dirichlet's formulas The first of these regularities stems from Peter Gustav Lejeune Dirichlet's work (in the 1830s) on the analytic formula for the class number of binary quadratic forms. Let q be a prime number, s a complex variable, and define a Dirichlet L-function as :L(s) = \sum_{n=1}^\infty\left(\frac{n}{q}\right)n^{-s}. Dirichlet showed that if q ≡ 3 (mod 4), then :L(1) = -\frac{\pi}{\sqrt q}\sum_{n=1}^{q-1} \frac{n}{q} \left(\frac{n}{q}\right) > 0. Therefore, in this case (prime q ≡ 3 (mod 4)), the sum of the quadratic residues minus the sum of the nonresidues in the range 1, 2, ..., q − 1 is a negative number. For example, modulo 11, :1, 2, 3, 4, 5, 6, 7, 8, 9, 10 (residues in bold) :1 + 4 + 9 + 5 + 3 = 22, 2 + 6 + 7 + 8 + 10 = 33, and the difference is −11. In fact the difference will always be an odd multiple of q if q > 3. In contrast, for prime q ≡ 1 (mod 4), the sum of the quadratic residues minus the sum of the nonresidues in the range 1, 2, ..., q − 1 is zero, implying that both sums equal \frac{q(q-1)}{4}. Dirichlet also proved that for prime q ≡ 3 (mod 4), :L(1) = \frac{\pi}{\left(2-\left(\frac{2}{q}\right)\right)\!\sqrt q}\sum_{n=1}^\frac{q-1}{2}\left(\frac{n}{q}\right) > 0. This implies that there are more quadratic residues than nonresidues among the numbers 1, 2, ..., (q − 1)/2. For example, modulo 11 there are four residues less than 6 (namely 1, 3, 4, and 5), but only one nonresidue (2). An intriguing fact about these two theorems is that all known proofs rely on analysis; no-one has ever published a simple or direct proof of either statement. Law of quadratic reciprocity If p and q are odd primes, then: ((p is a quadratic residue mod q) if and only if (q is a quadratic residue mod p)) if and only if (at least one of p and q is congruent to 1 mod 4). That is: : \left(\frac{p}{q}\right) \left(\frac{q}{p}\right) = (-1)^{\frac{p-1}{2} \cdot \frac{q-1}{2}} where \left(\frac{p}{q}\right) is the Legendre symbol. Thus, for numbers a and odd primes p that don't divide a: Pairs of residues and nonresidues Modulo a prime p, the number of pairs n, n + 1 where n R p and n + 1 R p, or n N p and n + 1 R p, etc., are almost equal. More precisely, let p be an odd prime. For i, j = 0, 1 define the sets :A_{ij}=\left\{k\in\{1,2,\dots,p-2\}: \left(\frac{k}{p}\right)=(-1)^i\land\left(\frac{k+1}{p}\right)=(-1)^j\right\}, and let :\alpha_{ij} = |A_{ij}|. That is, :α00 is the number of residues that are followed by a residue, :α01 is the number of residues that are followed by a nonresidue, :α10 is the number of nonresidues that are followed by a residue, and :α11 is the number of nonresidues that are followed by a nonresidue. Then if p ≡ 1 (mod 4) :\alpha_{00} = \frac{p-5}{4},\;\alpha_{01} =\alpha_{10} =\alpha_{11} = \frac{p-1}{4} and if p ≡ 3 (mod 4) :\alpha_{01} = \frac{p+1}{4},\;\alpha_{00} =\alpha_{10} =\alpha_{11} = \frac{p-3}{4}. For example: (residues in bold) Modulo 17 :1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 ::A00 = {1,8,15}, ::A01 = {2,4,9,13}, ::A10 = {3,7,12,14}, ::A11 = {5,6,10,11}. Modulo 19 :1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 ::A00 = {4,5,6,16}, ::A01 = {1,7,9,11,17}, ::A10 = {3,8,10,15}, ::A11 = {2,12,13,14}. Gauss (1828) introduced this sort of counting when he proved that if p ≡ 1 (mod 4) then x4 ≡ 2 (mod p) can be solved if and only if p = a2 + 64 b2. Pólya–Vinogradov inequality The values of (\tfrac{a}{p}) for consecutive values of a mimic a random variable like a coin flip. Specifically, Pólya and Vinogradov proved (independently) in 1918 that for any nonprincipal Dirichlet character χ(n) modulo q and any integers M and N, : \left|\sum_{n=M+1}^{M+N}\chi(n)\right| =O\left( \sqrt q \log q\right), in big O notation. Setting : \chi(n) = \left(\frac{n}{q}\right), this shows that the number of quadratic residues modulo q in any interval of length N is : \frac{1}{2}N + O(\sqrt q\log q). It is easy to prove that : \left| \sum_{n=M+1}^{M+N} \left( \frac{n}{q} \right) \right| In fact, : \left| \sum_{n=M+1}^{M+N} \left( \frac{n}{q} \right) \right| Montgomery and Vaughan improved this in 1977, showing that, if the generalized Riemann hypothesis is true then : \left|\sum_{n=M+1}^{M+N}\chi(n)\right|=O\left(\sqrt q \log \log q\right). This result cannot be substantially improved, for Schur had proved in 1918 that : \max_N \left|\sum_{n=1}^{N}\left(\frac{n}{q}\right)\right|>\frac{1}{2\pi}\sqrt q and Paley had proved in 1932 that : \max_N \left|\sum_{n=1}^{N}\left(\frac{d}{n}\right)\right|>\frac{1}{7}\sqrt d \log \log d for infinitely many d > 0. Least quadratic non-residue The least quadratic residue mod p is clearly 1. The question of the magnitude of the least quadratic non-residue n(p) is more subtle, but it is always prime, with 7 appearing for the first time at 71. The Pólya–Vinogradov inequality above gives . The best unconditional estimate is n(p) ≪ pθ for any , obtained by estimates of Burgess on character sums. Linnik showed that the number of p less than X such that n(p) > Xε is bounded by a constant depending on ε. The least quadratic non-residues mod p for odd primes p are: : 2, 2, 3, 2, 2, 3, 2, 5, 2, 3, 2, ... Quadratic excess Let p be an odd prime. The quadratic excess E(p) is the number of quadratic residues on the range (0,p/2) minus the number in the range (p/2,p) . For p congruent to 1 mod 4, the excess is zero, since −1 is a quadratic residue and the residues are symmetric under r ↔ p−r. For p congruent to 3 mod 4, the excess E is always positive. == Computational complexity ==

Two natural computational problems are: • Given a number a and a modulus n, determine whether a is a quadratic residue. • Calculate a square root of a (that is, a solution to x2 ≡ a (mod n)), assuming one exists. For prime moduli, both problems can be solved efficiently using the Tonelli-Shanks algorithm; the same is true for composite moduli whose prime factorization is known. In the case of a composite modulus with unknown prime factorization, the problem of identifying quadratic residues is known as the quadratic residuosity problem and is believed to be computationally difficult. Modulo a prime p, a quadratic residue a has 1 + (a|p) roots (i.e. zero if a N p, one if a ≡ 0 (mod p), or two if a R p and gcd(a,p) = 1.) In general if a composite modulus n is written as a product of powers of distinct primes, and there are n1 roots modulo the first one, n2 mod the second, ..., there will be n1n2... roots modulo n. The theoretical way solutions modulo the prime powers are combined to make solutions modulo n is called the Chinese remainder theorem; it can be implemented with an efficient algorithm. For example: : Solve x2 ≡ 6 (mod 15). :: x2 ≡ 6 (mod 3) has one solution, 0; x2 ≡ 6 (mod 5) has two, 1 and 4. :: and there are two solutions modulo 15, namely 6 and 9. : Solve x2 ≡ 4 (mod 15). :: x2 ≡ 4 (mod 3) has two solutions, 1 and 2; x2 ≡ 4 (mod 5) has two, 2 and 3. :: and there are four solutions modulo 15, namely 2, 7, 8, and 13. : Solve x2 ≡ 7 (mod 15). :: x2 ≡ 7 (mod 3) has two solutions, 1 and 2; x2 ≡ 7 (mod 5) has no solutions. :: and there are no solutions modulo 15. Prime or prime power modulus First off, if the modulus n is prime the Legendre symbol \left(\frac{a}{n}\right) can be quickly computed using a variation of Euclid's algorithm or the Euler's criterion. If it is −1 there is no solution. Secondly, assuming that \left(\frac{a}{n}\right)=1, if n ≡ 3 (mod 4), Lagrange found that the solutions are given by : x \equiv \pm\; a^{(n+1)/4} \pmod{n}, and Legendre found a similar solution if n ≡ 5 (mod 8): : x \equiv \begin{cases} \pm\; a^{(n+3)/8} \pmod{n}& \text{ if } a \text{ is a quartic residue modulo } n \\ \pm\; a^{(n+3)/8}2^{(n-1)/4} \pmod{n}& \text{ if } a \text{ is a quartic non-residue modulo } n \end{cases} For prime n ≡ 1 (mod 8), however, there is no known formula. Tonelli (in 1891) and Cipolla found efficient algorithms that work for all prime moduli. Both algorithms require finding a quadratic nonresidue modulo n, and there is no efficient deterministic algorithm known for doing that. But since half the numbers between 1 and n are nonresidues, picking numbers x at random and calculating the Legendre symbol \left(\frac{x}{n}\right) until a nonresidue is found will quickly produce one. A slight variant of this algorithm is the Tonelli–Shanks algorithm. If the modulus n is a prime power n = pe, a solution may be found modulo p and "lifted" to a solution modulo n using Hensel's lemma or an algorithm of Gauss. Determining whether a is a quadratic residue or nonresidue modulo n (denoted or ) can be done efficiently for prime n by computing the Legendre symbol. However, for composite n, this forms the quadratic residuosity problem, which is not known to be as hard as factorization, but is assumed to be quite hard. On the other hand, if we want to know if there is a solution for x less than some given limit c, this problem is NP-complete; however, this is a fixed-parameter tractable problem, where c is the parameter. In general, to determine if a is a quadratic residue modulo composite n, one can use the following theorem: Let , and . Then is solvable if and only if: • The Legendre symbol \left(\tfrac{a}{p}\right)=1 for all odd prime divisors p of n. • if n is divisible by 4 but not 8; or if n is divisible by 8. Note: This theorem essentially requires that the factorization of n is known. Also notice that if , then the congruence can be reduced to , but then this takes the problem away from quadratic residues (unless m is a square). == Number of quadratic residues ==

The list of the number of quadratic residues modulo n, for n = 1, 2, 3 ..., looks like: : 1, 2, 2, 2, 3, 4, 4, 3, 4, 6, 6, 4, 7, 8, 6, ... Formulas to count the number of squares modulo n are given by Stangl. Let s(n) give the number of squares modulo n. This is a multiplicative function, so it is completely characterized by its values at prime powers. Modulo powers of 2 • s(2^n) = \frac{2^{n-1}+4}{3} for even n • s(2^n) = \frac{2^{n-1}+5}{3} for odd n Modulo powers of odd primes • s(p) = \frac{p+1}{2} • s(p^2) = \frac{p^2-p+2}{2} • s(p^n) = \frac{p^{n+1}+p+2}{2(p+1)} for even n \geq 3 • s(p^n) = \frac{p^{n+1}+2p+1}{2(p+1)} for odd n \geq 3 == Applications of quadratic residues ==

Acoustics Sound diffusers have been based on number-theoretic concepts such as primitive roots and quadratic residues. Graph theory Paley graphs are dense undirected graphs, one for each prime p ≡ 1 (mod 4), that form an infinite family of conference graphs, which yield an infinite family of symmetric conference matrices. Paley digraphs are directed analogs of Paley graphs, one for each p ≡ 3 (mod 4), that yield antisymmetric conference matrices. The construction of these graphs uses quadratic residues. Cryptography The fact that finding a square root of a number modulo a large composite n is equivalent to factoring (which is widely believed to be a hard problem) has been used for constructing cryptographic schemes such as the Rabin signature and the oblivious transfer. The quadratic residuosity problem is the basis for the Goldwasser–Micali cryptosystem. The discrete logarithm is a similar problem that is also used in cryptography. Primality testing Euler's criterion is a formula for the Legendre symbol (a|p) where p is prime. If p is composite the formula may or may not compute (a|p) correctly. The Solovay–Strassen primality test for whether a given number n is prime or composite picks a random a and computes (a|n) using a modification of Euclid's algorithm, and also using Euler's criterion. If the results disagree, n is composite; if they agree, n may be composite or prime. For a composite n at least 1/2 the values of a in the range 2, 3, ..., n − 1 will return "n is composite"; for prime n none will. If, after using many different values of a, n has not been proved composite it is called a "probable prime". The Miller–Rabin primality test is based on the same principles. There is a deterministic version of it, but the proof that it works depends on the generalized Riemann hypothesis; the output from this test is "n is definitely composite" or "either n is prime or the GRH is false". If the second output ever occurs for a composite n, then the GRH would be false, which would have implications through many branches of mathematics. Integer factorization In § VI of the Disquisitiones Arithmeticae Gauss discusses two factoring algorithms that use quadratic residues and the law of quadratic reciprocity. Several modern factorization algorithms (including Dixon's algorithm, the continued fraction method, the quadratic sieve, and the number field sieve) generate small quadratic residues (modulo the number being factorized) in an attempt to find a congruence of squares which will yield a factorization. The number field sieve is the fastest general-purpose factorization algorithm known. == Table of quadratic residues ==

The following table lists the quadratic residues mod 1 to 75 (a means it is not coprime to n). (For the quadratic residues coprime to n, see , and for nonzero quadratic residues, see .) == See also ==