求解素数的n重境界（筛选法，试除法）

目录

试除法 

试除法代码块

对筛选法的解释

试除法与筛选法的区别 

筛选法代码块

---

试除法

假设要判断n是否为素数，试除法的操作为：用n去除 2 到 n - 1 ,如果发现有一个可以被n整除，那么n不是素数，如果这些数没有一个可以被n整除，那么n为素数。

在用试除法来判断时，我们所写出的代码也有几种境界。

>创建变量flag来辅助判断

>在创建变量flag的基础上过滤掉已知不是素数的偶数

>在前两种方法的基础上引入【数学定理】来极大缩小试除范围

【运用到的数学定理】

一个数若可以进行因数分解，那么分解时得到的两个数一定是一个小于等于sqrt(n)，一个大于等于sqrt(n)，所以对于每个数n，并不需要从2判断到n-1，遍历到sqrt(n)即可。因为若sqrt(n)左侧找不到约数，那么右侧也一定找不到约数。（sqrt是开平方）

注意：在使用库函数sqrt()时要引头文件<math.h>

试除法代码块

第一重境界：

#define _CRT_SECURE_NO_WARNINGS 1
#include<stdio.h>

//输出100到200之间的素数 
int main()
{
	int i = 0, count = 0;
	for (i = 100; i <= 200; i++)//产生100到200的数
	{
		int flag = 1;//辅助和优化判断过程
		int j = 0;
		for (j = 2; j < i; j++)
		{
			if (i % j == 0)
			{
				flag = 0;
				break;
			}
		}
		if (flag == 1)//出了内层for循环以后如果flag不被改为0，说明是素数
		{
			printf("%d ", i);
			count++;
		}

	}
	printf("\ncount = %d", count);

	return 0;
}

第二重境界：

#include<stdio.h>

//输出100到200之间的素数 
int main()
{
	int i = 0, count = 0;
	for (i = 101; i <= 200; i+=2)//过滤掉明知不可能是素数的的偶数
	{
		int flag = 1;//辅助和优化判断过程
		int j = 0;
		for (j = 2; j < i; j++)
		{
			if (i % j == 0)
			{
				flag = 0;
				break;
			}
		}
		if (flag == 1)//出了内层for循环以后如果flag不被改为0，说明是素数
		{
			printf("%d ", i);
			count++;
		}

	}
	printf("\ncount = %d", count);

	return 0;
}

第三重境界：

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

int main()
{
	int i = 0, count = 0;
	for (i = 101; i <= 200; i += 2)//过滤掉明知不可能是素数的的偶数
	{
		int flag = 1;//辅助和优化判断过程
		int j = 0;
		for (j = 2; j <= sqrt(i); j++)//试除范围缩小为2到sqrt(i)
		{
			if (i % j == 0)
			{
				flag = 0;
				break;
			}
		}
		if (flag == 1)//出了内层for循环以后如果flag不被改为0，说明是素数
		{
			printf("%d ", i);
			count++;
		}

	}
	printf("\ncount = %d", count);

	return 0;
}

对筛选法的解释

1既不是质数，也不是合数，因此首先就要把1划去，2是最小的素数，所以2不能筛去，但是区间内2的倍数要全部筛去，接着3是素数，3也不能筛去，同理3的倍数要全部筛去，然后留下5,7，筛去其全部倍数。

试除法与筛选法的区别 

筛选法和试除法其实有着本质上的区别，试除法是判断每一个数是不是素数来达到目的；而筛选法不是如此，筛选法是将不是素数的数全部去除，然后得到余下的数为素数 ，实现素数的查找。

筛选法代码块

第一重境界：

#include<stdio.h>

//输出100到200之间的素数
int main()
{
	int i = 0;
	int count1 = 0;//用于统计素数的个数
	int count2 = 0;//用于统计筛选的次数
	int arr[201] = { 0 };//创建数组
	for (i = 101; i <= 200; i++)//初始化数组
	{
		arr[i] = i;//下标与元素值相等
	}
	for (i = 101; i <= 200; i++)//外层for循环是为了产生下标
	{
		if (arr[i])//能进来说明是素数 注意，内层的for循环是嵌套在if语句里的
		{
			int j = 0;//j为数组元素下标，下面将通过j访问数组元素
			if (i >= 100 && i <= 200)//下标与元素值相等
			{
				printf("%d ", arr[i]);//打印素数
				count1++;//统计素数的个数
			}
			for (j = i + i; j <= 200; j += i)//找到arr[i]的倍数，因为下标与元素值相等，故素数 
											 //下标的倍数即素数本身的倍数
			{
				arr[j] = 0;//使其不满足if语句的判断条件
				count2++;//统计筛选的次数
			}
		}
	}
	printf("\ncount1 = %d\n", count1);
	printf("count2 = %d\n", count2);

	return 0;
}

显然，上述执行的结果是不对的。错误的原因在于筛选法适用于从一开始就筛选，即从2开始就筛选，上述代码是从101开始筛选，筛掉101的倍数，很显然，仅仅是101的2倍就已经超出了所要求查找的范围，因次导致错误筛选。

正确的代码如下：

#include<stdio.h>

//输出100到200之间的素数
int main()
{
	int i = 0;
	int count1 = 0;//用于统计素数的个数
	int count2 = 0;//用于统计筛选的次数
	int arr[201] = { 0 };//创建数组
	for (i = 2; i <= 200; i++)//初始化数组
	{
		arr[i] = i;//下标与元素值相等
	}
	for (i = 2; i <= 200; i++)//外层for循环是为了产生下标
	{
		if (arr[i])//能进来说明是素数 注意，内层的for循环是嵌套在if语句里的
		{
			int j = 0;//j为数组元素下标，下面将通过j访问数组元素
			if(i>=100 && i<=200)//下标与元素值相等
			{
				printf("%d ", arr[i]);//打印素数
				count1++;//统计素数的个数
			}
			for (j = i + i; j <= 200; j += i)//找到arr[i]的倍数，因为下标与元素值相等，故素数 
                                             //下标的倍数即素数本身的倍数
			{
				arr[j] = 0;//使其不满足if语句的判断条件
				count2++;//统计筛选的次数
			}
		}
	}
	printf("\ncount1 = %d\n", count1);
	printf("count2 = %d\n", count2);

	return 0;
}

在上述代码中，我们不难发现6既是2的倍数，又是3的倍数，因此6被重复筛选，为了提高效率，我们可以从素数的平方开始往后筛选，减少了被重复筛选的数字，尤其是在大范围筛选中得到体现，极大提高了程序的效率。请看代码，仅仅是在上述第一重境界上稍作修改。

#include<stdio.h>

//输出100到200之间的素数
int main()
{
	int i = 0;
	int count1 = 0;//用于统计素数的个数
	int count2 = 0;//用于统计筛选的次数
	int arr[201] = { 0 };//创建数组
	for (i = 2; i <= 200; i++)//初始化数组
	{
		arr[i] = i;//下标与元素值相等
	}
	for (i = 2; i <= 200; i++)//外层for循环是为了产生下标
	{
		if (arr[i])//能进来说明是素数 注意，内层的for循环是嵌套在if语句里的
		{
			int j = 0;//j为数组元素下标，下面将通过j访问数组元素
			if (i >= 100 && i <= 200)//下标与元素值相等
			{
				printf("%d ", arr[i]);//打印素数
				count1++;//统计素数的个数
			}
			for (j = i * i; j <= 200; j += i)//找到arr[i]的倍数，因为下标与元素值相等，故素数 
											 //下标的倍数即素数本身的倍数
			{
				arr[j] = 0;//使其不满足if语句的判断条件
				count2++;//统计筛选的次数
			}
		}
	}
	printf("\ncount1 = %d\n", count1);
	printf("count2 = %d\n", count2);

	return 0;
}

从323次筛选降到232次筛选，这里所求区间太短，还不是很明显，当求更长区间时，更能发现高效之处。