快速的素数算法(LogN的复杂度哦)---------1

如果要写一个素数算法,我想大家第一反应会是如下的代码:

bool	IsPrime(unsigned 	uValue){
	for(unsigned	int	i=	2; i< uValue; ++i){
		if(uValue%i==	0){
			return	false;
		}
	}
	return true;
}

        这个代码的正确性肯定是不需要否定的,但是,我们根据数学中的唯一分解定理可以推出,如果一个数N不是素数,那么一定存在x使得x|N(注意'|'读作整除),那么知道x和N/x都是N的约数,自然的有min(x, N/x)<= 根号(N)。所以这个代码可以往下优化如下的代码:

bool	IsPrime(unsigned 	uValue){
	unsigned	int	uMax= sqrt(uValue);	//在math.h里面声明了
	for(unsigned	int	i=	2; i<= uMax; ++i){
		if(uValue%i==	0){
			return	false;
		}
	}
	return true;
}

       这个代码基本上是大家都写的代码,其正确性也是勿用质疑的。然后,这个程序仍然不是最简的,事实上
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
是的,Miller-Rabin算法是一种常用的素数判断算法之一,其准确性较高,时间复杂度为 O(klog^3n),其中 k 为测试次数,n 为待判断的数。 Miller-Rabin算法的基本思想是:对于一个合数 n,有至少一半的 a 值,满足 a^(n-1) mod n != 1。因此,我们可以随机选取一些 a 值,进行计算,如果不等于 1,则 n 不是素数。如果等于 1,则 n 可能是素数,需要进行多次测试以提高准确性。 具体实现可以采用以下步骤: 1. 将待判断的数 n-1 分解成形如 (2^s) * d 的形式,其中 s 和 d 都为正整数,d 为奇数。 2. 随机选取一个整数 a,满足 1 < a < n。 3. 计算 b0 = a^d mod n,如果 b0 = 1 或 b0 = n-1,则 n 可能是素数,可以进行下一次测试。 4. 依次计算 b1 = b0^2 mod n,b2 = b1^2 mod n,...,bs = bs-1^2 mod n,如果某一步中 bi = 1 且 bi-1 != n-1,则 n 不是素数,可以结束测试。 5. 如果 bs != 1,则 n 不是素数,可以结束测试。否则,n 可能是素数,可以进行下一次测试。 需要注意的是,测试次数 k 要足够多,一般来说,10次左右的测试已经可以满足绝大多数情况。如果需要更高的准确性,可以增加测试次数。 以下是 Python 代码实现: ```python import random def is_prime(n, k=10): if n == 2 or n == 3: return True if n < 2 or n % 2 == 0: return False # 将 n-1 分解成 (2^s) * d 的形式 s, d = 0, n-1 while d % 2 == 0: s += 1 d //= 2 # 进行 k 次测试 for i in range(k): a = random.randint(2, n-2) x = pow(a, d, n) if x == 1 or x == n-1: continue for j in range(s-1): x = pow(x, 2, n) if x == n-1: break else: return False return True ``` 需要注意的是,该算法的时间复杂度为 O(klog^3n),其中 k 和 logn 都是常数,因此在实际应用中,可以快速判断较大的素数

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值