算法基本思想是找出(设置为假)所有能整除已知除数的那些位置,它们不是素数,则所剩下的位置都将是素数。为实现该算法,用循环一次检查每个位置(只检查2~sqrt(n)就行了),在那些位置设为真(还没有被有被素数的候选集合中排除的)的位置中找出该数字的所有倍数。当完成最外层循环时,就得到了筛子中的所有素数。
CODES:
#include<vector>
#include<iostream>
using namespace std;
int main()
{
//创建筛子,初始化为真
const int sievesize=100;
vector<bool> sieve(sievesize,true);
//现在查找为真的位置
for(int i=2;i*i<sievesize;i++)
if(sieve[i])
//找出倍数
for(int j=i+i;j<sievesize;j+=i)
sieve[j]=false;
//输出集合中剩下的值
for(int i=2;i<sievesize;i++)
if(sieve[i]) cout<<i<<"is prime\n";
return 0;
}