普通筛选:
#include<iostream>
#include<bits/stdc++.h>
using namespace std;
int main(){
int a[100];
int fag=0;//flg=0为素数,1为合数
int cnt=0;
for(int i=2;i<=100;i++){
fag=0;
for(int j=2;j<i;j++){
if((i%j)==0)fag=1;
}
if(fag==0)a[cnt++]=i;
}
for(int i=0;i<100;i++){
cout<<a[i]<<endl;
}
}
此方法虽然可以正确的筛选出素数,但筛选的时间复杂度为过高为n^2
若第二层循环枚举到√i,则时间复杂度为n√n依然过高。
埃式筛法
埃式筛法的关键在于提前声明一个标志全0数组fag
第一层循环每一个数时,把这个数的1~n倍(在给定范围内)对应的fag数组中的值改为1,而在判断素数时,若该位fag为0则为素数,为1则为合数。
#include<iostream>
#include<bits/stdc++.h>
using namespace std;
int main(){
int a[100]={0};
int fag[100]={0};//flg=0为素数,1为合数
int cnt=0;
for(int i=2;i<=100;i++)
{
if(fag[i]==0)a[cnt++]=i;
for(int j=i;j<=100;j+=i)
{
fag[j]=1;
}
}
for(int i=0;i<100;i++)
{
cout<<a[i]<<endl;
}
}
埃式筛法的时间复杂度为nlogn,依然不低。
欧拉筛法
#include<iostream>
#include<bits/stdc++.h>
using namespace std;
int main(){
int a[100]={0};
int fag[100]={0};//flg=0为素数,1为合数
int cnt=0;
for(int i=2;i<=100;i++)
{
if(fag[i]==0)
a[cnt++]=i;
for(int j=0;j<cnt;j++)
{
if(i*a[j]>=100)//判断是否越界
break;
fag[i*a[j]]=1;//筛选合数,并使fag该数变为1
if(i%a[j]==0)//欧拉筛选的灵魂所在
break;
}
}
for(int i=0;i<100;i++)
{
cout<<a[i]<<endl;
}
}
在埃式筛法中我们不难看出有很多合数被反复筛选了许多次,造成了浪费。
而欧拉筛法的关键在于第二层循环中的第二个break,如果这个数被之前的素数筛选过,则直接跳出循环,进行下一个数的运算。
而为什么只用判断是否被之前的素数筛选过,是因为任意一个合数都可以被某一个素数整除。
证明如下:
按照合数的定义,一个合数m 一定能被某个整数n整除,假设商是p,那么 m=np, n和p显然都比m小
如果n和p至少有一个是质数,那么结论就得到证明
如果n和p都不是质数,还是合数,那么n一定可以再分解成两个小于n的整数的乘积,n=n1p1
于是 m=n1p1p
如果n1,或者p1中有一个是质数,那么结论就得到证明了
如果n1,p1都不是质数,那么继续分解因数分下去。。。。
因为m是有限整数,,而分出的因数n1,p1,…等一个比一个小,最后总有一次得到一个因数是质数