素数求解的N种境界

本期介绍🍖
主要介绍：如何快速筛查出素数的各种方法👀。

---

---

1. 前言

  素数又称质数(prime number)，一个大于1的自然数，除了1和它本身外，不能被其他自然数整除的数被称为素数。例题：请编写代码找出1——100之间的素数，并打印在屏幕上。

---

2. 试除法

  试除法：通过不断试除1和本身之外的自然数，来判断该数是否为素数的方法（只要有一个自然是能够整除该数，则其就不是素数；若一个都不能整除该数，则其就是素数）。

2.1 试除2~(n-1)

#include<stdio.h>

int main()
{
	//1.产生1~100之间的数字，1不是素数忽略
	int i = 0;
	printf("打印1~100之间的素数：>\n");
	for (i = 2; i <= 100; i++)
	{
		//2.判断是否是素数
		int j = 0;
		int flag = 1;//表示是素数
		//产生除了1和它本身之外的数字
		for (j = 2; j <= i-1; j++)
		{
			//判断能否被整除
			if (i % j == 0)
			{
				flag = 0;//表示不是素数
				break;
			}
		}
		if (flag == 1)
			printf("%d ", i);
	}
	return 0;
}

---

2.2 排除2倍数的合数

  对上面实现求素数的代码进行优化，排除所有2的倍数。因为除了2之外所有的2的倍数都不可能是素数，没有必要再进行判断。代码如下。

#include<stdio.h>

int main()
{
	//1.产生1~100之间的数字，并排除所有2的倍数
	int i = 0;
	//不能忘记2也是素数
	printf("打印1~100之间的素数：>\n");
	printf("%d ", 2);
	for (i = 3; i <= 100; i += 2)
	{
		//2.判断是否是素数
		//产生2~i-1之间的数字
		int j = 0;
		int flag = 1;//表示是素数
		for (j = 2; j <= i - 1; j++)
		{
			//判断能否被整除
			if (i % j == 0)
			{
				flag = 0;//表示不是素数
				break;
			}
		}
		if (flag == 1)
			printf("%d ", i);
	}
	return 0;
}

---

2.3 试除2~sprt(n)

  如果一个数不是素数的话，它必然可以被因式分解，分解成两个数相乘。反过来，假如想证明这个数不是素数的话，是不是只要证明它可以被分解的其中一个数整除就可以，没有必要证明两个数都能整除。
  多分解几组自然数，大家会发现分解出来的两个因数必然是一大一小，且都相互趋近于被分解的那个数的算数平方根。假设被分解的数为n，sqrt(n)表示n的算数平方根，那么n分解出的两个因数必然会出现在以sqrt(n)为分界线的两边。左侧因数的范围：2 ~ sqrt(n)之间，右侧因数的范围：sqrt(n) ~ n-1之间。
  所以想判断n是否是素数，只需要看在2 ~ sqrt(n)这个范围内，是否有数可以整除n。如果可以被整除，那就不是素数；如果不能被整除，那就是素数。程序这样设计相较于之前能够减少素数一半的试除，大大优化了试除法的代码。代码如下。

#include<stdio.h>
#include<math.h>

int main()
{
	//1.产生1~100之间的数字,1不是素数忽略
	int i = 0;
	printf("打印1~100之间的素数：>\n");
	for (i = 2; i <= 100; i++)
	{
		//2.判断是否是素数
		//产生2~sqrt(i)之间的数字，sqrt是开平方的意思
		int j = 0;
		int flag = 1;//表示是素数
		for (j = 2; j <= sqrt(i); j++)
		{
			//判断能否被整除
			if (i % j == 0)
			{
				flag = 0;//表示不是素数
				break;
			}
		}
		if (flag == 1)
			printf("%d ", i);
	}
	return 0;
}

---

3. 筛选法

  筛选法的由来可以追溯到古希腊时期，那时有个人叫埃拉托斯特尼。他发现可以通过在涂蜡的木板记下数字，接着在蜡上以记一个点来作为删去一个数字，然后从最小的质数开始一个接着一个的划掉这些质数的倍数，最后那些没被点标记的数就是质数。
  通过该方法查找完质数后，涂蜡板全是密密麻麻小点，看上去就像一个筛子，所以就把着这种方法简称：筛选法。如下视频所示。

3.1 如何实现筛除法求素数

  要实现筛除法，首先需要考虑的是怎么记录一个数是否已经被筛除掉的状态。可以创建一个bool类型的数组，数组的下标就是需要被排查的数，下标所对应的数组元素就是记录该数的状态。如果没有被筛除状态为true，则表示是素数；反之被筛除了状态为false，则表示不是素数。
  从公认的最小质数2开始，保留2然后筛除2的所有倍数。由于2的倍数已经全被筛除，所以下个数必然是质数3，保留3然后筛除3的所有倍数。由于2、3的倍数已经全被筛除，所以下个数必然是质数5，保留5然后筛除5的所有倍数。由于2、3、5的倍数已经全被筛除，所以下个数必然是质数7，保留7然后筛除7的所有倍数……
  按照上面的步骤循环往复，筛除剩余的，也就是bool状态为true的数就是质数。代码如下。

#include<stdio.h>
#include<stdbool.h>

#define NUM 101

int main()
{
	bool arr[NUM];
	int num = 0;
	int prime = 0;
	//初始化数组，让数组每个元素都为true，表示为素数
	int i = 0;
	for (i = 2; i <= NUM - 1; i++)
		arr[i] = true;
	//找到素数，然后筛除该素数的倍数
	for (prime = 2; prime <= NUM - 1; prime++)//2是公认的最小素数
	{
		if (arr[prime] == true)
		{
			//筛除prime的倍数
			int num = 0;
			for (num = 2 * prime; num <= NUM - 1; num += prime)
			{
				arr[num] = false;//表示不是素数
			}
		}
	}
	//找出那些没有被筛除的数，就是素数
	printf("打印1~100之间的素数：>\n");
	for (i = 2; i <= NUM - 1; i++)
	{
		if (arr[i] == true)
			printf("%d ", i);
	}
	return 0;
}

---

3.2 优化代码

  一般在筛查质数的倍数，这个倍数会从2倍开始逐级增加。但是在使用过程中，会发现有些数字会重复筛查好多次。就譬如数字6，在筛查质数2的倍数时已经筛过了，但在筛查质数3的倍数时还会再筛一次，这就导致重复了。为了优化这种情况，可以从质数的平方处开始往后筛查，这样就减少了一部分的重复筛查。就譬如：筛查质数5的倍数，由于5的(2倍、3倍、4倍)已经在之前筛查过了，从5的自身倍数处（即：5*5）往后筛查即可。

#include<stdio.h>
#include<stdbool.h>

#define NUM 101

int main()
{
	bool arr[NUM];
	int num = 0;
	int prime = 0;
	//初始化数组，让数组每个元素都为true，表示为素数
	int i = 0;
	for (i = 2; i <= NUM - 1; i++)
		arr[i] = true;
	//找到素数，然后筛除该素数的倍数
	for (prime = 2; prime <= NUM - 1; prime++)//2是公认的最小素数
	{
		if (arr[prime] == true)
		{
			//筛除prime的倍数
			int num = 0;
			for (num = prime * prime; num <= NUM - 1; num += prime)
			{
				arr[num] = false;
			}
		}
	}
	//找出那些没有被筛除的数，就是素数
	printf("打印1~100之间的素数：>\n");
	for (i = 2; i <= NUM - 1; i++)
	{
		if (arr[i] == true)
			printf("%d ", i);
	}
	return 0;
}

---

4. 总结

  可以看出筛选法和试除法其实有着本质上的区别。试除法是判断每一个数是不是素数来达到目的，而筛选法却是将不是素数的数全部去除，然后得到余下的数来达到目的。当把找素数的范围提到1 ~ 10000以内，筛选法的执行次数远小于试除法的执行次数。所以可以看出，求素数的范围越大筛选法的优越性就越高。

---

这份博客👍如果对你有帮助，给博主一个免费的点赞以示鼓励欢迎各位🔎点赞👍评论收藏⭐️，谢谢！！！
如果有什么疑问或不同的见解，欢迎评论区留言欧👀。