HDU4548+素数

简单题。

 1 /*
 2 
 3 */
 4 #include<stdio.h>
 5 #include<string.h>
 6 #include<stdlib.h>
 7 #include<algorithm>
 8 #include<iostream>
 9 #include<queue>
10 #include<stack>
11 #include<math.h>
12 #include<map>
13 using namespace std;
14 const int maxn = 1000005;
15 int IsPrime[ maxn ],IsSpecialPrime[ maxn ];
16 int GetSum( int x ){
17     int ans = 0;
18     while( x ){
19         ans += x%10;
20         x/=10;
21     }
22     return ans;
23 }
24 
25 void init_prime(){
26     for( int i=1;i<maxn;i+=2 )
27         IsPrime[i] = 1;
28     for( int i=0;i<maxn;i+=2 )
29         IsPrime[i] = 0;
30     IsPrime[1] = 0;
31     IsPrime[2] = 1;
32     for( int i=3;i<maxn;i++ ){
33         if( IsPrime[i]==1 ){
34             int t,delta;
35             delta = i*2;
36             t = delta+i;
37             while( t<maxn ){
38                 IsPrime[t] = 0;
39                 t += delta;
40             }
41         }
42     }
43     //for( int i=maxn-1;;i-- ){
44         //if( IsPrime[i]==1 ){
45             //printf("max=%d\n",i);
46             //break;
47         //}
48     //}
49     //memset( IsSpecialPrime,0,sizeof(IsSpecialPrime) );
50     for( int i=0;i<maxn;i++ ){
51         if( IsPrime[i]==1&&IsPrime[GetSum(i)]==1 ){
52             IsSpecialPrime[i] = 1;
53         }
54         else{
55             IsSpecialPrime[i] = 0;
56         }
57     }
58     for( int i=1;i<maxn;i++ )
59         IsSpecialPrime[ i ] += IsSpecialPrime[ i-1 ];
60 }
61 int main(){
62     init_prime();
63     int ca;
64     scanf("%d",&ca);
65     int T = 1;
66     while( ca-- ){
67         printf("Case #%d: ",T++);
68         int L,R;
69         scanf("%d%d",&L,&R);
70         printf("%d\n",IsSpecialPrime[R]-IsSpecialPrime[max(0,L-1)]);
71     }
72     return 0;
73 }
View Code

 

posted @ 2013-08-06 10:51  xxx0624  阅读(192)  评论(0)  编辑  收藏  举报