poj 2262 Goldbach's Conjecture

#include <iostream>
using namespace std;
int prime[100000],end;
bool test(int a)
{
if(a%2==0)
return false;
for(int i=3;i*i<=a;i+=2)
{
if(a%i==0)
return false;
}
return true;
}
void init()
{
int j=0;
for(int i=3;i<=1000003;++i)
{
if(test(i))
prime[j++]=i;
}
end=j-1;
}
int find(int n)
{
int b=0,e=end,mid=(b+e)/2;
while(!(prime[mid]<n&&n<prime[mid+1]))
{
if(n>prime[mid+1])
{
b=mid+1;
mid=(b+e)/2;
}
else
{
e=mid;
mid=(b+e)/2;
}
}
return mid;
}
int main()
{
init();
int n,s,i,j,tag;
while(cin>>n&&n)
{
s=find(n);
tag=0;
for(i=0;i<=s;++i)
{
for(j=s;j>=0;--j)
{
if(prime[i]+prime[j]<=n)
break;
}
if(prime[i]+prime[j]==n)
{
tag=1;
printf("%d = %d + %d\n",n,prime[i],prime[j]);
break;
}
}
if(!tag)
printf("Goldbach's conjecture is wrong.\n");    //事实上不用打印这句，因为还未被证明Wrong
}
return 0;
}

//#include <iostream>        //ac，水过，不用打素数表，直接一个个地试验直到找出结果
//using namespace std;
//bool prime(int a)
//{
//    if(a%2==0)
//        return false;
//    for(int i=3;i*i<=a;i+=2)
//    {
//        if(a%i==0)
//            return false;
//    }
//    return true;
//}
//int main()
//{
//    int n,i;
//    while(cin>>n&&n)
//    {
//        for(i=3;i<n;i+=2)
//        {
//            if(prime(i)&&prime(n-i))
//                break;
//        }
//        printf("%d = %d + %d\n",n,i,n-i);
//    }
//    return 0;
//}