Java证明威尔逊定理

29 篇文章 8 订阅
3 篇文章 0 订阅

威尔逊定理是什么?

在初等数论中,威尔逊定理给出了判定一个自然数是否为素数的充分必要条件。即:当且仅当 p p p 为素数时: ( p − 1 ) ! ≡ − 1 ( m o d p ) (p - 1) ! ≡ -1 (mod p) (p1)!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) (p1)!62(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) (p1)!nab0(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(p1)!n(k2k)2nk20(modp)

必要性

p p p 是素数,取集合 A = 1 , 2 , 3 , . . . p − 1 A={1,2,3,...p -1} A=1,2,3,...p1 ; 则 A A A 构成模 p p p 乘法的简化剩余系,即任意 i ∈ A i∈A iA ,存在 j ∈ A j∈A jA ,使得 ( 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 ) x21(modp)
解得: x ≡ 1 ( m o d p ) x ≡ 1 ( mod p ) x1(modp) x ≡ p − 1 ( m o d p ) x ≡ p - 1 ( mod p ) xp1(modp)
其余两两配对;故而 ( p − 1 ) ! ≡ 1 ﹡ ( p − 1 ) ≡ − 1 ( m o d p ) ( p - 1 )! ≡ 1﹡( p -1 ) ≡ -1 ( mod p ) (p1)!1﹡(p1)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;
	}
}
  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

YoungGeeker

你的鼓励是我最大的动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值