PAT乙级 1007 素数对猜想

最新推荐文章于 2024-07-25 15:38:54 发布

Luuuuyang

最新推荐文章于 2024-07-25 15:38:54 发布

阅读量132

点赞数

分类专栏： PAT乙级素数

本文链接：https://blog.csdn.net/HPU_LY/article/details/86562822

版权

PAT乙级同时被 2 个专栏收录

27 篇文章 0 订阅

订阅专栏

素数

4 篇文章 0 订阅

订阅专栏

让我们定义dn为： $d_{n}=p_{n+1}-p_{n}$ ，其中 $p_{i}$ 是第个素数。显然有 $d_{1}=1$ ，且对于 n>1 有 $d_{n}$ 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

现给定任意正整数 $N(<10^{5})$ ，请计算不超过N的满足猜想的素数对的个数。

输入格式:

输入在一行给出正整数N。

输出格式:

在一行中输出不超过N的满足猜想的素数对的个数。

输入样例:

输出样例:

代码

#include<iostream>
#include<cstring>
using namespace std;
#define MAXN 100010
int prime[MAXN];
bool is_prime[MAXN];
int n=0;
void Prime(){
	int i,j;
	memset(is_prime,1,sizeof(is_prime));
	is_prime[0]=is_prime[1]=0;
	for(i=2;i<=MAXN;i++){
		if(is_prime[i]){
			prime[++n]=i;
			for(j=i+i;j<=MAXN;j+=i)
				is_prime[j]=0;
		}
	}
}
int ans[MAXN];
void Ans(){
	int i,j;
	ans[0]=ans[1]=ans[2]=0;
	for(i=3,j=2;i<MAXN;i++){
		if(i==prime[j]){
			if(prime[j]-prime[j-1]==2)
				ans[i]=ans[i-1]+1;
			else ans[i]=ans[i-1];
			j++;
		}
		else ans[i]=ans[i-1];
	}
}
int main(){
	Prime();
	Ans();
	int N;
	while(cin>>N){
		cout<<ans[N]<<endl;
	}
	return 0;
}