大數階乘

大數階乘。

階乘攀升特別快，在C中用int,double等，能計算的範圍相當有限；但我們對一個數的表示，按照通常的方法，例如123；如果以10爲基，則不同的只是每位所成基數的冪不同。

於是，想到使用一個數組保存數字，而基數可以任意。在計算中，應注意對基數的處理。

 

#include <stdio.h> #include <stdlib.h> #include <time.h> #define BUFFSIZE 50000 #define BASE 10000 int main(void) { int buff[BUFFSIZE]={0}; int N,i,k,j=0; int len=1; int add=0; FILE *fr,*fw; clock_t start,end; start=clock(); if((fr=fopen("input.txt","r"))==NULL) { perror("input.txt"); exit(EXIT_FAILURE); }//if if((fw=fopen("output.txt","w"))==NULL) { perror("input.txt"); exit(EXIT_FAILURE); }//if fscanf(fr,"%d",&N); N++; buff[0]=0; buff[1]=1; while(N-->1) { for(i=1;i<=len;i++) { k=buff[i]*N+add; add=k/BASE; buff[i]=k%BASE; } while(add>0){ buff[++len]=add%BASE; add/=BASE; } }//while len++; while(--len>0) fprintf(fw,"%d",buff[len]); fclose(fr); fclose(fw); end=clock(); printf("%lf",(((double)(end-start)/CLOCKS_PER_SEC)*1000.00)); return EXIT_SUCCESS; }

 

在測試中，我計算了10000的階乘；貌似基數越大，計算的精度越小；但基數如果取小，那麼計算的時間會增加。

在網上看到可以通過在C中嵌入彙編來提高速度，沒試過。