P3383 【模板】线性筛素数

题目背景
本题已更新,从判断素数改为了查询第 k 小的素数
提示:如果你使用 cin 来读入,建议使用 std::ios::sync_with_stdio(0) 来加速。

题目描述
如题,给定一个范围 n,有 q 个询问,每次输出第 k 小的素数。

输入格式
第一行包含两个正整数 n,q,分别表示查询的范围和查询的个数。

接下来 q 行每行一个正整数 k,表示查询第 k 小的素数。

输出格式
输出 q 行,每行一个正整数表示答案。

输入输出样例
输入 #1

100 5
1
2
3
4
5

输出 #1

2
3
5
7
11

说明/提示
【数据范围】

对于 100% 的数据,n = 10^8, 1≤q≤10^6,保证查询的素数不大于 n。

Data by NaCly_Fish.

//#include<bits/stdc++.h>
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int a[100000010];
int b[50000010],sum = 0;

void pri(int n){
    memset(a,1,sizeof(a));
    a[1] = 0;
    for(int i = 2;i <= n; i++) {
        if(a[i]) {
            b[++sum] = i;
        }
        for(int j = 1;j <= sum && i*b[j] <= n; j++) {
            a[i*b[j]] = 0;
            
            if(i%b[j] == 0){
                break;
            }
        }
    }
}

int main()
{
	int n, q;
	scanf("%d %d", &n, &q);
	pri(n);
	while (q--)
	{
		int k;
		scanf("%d", &k);
		printf("%d\n", b[k]);
	}
	return 0;
}

**

欧拉筛

**
算法原理
素数筛法的原理:把不是素数的筛掉,剩下的就都是素数。

欧拉线性筛的关键在于:每个合数只被它最大的非自身的因数筛掉。

例如:

12不会被4筛掉,而是被6筛掉;

45不会被9筛掉,而是被15筛掉。
模拟程序运行过程
我们以 12不被4筛掉而被6筛掉 为例:

当 i 循环执行到 i=4 时:

素数表prime中已有两个素数:2和3。

此时在check数组中已经被标记为1的数(合数)是:4、6、9。

因为i=4已经被标记为1,所以我们不将4添加进素数表prime。

此时执行 j 循环:

因为素数表prime中有2和3,所以预计将要被筛掉的数是4×2=8 和 4×3=12。

当4×2=8被筛掉以后,经过 imodprime[j]=0 判断可以知道4mod2=0, 此时应结束 j 循环,不再筛 4×3=12。

关键步骤解释
我们都知道12的最大因数(非自身,下同)是6,因此12应该被6筛掉,那么上述判断的原理是什么呢?

因为当出现 imodprime[j]=0 的情况时,意味着prime[j] 是 i 的一个因数(2是4的一个因数)。

既然4是12的因数,2又是4的因数,所以2也是12的因数。

由此可以知道,4不是12的最大因数,12的最大因数会在后面出现,即2所对应的6,12将会在i=6时被筛掉。

同理可以判断:9是45的因数,3又是9的因数,所以3也是45的因数,9不是45的最大因数,45的最大因数将会在后面出现,即3所对应的15,45将会在i=15时被筛掉。

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值