近似求PI

来源:互联网 发布:剑三女神花姐捏脸数据 编辑:程序博客网 时间:2024/05/13 06:08
5-10 近似求PI   (15分)

本题要求编写程序,根据下式求π\piπ的近似值,直到最后一项小于给定精度eps。

π2=1+1!3+2!3×5+3!3×5×7+⋯+i!3×5×⋯×(2×i+1)+⋯\frac{\pi}{2}=1+\frac{1!}{3}+\frac{2!}{3\times 5} + \frac{3!}{3\times 5\times 7}+ \cdots +\frac{i!}{3\times 5\times \cdots \times (2\times i+1)}+ \cdots2π=1+31!+3×52!+3×5×73!++3×5××(2×i+1)i!+

输入格式:

输入在一行中给出精度eps,可以使用以下语句来读输入:

scanf("%le", &eps);

输出格式:

在一行内,按照以下格式输出π\piπ的近似值(保留小数点后5位):

PI = 近似值

输入样例:

1E-5

输出样例:

PI = 3.14158

    分析“scanf("%le", &eps);这句话太重要了,告诉了我们输入的1E-5本身就是一个双浮点类型的数,直接输入即可,之后用两个函数计算判断出他的精度是否小于eps即可,注意的是为了保险,把所有的函数数据类型都设定成了double,int的话有很大的危险。

    #include <bits/stdc++.h>using namespace std;double jc(int q){    double ans=1;    for(int i=1; i<=q; i++)    {        ans*=i;    }    return ans;}double jc2(int q){    double ans=1;    for(int i=3; i<=2*q+1; i+=2)    {        ans*=i;    }    return ans;}int main(){    int q=1;    double eps;    scanf("%lf",&eps);    double ans=1.0;    double temp=1;    while(temp>=eps)    {        temp=jc(q)*1.0/jc2(q);;        ans+=temp;        q++;    }    printf("PI = %.5lf",ans*2);    return 0;}

      0 0
      原创粉丝点击