【数学相关】素数问题

素数是指只能被其自身和1整除的正整数。

素数判断：0,1,和负数都是非素数。从2到sqrt(n)依次判断能否被n整除，若存在某个数能够整除该数则该数不是素数；若这些数都不能整除它，则该数为素数。

【素数判定】

【题目描述】：给定一个数n，要求判断其是否为素数。

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

bool Judge(int x){
    if(x<2){
        return false;
    }
    int bound=sqrt(x);
    for(int i=2;i<=bound;i++){//从2到sqrt(n)依次判断能否被n整除
        if(x%i==0){//存在整数可以被n整除，n就不是质数
            return false;
        }
    }
    return true;
}

int main(){
    int n;
    while(cin>>n){
        if(Judge(n)){
            cout<<"YES"<<endl;
        }else{
            cout<<"NO"<<endl;
        }
    }
    return 0;
}

素数筛法

如何求一个范围内的所有素数呢？需要用到【素数筛法】。

【素数筛法】：遍历目标范围内的所有整数，若当前整数没有因为它是某个小于其的素数的倍数而被标记为非素数，则判定其为素数，并标记它的所有的倍数为非素数。然后继续遍历下一个数，直到遍历完目标范围内的所有整数。此时，所有未被标记成非素数的数即为要求的素数。

【素数】

【题目描述】:输入一个整数n(2<=n<=10000)。要求输出所有从1到这个整数之间个位为1的素数，若没有则输出-1。

【分析】在处理输入的数据之前，可以先利用素数筛法求出2到10000内的所有素数。当得到输入n之后，依次判断在1~n这个区间内的素数是否符合题目条件，若符合题目条件则输出该素数，否则继续下一个素数，直到将该区间内所有素数都判断完为止。

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

const int MAXN=1e5;
vector<int> prime;//保存质数的值
bool isPrime[MAXN];//判断是否是质数
//素数筛法
void Initial(){
    //初始认为都是质数
    for(int i=0;i<MAXN;i++){
        isPrime[i]=true;
    }
    //fill(isPrime,isPrime+MAXN,true);//填充函数
    isPrime[0]=false;
    isPrime[1]= false;
    for(int i=2;i<MAXN;i++){
        if(!isPrime[i]){//非质数就跳过
            continue;
        }
        prime.push_back(i);//把质数存入prime数组
        for(int j=i*i;j<MAXN;j+=i){//从i的平方开始标记尽可能避免重复工作
            isPrime[j]= false;//质数的整数倍都不是质数
        }
    }
    return;
}

int main(){
    Initial();
    int n;
    while(cin>>n){
        bool isOutput=false;//判断是否有输出
        for(int i=0;i<prime.size()&&prime[i]<n;i++){
            if(prime[i]%10==1){//判断该质数的个位是否为1
                isOutput=true;
                cout<<prime[i]<<" ";
            }
        }
        if(!isOutput){
            cout<<"-1";
        }
        cout<<endl;
    }
    return 0;
}

【输出第k个质数】

【分析】：利用素数筛法找出2-MAXN的所有质数并存入prime数组，输出第k个质数即输出prime[k-1]的值。

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

const int MAXN=1e5;
vector<int> prime;//保存质数的值
bool isPrime[MAXN];//判断是否是质数
//素数筛法
void Initial(){
    //初始认为都是质数
    for(int i=0;i<MAXN;i++){
        isPrime[i]=true;
    }
    //fill(isPrime,isPrime+MAXN,true);//填充函数
    isPrime[0]=false;
    isPrime[1]= false;
    for(int i=2;i<MAXN;i++){
        if(!isPrime[i]){//非质数就跳过
            continue;
        }
        prime.push_back(i);
        for(int j=i*i;j<MAXN;j+=i){
            isPrime[j]= false;//质数的整数倍都不是质数
        }
    }
    return;
}

int main(){
    Initial();
    int k;
    while(cin>>k){
        cout<<prime[k-1]<<endl;
    }
    return 0;
}