计算圆周率

根据下面关系式，求圆周率的值，直到最后一项的值小于给定阈值。

​2​​π​​=1+​3​​1​​+​3×5​​2!​​+​3×5×7​​3!​​+⋯+​3×5×7×⋯×(2n+1)​​n!​​+⋯

输入格式：

输入在一行中给出小于1的阈值。

输出格式：

在一行中输出满足阈值条件的近似圆周率，输出到小数点后6位。

输入样例：

0.01

输出样例：

3.132157

 

#include<iostream>
#include<iomanip>
using namespace std;

int main()
{
    float x,a[1000],sum,y;
    int i(1);
    a[0]=1;
    a[1]=1.0/3.0;
    sum=a[0]+a[1];
    cin>>x;
    while(a[i] >= x)
    {
        i++;
        a[i]=a[i-1]*i/(2*i+1);
        sum+=a[i];
    }

    cout<<fixed<<setprecision(6)<<2*sum<<endl;
    return 0;
}