《算法竞赛入门经典》4.13应用举例-孪生素数

//例题4-2
/*
*如果n和n+2都是素数，则称它们是孪生素数。输入m，输出两个数均不超过m的最大孪生素数。5<=m<=100。例如m=2时答案是17、19，m=1000时答案是881、883。
*/

//程序4-2 孪生素数（1）
#include <stdio.h>
/* do NOT use this if x is very large or small
*n太小时n=1会被错误地判断为素数，太大时i*i可能溢出
*/
//判断素数
int is_prime(int x)
{
    int i;
    for(i = 2; i*i <= x; i++)    //判断不超过sqrt(x)的整数i
        if(x % i ==0) return 0; //一旦发现有一个大于1的因子，立刻返回0（假）
    return 1;                    //最后返回1（真）
}

int main()
{
    int i, m;
    scanf("%d", &m);
    for(i = m-2; i >= 3; i--)
        if(is_prime(i) && is_prime(i+2))
        {
            printf("%d %d
", i, i+2);
            break;
        }
    return 0;
}

//程序4-2 孪生素数（2）
#include <stdio.h>
#include <math.h>
#include <assert.h>
//判断素数
int is_prime(int x)
{
    int i, m;
    assert(x >= 0);    //assert.h宏限制非法函数调用，当x>=0不成立时程序异常终止并给出提示消息
    if(x == 1) return 0;
    m = floor(sqrt(x) + 0.5);    //避免每次重复计算sqrt(x)和浮点误差
    for(i = 2; i*i <= x; i++)    //判断不超过sqrt(x)的整数i
        if(x % i ==0) return 0; //一旦发现有一个大于1的因子，立刻返回0（假）
    return 1;                    //最后返回1（真）
}

int main()
{
    int i, m;
    scanf("%d", &m);
    for(i = m-2; i >= 3; i--)
        if(is_prime(i) && is_prime(i+2))
        {
            printf("%d %d
", i, i+2);
            break;
        }
    return 0;
}