n的阶乘高精度算法解析.docVIP

n的阶乘高精度算法解析 在n的阶乘中，n大于16时n的阶乘得出的结果用长整形保存时总是溢出，这对于刚刚学习c的学生来说总是感觉很棘手，于是就上百度，在搜索栏打上“阶乘”两个字，在搜索结果里面，有“百度百科”的解释，打开网页，使劲往下翻，就看到了关于“n的阶乘高精度运算”的c代码： #include stdio.h #define N 5000 // modify it to hold larger number int main(void) { int n, i, j, s, up, f[N] = {0}; scanf(%d, n); for (i = 2, f[0] = 1; i = n; i++) for (j = up = 0; j N; j++) { s = f[j] * i + up; f[j] = s % 10; up = s / 10; } for (i = N-1; f[i] == 0; i--) ; for (; i = 0; i--) printf(%d, f[i]) ; printf(\n); return 0; } 若是感觉不是很好明白，可以改成以下代码，可能就会好明白一点了： #include stdio.h #define N 5000 // modify it to hold larger number int main(void) { int n, i, j, s, up, f[N] = {0}; scanf(%d, n); f[0] = 1; for (i = 2; i = n; i++)//外循环控制阶乘的数 { up=0; for (j=0; j N; j++) //内循环用于得出前一个数的阶乘结果乘以当前数i得出的阶乘结果 { s = f[j] * i + up; f[j] = s % 10; up = s / 10; } } for (i = N-1; f[i] == 0; i--) ; for (; i = 0; i--) printf(%d, f[i]) ; printf(\n); return 0; } 当我们第一次接触这个代码的时候，你会不会感到很无奈呢？是不是不明白为什么这么做啊？ 下面说一说我对这串代码的理解： 首先，要想要理解这串代码是怎么回事儿，最好的方法就是先找一个数带进去，找一张好的白纸，好好写一写，画一画，从0写到5。从更本质上来说，只把5带进去就行，因为 0！*1=1！； 1！*2=2!; 2!*3=3!; 3!*4=4!; 4!*5=5!; 但是3！=610;4!=24刚刚大于10,取5就更具有普遍性。 下面就讨论当n=5时的情况。 因为外循环是从i=2开始的所以我们也从i=2开始讨论。 内循环中对于for（j=0,up=0;jn;j++） 一：i=2时 1. j=0时：s=f[j]*i+up;-----s=f[0]*2+0=1*2+0=2； f[j]=s%10;-----f[0]=2%10=2; up=s/10;-----up=2/10=0; 2. j=1时：s=f[j]*i+up;-----s=f[1]*2+0=0*2+0=0; f[j]=s%10;-----f[1]=0%10=0; up=s/10;-----up=0/10=0; 当j大于等于2时的情况和“2.”时的情况相同，此时得到的结果为：f[0]=2;f[1]=0;f[2]=0;..... 二：i=3时 1.j=0时：s=f[j]*i+up;-----s=f[0]*3+0=2*3+0=6; f[j]=s%10;-----f[0]=6%10=6; up=s/10;-----up=6/10=0; 2.j=1时：s=f[j]*i+up;-----s=f[1]*3+0=0*3+0=0; f[j]=s%10;-----f[1]=0%10=0; up=s/10;-----up=0/10=0; 当j大于等于2时的情况和“2.”时的情况相同，此时得到的结果为：f[0]=6;f[1]=0;f[2]=0;...... 三：i=4时 1. j=0时：s=f[j]*i+up;-----s=f[0]*4+0=6*4+0=24; f[j]=s%10;-----f[0]=24%10=4; up=s/10;-----up=24/10=2; 2. j=1时:s=f[j]*i+up;-----s=f[1]*4+2=0*4+2=2; f[j]=s%10;-----f[1]=2%10=2; up=s/10;-----up=2/10=0; 3. j=2时：s=f[j]*i+up;-----s=f[2]*4+0=0*4+0=0; f[j]=s%10;-----f[2]=0%10=0; up=s/10;-----up=0/10=0;