埃拉托斯特尼(Eratosthene)筛法

 

下面是代码：

 

/* 功能：将1-n的素数输出
   使用了埃拉托斯特尼(Eratosthene)筛法
*/ 

#include <iostream>
#include <math.h>

using namespace std;

int prime[10000000];//存放需要筛选数
int sieve (int *prime,int n,int k)//n是1-n这个范围，k是当前把prime[k]作为筛子
{
    for (int i= k+1; i<=n ; i++)
    {
        //当前数prime[i]不为0且能整出筛子数，说明肯定是合数
        if (prime[i] && prime[i]%prime[k]== 0)
        {
            prime[i] = 0;//将合数赋值为0
        }
    }
    k++;
    //如果当前筛子数是0，则移到下一位
    while (!prime[k])
    {
        k++;
    }
    return k;
}

int main()
{
    int  n;
    while (cin >> n && n)
    {
    int k = 2;//从素数2开始
    for (int i = 1; i <= n; i++)
    {
        prime[i] = i;
    }

    while (prime[k]<= sqrt(n))
    {
        k = sieve(prime,n,k);//找到下一个筛子
    }
     for (int i = 2; i<= n; i++)
     {
         if (prime[i]!= 0)
         {
            cout << prime[i] <<" ";
         }
     }
    }
    return 0;
}

当然，这是一个非常初级的版本，如果你想找出你给定数组的素数，只要稍微改正即可。