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

原创 2013年12月03日 21:00:12
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;  
    }  


相关文章推荐

数论-hdu_4548_美素数

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4548 题目描述: Problem Description   小明对数的研究比较热爱,一...

hdu4548 美素数(C语言)

Problem Description   小明对数的研究比较热爱,一谈到数,脑子里就涌现出好多数的问题,今天,小明想考考你对素数的认识。   问题是这样的:一个十进制数,如果是素数,而且它的...

hdu4548 美素数

美素数 Problem Description   小明对数的研究比较热爱,一谈到数,脑子里就涌现出好多数的问题,今天,小明想考考你对素数的认识。   问题是这样的:一个十进制数,如果是素...

HDU 4548 美素数

链接 : http://acm.hdu.edu.cn/showproblem.php?pid=4548 美素数:本身是素数,并且各位数字的和也是素数 打出素数表,再对各位数字求和,访问...

HDU 4548.美素数【打表】【预处理】【8月11】

美素数 Problem Description   小明对数的研究比较热爱,一谈到数,脑子里就涌现出好多数的问题,今天,小明想考考你对素数的认识。   问题是这样的:一个十进制数,如果是素数,而...

HDU 4548美素数(简单题 保存结果时需要注意不要超时)

美素数 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Subm...
  • opm777
  • opm777
  • 2013年08月19日 21:18
  • 907

HDU 4548 美素数(打表呀,有个小技巧)

美素数 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Submis...

HDU 4548 美素数【打表法】

美素数 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Subm...

HDU4548 美素数

问题链接:HDU4548 美素数。基础练习题,用C语言编写程序。 问题简述:参见上述链接。 问题分析:美素数定义为一个数是素数并且其各位数字之和也是素数。 本题的目的是计算某个区间的美素数的个数。...

美素数 HDU 4548 TLE! = =|| 打表是一种态度。。。

 美素数 Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total ...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 4548 -------美素数
举报原因:
原因补充:

(最多只允许输入30个字)