C语言计算圆周率近似值的两种常用算法实现

周予安

发布时间：2026-10-16 | 3人浏览过

原创

利用C语言编程计算圆周率近似值，主要涉及级数展开与阶乘递推两种经典算法。本文将详细解析莱布尼茨公式（pi/4 = 1-1/3+1/5...）和阶乘关系式（pi/2 = 1+1/3+2!/(3×5)...）的实现原理，提供包含精度控制和变量处理的完整代码示例，解决从理论公式到实际代码转换中的关键问题。

基于莱布尼茨级数的实现

莱布尼茨级数是一种简单的无穷级数，其数学表达为：
π/4 = 1 - 1/3 + 1/5 - 1/7 + ...
要计算 π 的近似值，只需累加该级数直到最后一项的绝对值小于指定阈值（如 1e-4）。这种方法逻辑直观，适合初学者理解循环与浮点数运算。

第1步：准备条件与数据类型选择

由于涉及分数累加且需要控制小数精度，建议使用 float 或 double 类型存储结果。需注意循环终止条件为“最后一项绝对值”小于阈值，而非循环次数。

第2步：编写累加逻辑代码

以下是基于该公式的标准实现代码。程序通过变量 t 控制正负号交替，n 控制分母递增：

#include <stdio.h>
int main(void)
{
    float pi=0, t=1.0, n=1.0;
    // 当 1/n >= 1e-4 时继续累加
    while(1/n >= 1e-4) 
    {
        pi = pi + t * 1.0/n;
        // 符号取反：1 -> -1 -> 1
        t = -1 * t;
        // 分母每次增加2：1 -> 3 -> 5
        n += 2;
    }
    // 最后乘以4得到 pi
    pi = pi * 4;
    printf("n=%f, pi=%f\n", n, pi);
    return 0;
}

第3步：精度调整说明

代码中使用的阈值 1e-4 仅为示例。若需更高精度，可将阈值缩小至 1e-6 或更小，但会增加循环次数和计算时间。在实际工程中，应根据性能需求平衡精度与效率。

基于阶乘关系的进阶解法

在某些编程题目（如 LeetCode 风格题）中，可能会给出更复杂的递推公式。例如：
π/2 = 1 + 1/3 + 2!/(3×5) + 3!/(3×5×7) + ...
此类题目通常要求输入一个阈值，当某一项的值小于该阈值时停止计算，并输出最终结果。

第1步：分析数学规律

观察公式可知：

分子是阶乘序列：1!, 2!, 3!, ...
分母是奇数连乘：3, 3×5, 3×5×7, ...
每一项是分子除以分母，累加后需加1再乘以2得到 π

第2步：避免直接计算大阶乘

直接计算 n! 会导致数值溢出。更高效的方法是利用上一项的结果递推当前项。设上一项结果为 jie_guo_prev，则当前项 jie_guo = jie_guo_prev * i / (2*i + 1)。但在本题给出的参考实现中，为了清晰展示逻辑，分别维护了分子 fen_zi 和分母 fen_mu。

第3步：完整代码实现

以下代码实现了根据输入阈值动态计算 π 的功能，输出保留小数点后6位：

#include <stdio.h>
int main()
{
    double a; // 输入的阈值
    scanf("%lf", &a);
    
    int i = 1;
    double fen_zi = 1, fen_mu = 1; // 分子分母初值
    double jie_guo = 1, sum = 0;
    
    // 当当前项 jie_guo 仍大于等于阈值时，继续累加
    while(a < jie_guo){
        fen_zi *= i;           // 更新分子: i!
        fen_mu *= (2*i+1);     // 更新分母: 3*5*...*(2i+1)
        jie_guo = fen_zi / fen_mu; // 计算当前项值
        sum += jie_guo;        // 累加
        i += 1;
    }
    
    // 公式推导: sum + 1 对应 π/2，故 π = (sum + 1) * 2
    double s = (sum + 1) * 2;
    printf("%.6f", s);
    return 0;
}