7-183 计算圆周率

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

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

输入格式：

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

输出格式：

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

输入样例：

0.01

输出样例：

3.132157

代码长度限制

16 KB

时间限制

400 ms

内存限制

64 MB

#include <stdio.h>
double fact(int n);
int main(void)
//如果不需要从命令行中获取参数，请用int main(void)；否则请用int main(int argc, char *argv[])
{
	double eps;
	scanf("%le", &eps);

	int i;
	double sum=0.0,sub=1.0;
	for(i=0;;i++){
		sub*=2*i+1;
		sum+=fact(i)/sub;
		if(fact(i)/sub<eps){
			printf("%.6f",sum*2);	// \n
			break;
		}
	}
	
	return 0;
}

double fact(int n)
{
	int i;
	double factorial=1.0;
	if(n==0){	//为什么此处的if不能省略? 
		return 1.0;
	}
	for(i=2;i<=n;i++){
		factorial*=i;
	}
	return factorial;
}