根据下面关系式,求圆周率的值,直到最后一项的值小于给定阈值。
2π=1+31+3×52!+3×5×73!+⋯+3×5×7×⋯×(2n+1)n!+⋯
输入格式:
输入在一行中给出小于1的阈值。
输出格式:
在一行中输出满足阈值条件的近似圆周率,输出到小数点后6位。
输入样例:
0.01
输出样例:
3.132157
Note : ( 1 ) for the formula is Pi half, asking for is pi, and ( 2 ) observation of odd denominator to know the product, numerator is the factorial of an integer.
// date:2020/3/6
// author:xiezhg5
#include <stdio.h>
int main(void)
{
int i=1; //控制变量
double a,j=1.0,k=1.0,sum=1.0; //初始化赋初值
scanf("%lf",&a);
while((k/j)>a) //循环条件,最后一项大于阈值
{
j=j*(2*i+1); //分子为奇数成绩
k=k*i; //分母为阶乘
i++; //自增循环
sum=sum+(k/j); //计算圆周率的一半
}
printf("%.6lf",2.0*sum); //不要忘记乘以2
return 0;
}