题目描述
让我们定义 d n d_n dn 为: d n = p n + 1 − p n d_n=p_{n+1}−p_n dn=pn+1−pn,其中 p i p_i pi 是第 i 个素数。显然有 d 1 = 1 d_1=1 d1=1,且对于 n > 1 n>1 n>1有 d n d_n dn 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数 N ( < 1 0 5 ) N(<10^5) N(<105),请计算不超过 N N N 的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N。
输出格式:
在一行中输出不超过N的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
方法
解题思路
当我们要求 n 以内素数的个数时,我们采用素数筛法,先定义一个状态数组来标记每个数字的状态,false 表示 素数,true 表示合数。在遍历时如果数字 i 没有被标记为 true 就说明它是素数,并存入 primes 数组中,然后再循环把 i 的整数倍 (>=2) 标记为 true。
先用埃氏筛法把素数筛选出来存入 primes 数组中,然后遍历 primes,统计满足素数相邻且差值为 2 的对数。
代码
#include <iostream>
using namespace std;
const int N = 1e5;
int n;
bool st[N];
int primes[N];
int main()
{
cin >> n;
int k = 0;
for(int i = 2; i <= n; i++)
{
if(!st[i])
primes[k++] = i;
for(int j = 2 * i; j <= n; j += i)
st[j] = true;
}
int cnt = 0;
for(int i = 0; i < k - 1; i++)
if(primes[i + 1] - primes[i] == 2)
cnt++;
cout << cnt;
return 0;
}