让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N(<105),请计算不超过N的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N。
输出格式:
在一行中输出不超过N的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
算法:
用线性筛法求素数,然后遍历素数,符合条件则,res + 1,。
#include <iostream>
using namespace std;
const int N = 100010;
int primes[N], cnt;
bool st[N];
void get_primes(int n)
{
for(int i = 2; i <= n; i ++)
{
if(!st[i]) primes[cnt++] = i;
for(int j = 0; i * primes[j] <= n; j++)//给当前i乘上一个质因子
{
st[i * primes[j]] = true;
if(i % primes[j] == 0) break;
}
}
}
int main()
{
int n, res = 0;
cin >> n;
get_primes(n);
for(int i = 0; i < cnt; i++)
if(primes[i] + 2 == primes[i + 1])
res ++;
cout << res << endl;
return 0;
}