This theorem was stated by Bhāskara I (c. 600 - c. 680), Ibn al-Haytham (c.' 1000 AD), and John Wilson. Edward Waring announced the theorem in 1770, although neither he nor his student Wilson could prove it. Lagrange gave the first proof in 1771. There is evidence that Leibniz was also aware of the result a century earlier, but he never published it.
The following table shows the values of n from 2 to 30, (n-1)!, and the remainder when (n-1)! is divided by n. (The remainder when m is divided by n is written m mod n). The background color is pink for prime values of n, pale green for composite values.
Both of the proofs (for prime moduli) make use of the fact that the residue classes modulo a prime number are a field. See the article prime field for more details. Lagrange's theorem (in any field a polynomial of degree n has at most n roots) is needed for both proofs.
Composite modulus 
If n is composite it is divisible by some prime number q, where 2 ≤ q ≤ n − 2. If (n − 1)! were congruent to −1 (mod n) then it would also be congruent to −1 (mod q). But (n − 1)! ≡ 0 (mod q).
In fact, more is true. With the sole exception of 4, where 3! = 6 ≡ 2 (mod 4), if n is composite then (n − 1)! is congruent to 0 (mod n). The proof is divided into two cases: First, if n can be factored as the product of two unequal numbers, n = ab, where 2 ≤ a < b ≤ n − 2, then both a and b will appear in the product 1 × 2 × ... × (n − 1) = (n − 1)! and (n − 1)! will be divisible by n. If n has no such factorization, then it must be the square of some prime q, q > 2. But then 2q < q2 = n, both q and 2q will be factors of (n − 1)!, and again n divides (n − 1)!.
Prime modulus 
The result is trivial when p = 2, so assume p is an odd prime, p ≥ 3. Since the residue classes (mod p) are a field, every non-zero a has a unique multiplicative inverse, a−1. Lagrange's theorem implies that the only values of a for which a ≡ a−1 (mod p) are a ≡ ±1 (mod p) (because the congruence a2 ≡ 1 can have at most two roots (mod p).) Therefore, with the exception of ±1, the factors of (p − 1)! can be arranged in unequal pairs, where the product of each pair is ≡ 1 (mod p). This proves Wilson's theorem.
For example, if p = 11,
Prime modulus – another proof 
Again, the result is trivial for p = 2, so suppose p is an odd prime, p ≥ 3. Consider the polynomial
g has degree p − 1, leading term xp − 1, and constant term (p − 1)!. Its p − 1 roots are 1, 2, ..., p − 1.
h also has degree p − 1 and leading term xp − 1. Modulo p, Fermat's little theorem says it also has the same p − 1 roots, 1, 2, ..., p − 1.
f has degree at most p − 2 (since the leading terms cancel), and modulo p also has the p − 1 roots 1, 2, ..., p − 1. But Lagrange's theorem says it cannot have more than p − 2 roots. Therefore f must be identically zero (mod p), so its constant term (p − 1)! + 1 ≡ 0 (mod p). This is Wilson's theorem.
Wilson's theorem is useless as a primality test in practice, since computing (n − 1)! modulo n for large n is hard, and far easier primality tests are known (indeed, even trial division is considerably more efficient).
Using Wilson's Theorem, for any odd prime p = 2m + 1 we can rearrange the left hand side of
to obtain the equality
We can use this fact to prove part of a famous result: for any prime p such that p ≡ 1 (mod 4) the number (−1) is a square (quadratic residue) mod p. For suppose p = 4k + 1 for some integer k. Then we can take m = 2k above, and we conclude that
Wilson's theorem has been used to construct formulas for primes, but they are too slow to have practical value.
Gauss's generalization 
This further generalizes to the fact that in any finite abelian group, either the product of all elements is the identity, or there is precisely one element a of order 2 (but not both). In the latter case, the product of all elements equals a.
See also 
- O'Connor, John J.; Robertson, Edmund F., "Abu Ali al-Hasan ibn al-Haytham", MacTutor History of Mathematics archive, University of St Andrews.
- Edward Waring, Mediationes Algebraicae (Cambridge, England: 1770), page 218 (in Latin). In the third (1782) edition of Waring's Mediationes Algebraicae, Wilson's theorem appears as problem 5 on page 380. On that page, Waring states: "Hanc maxime elegantem primorum numerorum proprietatem invenit vir clarissimus, rerumque mathematicarum peritissimus Joannes Wilson Armiger." (A man most illustrious and most skilled in mathematics, Squire John Wilson, found this most elegant property of prime numbers.)
- Joseph Louis Lagrange, "Demonstration d'un théorème nouveau concernant les nombres premiers" (Proof of a new theorem concerning prime numbers), Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres (Berlin), vol. 2, pages 125–137 (1771).
- Giovanni Vacca (1899) "Sui manoscritti inediti di Leibniz" (On unpublished manuscripts of Leibniz), Bollettino di bibliografia e storia delle scienze matematiche ... (Bulletin of the bibliography and history of mathematics), vol. 2, pages 113-116; see page 114 (in Italian). Vacca quotes from Leibniz's mathematical manuscripts kept at the Royal Public Library in Hanover (Germany), vol. 3 B, bundle 11, page 10:
Original : Inoltre egli intravide anche il teorema di Wilson, come risulta dall'enunciato seguente:
"Productus continuorum usque ad numerum qui antepraecedit datum divisus per datum relinquit 1 (vel complementum ad unum?) si datus sit primitivus. Si datus sit derivativus relinquet numerum qui cum dato habeat communem mensuram unitate majorem."
Egli non giunse pero a dimostrarlo.
Translation : In addition, he [Leibniz] also glimpsed Wilson's theorem, as shown in the following statement:
"The product of all integers preceding the given integer, when divided by the given integer, leaves 1 (or the complement of 1?) if the given integer be prime. If the given integer be composite, it leaves a number which has a common factor with the given integer [which is] greater than one."
However, he didn't succeed in proving it.
- Landau, two proofs of thm. 78
- When n = 3, the only factors are ±1
- Gauss, DA, art. 78
- m = 1 and 2 have to be excluded because 1 ≡ −1 (mod 1 or 2).
The Disquisitiones Arithmeticae has been translated from Gauss's Ciceronian Latin into English and German. The German edition includes all of his papers on number theory: all the proofs of quadratic reciprocity, the determination of the sign of the Gauss sum, the investigations into biquadratic reciprocity, and unpublished notes.
- Gauss, Carl Friedrich; Clarke, Arthur A. (translator into English) (1986), Disquisitiones Arithemeticae (Second, corrected edition), New York: Springer, ISBN 0-387-96254-9
- Gauss, Carl Friedrich; Maser, H. (translator into German) (1965), Untersuchungen uber hohere Arithmetik (Disquisitiones Arithemeticae & other papers on number theory) (Second edition), New York: Chelsea, ISBN 0-8284-0191-8
- Landau, Edmund (1966), Elementary Number Theory, New York: Chelsea