ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

埃氏筛&欧拉筛~Biu~素数

2022-01-23 22:06:31  阅读:155  来源: 互联网

标签:prime 埃氏 int 合数 vis 素数 maxn Biu


两种方法筛素数

素数定义:大于0的数,除了1和他本身之外,没有其他数可以整除它。
最小的素数:2
合数定义:大于0的数,除了1和他本身外,还存在其他数可以整除它。
最小的合数:4

实际上合数和质数是相对立的。

埃氏筛:

先上代码:

#include<iostream>
#include<string.h>
using namespace std;
const int maxn=100;
int vis[maxn];
int main()
{
	memset(vis,1,sizeof(vis));//初始都设定为是素数
	vis[0]=vis[1]=0;//0 1 不是素数
	for(int i=2;i<=maxn;i++){
		if(vis[i]){
			for(int j=i*i;j<=maxn;j+=i){//这里为什么可以从i*i开始呢?后文有说明
				vis[j]=0;//所有是i的倍数的数都是合数,即不是素数
			}
		}
	}
	for(int i=0;i<=maxn;i++){
		if(vis[i]) cout<<i<<" ";
	}
	return 0;
}

说明:
为什么j可以从ii开始,而不是从i+i开始,是因为通过找规律可以知道,在i>2时, i(i-1)的数已经被筛完了,所以从i * i开始筛。
缺点:
埃氏筛的缺点很明显,例如12=26 =34,他奖会被筛多遍。这个时候欧拉筛就孕育而生,寻找最小质因子,每个合数只被最小质因子筛。

欧拉筛-euler

代码:

#include<bits/stdc++.h>
using namespace std;
const int maxn=100;
int vis[maxn],prime[maxn];

int main()
{
	memset(vis,1,sizeof(vis));
	memset(prime,0,sizeof(prime));
	vis[0]=vis[1]=0;
	for(int i=2;i<=maxn;i++){
		if(vis[i]){
			prime[0]++;//相当于一个记录素数的计数器
			prime[prime[0]]=i;//第i个数为prime
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<=maxn;j++){//只遍历素数表,也就满足了用最小素数筛的需求。
			vis[i*prime[j]]=0;
			if(i%prime[j]==0) break;//为什么这里需要break;思考一下
		}
	}
    for(int i=0;i<=maxn;i++){
        if(vis[i]) cout<<i<<" ";
    }
	return 0;
}

说明:
对于 i%prime[j] == 0 就break的解释 :当 i是prime[j]的倍数时,i = kprime[j],如果继续运算 j+1,i * prime[j+1] = prime[j] * k prime[j+1],这里prime[j]是最小的素因子,当i = k * prime[j+1]时会重复,所以才跳出循环。
举个例子 :i = 8 ,j = 1,prime[j] = 2,如果不跳出循环,prime[j+1] = 3,8 * 3 = 2 * 4 * 3 = 2 * 12,在i = 12时会计算。因为欧拉筛法的原理便是通过最小素因子来消除

完结!

标签:prime,埃氏,int,合数,vis,素数,maxn,Biu
来源: https://blog.csdn.net/CSDNqinzhike/article/details/122657591

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有