C语言计算不同参数几何级数出现超大值问题求助

1. 数学本质：发散几何级数的必然结果 #

几何级数求和公式适用于所有r≠1的场景，但当**|r|>1**时，随着项数n增大，r^(n+1)会呈指数级爆炸增长：

• 场景2：r=1.1，n=500，1.1^501本身就是天文数字，最终求和结果自然是超大值
• 场景3：r=2，n=100，2^101=25353012004564588032，乘以a=0.0001后与你输出的结果一致（浮点数精度导致末尾微小差异）

场景1的r=0.01<1，r^(n+1)趋近于0，级数收敛到a/(1-r)=2/0.99=2.020202...，所以结果正常。

2. 代码中的小瑕疵（不影响结果正确性，但可优化） #

• 冗余输出：场景2中printf("%.6Lf \n",formula);属于无效输出，此时formula未赋值，会打印随机值
• 数据类型不一致：场景3中a=0.0001、r=2.0未加L后缀，会被当作double类型赋值给long double变量，虽无严重影响，但建议统一用0.0001L、2.0L保持类型一致
• 无效循环：场景1中n=10000，但后面的项因r过小趋近于0，循环到后期属于无效计算

1. 明确结论 #

场景2和3的超大数值是正确的数学结果，并非代码逻辑错误。如果预期结果是“合理大小”，需检查场景参数是否符合实际需求（比如是否r的取值应小于1，或n的取值过大）。

2. 优化后的代码 #

以下代码修复了上述瑕疵，并添加了提前终止循环的逻辑（当项的绝对值小于精度阈值时停止，避免无效计算）：

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

int main(){
    // 几何级数计算脚本
    long double n = 0, a = 0, r = 0;
    long double summation = 0, formula = 0;
    const long double EPS = 1e-15; // 精度阈值，小于该值的项可忽略

    printf("==============================\n");
    printf("======== Switch Case =========\n");

    for(int i = 1; i <= 3; i++ ){
        switch (i){
            case 1:
                n = 10000L, a = 2.0L, r = 0.01L;
                summation = a;
                // 优化：项过小则提前终止循环
                for(int j = 1; j <= n; j++){
                    long double term = a * powl(r, j);
                    if(fabsl(term) < EPS) break;
                    summation += term;
                }
                formula = a * (1 - powl(r, n+1)) / (1 - r);
                printf("Case 1 summation = %.6Lf \n", summation);
                printf("Case 1 formula = %.6Lf \n", formula);
                printf("==============================\n");
                break break;
            case 2:
                n = 500L, a = 0.01L, r = 1.1L;
                summation = a;
                for(int j = 1; j <= n; j++){
                    summation += a * powl(r, j);
                }
                formula = a * (1 - powl(r, n+1)) / (1 - r);
                printf("Case 2 summation = %.6Lf \n", summation);
                printf("Case 2 formula = %.6Lf \n", formula);
                printf("==============================\n");
                break;
            case 3:
                n = 100L, a = 0.0001L, r = 2.0L;
                summation = a;
                for(int j = 1; j <= n; j++){
                    summation += a * powl(r, j);
                }
                formula = a * (1 - powl(r, n+1)) / (1 - r);
                printf("Case 3 summation = %.6Lf \n", summation);
                printf("Case 3 formula = %.6Lf \n", formula);
                printf("==============================\n");
                break;
        }
    }
    printf("=======End of Switch Case=======\n");
    return 0;
}

3. 额外提示 #

• 当|r|>1时，公式可等价转换为a*(r^(n+1)-1)/(r-1)，避免负数运算，结果一致
• long double的精度高于double，能处理更大范围的数值，但超过其表示极限时会输出inf（你的场景尚未达到该极限）

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