6-10 阶乘计算升级版 (数组模拟详细讲解)

本题要求实现一个打印非负整数阶乘的函数。

函数接口定义：

void Print_Factorial ( const int N );

其中N是用户传入的参数，其值不超过1000。如果N是非负整数，则该函数必须在一行中打印出N!的值，否则打印“Invalid input”。

裁判测试程序样例：

#include <stdio.h>

void Print_Factorial ( const int N );

int main()
{
    int N;
    
    scanf("%d", &N);
    Print_Factorial(N);
    return 0;
}

/* 你的代码将被嵌在这里 */

---

输入样例：

15

输出样例：

1307674368000

---

输入样例1:

1000

输出样例1:

402387260077093773543702433923003985719374864210714632543799910429938512398629020592044208486969404800479988610197196058
631666872994808558901323829669944590997424504087073759918823627727188732519779505950995276120874975462497043601418278094
646496291056393887437886487337119181045825783647849977012476632889835955735432513185323958463075557409114262417474349347
553428646576611667797396668820291207379143853719588249808126867838374559731746136085379534524221586593201928090878297308
431392844403281231558611036976801357304216168747609675871348312025478589320767169132448426236131412508780208000261683151
027341827977704784635868170164365024153691398281264810213092761244896359928705114964975419909342221566832572080821333186
116811553615836546984046708975602900950537616475847728421889679646244945160765353408198901385442487984959953319101723355
556602139450399736280750137837615307127761926849034352625200015888535147331611702103968175921510907788019393178114194545
257223865541461062892187960223838971476088506276862967146674697562911234082439208160153780889893964518263243671616762179
168909779911903754031274622289988005195444414282012187361745992642956581746628302955570299024324153181617210465832036786
906117260158783520751516284225540265170483304226143974286933061690897968482590125458327168226458066526769958652682272807
075781391858178889652208164348344825993266043367660176999612831860788386150279465955131156552036093988180612138558600301
435694527224206344631797460594682573103790084024432438465657245014402821885252470935190620929023136493273497565513958720
559654228749774011413346962715422845862377387538230483865688976461927383814900140767310446640259899490222221765904339901
886018566526485061799702356193897017860040811889729918311021171229845901641921068884387121855646124960798722908519296819
372388642614839657382291123125024186649353143970137428531926649875337218940694281434118520158014123344828015051399694290
153483077644569099073152433278288269864602789864321139083506217095002597389863554277196742822248757586765752344220207573
630569498825087968928162753848863396909959826280956121450994871701244516461260379029309120889086942028510640182154399457
156805941872748998094254742173582401063677404595741785160829230135358081840096996372524230560855903700624271243416909004
153690105933983835777939410970027753472000000000000000000000000000000000000000000000000000000000000000000000000000000000
000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
000000000000000000000000000000000000000000000000

---

由此显而易见,哪怕是最长的 unsigned long long 也无法满足存储需求 为此我们只能使用 "数组模拟法"  进行模拟运算 核心思想如下

num = 123456789          // 被模拟的数字                                           
nums[2022],length        // 模拟使用的数组,以及当前其长度
length = 0               // 记录nums模拟的长度 : length-1

while num>0:
    nums[length++] = num%10;
    num/=10;

执行完上述伪代码后 123456789被反向存入nums nums[0]=9,num[1]=[8]......nums[length-1] = 1 // nums[0] 是个位,nums[1]是十位 以此类推 -1 如果我们让其模拟的每个位数乘11 nums[0~length-1] *= 11; -2 经过步骤1的处理,此时nums每个位数的数字为 nums[0]=99,nums[1]=88,......nums[length-1] = 11 但是我们知道,我们用的可是 石 进制~ 这样显然是不合理的 so! z = 0~length : nums[z+1] = nums[z]/10 ,nums[z] %= 10 -3 步骤2处理后 nums[0]=9,nums[1]=(88+9)%10=7......nums[length-1]=(11+2)%10=3 ... ??? 诶,我们发现这样一乘之后,length的实际长度是应该要变大呀 那么我们是不是只要判断 if nums[length]>0: length+=1 就可以了? 上述样例似乎真的可以 但是 ....... 如果我要乘111的话呢 此时 nums[length-1] = (111+22)%10 那么 nums[length] = 0+13  由此可见 length 只加 1 肯定是不够的 那该怎么判断 length 加了多少 ? ... ... ... 思考完再看下面代码 下方代码对步骤 1与2 进行了合并 ~

---

void Print_Factorial ( const int N ) {
    if (N < 0) {
        printf("Invalid input");
        return;
    }
    int nums[10001] = {0,1}, length = 2; // 模拟数组

    for(int z=2;z<=N;z++){
        for(int z1=1;z1<=length;z1++)
        {
            nums[z1] = nums[z1]*z+nums[z1-1]/10;
            nums[z1-1]%=10;
        }
        while (nums[length]>0){
            nums[length+1] = nums[length]/10;
            nums[length]%=10;
            length++;
        }

    }

    for (int i = length-1; i >0 ; i--) printf("%d",nums[i]);

}

---