POJ 2909 / ZOJ 1657 : Goldbach\'s Conjecture - 筛法打素数表

来源:互联网 发布:义乌淘宝网店摄影培训 编辑:程序博客网 时间:2024/05/29 03:29
和前面的POJ 2262类似
输出及细节处理略有不同

题意:
给出一个偶数,要求在这个偶数的范围内,有几对素数和是等于这个偶数的

分析:
这次是验证哥德巴赫猜想的解有多少个,同样核心是打素数表  
思路和POJ 2262差不多,筛法打素数表 

注意素数对的不重复性! 

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int N=1000002;
bool b[N]={1,1,0};
int p[N/2],num=0;

void prime(){
int i,j,sq=(int)(sqrt(N*1.0));
for(i=1;i<sq;i++)
if(!b[i])
for(j=i*i;j<N;j+=i)
b[j]=1;
for(i=2;i<N;i++)
if(!b[i]) p[num++]=i;
}

int main(){
prime();
int n;
while(scanf("%d",&n)&&n){
int flag=0,count=0,i=0;
while(p[i]<n){
if(!b[n-p[i]]&&(n-p[i]==p[i])) count+=2;
else if(!b[n-p[i]]) count++;
i++;
}
printf("%d\n",count/2);
}
return 0;
}







0 0
原创粉丝点击