巴比伦平方根算法出现无限循环问题求助

你的代码在计算99999这类大数时陷入无限循环，核心原因是float类型的精度限制，以及你采用的绝对误差终止条件不适合大数场景。以下是具体的解决思路和修改方案：

问题根源 #

float类型仅包含23位尾数，对应约6-7位十进制有效数字。对于99999这类大数，其平方根约为316.227766，用float存储该近似值时，精度不足以让x*x - n的差值小于你设定的0.000001，导致循环永远无法终止。

解决方案 #

1. 改用相对误差终止条件 #

绝对误差阈值0.000001对大数不友好，换成相对误差判断更合理——检查近似值与真实值的相对偏差，而非绝对差值。比如可以直接判断迭代过程中x和y的相对变化：

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

float squareRoot(float n)
{
    float x = n;
    float y = 1;
    float e = 0.000001; /* 相对精度阈值 */

    // 检查x与y的相对偏差，替代原绝对误差条件
    while (fabs(x - y) > e * x) {
        x = (x + y) / 2;
        y = n / x;
    }

    return x;
}

int main()
{
    int n;
    printf("请输入要求平方根的数字\n");
    scanf("%d", &n);

    printf("%d的平方根是%.8f\n", n, squareRoot(n));
    return 0;
}

也可以用相对误差判断x*x与n的偏差：fabs(x*x - n) > e * n，效果类似。

2. 增加最大迭代次数限制 #

巴比伦法收敛速度极快，正常10次迭代内就能达到很高精度。给循环加一个迭代次数上限，即使精度没达标也强制退出，彻底避免死循环：

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

float squareRoot(float n)
{
    float x = n;
    float y = 1;
    float e = 0.000001;
    const int MAX_ITER = 100; // 设定最大迭代次数
    int iter_count = 0;

    while ((x*x - n) > e && iter_count < MAX_ITER) {
        x = (x + y) / 2;
        y = n / x;
        iter_count++;
    }

    if (iter_count == MAX_ITER) {
        printf("已达最大迭代次数，当前近似值为最优解\n");
    }

    return x;
}

int main()
{
    int n;
    printf("请输入要求平方根的数字\n");
    scanf("%d", &n);

    printf("%d的平方根是%.8f\n", n, squareRoot(n));
    return 0;
}

3. 优化初始近似值 #

原代码用x = n作为初始值，对大数来说收敛速度慢。可以改用更接近真实平方根的初始值，比如x = n / 2，或者如果作业允许，用x = sqrtf(n)直接获取一个接近的初始值，能大幅减少迭代次数。

内容的提问来源于stack exchange，提问作者Gustavo