gcd 素数埃式筛法欧式筛法（模板）

那就随便一点

已于 2022-08-27 10:01:34 修改

阅读量389

点赞数 1

分类专栏： # 基础算法模板文章标签： c++ 算法数据结构

于 2022-08-27 09:25:20 首次发布

本文链接：https://blog.csdn.net/m0_65508678/article/details/126554181

版权

基础算法模板专栏收录该内容

5 篇文章 1 订阅

订阅专栏

//求最大公约数。
int gcd(int a,int b){
	return !b ? a : gcd(b,a % b);//如果b = 0就返回 a (因为递推下去一定会使 b = 0，最后得到的a就是答案)，
}
//求最小公倍数。
int lcm(int a,int b){
	return a * b / gcd(a,b);
}
//判断素数。
int judge(int p){
	if(p == 1) return 0;
	for(int i = 2;i * i <=p;i ++ ){//平方根优化。
		if(p % i == 0) return 0;
	}
	return 1;
}
//埃氏筛法
const int N	 = 1e6 + 10;
int vis[N]{};//看是否是非素数。
int prime[N]{};
int tot = 0;
void prime_cal(int n){
	for(int i = 2;i <= n;i ++ ){
		if(!vis[i]){
			prime[ ++tot] = i;
		}
		for(int j =i + i;j<= n;j += i){
			vis[j] = 1;
		}
	}
}
//欧式筛法
int v[maxn], prime[maxn];  // v[i]是i的最小质因子；prime数组保存所有质数
void findPrime(int n)
{
	memset(v, 0, sizeof(v));  
	tot = 0;  // 质数数量
	for (int i = 2; i <= n; i ++)
	{
		if (!v[i]) {v[i] = i; prime[tot++] = i;}  // i是质数
		// 给当前的数i乘上一个质因子
		for (int j = 0; j < tot; j ++)
		{
			// 给当前的数i乘上一个质因子
			// i有比prime[j]更小的质因子，或者超出n的范围，停止循环
			if (prime[j] > v[i] || prime[j] > n/i) break;
			// prime[j]是合数i*prime[j]的最小质因子
			v[i*prime[j]] = prime[j];
		}
	}
}
另外：
#include<bits/stdc++.h>
using  namespace std;
//时间复杂度（n)
//适用低于10^7
const int MAX = 10000000;
void FindPrime(vector<int>& prime);
int main() {
	vector<int>prime;

	FindPrime(prime);
	for (auto i : prime)
		printf("%d ", i);
	return 0;
}
void FindPrime(vector<int>& prime) {
	array<int, MAX>flag;
	fill(flag.begin(), flag.end(), 0);
	for (int i = 2; i < MAX; i++) {
		if (flag[i]==0) {
			prime.push_back(i);
		}
		for (int j = 0; j < prime.size() && i * prime[j] < MAX; j++) {
			flag[i * prime[j]] = 1;
			if (i % prime[j] == 0) {
				break;
			}
		}
	}
}

埃式筛法

原理：素数的倍数肯定不是素数

算式

要得到自然数n以内的全部素数，必须把不大于

的所有素数的倍数剔除，剩下的就是素数。

给出要筛数值的范围n，找出以内的素数。先用2去筛，即把2留下，把2的倍数剔除掉；再用下一个质数，也就是3筛，把3留下，把3的倍数剔除掉；接下去用下一个质数5筛，把5留下，把5的倍数剔除掉；不断重复下去......。

步骤

详细列出算法如下：

列出2以后的所有序列：
- 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
标出序列中的第一个素数，也就是2，序列变成：
- 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
将剩下序列中，划掉2的倍数，序列变成：
- 2 3 5 7 9 11 13 15 17 19 21 23 25
如果这个序列中最大数小于最后一个标出的素数的平方，那么剩下的序列中所有的数都是素数，否则回到第二步。
本例中，因为25大于2的平方，我们返回第二步：
剩下的序列中第一个素数是3，将主序列中3的倍数划掉，主序列变成：
- 2 3 5 7 11 13 17 19 23 25
我们得到的素数有：2，3
25仍然大于3的平方，所以我们还要返回第二步：
序列中第一个素数是5，同样将序列中5的倍数划掉，主序列成了：
- 2 3 5 7 11 13 17 19 23
我们得到的素数有：2，3，5 。
因为23小于5的平方，跳出循环.