【总结】素数筛方法详解


前言

    在讲素数筛方法之前,先了解相关知识点:

  1. 素数: 在大于1的自然数中,除了1和它本身以外不再有其他因数的数。
  2. 对于一个合数x,必有一个范围在2~√x 的因数。

     素数筛方法主要有三种:

  1. 素数的判定(素数筛)
  2. 普通线性筛(埃氏筛法 / 埃拉托斯特尼(Eratosthenes)筛法)
  3. 优化后的线性筛(欧拉筛法 / 欧拉函数(Euler)筛)

一、素数的判定

时间复杂度:O(√n)

思路: 暴力枚举,很多情况下都会TLE。枚举2~~√n所有数,用n去试着除以,若有能整除的n为合数,若都不能整除,n就是质数了。

int isPrime(int x){
	//0和1既不是质数也不是合数 
	if(x<2)
		return 0;
	for(int i=2;i*i<=x;i++){
		if(!(x%i))
			return 0;
	}
	return 1;
}

二、埃氏筛法

时间复杂度:O(n loglogn )

思路: 0表示为不是素数,1表示是素数,先全部初始化为1,即先全部标记为素数。0和1均不是素数,然后从2开始,把2的倍数均标记为0,即不是素数,一直到大于n,然后从下一个素数3开始,进行同样的处理,到最后,数组中为1的即为素数。

代码1:

const int N=1e6;
bool vis[N];
void Erat_Prime(int n){
	memset(vis,1,sizeof(vis));
	vis[0]=vis[1]=0;	//0和1均不是素数
	for(int i=2;i*i<=n;i++){
		if(vis[i]){
			for(int j=i*i;j<=n;j+=i)
				vis[j]=0;
		}
	} 
}

代码2:

bool vis[maxn];
int prime[maxn],x;
void Erat_prime(int n) //埃氏筛
{
    for(int i=2;i<=n;i++)
    {
        if(!vis[i]) prime[x++]=i;
        for(int j=2;j*i<=n;j++)
        {
            vis[i*j]=true;
        }
    }
}

三、欧拉筛法

时间复杂度:O(n)

在用埃式筛法的同时,同一个数字也许会被筛选多次,比如6先被2筛选一次,再被3筛选一次,这样就浪费了很多不必要的时间。欧拉筛法就是在埃氏筛法的基础上,让每个合数只被它的最小质因子筛选一次,以达到不重复的目的。

用一条语句 if(i%prime[j]==0) break; 避免了重复筛选的发生。

const int N=1e5+10;
int vis[N];	//0表示素数,1表示非素数
int prime[N];	//只在这个函数有作用
void Euler_prime()  //欧拉筛法
{
    for(int i=2;i<=N;i++)
    {
        if(!vis[i]) prime[x++]=i;
        for(int j=0;j<x;j++)
        {
            if(i*prime[j]>N) break;
            vis[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}

四、例题

原题地址:【洛谷】P3383 【模板】线性筛素数

题目描述

如题,给定一个范围N,你需要处理M个某数字是否为质数的询问(每个数字均在范围1-N内)

输入输出格式

输入格式:

第一行包含两个正整数N、M,分别表示查询的范围和查询的个数。

接下来M行每行包含一个不小于1且不大于N的整数,即询问该数是否为质数。

输出格式:

输出包含M行,每行为Yes或No,即依次为每一个询问的结果。

输入输出样例

输入样例#1:

100 5
2
3
4
91
97

输出样例#1:

Yes
Yes
No
No
Yes

说明

时空限制:500ms 128M

数据规模:

对于30%的数据:N<=10000,M<=10000

对于100%的数据:N<=10000000,M<=100000

样例说明:

N=100,说明接下来的询问数均不大于100且不小于1。

所以2、3、97为质数,4、91非质数。

故依次输出Yes、Yes、No、No、Yes。

代码如下:

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std; 
const int N=1e7;
//n表示范围,m表示个数,数组a[i]表示输入的数、数组b[i],0表示不是为素数,1表示是素数 
int n,m,x,a[N],b[N];	
int main() 
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)	//输入m个数 
		scanf("%d",&a[i]);
	memset(b,1,sizeof(b));	//初始化,把所有的数标记为是素数 
	b[0]=b[1]=0;	//0和1不是素数 
	for(int i=2;i<=n;i++)	//从2开始遍历 
	{
		if(b[i])	//如果标记为素数,即未改变初始化的值 
		{
			for(int j=i+i;j<=n;j+=i)	//把所有该数的倍数标记为0,即不是素数 
				b[j]=0;
		}
	}
	for(int i=1;i<=m;i++)
	{
		if(b[a[i]])	//如果输入的数是素数,输出Yes 
			printf("Yes\n");
		else	//否则输出No 
			printf("No\n");
	}
	return 0;
}
  • 3
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
ExecutorService是Java中的一个接口,它是Executor的子类接口,提供了一些额外的方法来管理任务的生命周期和获取异步任务的执行结果。ExecutorService,我们可以更方便地线程池来执行任务。 ExecutorService的一些常用方法包括: 1. submit方法:用于提交一个Runnable或Callable的任务,并返回一个Future对象。对于Runnable任务,Future的get方法将返回null;对于Callable任务,Future的get方法将返回任务的执行结果。 2. invokeAny方法:用于提交一组Callable任务,并返回其中一个任务的执行结果。它会阻塞直到有一个任务完成并返回结果,忽略其他任务的执行结果。 3. invokeAll方法:用于提交一组Callable任务,并返回所有任务的执行结果。它会阻塞直到所有任务都完成并返回结果。 4. awaitTermination方法:用于等待所有任务在指定的时间内完成。如果所有任务在超时时间内完成,则返回true;否则返回false。 5. isShutdown方法:用于判断ExecutorService是否已经关闭。 6. shutdown方法:用于平滑地关闭ExecutorService。它将停止接收新的任务,并等待已经提交的任务完成。 7. shutdownNow方法:用于立即关闭ExecutorService。它将尝试取消所有正在执行的任务,并返回等待执行的任务列表。 8. isTerminated方法:用于判断所有任务是否已经完成。 这些方法可以帮助我们更灵活地管理和控制线程池中的任务的执行。比如,通过submit方法可以提交任务并获取任务的执行结果;通过invokeAny方法可以获取一组任务中最先完成的任务的执行结果;通过shutdown方法可以平滑地关闭线程池等。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [Java高并发编程中Executor、ExecutorService的使用及详细介绍-刘宇](https://blog.csdn.net/liuyu973971883/article/details/107947496)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* [java ExecutorService使用方法详解](https://download.csdn.net/download/weixin_38689477/12790401)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值