根据下面关系式,求圆周率的值,直到最后一项的值小于给定阈值。
2π=1+31+3×52!+3×5×73!+⋯+3×5×7×⋯×(2n+1)n!+⋯
输入格式:
输入在一行中给出小于1的阈值。
输出格式:
在一行中输出满足阈值条件的近似圆周率,输出到小数点后6位。
输入样例:
0.01
输出样例:
3.132157
a=input()
a=float(a)
b=1
nm=1
n=1
m=1
sum=0
if(a<1):
while(a<b):
if(nm>1):
n=n*(nm-1)
m=m*(2*nm-1)
else:
n=1
m=1
sum+=b
nm+=1
b=n/m
print("{:.6f}".format((sum-1+b)*2))