大数阶乘

大数阶乘。

阶乘攀升特别快，在C中用int,double等，能计算的范围相当有限；但我们对一个数的表示，按照通常的方法，例如123；如果以10为基，则不同的只是每位所成基数的幂不同。

于是，想到使用一个数组保存数字，而基数可以任意。在计算中，应注意对基数的处理。

 

#include <stdio.h>#include <stdlib.h>#include <time.h>#define BUFFSIZE 50000#define BASE 10000intmain(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);}//ifif((fw=fopen("output.txt","w"))==NULL){perror("input.txt");exit(EXIT_FAILURE);}//iffscanf(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;}}//whilelen++;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中嵌入汇编来提高速度，没试过。