python acm 素数个数_湘潭大学OJ-1098求区间内素数个数问题

求区间内素数个数问题

题目描述

Description

给定两个非负整数a，b，其中0<= a,b<=1,000,000，请计算这两个数之间有多少个素数。限制：Time Limit : 1000 MS Memory Limit : 65536 KB

输入

第一行是一个整数K(1<=K<=1000)，表示有多少个样例，每个样例占一行，是两个整数a和b，每个整数之间用一个空格隔开。

输出

每行输出一个样例的结果。

Sample Input

2

2 3

17 19

Sample Output

2

2

解题思路

这个题目思路很简单，就是判断素数，然后计算个数。

主要的问题是如何判断某个数是否是素数，以及采用什么方式来计算个数

首先这个题目中a,b的值可能会很大，所以常规的求素数的方法是不可取的。相对简单的一个方法就是在循环外面进行打表，之后只要看对应下标的值来判断该数是否为素数就好了。现有一个数组noprime，下标 i 对应的noprime[ i ] 为1，则表示 i 为非素数，对应的noprime[ i ] 为0，则表示 i 为素数。(若一个数为某个数的倍数，那么这个数一定是非素数)

int noprime[1000050];

noprime[0] = 1;

noprime[1] = 1;

for(int i = 2; i<= 1000000; i++) {

if(!noprime[i]) {

for(int j = 2*i; j<= 1000000; j+=i) {

noprime[j] = 1;

}

}

}

代码如下

第一种方法

打表后直接循环a到b，判断每一个数是否为素数。

#include

#include

using namespace std;

int noprime[1000050];

int main() {

noprime[0] = 1;

noprime[1] = 1; //1表示非素数，0表示素数

for(int i = 2; i<= 1000000; i++) {

if(!noprime[i]) {

for(int j = 2*i; j<= 1000000; j+=i) {

noprime[j] = 1;

}

}

}

int k;

scanf("%d",&k);

while(k--) {

int a,b;

scanf("%d%d",&a,&b);

if(a>b) swap(a,b);

int num = 0;

for(int i = a; i<= b; i++) { //直接在a和b范围内数素数的个数

if(!noprime[i]) num++;

}

printf("%d\n",num);

}

return 0;

}

第二种方法

打表后，采用前缀和算法直接得出素数的个数。用一个数组s[ ]记录到该数为止前面一共有多少个素数，那么若要得出a 到 b的素数个数，可以直接用s[b]-s[a-1]得出。

比如说现在判断0到7的素数个数：

0为非素数，初始化s[0]=0;

1为非素数，所以到1为止一共有0个素数，就是在s[0]的基础上加0；

2为素数，到2为止一共有1个素数，就是在s[1]的基础上加上1；

3为素数，到3为止一共有2两个素数，就是在s[2]的基础上加上1；

…

所以我们注意到我们只需要考察当前这个数是否为素数，若为素数，就需要在前面得出的个数上加1，若不是素数，则仍保持为之前得出的数字。这样s[ i ]表示的就是到目前 i 这个位置为止，素数的个数。

int s[1000050];

s[0] = 0;

for(int i = 1; i<= 1000000; i++) {

if(noprime[i]==1) s[i] = s[i-1];

else if(noprime[i]==0) s[i]=s[i-1]+1;

}

总的代码如下。

#include

#include

using namespace std;

int s[1000050];

int noprime[1000050];

int main() {

noprime[0] = 1;

noprime[1] = 1;

for(int i = 2; i<= 1000000; i++) {

if(!noprime[i]) {

for(int j = 2*i; j<= 1000000; j+=i) {

noprime[j] = 1;

}

}

}

s[0] = 0;

for(int i = 1; i<= 1000000; i++) {

if(noprime[i]==1) s[i] = s[i-1];

else if(noprime[i]==0) s[i]=s[i-1]+1;

}

int k;

scanf("%d",&k);

while(k--) {

int a,b;

scanf("%d%d",&a,&b);

if(a>b) swap(a,b);

if(a==0) printf("%d\n",s[b]);

else printf("%d\n",s[b]-s[a-1]);

}

}

标签：noprime,OJ,int,个数,1098,素数,1000050,include

来源： https://blog.csdn.net/weixin_43966538/article/details/106255569