#include "Stdio.h"
#include "Conio.h"
int main(void)
{
    int nums[10000]; /*nums數組越大,且改變為長整型,可計算更大的階乘,但改變數組類型會變慢*/
    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; /*這里有改進的余地*/
    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); /*進位計算*/
            nums[j]=nums[j]%10;
        }
        while(nums[k]>=10)
        {
            nums[k+1]=nums[k+1]+(int)(nums[k]/10); /*把結果數字的最高位分別放進nums數組中對應的位中,也可改進*/
            nums[k]=nums[k]%10;
            k++;
        }
    }
    printf("\nn!=");
    for(i=k;i>=0;i--) printf("%d",nums[i]); /*輸出結果*/
    printf("\nk=%d",k+1);  /*k+1是最后結果的位數*/
    getch();
    return 0;
}

這個小程序根據整數乘法的原理，可突破計算機本身及編程語言存儲類型大小限制，計算高精度的階乘,如果nums數組大小為10000，大概可以計算3000的階乘，即3000!，3000!計算的最后結果數字是9155位，其實這個程序還有很多可以改進的地方，初用C語言，沒有把C語言的一些優勢發揮出來。還有，我在VB里寫好后運行測試發現和C語言比，在這算法一樣的情況下，效率是天差地別的！
本程序在TC2.0下編譯通過！
（還有，標顏色標得我手軟了！）

相關文章