LightOJ 1197（大数范围筛区间内素数个数）

最新推荐文章于 2020-03-01 08:52:08 发布

tick_tokc97

最新推荐文章于 2020-03-01 08:52:08 发布

阅读量650

点赞数

分类专栏： ACM-数论文章标签：数论 ACM 素数

本文链接：https://blog.csdn.net/tick_tock97/article/details/71104203

版权

ACM-数论专栏收录该内容

23 篇文章 0 订阅

订阅专栏

对于大范围内的素数个数，可以通过模拟筛法筛区间内素数个数，先打表，再根据表进行筛素数操作

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cassert>
#include <iostream>
#include <string>
#include <vector>
using namespace std;
typedef long long ll;

const int MAXN = 1000000+10;

bool is_prime[MAXN];//true的为合数，false的为素数
int prime[MAXN/10];//存素数
bool num[MAXN];
int tol;
void find_prime()
{
    memset(is_prime,false,sizeof(is_prime));
    tol=0;
    is_prime[1]=1;
    for(int i=2; i<MAXN; i++)
    {
        if(!is_prime[i])
        {
            prime[tol++]=i;
            for(int j=i*2; j<MAXN; j+=i)
                is_prime[j]=1;
        }
    }
}


int main()
{
//    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
    int T;
    cin>>T;
    find_prime();
    for(int cas=1; cas<=T; cas++)
    {
        int a,b;
        scanf("%d %d",&a,&b);
        int n=b-a;
        memset(num,false,sizeof(num));
        for(int i=0;i<tol&&prime[i]*prime[i]<=b;i++)
        {
            int j=0;
            if(a%prime[i]!=0)
                j=prime[i]-a%prime[i];//找到第一个需要筛掉的数(j+a)%prime[i]==0
            if(a<=prime[i])
                j+=prime[i];//(j+a)/prime[i]==1,则（j+a）是素数，向下推一个
            for(;j<=n;j+=prime[i])
                num[j]=true;
        }
        int ans=0;
        for(int j=0;j<=n;j++)
        {
            if(!num[j]) ans++;
        }
        if(a==1) ans--;
        printf("Case %d: %d\n",cas,ans);
    }
    return 0;
}