质数(素数)的几种求法

关注

原创最新推荐文章于 2025-05-05 21:45:56 发布 8.3k 阅读 AI 写同款· GEO 优化 ›

质数（素数）的几种判断方法如果使用素筛法时，参见代码段6的init_prime() 函数，当i = 2，3，5时，30都会被标记一次，但是线筛法的思想是，找到最小素因子2所对应的因子15，用15*2 标记出30，因为15是30的最大因子，是30的对后一个因式，是i值循环加加的过程中，最后一次得到30的机会，用来标记30，30的有且一次的标记机会。当前质数为2时，标记2的倍数，6被标记了一次，当执行到质数3时，6作为3的倍数，又被标记了一次，所以，6被标记了2次。30存在的质数因子为：2，3，5，所以，它将被标记3次。阅读详情

质数就是只能被1和它本身整除的数。求质数在程序设计和算法中很常见，尤其是在密码学中经常用到质数。最通常的想法就是依据定义来求质数

如求100以内的质数，我们会写出如下的算法：

int prime_1(int num)
{
	int cur=2;
	int count=0;   //记录质数个数
	while (cur<num)
	{
		int factor=2;  //每次循环从2开始判断
		while (factor<cur)
		{
			if (cur%factor==0)
			{
				break;
			}
			factor++;
		}
		if (factor==cur)
		{
			cout<<cur<<" ";
			count++;
		}
		cur++;
	}
	return count;
}

对于每一个质数都要除到比它小的那个数才能确定下来，这样会做许多无用功。比如11，比如从2到10判断9次才能确定下来。那能不能减少判断次数呢，许多人就会想到只要判断到n/2即可，后边的数不用判断。是的，只要不能被n/2以内的数整除，n就是质数。我们再来想想，一个数n能被另一个数m整除，整除结果为k。可以表示为n=m*k，m和k就称作m的因子，m和k的大小有三种情况，m<k, m=k, m>k，我们只需判断到m=k=sqrt(n)为止，后边的就无需判断，因为如果有大于sqrt(n)的因子那么另一个因子一定小于sqrt(n),小于sqrt(n)的情况已经判断过了，所以不用再次判断了，这样就能减少判断次数。

int prime_2(int num)
{
	int count=0;
	int cur=2;
	while (cur<num)
	{
		int maxcmp=(int)sqrt(cur);
		int factor=2;
		bool flag=true;  //设定标志位，如果不是质数置为false
		while(factor<=maxcmp)
		{
			if (cur%factor==0)
			{
				flag=false;
				break;
			}
			factor++;
		}
		if (flag)
		{
			cout<<cur<<" ";
			count++;
		}
		cur++;
	}
	return count;
}

利用一个性质：一个合数(非质数)都能分解为若干个比它小的质数的乘积。12=2*2*3, 15=3*5,如果我们把以前的质数保存起来，只需判断是否能被比它小的质数整除即可，结合上面的分析只需判断能否被不大于sqrt(n)的质数整除即可，按照这个思路，写出如下代码：

int prime_3(int num)
{
	int count=0;
	int *result=new int[num];    //记录找到的质数
	result[0]=2;
	int cur=3;
	while(cur<num)
	{
		bool flag=true;  //设定标志位，如果不是质数置为false
		int maxcmp=(int)sqrt(cur);
		int index=0;
		while(result[index]<=maxcmp)
		{
			if (cur%result[index]==0) //只判断能否被之前保存的质数整除
			{
				flag=false;
				break;
			}
			index++;
		}
		if (flag)
		{
			count++;
			result[count]=cur;
		}
		cur++;
	}
	for (int index=0;index<count+1;index++)
	{
		cout<<result[index]<<" ";
	}
	delete []result;
	return count;
}