HDOJ 4548 美素数

本文介绍了一种利用素数筛算法结合数字之和判断的方法,通过两层筛选来快速找出特定范围内的素数。首先使用埃拉托斯特尼筛法获取基本的素数列表,接着对每个数字进行其各位数字之和是否为素数的判断,最终实现高效求解。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<iostream>
using namespace std;
const int L=1000005,inf=1<<30,maxn=1005;
int prime[L];
bool is[L];
void getPrime()
{
    fill(is,is+L,1);
    is[1]=0;
    int np=0;
    for(int i=2;i<L;i++)
        if(is[i])
        {
            prime[++np]=i;
            for(int j=2*i;j<L;j+=i) is[j]=0;
        }
}
void Get()
{
    for(int i=1;i<=1000000;i++)
    {
        if(!is[i])
            prime[i]=0;
        else
        {
            int t=i,s=0;
            while(t)
            {
                s+=t%10;
                t/=10;
            }
            if(is[s])
                prime[i]=1;
            else prime[i]=0;
        }
    }
}
int b[1000005]={0};
int main()
{
    getPrime();
    Get();
     for(int i=1;i<=1000010;i++)
    {
               if(prime[i])
                      b[i]=b[i-1]+1;
               else
                      b[i]=b[i-1];
    }

    int t,s,i,l,r;
   // scanf("%d",&t);
    int cnt=0;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&l,&r);

        printf("Case #%d: %d\n",++cnt,b[r]-b[l-1]);

    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值