7-15 计算圆周率 (15 分)

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

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

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

输入样例：

0.01

结尾无空行
输出样例：

3.132157

第一次写的代码可以运行，但是无法通过，提示直接计算阶乘会超出长整型

#include <stdio.h>
int funup(int n);
int fundown(int n);
int main ()
{
    int i = 0;
    double prn;
    double result = 0,iteam = 1;
    scanf("%lf",&prn);
    while(iteam>prn)
    {
        iteam = 1.0*funup(i)/fundown(i);
        result += iteam;
        i++;
    }
    printf("%f",2*result);
    return 0;
}
int funup(int n)
{
    int result = 1;
    if(n == 0)
    {
        result = 1;
    }
    else
    {
        for(int i = 1;i<=n;i++)
        {
            result *= i;
        }
    }
    return result;
}
int fundown(int n)
{
    int result = 1;
    for(int i =0;i<=n;i++)
    {
        result *= (2*i+1);
    }
    return result;
}

后来参考了别人的代码写的就通过了

#include <stdio.h>

int main()
{
    double up = 1,down = 1,item = 1,t = 0,sum = 1;
    int i = 1;
    scanf("%lf",&t);
    do
    {
        up *= i;
        down *= (2*i+1);
        item = 1.0*up/down;
        sum += item;
        i++;
    }while(item>t);
    printf("%f",sum*2);
    return 0;
}

有点不明白，为什么第一次单独计算分子的阶乘就说超出范围了？第二次“up *= i;”不也是算出来阶乘吗？

还有用递归的方法，更巧妙了

#include<bits/stdc++.h>
using namespace std;
double getnum(int x){
	if(x == 1) return 1;
	else return getnum(x - 1) * x;
}
double getnnum(int x){
	if(x == 1) return 1; 
	else return getnnum(x - 1) * (2 * x - 1);
}
int main(){
	double pi = 1, n = 1, m, up = 1, down = 1;
	cin >> m;
	while(up / down >= m){
		up = getnum(n);
		down = getnnum(n + 1);
		pi += up / down;
		n++;
	}
	printf("%.6lf", (double)pi * 2);
	return 0;
}

代码内容参考了这两篇文章
https://blog.csdn.net/weixin_45031646/article/details/90147530
https://blog.csdn.net/qiguanjiezl/article/details/102737041