威尔逊定理是什么?
在初等数论中,威尔逊定理给出了判定一个自然数是否为素数的充分必要条件。即:当且仅当 p p p 为素数时: ( p − 1 ) ! ≡ − 1 ( m o d p ) (p - 1) ! ≡ -1 (mod p) (p−1)!≡−1(modp) ,但是由于阶乘是呈爆炸增长的,其结论对于实际操作意义不大,但借助计算机的运算能力有广泛的应用,也可以辅助数学推导。
历史沿革
该定理是以英格兰数学家爱德华·华林的学生约翰·威尔逊命名的,尽管这对师生都未能给出证明。华林于1770年提出该定理,1773年由拉格朗日首次证明。
数学方法证明
充分性
如果 p p p 不是素数,当 p = 4 p = 4 p=4 时,显然 ( p − 1 ) ! ≡ 6 ≡ 2 ( m o d p ) (p-1)!≡6≡2(mod p) (p−1)!≡6≡2(modp), 当 p > 4 p>4 p>4 时,若 p p p 不是完全平方数,则存在两个不等的因数 a a a 、 b b b 使得 a b = p ab=p ab=p ,则 ( p − 1 ) ! ≡ n a b ≡ 0 ( m o d p ) (p-1)!≡nab≡0(mod p) (p−1)!≡nab≡0(modp);若 p p p 是完全平方数即 p = k 2 p=k^2 p=k2,因为 p > 4 p>4 p>4 ,所以 k > 2 k>2 k>2 , ( k , 2 k ) < p , ( p − 1 ) ! ≡ n ( k ∗ 2 k ) ≡ 2 n k 2 ≡ 0 ( m o d p ) (k,2k)<p,(p-1)!≡n(k*2k)≡2nk^2≡0(mod p) (k,2k)<p,(p−1)!≡n(k∗2k)≡2nk2≡0(modp)。
必要性
若
p
p
p 是素数,取集合
A
=
1
,
2
,
3
,
.
.
.
p
−
1
A={1,2,3,...p -1}
A=1,2,3,...p−1 ; 则
A
A
A 构成模
p
p
p 乘法的简化剩余系,即任意
i
∈
A
i∈A
i∈A ,存在
j
∈
A
j∈A
j∈A ,使得
(
i
j
)
≡
1
(
m
o
d
p
)
( i j ) ≡ 1 ( mod p )
(ij)≡1(modp) 。
那么A中的元素是不是恰好两两配对呢? 不一定,但只需考虑这种情况:
x
2
≡
1
(
m
o
d
p
)
x^2 ≡ 1 ( mod p )
x2≡1(modp)
解得:
x
≡
1
(
m
o
d
p
)
x ≡ 1 ( mod p )
x≡1(modp) 或
x
≡
p
−
1
(
m
o
d
p
)
x ≡ p - 1 ( mod p )
x≡p−1(modp)
其余两两配对;故而
(
p
−
1
)
!
≡
1
﹡
(
p
−
1
)
≡
−
1
(
m
o
d
p
)
( p - 1 )! ≡ 1﹡( p -1 ) ≡ -1 ( mod p )
(p−1)!≡1﹡(p−1)≡−1(modp) 。
Java方法证明
import java.math.BigInteger;
import java.util.Scanner;
public class WilsonTheorem {
public static void main(String[] args) {
System.out.println("请输入一个质数:");
Scanner in = new Scanner(System.in);
int p;
while (true) {
try {
p = in.nextInt();
if (!judge(p))
throw new Exception();
BigInteger factorial = new BigInteger("1");
for (int i = 2; i < p; i++) {
factorial = factorial.multiply(new BigInteger(i + ""));
}
factorial = factorial.add(new BigInteger("1"));
System.out.println("(p-1)! + 1 = " + factorial);
BigInteger result = factorial.divide(new BigInteger(p + ""));
System.out.println("(p-1)! + 1除以p = " + result);
} catch (Exception e) {
System.out.println("输入错误,请重新输入!");
continue;
}
}
}
/**
* 判断一个是否是质数的方法
*
* @param num
* @return
*/
public static boolean judge(int num) {
int i;
for (i = 2; i <= num / 2; i++) {
if (num % i == 0) break;
}
return i > num / 2 ? true : false;
}
}