根据下面关系式,求圆周率的值,直到最后一项的值小于给定阈值。
2π=1+31+3×52!+3×5×73!+⋯+3×5×7×⋯×(2n+1)n!+⋯
输入格式:
输入在一行中给出小于1的阈值。
输出格式:
在一行中输出满足阈值条件的近似圆周率,输出到小数点后6位。
输入样例:
0.01
输出样例:
3.132157
#include<stdio.h>
double my(int x)
{
if (x == 0)
return 1;
return x * my(x - 1);
}
double My(int x)
{
if (x == 0)
return 1;
return (2 * x + 1) * My(x - 1);
}
int main()
{
double k = 1;
double n = 0;
scanf("%lf", &n);
double pi = 0;
int i = 0;
double sum = 1.0;
while (k>n)
{
i++;
k = my(i) / My(i);
sum += k;
}
printf("%lf", 2*sum);
return 0;
}