判断素数的几种求法

总结一下判断素数的几种求法

素数的概念

素数及除了1和本身之外，不能被其他数整除的一类数，1既不是素数也不是合数

素数的判断

判断n是否能被2，3，…n-1中的一个整除，只有全部都不能被整除才能判断为，而只要其中一个可以被整除，就可以判断为非素数。

由于该算法的复杂度为o(n)，而且求素数通常为整个程序中的一部分，故考虑进行优化

设数n在2-n-1中存在n的约数k，那么由于k*(n/k) ==n，n/k也为其中的一个约数，且k与n/k中一定满足其中一个小于sqrt(n)，另一个大于sqrt(n)（该式子表示根号的意思）。所以只需要判断n能否被2,3…sqrt(n)中的一个整除，即可判断n是否为素数，而该算法的时间复杂度变为O(sqrt(n))

代码如下：

bool isPrime(int n){
	if(n <= 1) return false; //特判
	int sqr = (int)sqrt(1.0 * n); //根号n
	for(int i = 2, i <= sqr; i++){
		if(n % i == 0) return false;
	}
	return true;
}

这个方法如果不会超过int型变量的范围上界的话可以用以下更加简洁的方法

bool isPrime(int n){
	if(n <= 1) return false;
	for（int i = 2; i * i < n; i++）{
		if(n % i == 0) return false;
	}
	return true;
}

这样写可能会导致n接近int型变量的范围上界时导致i*i溢出

素数表的获取

素数表的获取即判断1-n内的数是否为素数，如果是则加入素数表中

代码如下：

//获取1-100的素数表
#include<stdio.h>
#include<math.h>
bool isPrime(int n){
	if(n <= 1) return false;
	for（int i = 2; i * i < n; i++）{
		if(n % i == 0) return false;
	}
	return true;
}
const int maxn = 101;//表长
int prime[maxn], pNum = 0;//prime数组存放所有素数，pNum为素数个数
bool p[maxn] = {0}; //p[i] == true 表示i为素数
void Find_Prime(){
	for(int i = 1; i < maxn; i++){
		if(isPrime(i) == true){
			prime[pNum++] = i;
			p[i] = true;
		}
	}
}

int main(){
	Find_Prime();
	for(int i = 0; i < pNum; i++){
		printf("%d", prime[i]);
	}
	return 0;
}

上面的算法对于n在10^5内都是可以承受的，但是如果出现更大范围的素数表，o(n*sqrt(n))就可能会过于复杂，因此有一个更加高效的算法，叫做“筛法”

筛法

引用《算法笔记》书中对该方法的讲解

代码如下：

#include<stdio.h>
const int maxn = 101;
int prime[maxn], pNum = 0;
bool p[maxn] = {0};
void Find_Prime(){
	for(int i = 2; i < maxn; i++){
		if(p[i] == false){
			prime[pNum++] = i;
			for(int j = i + i, j < maxn; j += i){
				p[j] = true;
			}
		}
	}
}

int main(){
	Find_Prime();
	for(int i = 0; i < pNum; i++){
		printf("%d", prime[i]);
	}
	return 0;
}

此种算法的时间复杂度就只有O(n * loglogn)