【学习笔记】P3383 【模板】线性筛素数

题目描述

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

输入格式

第一行包含两个正整数 n,q，分别表示查询的范围和查询的个数。
接下来 q 行每行一个正整数 k，表示查询第 k 小的素数

输出格式

输出 q行，每行一个正整数表示答案。

原题点击欧拉筛跳转

入门——暴力枚举

从2到n遍历每一个数，对每个数k进行检测，看它是否为素数，符合要求的数放入prime数组。检测方法：j从1到根号k，检测k mod j是否能被整除，如果能，则不是素数。

复杂度分析：O(n^(3/2))

外层循环n次，内层循环根号n次。

缺点分析：极慢

进阶——埃氏筛

核心思路：打标记。（标记范围内所有素数的倍数）
过程：遍历从2到n，对每个k，检测这个数是否被标记过，如果没有被标记过，则放入prime数组，并且将n以内所有k的倍数全部打上标记；如果这个数已经被标记过，则跳过它。

复杂度分析 O(nlogn)

易有1+1/2+1/3……+1/n=lnn
筛选n以内每个素数的倍数，标记的次数差不多是1/2n+1/3n+1/5n+……
=(1/2+1/3+1/5+……)*n约等于nlnn（注：推导过程并不严谨，若有不足之处欢迎指出）

缺点分析：一个数可能会被重复标记（如 12 会被2，3重复标记）

优化——欧拉筛

核心思路：在埃氏筛的基础上，保证每个合数只被标记一次（一个合数只能被它最大的因数标记。如12，只能由6来标记）

复杂度分析 O(n)

下面附上AC代码（如果过不了就用scanf和printf）

#include<iostream>
#include<vector>
using namespace std;

const int maxn=1e6+10;
vector<int>prime;
int vis[(int)1e8+10];

int main(){    
    std::ios::sync_with_stdio(0);//cin cout优化    
    std::cout.tie(0);
    std::cin.tie(0);
            
    int n,q;    
    cin>>n>>q;
    for(int i=2;i<=n;i++){//欧拉筛        
        if(!vis[i]) prime.push_back(i);
        for(int j=0;j<prime.size()&&prime[j]*i<=n;j++){            	 
            vis[i*prime[j]]=1;
	    if(i%prime[j]==0)break;        
	}   
    }
    
    for(int i=1;i<=q;i++){        
        int k;        
        cin>>k;        
        cout<<prime[k-1]<<endl;    
    }    
    return 0;
}