double factorial(int x) {
if (x >= 1)
return x * factorial(x - 1);
else
return 1;
}
当我尝试得到 1000000 的阶乘时，它会导致分段错误。请任何人都可以帮助我吗？？
【问题讨论】：
标签： c factorial abstract-data-type
添加一个答案，它不会给出该阶乘的精确值，而是给出该结果数的10的幂。
为此，让我们使用 y = log10(x) 给出 x = 10y >y。
例如 log10(1000000) 是 6（6 个零，106）。 log10(9500) 约为 3.98，因为我们接近 104。
那么log(a * b) 是log(a) + log(b)，这在我们的例子中非常有用，因为fact(x) 是x*(x-1)*(x-2)*...*1)。
现在正如@dbush 所说，您不能递归循环 1M 次，因为保存返回地址（和局部变量）的堆栈将很快溢出。所以让我们做一个循环来代替
double d,pow10 = 0.0;
for(d=1 ; d<=1000000 ; d++) pow10 += log10(d);
printf("Result is 10^ %.02lf\n", pow10);
做到了，结果是~105565708.92，一个有~5565709位数的数字...
一个不需要循环的封闭公式，可以提供出色的阶乘近似。再次计算将达到限制，所以做 log10，有 log(ab) = b*log(a)
1000000*log(1000000) - 1000000*log(e) + log(sqrt((2000000+1/3)*PI))
也给出 ~5565708.92，阶乘为 ~105565708.92
【讨论】：
您不能让函数递归调用自身 1000000 次。这很容易溢出堆栈，导致段错误。
您需要切换到迭代解决方案：
double factorial(int x)
{
double result = 1;
while (x >= 1) {
result *= x;
x--;
}
return result;
}
另请注意，您没有足够的精度来存储此值的所有数字。为此，您需要使用像 GMP 这样的大数字库。
【讨论】：
double 是否有意义，或者例如unsigned long long比double好吗？
unsigned long long 可能有 64。double 的指数范围仅上升到 1024。即使 IEEE 四元精度也仅上升到 16384。没有办法表示 1000000 的近似值！没有具有巨大指数范围的自定义软浮点类型，并且没有很多 GB 内存就无法表示确切的值。
可以迭代计算任意大的阶乘，而无需使用库来处理大量数字。 我用 Java 编写了一个解决方案（针对 ProjectEuler 问题），我已经针对谷歌验证了高达 1000！，并运行了 5000！在 5 年前的桌面上不到一分钟（可能很慢，因为它正在控制台上显示输出。）。 不到 100 行代码，包括 cmets 和打印输出。 （不共享代码，因为这是 ProjectEuler 协议的一部分）。
【讨论】：