C语言----大整数阶乘代码详解

输入：一个整数n

输出：n!

注意：n! = n × (n-1)!，并且 0! = 1

说明：我们知道阶乘就是从1开始一直乘到它本身的一串式子，小的阶乘我们还可以用递归或者循环解决，如下列程序：

【1】递归实现

#include <stdio.h>

// 递归实现阶乘函数
int factorial_recursive(int n) {
    if (n == 0) {
        return 1; // 0 的阶乘定义为 1
    } else {
        return n * factorial_recursive(n - 1); // n! = n * (n-1)!
    }
}

int main() {
    int n;
    scanf("%d", &n);

    if (number < 0) {
        printf("阶乘只适用于非负整数。\n");
    } else {
        int res = factorial_recursive(number);
        printf("%d\n", res);
    }

    return 0;
}

终端输出结果：

【2】循环实现

#include <stdio.h>

// 循环实现阶乘函数
int factorial_iterative(int n) {
    int result = 1; // 使用 int 处理较小阶乘
    for (int i = 1; i <= n; i++) {
        result *= i; // result = result * i
    }
    return result;
}

int main() {
    int number;
    scanf("%d", &number);

    if (number < 0) {
        printf("阶乘只适用于非负整数。\n");
    } else {
        unsigned long long result = factorial_iterative(number);
        printf("%d\n", result);
    }

    return 0;
}

其实熟悉的小伙伴就知道一切循环都可以用递归来表示哈哈

但上面的只是处理较小阶乘的程序时可以用，我们发现，阶乘的增长速度非常快，这个速度远远高于指数级增长，就像n到9的时候，9！=362880.

然而我们的整型或者浮点型都是比较小的，远远不能实现大阶乘计算，所以这里我们运用了数组的结构老保存每一位阶乘结果，是不是很神奇。

程序如下：

#include <stdio.h>

int main() {
    int a[10005], i, j, l = 0, s, n, cnt = 0; //数组a可以保存最大位数为10005的阶乘结果
    scanf("%d", &n);
    a[0] = 1;  //阶乘都是从1开始相乘

    for(i = 1; i <= n; i++) {  //循环模拟阶乘的计算过程
        s = 0; //存储计算的中间结果
        for(j = 0; j <= l; j++) {
            s += a[j]*i;
            a[j] = s%10;
            s /= 10;
        }
        while(s != 0) {
            l++;  //进一位
            a[l] = s%10;
            s /= 10;
        }
    }
    for(i = l; i >= 0; i--) { //倒序输出结果
        cnt++;
        printf("%d", a[i]);
        if(cnt%5 == 0) printf(" ");
        if(cnt%25 == 0) printf("\n");
    }

    return 0;
}

终端输出结果：

简化一下：

#include <stdio.h>

int main() {
    int n, a[10005], l = 0, s;
    scanf(" %d", &n);
    a[0] = 1;

    for(int i = 1; i <= n; i++) {
        s = 0;
        for(int j = 0; j <= l; j++, s /= 10) {
            s += a[j]*i;
            a[j] = s%10;
        }
        for(; s; s /= 10) a[++l] = s%10;
    }

    for(int i = l, cnt = 1; i >= 0; i--, cnt++) {
        printf("%d", a[i]);
        if(cnt%5 == 0) printf(" ");
        if(cnt%25 == 0) printf("\n");
    }
    return 0;
}

终端结果：