poj 2262 Goldbach's Conjecture 【素数筛】

这题必然的筛选法，出现了2个问题：1.开始开了一个 result 数组（全局变量），想把素数挨个存进来，虽然还计数估算了的，一直的runtime error ， 后来发现是多此一举，去掉之后就变成wrong answer，看来可以编译了。这么说来，一百万对于两个大数组还是有点吃不消的  2.筛的时候一定要筛完整，for(j=2*i;j<MAXN;j+=i)prime[j]=1;这里开始用的

j<(MAXN/i)，明显就错了。

另外我很好奇题目中说的"Goldbach's conjecture is wrong."，很明显不会wrong的，wrong了科学家们还继续证明干嘛，可是我AC以后看网上的结题报告居然还真有把wrong考虑进去的人。。。

2个错误都存在的代码

#include <stdio.h>#define MAXN 1000002#define PrimeNum 671000int prime[MAXN];//用筛法求素数，1代表不是素数（被筛掉）int result[PrimeNum];int main(){//先打出素数表prime[0]=prime[1]=1;    //开始去掉prime[0]和prime[1]int i,j;for(i=2;i<MAXN;i++){if(prime[i]==0)//如果prime[i]是素数就把他的倍数都筛掉{for(j=2*i;j<(MAXN/i);j+=i)prime[j]=1;}}//test 1/*int count=0;for(i=2;i<MAXN;i++)if(prime[i]==0)count++;printf("一百万以内的素数有 %d 个\n",count);*///test 2/*for(i=2;i<100;i++)if(prime[i]==0)printf("%d ",i);*/int count=0;for(i=3;i<100;i++)if(prime[i]==0){result[count]=i;count++;}/*for(i=0;i<100;i++)printf("%d ",result[i]);*/int n;while(scanf("%d",&n)){if(n==0)break;for(i=0; ;i++)if(prime[n-result[i]]==0){printf("%d = %d + %d\n",n,result[i],n-result[i]);break;}}return 0;}

正确的代码

#include <stdio.h>#define MAXN 1000002int prime[MAXN];//用筛法求素数，1代表不是素数（被筛掉）int main(){//先打出素数表prime[0]=prime[1]=1;    //开始去掉prime[0]和prime[1]int i,j;for(i=2;i<MAXN;i++){if(prime[i]==0)//如果prime[i]是素数就把他的倍数都筛掉{for(j=2*i;j<MAXN;j+=i)prime[j]=1;}}int n;while(scanf("%d",&n)){if(n==0)break;for(i=3; ;i++)if(prime[i]==0 && prime[n-i]==0){printf("%d = %d + %d\n",n,i,n-i);break;}}return 0;}