题目背景
本题已更新,从判断素数改为了查询第 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时被筛掉。