poj_2262 Goldbach's Conjecture

//C++代码
#include<iostream>
using namespace std;

bool a[1000001];
int p[300001];

void prim(){  //素数筛法，筛出1-1000000里的素数
int i,j,num=0;
for(i=0;i<=1000000;i++) a[i]=true;
for(i=2;i*i<=1000000;i++){
if(a[i]){
j=i+i;
while(j<=1000000){
a[j]=false;
j+=i;
}
}
}
for(i=2;i<=1000000;i++){
if(a[i]) p[num++]=i;
}
}

int main(){
prim();
int n;
while(cin>>n,n){
for(int i=0;;i++){
if(a[n-p[i]]){
cout<<n<<" = "<<p[i]<<" + "<<n-p[i]<<endl;
break;
}
}
}
return 0;
}