hdu 4548 初始化+二分 *

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

链接：点我

先找出所有美素数，然后二分找位置，我这里直接用stl了，位置相减即可，注意区间端点的判断

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 #include<map>
 using namespace std;
 #define MOD 1000000007
 #define pb(a) push_back(a)
 const int INF=0x3f3f3f3f;
 const double eps=1e-;
 typedef long long ll;
 #define cl(a) memset(a,0,sizeof(a))
 #define ts printf("*****\n");
 const int MAXN=;
 int a[];
 int n,m,tt,cnt;
 bool notprime[MAXN];//值为false表示素数，值为true表示非素数
 void init()
 {
     memset(notprime,false,sizeof(notprime));
     notprime[]=notprime[]=true;
     for(int i=;i<MAXN;i++)
     if(!notprime[i])
     {
         if(i>MAXN/i)continue;//防止后面i*i溢出(或者i,j用long long)
         //直接从i*i开始就可以，小于i倍的已经筛选过了,注意是j+=i
         for(int j=i*i;j<MAXN;j+=i)
         notprime[j]=true;
     }
 }
 bool fun(int x)
 {
     int w=x;
     int sum=;
     while(w)
     {
         sum+=(w%);
         w/=;
     }
     if(!notprime[x]&&!notprime[sum])    return ;
     else return ;
 }
 int main()
 {
     int i,j,k;
     #ifndef ONLINE_JUDGE
     freopen("1.in","r",stdin);
     #endif
     a[]=;
     int tot=;
     init();
     for(i=;i<=;i++)
     {
         if(fun(i))
         {
             a[tot++]=i;
         }
     }
     int tt;
     scanf("%d",&tt);
     int ca=;
     while(tt--)
     {
         ca++;
         int aa,bb;
         scanf("%d%d",&aa,&bb);
         if(aa==bb)
         {
             if(fun(aa))
             {
                 printf("Case #%d: 1\n",ca);
             }
             else
             {
                 printf("Case #%d: 0\n",ca);
             }
             continue;
         }
         int pos1,pos2;
         int f1=fun(aa);
         int f2=fun(bb);
         pos1=lower_bound(a,a+,bb)-a;
         pos2=lower_bound(a,a+,aa)-a;
         //printf("%d %d\n",pos1,pos2);
         int q;
         if((f1==&&f2==)||(f1==&&f2==))
         {
             q=;
         }
         else q=;
         printf("Case #%d: %d\n",ca,pos1-pos2+q);
     }
 }