筛法求素数

中心思想：把从1开始的、某一范围内的正整数从小到大顺序排列， 1不是素数，首先把它筛掉。剩下的数中选择最小的数是素数，然后去掉它的倍数。依次类推，直到筛子为空时结束。

举例： 从1-1000中求素数。

            首先设置一个标记数组，把除了0和1之外的数字都设置成true。

            接着从2开始，凡是2的倍数且小于1000的数字，全部标记成false。然后选择下一个素数3。

            然后从3开始，凡是3的倍数且小于1000的数字，全部标记成flase。然后选择下一个素数5。

            ......

            直到开始的数字大于1000时，循环结束。

            最后剩下来被标记成ture的数字就是1-1000之间的素数了。

C++代码：

//没有优化后的筛法
#include"cstdio"  
#include"cstring"  
using namespace std;
#define MAX 100000//求MAX范围内的素数  
long long su[MAX], cnt;
bool isprime[MAX];
void prime()
{
	cnt = 1;//cnt用来计数
	memset(isprime, 1, sizeof(isprime));//初始化认为所有数都为素数  
	isprime[0] = isprime[1] = 0;//0和1不是素数  
	for (long long i = 2; i <= MAX; i++)//
	{
		if (isprime[i])
			su[cnt++] = i;//保存素数i  
		for (long long j = 2*i; j<MAX;j+=i)
		{
			isprime[j] = 0; //将i的倍数都置为flase
		}
	}
}
int main()
{
	prime();
	for (long long i = 1; i<cnt; i++)
		printf("%d  ", su[i]);
	return 0;
}

优化后的筛法：

//优化后的筛法，复杂度接近于O(n)
#include"cstdio"  
#include"cstring"  
using namespace std;
#define MAX 100000//求MAX范围内的素数  
long long su[MAX], cnt;
bool isprime[MAX];
void prime()
{
	cnt = 1;
	memset(isprime, 1, sizeof(isprime));//初始化认为所有数都为素数  
	isprime[0] = isprime[1] = 0;//0和1不是素数  
	for (long long i = 2; i <= MAX; i++)
	{
		if (isprime[i])
			su[cnt++] = i;//保存素数i  
		for (long long j = 1; j<cnt&&su[j] * i<MAX; j++)
		{
			isprime[su[j] * i] = 0;//筛掉小于等于i的素数和i的积构成的合数  
		}
	}
}
int main()
{
	prime();
	for (long long i = 1; i<cnt; i++)
		printf("%d  ", su[i]);
	return 0;
}