HDU 4548 美素数 （素数打表+前缀和）

[align=center][/align]

美素数

[align=center]Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 7548    Accepted Submission(s): 2661
[/align]

[align=left]Problem Description[/align]
　　小明对数的研究比较热爱，一谈到数，脑子里就涌现出好多数的问题，今天，小明想考考你对素数的认识。

　　问题是这样的：一个十进制数，如果是素数，而且它的各位数字和也是素数，则称之为“美素数”，如29，本身是素数，而且2+9 = 11也是素数，所以它是美素数。

　　给定一个区间，你能计算出这个区间内有多少个美素数吗？
 

[align=left]Input[/align]
第一行输入一个正整数T，表示总共有T组数据(T <= 10000)。

接下来共T行，每行输入两个整数L，R(1<= L <= R <= 1000000)，表示区间的左值和右值。
 

[align=left]Output[/align]
对于每组数据，先输出Case数，然后输出区间内美素数的个数（包括端点值L,R）。

每组数据占一行，具体输出格式参见样例。
 

[align=left]Sample Input[/align]

3
1 100
2 2
3 19

 

[align=left]Sample Output[/align]

Case #1: 14
Case #2: 1
Case #3: 4

 

[align=left]Source[/align]
2013金山西山居创意游戏程序挑战赛——初赛（2）

题目链接：HDU 4548 美素数
[align=center]
[/align]
[align=center]素数打表+前缀和
[/align]


代码：

#include<cstdio>
#include<cstring>
using namespace std;
const int M=1e6+1;
bool num[M];//素数打表
int  ans[M];//前缀和，ans[i]表示从1到i 共有多少个美素数

void init_prime()
{
memset(num,1,sizeof(num));
num[1]=0;
for(int i=2;i<=M;i++)
{
if(num[i])
{
for(int j=2*i;j<=M;j+=i)
num[j]=0;
}
}
}
void init_ans()
{
memset(ans,0,sizeof(ans));
for(int i=2;i<=M;i++)
{
if(num[i])
{
int s=0;
int temp=i;
//s为i各位数字和
while(temp>0)
{
s+=temp%10;
temp/=10;
}

if(num[s])
{
ans[i]=ans[i-1]+1;
}
else
{
ans[i]=ans[i-1];
}
}
else ans[i]=ans[i-1];
}

}

int main()
{
int n,Case=0;
init_prime();
init_ans();
scanf("%d",&n);
while(n--)
{
int L,R;
scanf("%d%d",&L,&R);
//ans[R]-ans[L-1]计算区间 L-R内美素数的个数
printf("Case #%d: %d\n",++Case,ans[R]-ans[L-1]);
}
return 0;
}