Sieve of Eratosthenes
使用埃拉托斯特尼筛选法计算小于100000的素数。
埃拉托斯特尼筛选法是最为知名的产生素数的筛选法,适用于产生最小的N个素数。
该方法的唯一缺点是使用的存储空间大,可以进一步改进。
另外,该算法也不适用于计算某个范围内的全部素数。
C++版使用的标志是布尔标志,比起C语言版(用整数数组作为标识太浪费空间;用位运算逻辑太复杂,易出错),使用的空间上改进相当多,并且逻辑也十分简单。
这是一个基础程序,用到的时候可以拷贝修改加以利用。
程序如下:
/*
* 筛选法计算小于100000的素数
* 筛选法求最小的素数序列,原本是不朽的数学家们手工计算的方法,
*/
#include <iostream>
#include <cmath>
using namespace std;
const int MAXN = 100000;
bool sieveflag[MAXN+1] = {false, false, true};
void esieve(bool sflag[], int n)
{
// 初始化
for(int i=3; i<=n; i++) {
sflag[i++] = true;
sflag[i] = false;
}
// 筛选
int max = sqrt(n);
for(int i=3; i<=max; i++){
if(sflag[i]) {
for(int j=i+i; j <= n; j+=i)
sflag[j] = false;
}
}
}
int main()
{
esieve(sieveflag, MAXN);
for(int i=2, num=0; i<=MAXN; i++){
if(sieveflag[i]){
++num;
cout << num << ":" << i <<endl;
}
}
return 0;
}