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[],b[];
 int N=;
　　　　//打素数表
 void set_prime()
 {
     prime[]=;      //这里是标记多余的一个数为1
     prime[]=;      //因为1不是素数所以要标记为1
     for(int i=; 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]=;
         }
     }
 }
 
　　　//打美素数的个数表..
 void set_ansbiao()
 {
     int time=,k;
     for(int i=; i<N; i++)
     {
         if(!prime[i])   //是素数
         {
             int sum=;
             k=i;
             while(k)
             {
                 sum+=k%;
                 k=k/;
             }
             if(!prime[sum])   //它的和是素数...
                 time++;       // 个数加1
         }
         b[i]=time;         //个数用数组储存下来
     }
 }

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