美素数(素数表与线段树查询)

Problem Description
  小明对数的研究比较热爱,一谈到数,脑子里就涌现出好多数的问题,今天,小明想考考你对素数的认识。
  问题是这样的:一个十进制数,如果是素数,而且它的各位数字和也是素数,则称之为“美素数”,如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

//线段树就有查询与更新的功能,首先说说建树,如果该数是美素数,那么该数对应的节点值为1,否则为0,输出结果即(l,r)区间中1的个数。学了线段树得用上啊,不能光理解,要懂得应用。

#include<stdio.h>
#include<string.h>
#define Max 1000005
int vis[Max];
struct 
{
	int left,right,data;
}tree[4*Max];
int get_num(int x)
{
	int sum=0;
	while(x)
	{
		sum+=(x%10);
		x/=10;
	}
	return sum;
}
void make_prime()
{
	int i,j;
	memset(vis,0,sizeof(vis));
	vis[1]=1;
	for(i=2;i<Max;i++)
		if(!vis[i])
			for(j=i+i;j<Max;j+=i)
				vis[j]=1;
}
void build(int l,int r,int nd)
{
	tree[nd].left=l,tree[nd].right=r,tree[nd].data=0;
	if(!vis[l]&&!vis[get_num(l)]) tree[nd].data=1;
	if(l==r) return;
	int mid=(l+r)/2;
	build(l,mid,2*nd+1);
	build(mid+1,r,2*nd+2);
	tree[nd].data=(tree[2*nd+1].data+tree[2*nd+2].data);
}
int query(int l,int r,int nd)
{
	if(tree[nd].left==l&&tree[nd].right==r) return tree[nd].data;
	int mid=(tree[nd].left+tree[nd].right)/2;
	if(l<=mid)
	{
		if(r>=mid+1) return query(l,mid,2*nd+1)+query(mid+1,r,2*nd+2);
		else return query(l,r,2*nd+1);
	}
	else return query(l,r,2*nd+2);
}
int main()
{
	//freopen("b.txt","r",stdin);
	make_prime();
	build(1,Max-1,0);
	int t,l,r,cas=0;
	scanf("%d",&t);
	while(t--)
	{
		printf("Case #%d: ",++cas);
		scanf("%d %d",&l,&r);
		printf("%d\n",query(l,r,0));
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值