Digit Primes

最新推荐文章于 2023-06-11 17:34:25 发布

原创最新推荐文章于 2023-06-11 17:34:25 发布 · 577 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#树状数组 #素数筛

数论同时被 2 个专栏收录

9 篇文章

订阅专栏

树状数组

6 篇文章

订阅专栏

链接：http://acm.hust.edu.cn/vjudge/problem/25171/origin

题目：

A prime number is a positive number, which is divisible by exactly two different integers. A digit prime

is a prime number whose sum of digits is also prime. For example the prime number 41 is a digit prime

because 4 + 1 = 5 and 5 is a prime number. 17 is not a digit prime because 1 + 7 = 8, and 8 is not

a prime number. In this problem your job is to nd out the number of digit primes within a certain

range less than 1000000.

题意：找到所给区间里数字素数的个数，所谓数字素数就是把所有数位上的数字加到一起的和也是素数。

分析：首先数据规定了不高于1000000，所以基本可以用于处理打表来实现，打表的方法是素数筛，之后二分去数应该是可以的，这里打表之后使用树状数组记录的，树状数组是一种查询速度非常快的记录工具，之后的文章里会仔细讲解。

题解：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <string>
#include <cstring>
#include <functional>
#include <cmath>
#include <cctype>
#include <cfloat>
#include <climits>
#include <complex>
#include <deque>
#include <list>
#include <set>
#include <utility>
using namespace std;

bool vi[1000010];
int pri[1000010];
int s[1000010];
void add(int x)
{
	while(x<=1000000){
		s[x]++;
		x+=x&(-x);
	}
}
int sum(int x)
{
	int ans=0;
	while(x>0){
		ans+=s[x];
		x-=x&(-x);
	}
	return ans;
}
int get(int x)
{
	return sum(x)-sum(x-1);
}
void getprime()
{
	int num=0;
	for(int i=2;i<=1000000;i++){//素数筛
		if(!vi[i]){
			pri[num++]=i;
			for(int j=i+i;j<=1000000;j+=i)
				vi[j]=true;
		}
	}
	for(int i=0;i<num;i++){
		int check=0,temp=pri[i];
		while(temp){
			check+=temp%10;
			temp/=10;
		}
		if(!vi[check])
			add(pri[i]);
	}
}
int main()
{
	getprime();
	int n;
	scanf("%d",&n);
	while(n--){
		int a,b;
		scanf("%d %d",&a,&b);
		printf("%d\n",sum(b)-sum(a-1));
	}
	return 0;
}