题目要求
让我们定义 d n d_n dn为: d n d_n dn= p n + 1 p_{n+1} pn+1 - p n p_n pn,其中 p i p_i pi是第i个素数。显然有 d 1 d_1 d1=1,且对于n>1有 d n d_n dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N(< 1 0 5 10^5 105),请计算不超过N的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N。
输出格式:
在一行中输出不超过N的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
代码
#include <iostream>
#include<math.h>
using namespace std;
int isPrime(int n)
{
for (int i=2;i<=(int)sqrt(n);i++)
if(n % i == 0) return 0;
return 1;
}
int main()
{
int n,num=0,count=0;
int x[10000];
cin >> n;
for(int i = 2;i<=n;i++)
{
if(isPrime(i)) x[num++]=i;
}
for(int i=1;i<num;i++)
{
if(x[i]-x[i-1] == 2) count++;
}
cout << count;
}
思路
题目很简单,第一步先找出输入以内的所有素数,再对相邻的素数求差,如果满足差为2,则输出。
- 那么唯一的难点可能还是在素数判断上:
int isPrime(int n)
{
for (int i=2;i<=(int)sqrt(n);i++)
if(n % i == 0) return 0;
return 1;
}
如上定义一个函数来判断是否为素数,即遍历从2到sqrt(n)
,如果可以整除,则不是素数。
- 将相应的素数放到一个数组里面
for(int i = 2;i<=n;i++)
{
if(isPrime(i)) x[num++]=i;
}
- 对相邻两者求差,为2则输出
for(int i=1;i<num;i++)
{
if(x[i]-x[i-1] == 2) count++;
}
如此,完成!