让我们定义dn为:dn=pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N
(<105),请计算不超过N
的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N
。
输出格式:
在一行中输出不超过N
的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
/*
题目初解:
题目要从2开始查找直到N的素数,不妨将N设置到最大,题目要求素数
对之间的差值dn为2。逐个素数查找相减就可以,考虑到素数除2
以为只能是奇数,那么可以间隔2进行寻找。题目寻找的只是个数,不需要
对素数进行存储,因此我们只需要知道当前素数和下一个素数即可,然后计
算差值,满足dn==2则计数器加1即可。
补充:
另外考虑一点,就是素数检验的边界,因子的大小,只需要在算术平方根以内就可了
因为比算术平方根大的的因子一定存在着一个比算数平方根小的因子。因子一定是成对
出现的,除了算术平方根;
比如100的算术平方根是10,2*50==100,虽然边界是,但是2对应着比10大
的因子50,因此如果只要在取平方根之前的数中没有因子,那么就一定是素数了
这样可以减少很多次循环,避免超时。
*/
#include<iostream>
#include<cmath>
using namespace std;
int main(){
//输入N;
int N;
cin>>N;
int x; //定义N中循环变量;
int i; //定义整除变量;
int fence=0; //每一个整数的整除边界;
int Pm=1,Pn=0; //相邻素数;
int count=0; //素数对计数器;
int flag=1; //素数标志,若为1则是素数,为0则不是;
for(int x=2;x<=N;x++){ //逐个整数处理;
//判断当前整数是否为素数;
//优化素数判断条件;
fence=(int)sqrt((double)x); //这里注意sqrt函数的参数和返回值都是double类型。
for(int i=2;i<=fence;i++){
if(x%i==0){
flag=0; //可以整除表示有因子,即不是素数;
}
}
if(flag==1){
if(x!=2){ //这里排除了2这种情况,因为我Pn设置的0,不排除会多出一对2-0=0;
Pm=x;
if(Pm-Pn==2){
count++;
}
Pn=Pm; //将素数检测调至当前素数;
}
}
else{
flag=1; //若上一个数不是素数,则需要将标志重置;
}
}
//输出;
cout<<count<<endl;;
}
/*
解题反思:
第一次做这个题的时候大数运行超时,具体原因不确定,但是在CSDN上看到
素数判断条件的优化,与之相较的是我自己的代码并没有进行判断优化,而是直接
从2开始循环增大除数到x-1,耗费较多循环次数导致运行时间很长,并且我采用的
是穷举列出素数表,再进行相邻的相减判断差值,但这里素数对差值的判断是必须的。
个人认为此处没有太大优化价值,核心点在于判断条件优化后大幅度减少循环次数,
从而减少运行时间,还有编写出错时一点:没有写x不是素数时标志重置,导致flag
变0后一直为0,运行结果不正确。然后则是sqrt(2)这里,与其他素数判断不同,这里
没有进入判断的循环中。
*/