基于Fibonacci数列计算黄金比例的代码报错，如何完美运行？

原代码存在的问题 #

• 参数类型不合理：fibonacci函数接收double类型参数，但实际调用时传入的是整数，递归到非整数n时会出现逻辑混乱，甚至导致栈溢出。
• 整数溢出：fibonacci返回int类型，斐波那契数列增长极快，n≥47时数值会超过32位int的最大值，变为负数或错误值，后续计算比例时会出现除以零或异常结果，这就是你看到“integral goes to zero”错误的核心原因。
• 递归效率极低：递归计算斐波那契的时间复杂度为O(2^n)，n稍大就会变得极慢，甚至触发栈溢出。
• 循环边界错误：main中i从1开始，第一次循环会调用fibonacci_golden_ratio(0)，此时会计算fibonacci(0)/fibonacci(-1)，属于非法操作，导致逻辑混乱。
• 非标准库依赖：conio.h是Windows专属的非标准库，在GCC等跨平台编译器中无法使用，getch()也需要替换。

修改后的代码 #

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

// 迭代计算第n项斐波那契数，用double避免溢出
double fibonacci(int n) {
    if (n == 0) return 0.0;
    if (n == 1 || n == 2) return 1.0;
    
    double prev_prev = 1.0; // F(n-2)
    double prev = 1.0;      // F(n-1)
    double current;
    
    for (int i = 3; i <= n; i++) {
        current = prev_prev + prev;
        prev_prev = prev;
        prev = current;
    }
    return prev;
}

double fibonacci_golden_ratio(int n) {
    // 确保n至少为2，避免除以零
    if (n < 2) {
        printf("错误：n必须大于等于2才能计算黄金比例\n");
        return -1.0;
    }
    return fibonacci(n) / fibonacci(n - 1);
}

int main() {
    int i;
    double prev_phi = 0.0;
    double current_phi;
    
    // 从i=2开始，因为F(2)/F(1)是第一个有效比例
    for (i = 2; ; i++) {
        current_phi = fibonacci_golden_ratio(i);
        printf("F(%d) = %.0f , phi = %.5lf\n", i, fibonacci(i), current_phi);
        
        // 当两次比例的差值小于0.0001时停止迭代
        if (fabs(current_phi - prev_phi) < 0.0001) {
            break;
        }
        prev_phi = current_phi;
    }
    
    printf("\n最终逼近的黄金比例为: %.5lf\n", current_phi);
    getchar();
    return 0;
}

关键修改说明 #

• 改用迭代实现斐波那契：迭代的时间复杂度为O(n)，效率远高于递归，同时避免栈溢出风险。
• 使用double存储斐波那契数：double可以精确表示2^53以内的整数，足够计算到很大的n值，彻底解决整数溢出问题。
• 修正参数类型：fibonacci函数接收int类型参数，符合计算整数项斐波那契数的需求。
• 调整循环起始点：从i=2开始，避免计算无效的n=0或n=1的比例，确保每次计算的都是有效的F(n)/F(n-1)。
• 增加边界检查：在fibonacci_golden_ratio中检查n是否≥2，防止除以零的错误。
• 移除非标准库：去掉conio.h，用标准库的getchar()代替getch()，提升代码的跨平台兼容性。
• 优化比例计算逻辑：循环中只计算一次当前比例，避免重复调用函数，提升运行效率。

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