下面是代码:
/* 功能:将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;
}
当然,这是一个非常初级的版本,如果你想找出你给定数组的素数,只要稍微改正即可。