高精度階乘小程序

#include "Stdio.h"
#include "Conio.h"
int main(void)
{
int nums[10000]; /*nums數組越大,且改變?yōu)殚L(cháng)整型,可計算更大的階乘,但改變數組類(lèi)型會(huì )變慢*/
int k,j;
int n,i;
printf("Please input a integer n:");
scanf("%d",n);
for(i=1;i=n-1;i++) nums[i]=0; /*這里有改進(jìn)的余地*/
nums[0]=1;
k=0;
for(i=2;i=n;i++)
{
for(j=0;j=k;j++) nums[j]=nums[j]*i;
for(j=0;j=k-1;j++)
{
nums[j+1]=nums[j+1]+(int)(nums[j]/10); /*進(jìn)位計算*/
nums[j]=nums[j]%10;
}
while(nums[k]>=10)
{
nums[k+1]=nums[k+1]+(int)(nums[k]/10); /*把結果數字的最高位分別放進(jìn)nums數組中對應的位中,也可改進(jìn)*/
nums[k]=nums[k]%10;
k++;
}
}
printf("n!=");
for(i=k;i>=0;i--) printf("%d",nums[i]); /*輸出結果*/
printf("k=%d",k+1); /*k+1是最后結果的位數*/
getch();
return 0;
}

這個(gè)小程序根據整數乘法的原理，可突破計算機本身及編程語(yǔ)言存儲類(lèi)型大小限制，計算高精度的階乘,如果nums數組大小為10000，大概可以計算3000的階乘，即3000!，3000!計算的最后結果數字是9155位，其實(shí)這個(gè)程序還有很多可以改進(jìn)的地方，初用C語(yǔ)言，沒(méi)有把C語(yǔ)言的一些優(yōu)勢發(fā)揮出來(lái)。還有，我在VB里寫(xiě)好后運行測試發(fā)現和C語(yǔ)言比，在這算法一樣的情況下，效率是天差地別的！
本程序在TC2.0下編譯通過(guò)！