通过质数筛筛取素数是较常用的方法,而欧拉筛是一种效率较高的算法。
下面介绍这种筛如何实现。
#include <stdio.h>
#include <stdbool.h>
bool flag[100];//标记素数
int prime[100];//储存素数
int count;
void primes(int n){
int i=0;
for(i=2;i<=n;i++){
if(!flag[i])
prime[++count]=i;
int j=0;
for(j=1;j<=count&&prime[j]<=n/i;j++){
flag[i*prime[j]]=true;
if(i%prime[j]==0)//注意
break;
}
}
}
int main(){
int n=0;
scanf("%d",&n);
primes(n);
int i=1;
while(prime[i]){
printf("%d\t",prime[i]);
i++;
}
return 0;
}
(这里解释一下,由于在c语言的标准库中没有bool类型,所以需要使用<stdbool.h>这个头文件。)
欧拉筛是通过运用合数可以分解成两个数的乘积,如16,分解为1*16,2*8,4*4的乘积,其中存在一种情况为一个质数(最小质因子)与另一个数的乘积,比如这里的例子是2*8,2为最小质因子。因此可以用这种方法将合数筛出,剩余的便为质数。
代码中标注的注意位置的代码是为了避免多次对一个数的筛取。比如,当i=10进行筛数时,此时prime[]中已存有{2,3,5,7,9},则能筛掉的数有{20,30,50,70,90},但是它们的最小质因子为2,也就是说它们可以在i=10,i=15,i=25,i=35,i=45时用最小质因子2进行筛掉,为避免重复筛取它们,现在仅先将20筛出,然后break出来,后面几个数在后续运行中逐个筛掉。