近似求PI

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)}+ \cdots​2​​π​​=1+​3​​1!​​+​3×5​​2!​​+​3×5×7​​3!​​+⋯+​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;
}