大数阶乘

大数阶乘。

阶乘攀升特别快，在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中嵌入汇编来提高速度，没试过。