HDU 4548 美素数

Description

　　小明对数的研究比较热爱，一谈到数，脑子里就涌现出好多数的问题，今天，小明想考考你对素数的认识。  　　问题是这样的：一个十进制数，如果是素数，而且它的各位数字和也是素数，则称之为“美素数”，如29，本身是素数，而且2+9 = 11也是素数，所以它是美素数。  　　给定一个区间，你能计算出这个区间内有多少个美素数吗？

 

Input

第一行输入一个正整数T，表示总共有T组数据(T <= 10000)。  接下来共T行，每行输入两个整数L，R(1<= L <= R <= 1000000)，表示区间的左值和右值。

 

Output

对于每组数据，先输出Case数，然后输出区间内美素数的个数（包括端点值L,R）。  每组数据占一行，具体输出格式参见样例。

 

Sample Input

3

1 100

2 2

3 19

 

Sample Output

Case #1: 14

Case #2: 1

Case #3: 4

 

 

解题思路：为了不超时，就要打素数表了。这里还要打一个美素数的个数表.. 

 

 

 

代码如下：（详情看注释）

 

 

#include <stdio.h>
#include <math.h>
#include <string.h>
int prime[1000005],b[1000005];
int N=1000000;
　　　　//打素数表
void set_prime()
{
    prime[0]=1;      //这里是标记多余的一个数为1
    prime[1]=1;      //因为1不是素数所以要标记为1
    for(int i=2; i<=sqrt(N); i++)
    {
        if(!prime[i])
        {
            for(int j=i*i; j<N; j+=i)   //这里表示i的倍数不是素数。（用j+j更加直观一些，用j*j更省时，但是为什么可以用j*j，我也一时没反应过来，毕竟是小白.. (⊙﹏⊙)b）
                prime[j]=1;
        }
    }
}

　　　//打美素数的个数表..
void set_ansbiao()
{
    int time=0,k;
    for(int i=0; i<N; i++)
    {
        if(!prime[i])   //是素数
        {
            int sum=0;
            k=i;
            while(k)
            {
                sum+=k%10;
                k=k/10;
            }
            if(!prime[sum])   //它的和是素数...
                time++;       // 个数加1
        }
        b[i]=time;         //个数用数组储存下来
    }
}

int main()
{
    set_prime();
    set_ansbiao();
    int T,cas=0;
    scanf("%d",&T);
    while(T--)
    {
        int L,R;
        scanf("%d%d",&L,&R);
        cas++;
        printf("Case #%d: %d
",cas,b[R]-b[L-1]);    // 输出个数....
    }
    return 0;
}