关闭

HDU 4548 -------美素数

177人阅读 评论(0) 收藏 举报
分类:
http://acm.hdu.edu.cn/showproblem.php?pid=4548
美素数是指所有位累加和仍是一个素数的素数。
所以,可以先打表求出1——1000000之间所有的素数。1——1000000之间的所有素数个数为78498个。
然后在这么多个素数中进行判断,将是美素数的数再存到一个数组里面,1——1000000之间的美素数有30123个。
以后每次区间询问的时候,就可以在该数组进行遍历查找。


代码:
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <algorithm>
#include <math.h>

using namespace std;

#define N 1000001

bool visit[N];
int prime[N];

int perfect_prime[N];

int num = 0;
void make_prime()
{
    memset(visit,true,sizeof(visit));

    for(int i=2; i<=N; i++)
    {
        if(visit[i])
        {
            num++;
            prime[num] = i;
        }
        for(int j=1; (j<=num)&&(i*prime[j] <= N); j++)
        {
            visit[i*prime[j]] = false;
            if(!(i%prime[j])) break;
        }
    }
}

bool judge_prime(int n)
{
    if(n == 2) return true;
    int k = (int)sqrt(double(n));
    for(int i=1; i<=k; i++)
    {
        if(!(n%prime[i])) return false;
    }
    return true;
}

bool bit_prime(int n)
{
    int ans = 0;
    while(n)
    {
        ans+= n%10;
        n/=10;
    }
    return judge_prime(ans);
}

int init()
{
    make_prime();
    int cnt=0;
    for(int i=1; i<=num; i++)
    {
        if(bit_prime(prime[i]))
        {
             perfect_prime[cnt++] = prime[i];
        }
    }
    return cnt;
}
int main()
{
    int t;
    int n,m;
    int _case = 0;
    int ret = init();
    //printf("%d  -- %d\n",num,ret);
    scanf("%d%*c",&t);
    while(t--)
    {
        _case++;
        int ans = 0;
        scanf("%d%d",&n,&m);
        for(int i=0; i<ret; i++)
        if(perfect_prime[i]>=n && perfect_prime[i]<=m) ans++;
        printf("Case #%d: %d\n",_case,ans);
    }
    return 0;
}


0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:376次
    • 积分:22
    • 等级:
    • 排名:千里之外
    • 原创:2篇
    • 转载:0篇
    • 译文:0篇
    • 评论:0条
    文章存档