codevs 3223 素数密度(埃氏筛法)

题目描述 Description
给定区间[L, R](L <= R <= 2147483647,R-L <= 1000000),请计算区间中素数的个数。

输入描述 Input Description
两个数L和R

输出描述 Output Description
一行,区间中素数的个数

样例输入 Sample Input
2 11

样例输出 Sample Output
5

数据范围及提示 Data Size & Hint
详见试题

题解:这道题我们可以用埃氏筛法来做,我们先筛出[1,r]区间内的素数。然后用这些素数筛去[l,r]的非素数,剩下的就是素数个数。这道题的细节处理还是比较多的,所以打代码的时候要仔细。

代码如下

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
#define int long long
//这道题会爆int所以我在这里把int换成long long
//define是纯文本替换,不要问我为什么不手动改╭(╯^╰)╮ 
bool used[2000000],used2[2000000];
int pri[2000000];
#undef int
int main()//主函数不能改成long long啊! 
{
    #define int long long
    int l,r,p=0,sum=0;
    scanf("%lld%lld",&l,&r);
    for(int i=2;i<=sqrt(r);i++)//埃氏筛找[1,r]的素数 
        if(!used[i])
        {
            pri[++p]=i;
            for(int j=i*i;j<=sqrt(r);j+=i)
                used[j]=1;
        }
    for(int i=1;i<=p;i++)
    {
        int ll=l/pri[i];
        ll=ll*pri[i]+pri[i];//找大于l的能被筛去的最小数 
        if(l%pri[i]==0) ll-=pri[i];//判断一下l本身能被筛去的情况 
        for(int j=ll;j<=r;j+=pri[i])
            used2[j-l+1]=1;//把[l,r]所有的非素数筛去 
    }
    for(int i=l;i<=r;i++)
        if(!used2[i-l+1]) sum++;
    printf("%lld\n",sum);//输出素数个数 
    return 0;
}
/*
2146481947 2147481947
46609
*/
发布了81 篇原创文章 · 获赞 2 · 访问量 4万+

猜你喜欢

转载自blog.csdn.net/Loi_YZS/article/details/53039602
今日推荐