POJ2262 素数筛

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
const int MAX=1e7;

bool isprime[MAX];
int num_prime,prime[MAX];

void get_prime()
{
memset(isprime,true,sizeof(isprime));
long long i,j;
num_prime=0;
for(i=2;i<MAX;i++)
{
if(isprime[i]==true)
{
for(j=i*i;j<MAX;j+=i)
{
isprime[j]=false;
}
prime[num_prime++]=i;
}
}

}

int main()
{
get_prime();
int n;
while(~scanf("%d",&n)&&n)
{
int i;
for( i=0;prime[i]<=n/2&&i<num_prime;i++)
{
if(isprime[n-prime[i]])
break;
}

printf("%d = %d + %d\n",n,prime[i],n-prime[i]);
}
return 0;
}