專注電子技術(shù)學(xué)習(xí)與研究
當(dāng)前位置:單片機(jī)教程網(wǎng) >> MCU設(shè)計(jì)實(shí)例 >> 瀏覽文章

高精度階乘小程序

作者:佚名   來(lái)源:本站原創(chuàng)   點(diǎn)擊數(shù):  更新時(shí)間:2009年12月04日   【字體:

#include "Stdio.h"
#include "Conio.h"
int main(void)
{
    int nums[10000]; /*nums數(shù)組越大,且改變?yōu)殚L(zhǎng)整型,可計(jì)算更大的階乘,但改變數(shù)組類型會(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)位計(jì)算*/
            nums[j]=nums[j]%10;
        }
        while(nums[k]>=10)
        {
            nums[k+1]=nums[k+1]+(int)(nums[k]/10); /*把結(jié)果數(shù)字的最高位分別放進(jìn)nums數(shù)組中對(duì)應(yīng)的位中,也可改進(jìn)*/
            nums[k]=nums[k]%10;
            k++;
        }
    }
    printf("\nn!=");
    for(i=k;i>=0;i--) printf("%d",nums[i]); /*輸出結(jié)果*/
    printf("\nk=%d",k+1);  /*k+1是最后結(jié)果的位數(shù)*/
    getch();
    return 0;
}

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

關(guān)閉窗口

相關(guān)文章