PTA 习题6-5 使用函数验证哥德巴赫猜想 (20分)

习题6-5 使用函数验证哥德巴赫猜想 (20分)

本题要求实现一个判断素数的简单函数，并利用该函数验证哥德巴赫猜想：任何一个不小于6的偶数均可表示为两个奇素数之和。素数就是只能被1和自身整除的正整数。注意：1不是素数，2是素数。

函数接口定义：

int prime( int p );
void Goldbach( int n );

其中函数prime当用户传入参数p为素数时返回1，否则返回0；函数Goldbach按照格式“n=p+q”输出n的素数分解，其中p≤q均为素数。又因为这样的分解不唯一（例如24可以分解为5+19，还可以分解为7+17），要求必须输出所有解中p最小的解。

裁判测试程序样例：

#include <stdio.h>
#include <math.h>

int prime( int p );
void Goldbach( int n );

int main()
{
    int m, n, i, cnt;

    scanf("%d %d", &m, &n);
    if ( prime(m) != 0 ) printf("%d is a prime number\n", m);
    if ( m < 6 ) m = 6;
    if ( m%2 ) m++;
    cnt = 0;
    for( i=m; i<=n; i+=2 ) {
        Goldbach(i);
        cnt++;
        if ( cnt%5 ) printf(", ");
        else printf("\n");
    }

    return 0;
}

/* 你的代码将被嵌在这里 */

输入样例：

89 100

输出样例：

89 is a prime number
90=7+83, 92=3+89, 94=5+89, 96=7+89, 98=19+79
100=3+97, 

代码：

int prime( int p )//此处我直接照搬原来的代码了。
{
    int i,half,flag=1;
    if(p<=1)return 0;
    if(p==2)return 1;
        half=sqrt(p);//这里省事的话可以让half=p/2，不过数学意义上的临界点应该是开方。
        for(i=2;i<=half;i++)
        {
            if(p%i==0)flag=0;//flag是我们定义的一个标记变量，如果能被某个>1的数整除（不是素数）flag的值就会变成0
        }
    return flag;
}
void Goldbach( int n )//
{
    int Now=3,Another=n-Now;//奇素数，初始值给个3没毛病
    if(prime(Another))printf("%d=%d+%d",n,Now,Another);//这里是之前说过的利用1为真
    else
    {
        for(;1;)//这是一个死循环，靠break强制跳出。
        {
            for(;!(prime(++Now)););//利用1为真和++i先自增后参与运算的原理（前面几题有解释过用法），实现把Now替换为下一个质数。
            Another=n-Now;
            if(prime(Another))
            {
                printf("%d=%d+%d",n,Now,Another);
                break;
            }
        }
        
    }
    return 0;
}

一些解释：

   代码里用了一些花里胡哨的判断，前面的博客中有解释的。