hdu 4548 素数打表

链接：http://acm.hdu.edu.cn/showproblem.php?pid=4548
方法：筛法打表，再将是素数的中判断是不是美素数。比较裸。但是就是一开始这样做会超时。重新用一个数组来保存到当前数字时已经出现了几个美素数，最后将区间两端所对应的相互减一下就可以了。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#define M 1000001
bool is_prime[M]; //素数表
bool is[M]; //美素数表
int dp[M];
bool judge(char ss[])
{
    int sum = 0;
    for(int i = 0;i < strlen(ss);i++)
    {
        sum += ss[i]-'0';
    }
    return is_prime[sum];
}
int main()
{
    fill(is_prime,is_prime+M,true);
    for(int i = 2;i <= M;i++)  //筛法打表
    {
        if(is_prime[i])
        {
            for(int j = 2*i;j <= M;j += i)
                is_prime[j] = false;
        }
    }
    is_prime[1] = false;
    for(int i = 1;i <= M;i++)
    {
        char s[100];
        sprintf(s,"%d",i);
        if(is_prime[i] && judge(s))
            is[i] = true;
    }
    fill(dp,dp+M,0);
    for(int i = 1;i <= M;i++)  //运用一个dp[i]记录下到i这个数为止有几个素数(包括i)
    {
        if(is[i])
            dp[i] = dp[i-1]+1;
        else dp[i] = dp[i-1];
    }
    int t;
    while(scanf("%d",&t)==1)
    {
        for(int i = 1;i <= t;i++)
        {
            int a,b;
            scanf("%d %d",&a,&b);
            int ans;
            if(!is[a])
                ans = dp[b]-dp[a];
            else ans = dp[b]-dp[a]+1;
            printf("Case #%d: %d\n",i,ans); //也可以直接输出dp[b]-dp[a-1] 因为如果a不符合就是对的，符合的话刚好不会减掉a这个
        }
    }
    return 0;
}