PAT 乙级 1007. 素数对猜想 (20) c++ 筛选法求素数

来源:互联网 发布:sop8单片机的脚位介绍 编辑:程序博客网 时间:2024/05/29 19:40

PAT 乙级 1007. 素数对猜想 (20) c++ 筛选法求素数

让我们定义 dn 为:dn = pn+1 - pn,其中 pi 是第i个素数。显然有 d1=1 且对于n>1有 dn 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

现给定任意正整数N (< 105),请计算不超过N的满足猜想的素数对的个数。

输入格式:每个测试输入包含1个测试用例,给出正整数N。

输出格式:每个测试用例的输出占一行,不超过N的满足猜想的素数对的个数。

输入样例:
20
输出样例:
4
思路分析: 首先这道题如果用普通的求素数的方法写可能会超时,毕竟10的五次方了 然后还要去循环判断每一个数,我也是在网上看了别人的代码才知道有筛选法求素数的。筛选法求素数的基本思想就是把从1开始的、某一范围内的正整数从小到大顺序排列, 1不是素数,首先把它筛掉。剩下的数中选择最小的数是素数,然后去掉它的倍数。依次类推,直到筛子为空时结束。。 最根本的思想可以说是不难,但实现起来。。其实实现起来也不难。

代码实现:

#include<iostream>using namespace std;bool isprime [100001]; // bool数组标志1-100000的数据是否为素数int count = 0,tmp =2;void   getprime(int n){   //运用筛选法获得n以内的全部素数    for(int i=2;i<=n;i++){  //初始化为true即默认都是素数        isprime[i]=true;    }    for(int j=2;j<=n;j++){   //依次遍历后面的所有标志为true的数据    if(isprime[j]==true){    //从第一个素数2开始判断 设置2的倍数(非素数)为false 因为素数的倍数肯定不是素数他至少有一个因子是该素数本身        for(int m=2;j*m<=n;++m){//  m为倍数 j为素数 去掉j的m倍即为去掉非素数            isprime[j*m]=false;        }    }}for(int k=2;k<=n;k++){        if(isprime[k] == true){   //如果这个数据是素数             if(k - tmp == 2){  //如果这个数减去上一个相邻的素数结果为2 即素数对数加1                count ++;                tmp = k;    //更新tmp为本次的素数            }            else{                tmp = k;            }        }    }cout<<count;}int main(){    int N;    scanf("%d",&N);    getprime(N);    return 0;}