让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N(<10次方5),请计算不超过N的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N。
输出格式:
在一行中输出不超过N的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
#include<stdio.h>
#include<math.h>
const int maxn = 100005;
int prime[maxn], pNum = 0;
bool p[maxn] = {0}; //false表示是素数,true为合数
//传统做法,以后只记埃氏筛法
/*bool isPrime(int n) {
if(n <= 1) {
return false;
}
int sqr = (int) sqrt (1.0 * n);
for(int i = 2; i < sqr; ++i) {
if(n % i == 0) {
return false;
}
}
return true;
}
void Find_prime_2(int n) {
for(int i = 0; i < n; ++i) {
if(isPrime(i) == true) {
prime[pNum++] = i;
p[i] = true; //这个true是素数
}
}
}*/
void Find_prime(int n) {
for(int i = 2; i <= n; ++i) { //埃氏筛选法
if(p[i] == false) {
prime[pNum++] = i;
for(int j = i + i; j <= n; j += i) { //剔除i的倍数
p[j] = true;
}
}
}
}
int main() {
int n;
scanf("%d", &n);
Find_prime(n);
int count = 0; //用于计数素数对
for(int i = 1; i < pNum; ++i) {
if(prime[i] - prime[i - 1] == 2) {
count++;
}
}
printf("%d", count);
} ```
return 0;