hdu4548 美素数

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4548

#include <stdio.h>  
#include <string.h>  
#define MAXN 1000000 
#define SIZE 1000002  

int prime[SIZE],valid[SIZE],ans[SIZE];  

void getPrime()//素数筛选O(N)  
{  
	int i,j,tot=0;  
	memset(valid,1,sizeof(valid));//valid[i]=1是素数  
	for(i=2;i<=MAXN;++i)  
	{  
		if(valid[i])  
		{  
			tot++;  
			prime[tot]=i;  
		}  
		for(j=1;((j<=tot)&&(i*prime[j]<=MAXN));++j)  
		{  
			valid[i*prime[j]]=0;  
			if(i%prime[j]==0)  
				break;  
		}  
	}  
	valid[0]=valid[1]=0;
}  

int getSum(int n)
{//求n各位数字和
	int sum=0;
	while(n)
	{
		sum+=n%10;
		n/=10;
	}
	return sum;
}

void getAns()
{
	int i;
	ans[0]=ans[1]=0;
	for(i=2;i<=MAXN;++i)
	{
		if(valid[i]&&valid[getSum(i)])
			ans[i]=ans[i-1]+1;
		else
			ans[i]=ans[i-1];
	}
}

int main()
{
	int L,R,test,cases=1;
	getPrime();
	getAns();
	scanf("%d",&test);
	while(test--)
	{
		scanf("%d %d",&L,&R);
		printf("Case #%d: %d\n",cases++,ans[R]-ans[L-1]);
	}
	return 0;
}