今天的小明爱质数

最新推荐文章于 2024-08-01 13:27:05 发布

JackeyHook

最新推荐文章于 2024-08-01 13:27:05 发布

阅读量409

点赞数 2

分类专栏： C++ 文章标签：算法 c++

本文链接：https://blog.csdn.net/oTnLg/article/details/127592579

版权

C++ 专栏收录该内容

4 篇文章 0 订阅

订阅专栏

一名在校大学生，因为自己很菜，所以就更想提高自己，所以就在练习算法。

今天看到的一道题，具体内容是这样子的

一开始，看见题目就觉得很简单，写个判断质数的函数在来几个for循环就可以写完了，觉得so easy，写完之后，自测运行能通过，一提交就是超时了。然后实在是没有思路，去网上找相关的题目来看看。在c语言网站找到了一个类似的题目。大家练题的时候也可以去这个网站上面看看。Dotcpp编程（C语言网） - 实用的编程学练平台

找到一个类似的题目

#include <iostream>
#include <cmath>
#include <cstring>
using namespace std;
int main(){
    bool isprimeA[31623] = {};
    bool isprimeB[1000001] = {};
    for (int i = 2; i <= sqrt(31622); ++i){从区间[2:sqrt(sqrt(1000000000))]中选择质数
        if (!isprimeA[i]){
            for (int j = i*i; j <= 31622; j += i)//筛选区间[2:sqrt(1000000000)]里的合数
                isprimeA[j] = true;
        }
    }
    int L, R;
    while (cin >> L >> R){
        for (int i = 2; i <= sqrt(R); ++i){//从区间[2:sqrt(R)]中选择质数
            if (!isprimeA[i]){
                int j = L / i + (bool)(L%i);//选择大于等于L并且可以被i整除的第一个合数
                for (j = j > 2 ? j*i : 2 * i; j <= R; j += i)//筛选区间[L:R]里的合数
                    isprimeB[j - L] = true;
            }
        }
        int count = 0;
        for (int i = 0; i <= R-L; ++i){//统计区间[L:R]里面的质数
            if (!isprimeB[i])
                count++;
        }
        cout << count << endl;
        memset(isprimeB, 0, 1000001);//将区间[L:R]清零，以便下一次使用
    }
    return 0;
}

具体思路就是：

数据比较大，所以需要使用埃氏筛法来做这道题目**

埃氏筛法：将所有质数的倍数筛选出来，那么剩下的就是质数。

数据规模约定(2<=L<=R<=1000000000,R-L<=1000000)，所以我们只需要计算出sqrt(1000000000)以内的所有质数，然后再从区间[L:R]中将这些质数的倍数筛去，那么区间[L:R]中剩下的数的个数，就是区间[L:R]中质数的个数。

可能有的同志要问：为什么我们只要计算出sqrt(1000000000)以内的所有质数呢？（sqrt(1000000000) == 31622）

答：因为一个合数一定会分两个正整数相乘的形式(num=x*y),假设x和y都大于sqrt(num),这是不成立的，所以x和y中必然有一个是小于等于sqrt(num),所以我们只需要得到sqrt(num)以内的质数，就可以判断num是否为合数了。

我就理解了一会儿，就用埃氏筛法来做这道题，还是超时了，又来有了解到有比埃氏筛法在时间复杂度上更快的欧拉筛法，又去学习了一下，其实学了有点久。欧拉筛法和埃氏筛法的不同点在于埃氏筛法一个数可能要筛很多次，在时间上面就会浪费，对于欧拉筛法呢，一个数只用自己的最小质因子来筛，所以时间复杂度上面又会少一点。

后来用了欧拉筛法还是没有过，去问了一下自己的学长，这道题需要欧拉加上前缀和一起做才可以通过，有不得不去学习了一会儿前缀和，有兴趣了解的同学可以在B站上面找到相应的资源。

最后经过我长时间的努力，还是做出来了这道题。

#include<iostream>
#include<algorithm>
using namespace std;
const int N = 15000000;
int prime[N], cnt,arr[N]={0};
bool st[N];
int main()
{
    int sum[N];
    for (int i = 2; i < N; i++)
        {
            if (st[i]==false) {
                prime[cnt++] = i;
                arr[i] = 1;
            }
            for (int j = 0; prime[j] < N / i; j++) {
                st[prime[j] * i] = true;
                if (i % prime[j] == 0) break;
            }
        }
    for (int i = 0; i < N; i++) {
        sum[i] = sum[i - 1] + arr[i];
    }
    int m;
    scanf("%d", &m);
    while (m--) {
        int L, R;
        scanf("%d%d", &L, &R);
        printf("%d\n",sum[R]-sum[L-1]);
    }
    return 0;
}

这是牛客网上面的题，到了最后，答案也没有错了，还是给我报错，说是数组越界了。后来看了csdn上面有相关的问题，才知道是因为数组设置小了。所以报错了。

这道题的思路就是用欧拉筛发先把范围里面的质数筛出来，就不要用while循环去写一次次的欧拉筛，而是用前缀和得出想要的答案。