【C++算法模板】数论:欧拉筛,线性查找质数的算法

文章目录

1)传统找质数的方法(优化筛选次数)

c 复制代码
bool isPrime(int num) {
    for(int i=2;i<=sqrt(num)) {
        if(num%i==0)
            return false;
    }
    return true;
}
  • 如果要找从 [ 1 , 1 e 6 ] [1,1e6] [1,1e6] 中的所有质数,时间复杂度很高

2)欧拉筛

  • 算法思想:遍历到 2 2 2 的时候,筛掉范围内所有 2 2 2 的倍数(因为除了 1 1 1 和自身以外,一定能被 2 2 2 整除),到 3 3 3 的时候,筛掉所有 3 3 3 的倍数···
  • 注意:如果计算 [ l , r ] [l,r] [l,r] 之间出现的质数的个数?可以用前缀和的思想;当 n n n 过大时, i × i i×i i×i 容易出现数组越界的错误,即可能 R u n t i m e E r r o r RuntimeError RuntimeError ,此时要将线性筛中第二个 f o r for for 中的 j = i × i j=i×i j=i×i 改为 j = i + i j=i+i j=i+i
c 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=1e5;
int f[N]; // 下标为i时,记录1~i出现的所有质数的数量
bool vis[N]; // 是否已经访问过
int p[N]; // 用于存放素数
int idx,n; // idx是存放素数的遍历因子
void get_primes(int n)
{
	f[1]=0;
	for(int i=2;i<=n;i++)
	{
		// 如果vis[i]为false才需要遍历
		if(!vis[i]){
			f[i]=f[i-1]+1; // 计算前缀和
			p[++idx]=i; // 是素数,存起来
            // 如果出现RuntimeError,将j=i+i,这样就不会数组越界
			for(int j=i*i;j<=n;j+=i) // 将素数i的倍数全部标记为合数,则无需遍历
				vis[j]=true;
		} else 
			f[i]=f[i-1]; // 向下传递素数个数
	}
}

int main()
{
	// 如:输入1000即打印0~1000以内的素数
	cin>>n;  
	int l,r; // 左右区间
	cin>>l>>r;
	get_primes(n);
	cout<<"1~"<<n<<"之间的素数分别是:";
	for(int i=1;i<=idx;i++) {
		cout<<p[i];
		if(i!=idx)
			cout<<",";
	}
	cout<<endl;
	cout<<l<<"~"<<r<<"所出现的素数个数为:"<<f[r]-f[l-1]<<endl;
	return 0;
}
  • 最简单的模板
c 复制代码
const int N=5e4+5; // 注意这里没有开到2^9,只要比sqrt(2^9)大即可
int primes[N],cnt;
bool st[N];
int ans[N],len;

// 线性筛模板
void get_primes(int n) {
	for(int i=2;i<=n;i++) {
		if(!st[i]) primes[cnt++]=i;
		for(int j=0;primes[j]*i<=n;j++) {
			st[primes[j]*i]=true;
			if(i%primes[j]==0) break;
		}
	}
}

// 为什么特判x>=N的情况?因为为了节省内存或者没必要,本题中只要保证MAXN比sqrt(理论最大值)大即可
bool is_prime(int x) {
	if(x<N) return !st[x];
	for(int i=0;primes[i]<=x/primes[i];i++)
		if(x%primes[i]==0)
			return false;
	return true;
}
相关推荐
ZhiqianXia1 天前
C++ 常见代码异味(Code Smells)
c++
编码追梦人1 天前
基于 STM32 的智能语音唤醒与关键词识别系统设计 —— 从硬件集成到算法实现
stm32·算法·struts
循着风1 天前
二叉树的多种遍历方式
数据结构·算法
老猿讲编程1 天前
C++中的奇异递归模板模式CRTP
开发语言·c++
.格子衫.1 天前
022数据结构之树状数组——算法备赛
数据结构·算法·1024程序员节
黑科技Python1 天前
生活中的“小智慧”——认识算法
学习·算法·生活
Yupureki1 天前
从零开始的C++学习生活 16:C++11新特性全解析
c语言·数据结构·c++·学习·visual studio
紫荆鱼1 天前
设计模式-迭代器模式(Iterator)
c++·后端·设计模式·迭代器模式
sali-tec1 天前
C# 基于halcon的视觉工作流-章52-生成标定板
开发语言·图像处理·人工智能·算法·计算机视觉
IT古董1 天前
【第五章:计算机视觉-项目实战之推荐/广告系统】2.粗排算法-(4)粗排算法模型多目标算法(Multi Task Learning)及目标融合
人工智能·算法·1024程序员节