1149:哥德巴赫猜想

Description

哥德巴赫（Goldbach ]C.，1690.3.18~1764.11.20）是德国数学家；出生于格奥尼格斯别尔格（现名加里宁城）；曾在英国牛津大学学习；原学法学，由于在欧洲各国访问期间结识了贝努利家族,所以对数学研究产生了兴趣；曾担任中学教师。

1725年，到了俄国，同年被选为彼得堡科学院院士；1725年~1740年担任彼得堡科学院会议秘书；1742年，移居莫斯科，并在俄国外交部任职。 1742年，哥德巴赫在教学中发现，每个不小于6的偶数都是两个素数（只能被1和它本身整除的数）之和。如6＝3＋3，14＝3＋11等等。公元1742年6月7日哥德巴赫写信给当时的大数学家欧拉，欧拉在6月30日给他的回信中说，他相信这个猜想是正确的，但他不能证明。叙述如此简单的问题，连欧拉这样首屈一指的数学家都不能证明，这个猜想便引起了许多数学家的注意。从哥德巴赫提出这个猜想至今，许多数学家都不断努力想攻克它，但都没有成功。 我们不需要你去证明哥德巴赫猜想。 如果哥德巴赫猜想是正确的，一个（不小于6的）偶数，都是两个素数之和。那么这个偶数能被至少一个素数对表示，如14，即可以表示为14＝3＋11，也可以表示为14=7+7。不同的偶数对应的素数对的数目是不一样的，如偶数6，就只能表示为6=3+3。对于每个给定的偶数，我们希望知道有多少素数对的和等于该偶数。

Input

有多组测试数据。每组测试数据占一行，包含唯一的一个正偶数n.(6 <= n <= 224,)。 输出以EOF结束

Output

对于每个输入的偶数，输出一行包含唯一的一个整数：表示有多少个素数对的和是输入的偶数。

Sample Input

6
14

Sample Output

1
2

        补充内容：用Eratosthenes筛法 ( 这个人是谁，名字这么拗口 ) 构造1-n的素数表。  

        筛法的思想非常简单，对于不超过n的每个非负整数p，删除2p，3p，4p，···，当处理完所有数之后，还没有被删除的就是素数。如果用vis[i]表示i已经被删除，筛法的代码可以写成这样:

memset(vis ,0 ,sizeof(vis));
for (int i = 2 ; i <= n ; i++)
    for (int j = i*2 ; j<= n ; j += i) 
        vis[j] = 1;

尽管可以继续改进，但这份代码已经相当高效了。为什么？因为给定外层循环变量i，内层循环的次数是( (n/i)的下界 - 1 )，这个数小于n/i。 这样，循环的总次数小于

n/2 + n/3 + ··· + n/n = O(n*logn)。这个结论来源于欧拉在 1734 年得到的结果：1 + 1/2 + 1/3 + ··· + 1/n =ln(n+1) + γ ，其中欧拉常数 γ=0.577218，这样低的时间复杂度已足以完成题目所需的规定时限。

下面我们来改进这份代码。首先，在“对于不超过n的每个非负整数p”中，p可以限定为素数——只需在第二重循环前加一个判断 if(!vis[i])即可。另外，内层循环也不必从i*2开始——它已经在i=2时被筛掉了。改进后的代码如下：

int m = sqrt(n+0.5);
int c = 0;
memset(vis, 0,sizeof(vis));
for (int i = 2 ; i <= m ; i++) if(!vis[i])
{
    prime[c++] = i;
    for (int j = i*i ; i<= n ; j += i)
    vis[j]=1;
}

注意，上面的代码增加了一项新功能：把素数保存在素数数组里，其中c为素数总数，prime[i]是编号为i的素数(从0开始编号。)。

#include <cstdio>
#define N 8388608
bool vis[N]={false};
int prime[1100000] , Nb_Prime;
int sieve()
{
    int i, j, k;
    for(i=0;i<2046;++i)
        if(!vis[i])
            for(k=(i<<1)+3,j=k*i+i+k;j<N;j+=k)
            vis[j]=true;
            for(i=j=0;i<N;++i)
                if(!vis[i])
                    prime[j++]=(i<<1)+3;
    return j;
}
int main()
{
    int i,n,sum,hn;
    Nb_Prime=sieve();
    while(scanf("%d",&n)!=EOF)
    {
        hn=n>>1;
        for(i=sum=0;prime[i]<=hn;++i)
        if(!(vis[(n-prime[i]-3)>>1]))
        sum++;
        printf("%d\n",sum);
    }
    return 0;
}