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