让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N
(<10^5),请计算不超过N
的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N
。
输出格式:
在一行中输出不超过N
的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
代码:
#include<iostream>
using namespace std;
int main()
{
int n,count=0; //count用于计数
cin >> n;
int i,j,B[100001];
bool A[100001]; //因为测试数据最大为99999,故需要把用作标记的 bool 数组的长度设置很大
for (i = 3; i<=n; i++)
{
if (i%2)
A[i] = true;
else
A[i] = false; //所有大于2的偶数都不是素数
}
A[2] = true;
for (i = 3; i*i <= n; i+=2)
{
if (A[i])
{
for (j = i + i; j <=n; j += i)
{
A[j] = false; //若A[i]为素数,则其倍数一定不是素数
}
}
}
for (i = 2,j=0; i <=n; i++)
{
if (A[i])
{
B[j] = i; //将素数记录到数组B[]中
j++;
}
}
for (i = 1; i < j; i++)
{
if (B[i] - B[i - 1] == 2) //判断相邻素数的差是否为2
{
count++;
}
}
cout << count;
return 0;
}
注:需要注意的地方已添加注释;
在这里引用一位博友的方法,他的方法更好一些,也巧妙一些,不需要用到数组;
博友的地址:https://blog.csdn.net/whl_program/article/details/76601993;(代码如下,在其基础上添加了一些注释,方便理解)
#include <iostream>
using namespace std;
int main()
{
int n;
cin >> n;
int count = 0;
int a = 2, b = 3; //已知2、3是素数,但不是满足题意的相邻素数
for (int i = 5; i <= n; i++) {
bool temp = true;
for (int j = 2; j*j <= i; j++) //j*j很重要,减少运行时间
{
if (i%j == 0) {
temp = false;
break;
}
}
if (temp) { //判断为素数
a = b; //通过赋值得到相邻的素数
b = i;
if ((b - a) == 2) //判断是否满足题意
count++;
}
}
cout << count;
return 0;
}
希望大家共同学习,共同进步!