素数筛法

一、素数

素数筛选，判断<MAXN的数是不是素数，模板by kuangbin

//判断小于MAXN的数是否是素数
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
const int MAXN=1e5+5;
bool notprime[MAXN];            //表，false是素数，true不是素数
void isprimer(){
    memset(notprime,false,sizeof(notprime));
    notprime[0]=notprime[1]=true;
    for(int i=2;i<MAXN;i++){
        if(!notprime[i]){
            if(i>MAXN/i)    continue;       //防止后面i*i溢出（或者i,j用Long long）
            //直接从i*i开始就可以，小于i倍的已经筛过了，注意是j+=i;
            for(int j=i*i;j<MAXN;j+=i){
                notprime[j]=true;
            }
        }
    }
}
int main(){
    int n;
    isprimer();
    cin>>n;
    for(int i=0;i<n;i++){
        if(notprime[i]==false)
            cout<<i<<" ";
    }
    return 0;
}

二、PAT乙级：1013 数素数 

/*
令 Pi表示第 i 个素数。现任给两个正整数 M≤N≤10^4,请输出 PM到 PN的所有素数。

输入格式：
输入在一行中给出 M 和 N，其间以空格分隔。

输出格式：
输出从 PM到 PN的所有素数，每 10 个数字占 1 行，其间以空格分隔，但行末不得有多余空格。

输入样例：
5 27

输出样例：
11 13 17 19 23 29 31 37 41 43
47 53 59 61 67 71 73 79 83 89
97 101 103
*/
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
using namespace std;
const int MAXN=1e7+10;
int primer[MAXN],num=0;
bool notPrimer[MAXN];
void isPrimer(){
    memset(notPrimer,false,sizeof(notPrimer));
    notPrimer[0]=notPrimer[1]=true;
    for(int i=2;i<MAXN;i++){
        if(!notPrimer[i]){
            primer[num++]=i;
            if(i>MAXN/i)    continue;
            for(int j=i*i;j<MAXN;j+=i)
                notPrimer[j]=true;
        }
    }
}
int main(){
    int n,m,ans=0;
    isPrimer();
    cin>>m>>n;
    for(int i=m;i<=n;i++){
        cout<<primer[i-1];
        ans++;
        if(ans%10!=0 &&i<n)
            cout<<" ";  
        else{
            cout<<endl;
        }                    
    }
    return 0;
}