米勒-拉宾算法

米勒-拉宾算法:快速判断一个数是不是素数

需要用到的定理:

最小费马定理:如果n是素数,则(a ^ (n - 1)) % n恒等于1。

快速模取幂

米勒-拉宾算法就是结合上面两种,通过不断判断fmod(a, n - 1, n)的值是否为1来判断。这是一个概率算法,如果为1,不一定为素数,不为1,则必定是合数。循环判断多次就会让概率变得极为的小。

算法模板:

#include <iostream>

using namespace std;

int fmod(int a, int b, int c)//快速模取幂
{
	if(b == 1) return a;
	int t = fmod(a, b / 2, c);
	t = (t * t) % c;
	if(b & 1) t = (t * a) % c;
	return t;
}

int check(int n)//米勒-拉宾算法
{
	srand(time(0));
	for(int i = 0; i < 100; ++ i)
	{
		if(fmod(rand() % (n - 1) + 1, n - 1, n) != 1)//a的值需要变化,所以用到随机函数;a的取值为[1,n-1],因为这样就会保证如果n为素数,则n与a必互质
			return 0;
	}
	return 1;
}

int main()
{
	int n;
	while(scanf("%d", &n) != EOF)
	{
		if(n < 2)
		{
			printf("非素数\n");
			continue;
		}
		if(check(n))
			printf("素数\n");
		else
			printf("非素数\n");
	}
	return 0;
}


  • 0
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
费马素性检验算法米勒-拉宾素性检验算法和Solovay-Stassen素性检验算法都是判断一个数是否为素数的经典算法,它们的区别和联系如下: 1. 费马素性检验算法: 费马素性检验算法是最早提出的素性检验算法之一,它基于费马小定理:如果p是质数,a是不被p整除的任意整数,则a^(p-1) ≡ 1 (mod p)。因此,如果对于给定的a和p,a^(p-1) mod p = 1,则p可能是质数,否则p一定是合数。 2. 米勒-拉宾素性检验算法米勒-拉宾素性检验算法是一种基于费马小定理的改进算法,它随机选取若干个数a,对每个a进行费马小定理的检验,如果某个a使得a^(p-1) mod p ≠ 1且a^((p-1)/2) mod p ≠ -1,则p一定是合数,否则p可能是质数。米勒-拉宾素性检验算法的可靠性与选取的a有关,但是经过多次检验,误判的概率可以降低到极小。 3. Solovay-Stassen素性检验算法: Solovay-Stassen素性检验算法是一种基于欧拉准则的改进算法,它首先利用欧拉准则判断给定的数n是否为合数,如果n是合数,则返回false;如果n是质数,则随机选取若干个数a,并计算它们的雅可比符号,如果所有的雅可比符号都满足a^((n-1)/2) ≡ Jacobi(a, n) (mod n),则n可能是质数,否则n一定是合数。 这三种算法都是基于数学定理来判断一个数是否为素数。费马素性检验算法最简单,但误判率较高,适合于小范围的素数判断;米勒-拉宾素性检验算法相对来说误判率较低,适合于大范围的素数判断;Solovay-Stassen素性检验算法也是误判率较低的算法,但比米勒-拉宾素性检验算法更加复杂。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值