素数

素数又称质数，是指除1和本身之外，不能被其他数整除的一类数。即对给定的整数你，如果对任意的正整数a（1<a<n)，都有n%a!=0成立，那么称n是素数，否则为合数。其中，1既不是素数，也不是合数。

1.素数的判断

从定义得，整数n是否为素数，需要判断n能否被2，3，4…n-1中的一个整除，这时候有算法的复杂度为O(n)，但是有个证明：

设n的约数为k，即n%k=0；
那么k *（n / k）== n, 可知 n/k也是n的一个约数，
所以 k 与 k / n中一定满足：
其中一个小于等于sqrt（n），另一个大于等于sqrt（n）

sqrt（n）为根号n。
这告诉我们，现在只需判断n能否被2，3，…,sqrt(n)中的一个整除（sqrt(n)向下取整），即可判定n是否为素数。改算法复杂度为O(sqrt（n）)。

2.素数表的获取

现在输出100以内的素数，代码如下：

#include <cstdio>
#include <math.h>
const int maxn = 101;
int prime[maxn], pNum = 0;     //prime存储数字 , pNum为素数个数 
bool p[maxn] = {0};			   //p[i] == true表示i是素数； 

bool isPrime(int n){ //判断n是否为素数 
	if(n <= 1) return false;  // 特判
	int sqr = (int)sqrt(1.0 * n);//需要添加<math.h>头文件。由于sqrt参数要求是浮点数，需要在n前面乘以1.0，并强制类型转换
	for(int i = 2; i <= sqr; i++){
		if(n % i == 0) return false;
	}
	return true;
}
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;
}

这种方法是将素数保存在prime[pNum++]中，p[i]可以表示i是否为素数。
此方法枚举部分的复杂度为O(n),判断素数的复杂度为O(√n)，总的复杂度为O(n√n)。但是对n不超过10^5的大小没有问题。

埃氏筛法
此为更高效的算法，时间复杂度为O(nloglogn),可以适用于更大的范围。
如果求1~15中的素数
①.2是素数，去掉2的倍
留下： 2 3 5 7 9 11 13 15
②.3没有被删去，是素数，去掉3的倍数
留下：2 3 5 7 11 13 15
③.4被①删去了，因此4不是素数
④.5没有被删去，是素数，删去所有5的倍数
留下：2 3 5 7 11 13
⑤.6已经被删去了，6不是素数
。。。。
直至得到：2 3 5 7 11 13。
代码如下：

#include <cstdio>
const int maxn = 101; //求取100以内的所有素数 
int prime[maxn], pNum = 0;//prime存储数字 , pNum为素数个数
bool p[maxn] = {0};//如果i为素数，p[i]为false，否则，为true

void Find_Prime(){
	for(int i = 2; i < maxn; i++){
		if(p[i] == false){ //如果为素数
			prime[pNum++] = i;//把素数i存到prime数组中
			for(int j = i + i; j < maxn; j += i){//筛去所有i的倍数
				p[j] = true;
			}
		}
	}
}
int main(){
	Find_Prime();
	for(int i = 0; i < pNum; i++){
		printf("%d ", prime[i]); 
	}
	return 0;
} 

3.真实例题

令Pi表示第i个素数。现任给两个正整数M <= N <= 104，请输出PM到PN的所有素数。

输入格式：

输入在一行中给出M和N，其间以空格分隔。

输出格式：

输出从PM到PN的所有素数，每10个数字占1行，其间以空格分隔，但行末不得有多余空格。

输入样例：

5 27

输出样例：

11 13 17 19 23 29 31 37 41 43
47 53 59 61 67 71 73 79 83 89
97 101 103

代码如下：

#include <cstdio>
const int maxn = 10000001; 
int prime[maxn], num = 0;
bool p[maxn] = {0};

void Find_prime(int n){
	for(int i = 2; i< maxn; i++){
		if(p[i] == false){
			prime[num++] = i;
			if(num >= n) break;
			for(int j = i + i; j < maxn; j += i){
				p[j] = true;
			}
		}
	}

}
int main(){
	int m, n, count = 0;
	scanf("%d%d", &m, &n);
	Find_prime(n);
	for(int i = m; i <= n; i++ ){
		printf("%d", prime[i - 1]);  //从下标0开始 
		count++;
		if(count % 10 != 0 && i < n) printf(" "); 
		else printf("\n");
	}
	return 0;
}