From: http://blog.csdn.net/dinosoft/article/details/5829550

===========================================================================

以上完整源码(a.cpp)

/*
功能: 	求[0, 20000)间的所有素数, 假设内存空间足够
环境:	Linux C++
编译:	g++ -o a a.cpp -Wall -Os
总结:	这两种算法在性能上差距不是很大, 当N个数较少时, 还是"一般线性筛选法"速度更快,
但是当N较大时, 快速线性筛选法的优势更加明显
*/
#include <stdio.h>
#include <string.h>
#define N 20000
// 一般线性筛选法: 会出现重复筛选同一个数
void make_prime(int primes[], int& cnt)
{
bool bPrime[N];		// 质数标志数组
cnt = 0;			// 素数个数
memset(bPrime, true, sizeof(bPrime));// 假设全是素数
bPrime[0] = false;					// 0: 非素数
bPrime[1] = false;					// 1: 非素数
for (int i = 2; i < N; i++)
{
if (bPrime[i])					// i是素数
{
primes[cnt++] = i;			// 将素数i保存到bPrimes[]中
// 作筛选: i的倍数都不是素数
for (int k = i * i; k < N; k += i)	// 将素数i的倍数全置为非素数标志
bPrime[k] = false;
}
}
}
// 快速线性筛选法: 不会出现重复筛选同一个数
void getPrimes(int primes[], int& cnt)
{
bool bPrime[N];		// 素数标志数组
cnt = 0;		// 素数个数
memset(bPrime, true, sizeof(bPrime));	// 假设全部为素数
bPrime[0] = false;	// 0: 非素数
bPrime[1] = false;	// 1: 非素数
for(int i = 2; i < N; i++)
{
if(bPrime[i])	// i是素数
primes[cnt++] = i;			// 保存素数i
// 作筛选: i与素数的乘积都不是素数
for(int j = 0; j < cnt && i * primes[j] < N; j++)
{
bPrime[i * primes[j]] = false;	// 置非素数标志
if(i % primes[j] == 0)		// i为素数的倍数
break;
}
}
}
int main()
{
int primes[N];		// 保存所有素数
int cnt = 0;		// 素数个数
#if 1
make_prime(primes, cnt);	// 调用一般线性筛选法
#else
getPrimes(primes, cnt);		// 调用快速线性筛选法
#endif
for(int i = 0; i < cnt; i++)
printf("primes[%d] = %d\n", i, primes[i]);
printf("\n素数个数cnt=%d\n", cnt);
return 0;
}