筛选素数的更高级做法(埃式筛法、欧拉筛法)

三种筛选0-100间素数的方法

普通筛选:

#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=n1
p1
于是 m=n1p1p
如果n1,或者p1中有一个是质数,那么结论就得到证明了
如果n1,p1都不是质数,那么继续分解因数分下去。。。。
因为m是有限整数,,而分出的因数n1,p1,…等一个比一个小,最后总有一次得到一个因数是质数

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值