蓝桥杯 算法提高 找素数 【思维找素数 + 筛法】

算法提高 找素数

时间限制：1.0s 内存限制：256.0MB

问题描述

　　给定区间[L, R] ， 请计算区间中素数的个数。

输入格式

　　两个数L和R。

输出格式

　　一行，区间中素数的个数。

样例输入

2 11

样例输出

5

数据规模和约定

　　2 <= L <= R <= 2147483647 R-L <= 1000000

题意： 略

分析： 知道一般的素数筛法的应该都知道原理，这里其实就是模拟一下那里的筛法，简叙如下:

• 1.预处理出[0,2147483647]内的所有素数
• 2.用[0,r - l]来表示[l,r]，找到对应关系
• 3.利用筛法的思想进行筛数，比如第一个素数 2，[l,l+2]其中肯定有一个是2的倍数，找到起点，然后往后扫一遍筛掉即可

• 4.考虑3的做法是否合理，因为 r - l的范围是1e6，时限上过不去的，考虑优化，我们可以想到其实我们可以在O(1)的时间内找到起点，只有几种情况:

  • <1> 当素数p >= l， 我们只需将起点放在 p - l处即可
  • <2> 否则，考虑l % p是否为0，如若为0，则就在0点，否则在 p - l % p处
  • <3> 如果当前起点 + l为素数时，就要往后走p个单位

  ●5.O(n)统计可以

参考代码

#include <bits/stdc++.h>

using namespace std;

const int maxn = 1e6 + 7;

int m = (int) sqrt(2147483647+0.5) + 10;
bool isp[maxn];
int p[maxn/10];
int len;

void init() {
    isp[0] = isp[1] = true;
    for (int i = 2; i < m; i++) {
        if (!isp[i]) p[++len] = i;
        for (int j = 1; j <= len && p[j]*i < m; j++) {
            isp[i*p[j]] = true;
            if(i % p[j] == 0) break;
        }
    }
}

int main() {
    init();
    memset(isp,false,sizeof isp);
    int l,r;cin>>l>>r;
    int Len = r - l + 1;
    for (int i = 1; i <= len; i++) {
        int L;
        if(p[i] >= l) L = p[i] - l;
        else {
            if(l % p[i] == 0) L = 0;
            else L = p[i] - l%p[i];
        }
        if(l + L == p[i]) L += p[i];
        for (int j = L; j <= Len; j += p[i]) {
            isp[j] = true;
        }
    }
    int cnt = 0;
    for(int i = 0; i < Len; i++) {
        if (!isp[i]) cnt++;
    }
    if(l == 1) cnt--;
    cout<<cnt<<endl;
    return 0;
}

• 如有错误或遗漏，请私聊下UP，thx