POJ-2262-Goldbach's Conjecture 解题报告

       验证哥德巴赫猜想。题意：任何一个大于4的偶数都能写成两个奇素数的和。现在给你一个大于等于6的数n，请你分解成n = a + b的形式，并且a和b都是奇素数，如果有多组分解形式，要求输出b - a最大的那种。如果不能够分解，输出指定的信息。

       我的解题思路：首先，根据计算机科学家运行的结果，在ACM的世界里（64位整型范围内），哥德巴赫猜想是成立的，所以一定能够分解。因此将范围内的素数筛出来，然后进行枚举判断就是了。

       我的解题代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>

using namespace std;

const int N = 1000001;

bool isprime[N];
int primes[N], pn;
int n;

void InitRead();

void DataProcess();

void FastSieve(int maxn);

int main()
{
    InitRead();
    while (~scanf("%d", &n))
    {
        if (n == 0) break;
        DataProcess();
    }
    return 0;
}

void InitRead()
{
    memset(isprime, true, sizeof(isprime));
    isprime[0] = isprime[1] = false;
    pn = 0;
    FastSieve(N-1);
    return;
}

void DataProcess()
{
    for (int i=0; i<pn; ++i)
    {
        if (isprime[n - primes[i]])
        {
            printf("%d = %d + %d\n", n, primes[i], n - primes[i]);
            break;
        }
    }
    return;
}

void FastSieve(int maxn)
{
    for (int i=2; i<=maxn; ++i)
    {
        if (isprime[i]) primes[pn++] = i;
        for (int j=0; j<pn; ++j)
        {
            if (i * primes[j] > maxn) break;
            isprime[i * primes[j]] = false;
            if (i % primes[j] == 0) break;
        }
    }
    return;
}