求e的值，gcc运行结果为inf？

GeBron

用公式e = 1 + 1/1! + 1/2! + 1/3! + 1/4! + ...

编写的程序如下：

1. #include <stdio.h>
   
2. 
   
3. int main(void)
   
4. {
   
5.     long fact(int n);
   
6.     int n;
   
7.     double i,sum;
   
8.     printf("enter n:");
   
9.     scanf("%d",&n);
   
10.     for (i = 0.0,sum = 0.0; i <= n; i++)
    
11.     {
    
12.         sum = sum + 1.0/fact(i);
    
13.     }
    
14.     printf("%lf\n",sum);
    
15. }
    
16. 
    
17. long fact(int n)
    
18. {
    
19.     long fac;
    
20.     if (n == 0 || n == 1)
    
21.         fac = 1;
    
22.     else
    
23.         fac = fact(n - 1) * n;
    
24.     return fac;
    
25. }

复制代码

当n = 5时，运算结果为：2.716667

当n = 50时，就显示为：inf

不知道具体什么原因，愿指教。

先谢谢了~

KBTiller

明显的一个错误是溢出
long fact(int n)

fact(50)

溢出是未定义行为
产生任何结果都不奇怪

KBTiller

与结果无关
for (i = 0.0,sum = 0.0; i <= n; i++)
这个循环写得也不符合规范：
用double类型变量做循环变量
循环变量应该是整数类型

KBTiller

不规范的地方还有
类型不一致

fact(i);     //double i,sum;

long fact(int n)

x5miao

还有n的类型也不应该是有符号整形，不然如果用户输入负数，lz你的程序就悲剧了。。。。

x5miao

而且每次计算i的阶乘都得从1到i递归一次效率多低啊，要是每次计算i的阶乘以后能保存下来当计算（i+1）的时候再直接采用不更好啊

GeBron

回复 2# KBTiller


是这个求阶乘的函数定义有错吗？

GeBron

回复 3# KBTiller


    刚开始学C语言，是有不少不规范的。谢谢指出。

GeBron

回复 6# x5miao


    好的想法。谢谢。

KBTiller

回复 7# GeBron


    两说。在一定范围内是正确的，一旦结果超过了long的表示范围，就是错的