素数筛法

筛法是一种简单检定素数的算法。据说是古希腊的埃拉托斯特尼（Eratosthenes，约公元前274～194年）发明的，又称埃拉托斯特尼筛法（sieve of Eratosthenes）。

具体做法是：给出要筛数值的范围

   

n，找出

   

n以内的素数p 1，p 2，p 3，......，p k。先用2去筛，即把2留下，把2的倍数剔除掉；再用下一个素数，也就是3筛，把3留下，把3的倍数剔除掉；接下去用下一个素数5筛，把5留下，把5的 倍数剔除掉；不断重复下去......。

因为 希腊人是把数写在涂腊的板上，每要划去一个数，就在上面记以小点，寻求 质数的工作完毕后，这许多小点就像一个筛子，所以就把 埃拉托斯特尼的方法叫做“ 埃拉托斯特尼筛法”，简称“筛法”。

详细列出算法如下：

1. 列出2以后的所有序列：
   • 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
2. 标出序列中的第一个素数，也就是2，序列变成：
   • 23 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
3. 将剩下序列中，划摽2的倍数（用删除线标出），序列变成：
   • 2 3 45 67  89 1011  12 13  14 15  16 17  18 19  20 21  22 23  24 25
4. 如果现在这个序列中最大数小于最后一个标出的素数的平方，那么剩下的序列中所有的数都是素数，否则回到第二步。
5. 本例中，因为25大于2的平方，我们返回第二步：
6. 剩下的序列中第一个素数是3，将主序列中3的倍数划出（删除线 ），主序列变成：
   • 2 3 4  6 7 8  9 10 11  12 13 14  15 16 17  18 19 20  21 22 23  24 25
7. 我们得到的素数有：2，3
8. 25仍然大于3的平方，所以我们还要返回第二步：
9. 现在序列中第一个素数是5，同样将序列中5的倍数划出，主序列成了：
   • 2 3 4 5 6 7 8 9  10 11 12 13 14  15 16 17 18 19  20 21 22 23 24  25
10. 我们得到的素数有：2 3 5 。
11. 因为25等于5的平方，跳出循环.

结论：去掉划线的数字，2到25之间的素数是：2 3 5 7 11 13 17 19 23。

具体就是从２开始，筛掉２的倍数，接着３，筛掉３的倍数，接着４，由于４是２的倍数，所以直接跳过，接着５，筛掉５的倍数，以此类推，最后得出２－ｎ区间内的所有素数。

#include<iostream>
using namespace std;
#include<cstring> 
#include<cstdio> 
#include<cmath>
const int MAXN = 2000001;
bool a[MAXN];
void get_Prime()
{
	memset(a, true, sizeof(a));
	a[0] = a[1] = false;
	for(int i = 2; i <= sqrt((double)MAXN); ++i)
	{
		if(a[i])
		for(int j = i*i; j <＝ MAXN; j += i)
		{
			
			a[j] = false;
		}
	}
}
int main()
{
	int m;
	get_Prime();
	while(cin >> m, m)
	{
		for(int i = 1; i <= m; ++i)
		{
			if(a[i])
			{
				cout << i << ' ';
			}
		}
		putchar('\n');
		cout << endl;
	}
	return 0;
}