数学概念——H 最美素数

#include <cstdio>
#include <cstring>
using namespace std;
bool prime[1000010];
bool beauty[1000010];
int dp[1000010];
void findprime()
{
    memset(prime,true,sizeof(prime));
    memset(beauty,true,sizeof(beauty));
    prime[0]=prime[1]=false;
    beauty[0]=beauty[1]=false;
    for(int i=2;i<=1000;i++)
    {
        for(int j=2*i;j<=1000000;j+=i)
            {
                prime[j]=false;
                beauty[j]=false;
            }
    }
}
void beautyprime()
{    int sum,k;
     for(int i=2;i<=1000000;i++)
    {
        sum=0;k=i;
        if(prime[i]==true)
        {
            while(k!=0)
            {
                sum+=k%10;
                k/=10;
            }
            if(prime[sum]!=true)
                beauty[i]=false;
        }
    }
}
void DP()
{
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=1000000;i++)
        if(beauty[i]==true)
           dp[i]=dp[i-1]+1;
        else
           dp[i]=dp[i-1];
}
int main()
{
    int a,b,T,Case=0;
    findprime();
    beautyprime();
    DP();
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&a,&b);
        printf("Case #%d: %d\n",++Case,dp[b]-dp[a-1]);
    }
        return 0;
}