PI的求法

1671年(J.Gregory)

π/4=1-1/3+1/5-1/7+1/9-1/11+....


#include <stdio.h>
#include <string.h>
#include <stdlib.h>

int main()
{
    double pi = 0.0;
    int flag, i;
    int n;

    while(~scanf("%d", &n))
    {
        pi = 1.0;
        flag = 0;
        for( i = 3; i < n; i += 2 )
        {
            pi += (flag == 0 ? -1.0 : 1.0)/i;
            flag = !flag;
        }

        printf("%lf\n", 4.0*pi);
    }

    return 0;
}



欧拉(Leonhard Euler)

π^2/8=1+1/3^2+1/5^2+1/7^2+.....


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

int main()
{
    double pi = 0.0;
    int n, i;

    while(~scanf("%d", &n))
    {
        pi = 1.0;
        for( i = 3; i <= n; i += 2 )
        {
            pi += 1.0/(i*i);
        }
        pi *= 8;
        pi = sqrt( pi );
        printf("%lf\n", pi);
    }

    return 0;
}

在Python中,有多种公式法可用于计算圆周率&pi;,以下是几种常见的方法及其代码示例: ### 公式法1:莱布尼兹公式 莱布尼兹公式为:$ \pi = 4\times(1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + ...) $,该公式简单而优美,但收敛速度较慢。 ```python s = 0 n = 1000000 # 迭代次数 for i in range(n): term = (-1) ** i / (2 * i + 1) s += term pi = 4 * s print(f"使用莱布尼兹公式计算的圆周率近似值: {pi}") ``` ### 公式法2:示例代码中的公式法 根据引用[1]中的公式法代码,其通过特定的公式循环计算圆周率。 ```python i = 1 j = 1 s = 0 print('******圆周率公式法******') for i in range(1, 100): print(f'循环{i}次,第一次计算{j=}') x = (1 / (2 ** j) + 1 / (3 ** j)) * (1 / j) print(f'{x=}') j = j + 2 print(f'循环{i}次,第二次计算{j=}') y = -((1 / (2 ** j) + 1 / (3 ** j)) * (1 / j)) print(f'{y=}') j = j + 2 print(f'循环{i}次计算完成后,{j=}') z = x + y s = s + z print(f'{s=}') print('————————————————————') print(f'圆周率的计算结果如下:{4 * s}') ``` ### 公式法3:BBP公式 BBP公式由David Bailey、Peter Borwein和Simon Plouffe于1995年共同发表,它打破了传统的圆周率算法,可以计算圆周率的任意第n位,而不用计算前面的n - 1位,为圆周率的分布式计算提供了可行性[^3]。 ```python def bbp(n): pi = 0 for k in range(n): pi += (1 / 16 ** k) * ( 4 / (8 * k + 1) - 2 / (8 * k + 4) - 1 / (8 * k + 5) - 1 / (8 * k + 6) ) return pi n = 100 # 迭代次数 result = bbp(n) print(f"使用BBP公式计算的圆周率近似值: {result}") ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值