找素数 埃氏筛法详解 pat B1013数素数

题目在此

思路

本题主要采用的“埃氏筛法”是众多筛法中最简单且最容易理解的一种，即Eratosthenes筛法。更优的欧拉筛法可以达到O(n)的时间复杂度，此处不与赘述。
introduction：
素数筛法的关键就在一个“筛”字。算法从小到大遍历所有数，对每一个素数，筛去它的所有倍数，剩下的就都是素数了。可能有读者问，一开始并不知道哪些是素数，何来的“对每一个素数”呢？

下面看一个例子：求1~15中所有的素数。

1. 2是素数（唯一需要事先确定），因此筛去所有2的倍数，即4、6、8、10、12、14。

   2 3 4 5 6 7 8 9 10 11 12 13 14 15

2. 3没有被前面的步骤筛去，因此3是素数，筛去所有3的倍数，即6，9，12，15

   2 3 4 5 6 7 8 9 10 11 12 13 14 15

3. 4不在了，看下一个

4. 5还在，筛去10，15。

   2 3 4 5 6 7 8 9 10 11 12 13 14 15

5. 6已经被筛了，看下一个。

6. 7还在，筛去14。
   2 3 4 5 6 7 8 9 10 11 12 13 14 15

…以此类推最后就剩下了2，3，5，7，11，13，筛选完成。

上面的例子就可以发现，当从小到大达到某个数a时，如果a没有被前面步骤的数筛去，那么a一定是素数。这是因为，如果a不是素数，那么a一定有小于a的素因子，这样在之前的步骤中a一定会被筛掉，所以，如果当遍历到a时还没有被筛掉，那么a一定是素数。

至于“筛”这个动作的实现，可以使用一个bool类型的数组p来标记，如果a被筛掉，那么p[a]为true；否则，p[a]为false。在程序开始时可以初始化p数组全为false。其算法的复杂度为O（nloglogn）

好了，埃氏筛法介绍完，该说题目思路了。题目就是要输出从2开始的第n个素数到第m个素数。注意点：

1. 用筛法，或者非筛法都可以解决这个题，在算法中只需要添加一句控制素数个数的语句：

   if(pnum >= n) break;
   

2. 小技巧安利：由于空格在测试的时候看不出来，因此如果提交返回“格式错误”，可以把程序中的空格改成其他符号（比如#）来输出，看看是哪里多了个空格。

3. 考虑到不知道第10⁴个素数有多大，不妨将测试上限设置得大一些，反正在素数个数超过n时是会中断的，不影响时间复杂度。当然也可以先用程序测试下第10⁴个素数是多少，然后再用这个数作为上限。

4. 本题在素数表生成过程中其实就可以直接输出，不过看起来会显得比较冗杂，因此还是应先生成完整素数表，然后再按照格式输出。

5. find_prime()函数中要记得是i < maxn而不是i <= maxn,否则程序一运行就会崩溃，在main函数中要记得调用find_prime()函数，否则就不会出结果。

AC代码：

#include <cstdio>
const int maxn = 1000001;
bool p[maxn]={0};
int pnum =0,prime[maxn];
void Find_Prime(int n){
	for(int i=2;i<maxn;i++){
		if(p[i]==false)
		{
			prime[pnum++] = i;
			if(pnum>=n) break;
			for(int j = i+i;j<maxn;j+=i){
				p[j]=true;
			}
		}
	}
}
int main(){
	int count=0,m,n;
	scanf("%d %d",&m,&n);
	Find_Prime(n);
	for(int i=m;i<=n;i++){
		printf("%d",prime[i-1]);
		count++;
		if(count%10!=0 && i<n) printf(" ");
		else printf("\n");
	}
	return 0;
}