选美素数~~~~~~~(⊙o⊙)

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

思路：

打表。质数的倍数一定是合数，利用这点我们就可以把合数都标记出来

 1 #include<iostream>
 2 #include<cstring>
 3 using namespace std;
 4 const int M = 1000010;
 5 int isprime[M];
 6 int isbuprime_sum[M];
 7 void prime()
 8 {
 9     memset(isprime, 0, sizeof(isprime));
10     memset(isbuprime_sum, 0, sizeof(isbuprime_sum));
11     isprime[0] = 1;
12     isprime[1] = 1;                     //初始化
13     for (int i = 2; i < M; i++)
14     {
15         if (isprime[i])                        //判断，是合数就直接跳过
16             continue;
17         for (int j = i+i; j < M; j += i)       //先把素数的倍数（就是合数）都标记好
18         {
19             isprime[j] = 1;
20         }
21     }
22 }
23
24 int sumsum(int k)
25 {
26     int sum = 0;
27     while (k)
28     {
29         sum += k % 10;                       //各个位数的和，以便后面判断
30         k /= 10;
31     }
32     return sum;
33 }
34
35 int main()
36 {
37     prime();
38     isbuprime_sum[0] = isbuprime_sum[1] = 0;
39     for (int i = 2; i < M; i++)
40     {
41         if (!isprime[i] && !isprime[sumsum(i)])          //之前素数被标记为0了
42            isbuprime_sum[i] = isbuprime_sum[i - 1] + 1; //是美素数就计一次数
43         else
44              isbuprime_sum[i] = isbuprime_sum[i - 1];   //不是美素数，当前计数不变
45     }
46     int T,ans=1;
47     cin >> T;
48     while (T--)
49     {
50         int R, L;
51         cin >> L >> R;
52         cout << "Case #" << ans++ << ": " << isbuprime_sum[R] - isbuprime_sum[L - 1] << endl; //输出L-1因为包括左右端点
53
54     }
55
56     return 0;
57 }

心得：

再get一招，可以这样找素数，简直棒棒的，哈哈，集训太快了，转眼就快结束了。还真有点舍不得就回家了呢，说多了，学到知识更重要，加油加油啦~~~~~~~?(^∇^*)