HDU--4548--美素数

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

Input

第一行输入一个正整数T,表示总共有T组数据(T <= 10000)。
接下来共T行,每行输入两个整数L,R(1<= L <= R <= 1000000),表示区间的左值和右值。

Output

对于每组数据,先输出Case数,然后输出区间内美素数的个数(包括端点值L,R)。
每组数据占一行,具体输出格式参见样例。

Sample Input
3
1 100
2 2
3 19
Sample Output
Case #1: 14
Case #2: 1
Case #3: 4
#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;
const int maxa=1000000+10;
int	n,m;
int isprime[maxa];bool prime(int n){
	if(n==1)	return false;
	for(int i=2;i<=sqrt(n);i++)
		if(n%i==0)	return false;
	return true;
}bool mprime(int n){//1000000
	if(n<10)	return true;//1
	else if(n<100)	return prime(n/10+n%10);//12
	else if(n<1000)	return prime(n%10+n/10%10+n/100);//123
	else if(n<10000) return prime(n%10+n/10%10+n/100%10+n/1000);//1234
	else if(n<100000) return prime(n%10+n/10%10+n/100%10+n/1000%10+n/10000);//12345
	else if(n<1000000) return prime(n%10+n/10%10+n/100%10+n/1000%10+n/10000%10+n/100000);//123456
}
void init(){
	memset(isprime,0,sizeof isprime);
	for(int i=1;i<=maxa;i++){
		if(prime(i)&&mprime(i))
			isprime[i]=isprime[i-1]+1;
		else isprime[i]=isprime[i-1];
	}
}
int main(){
	init();
	int t,k=1;
	scanf("%d",&t);
	while(t--){
		scanf("%d%d",&n,&m);
		printf("Case #%d: %d\n",k++,isprime[m]-isprime[n-1]);
	}
	return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值