//大数阶乘// 07.4.27 #include<stdio.h>#include<time.h> int main(){    double duration;     int a[100000]={0};    int N,i,c=1,j;     a[1]=1;    clock_t start, finish;     printf("大数阶乘,请输入一个数\n");    scanf("%d",&N);    start = clock();    for (j=1;j<=N;j++)    {        for (i=1;i<=c;i++)        {            a[i]=a[i]*j;        }         for (i=1;i<=c;i++)          {            /*不能写成(i=c;i>0;i--)因为这样将漏掉低位的进位,              例如a[3]=7*7,a[2]=3*3,a[1]=2*5,当i减到1时a[1]=10需要向a[2]进一位，而当仅为后              a[2]=10又该向a[3]进位，但此时循环已是最后一次循环，故不能实现向a[3]进位              a[4]=4 a[3]=9 a[2]=10 a[1]=0            */            if (a[i]>=10)            {                 a[i+1]=a[i+1]+a[i]/10;                a[i]=a[i]%10;                if(a[c+1]!=0)                 c++;            }        }     //若将if(a[c+1]!=0)语句放于此处(for (i=1;i<=c;i++))之外则造成a[c+1]不是一位数,     //因为a[c+1]这个元素将不止一位数，当进行例如1000!的计算时     //会由于a[c+1]数值的溢出而造成计算不准确。    }    printf("result:\n");     for(i=c;i>0;i--) {  printf("%d",a[i]);    }    finish = clock();    printf("\n一共%d位\n",c);     duration = (double)(finish - start)/CLOCKS_PER_SEC ;    printf( "%f ms\n", duration );     getchar();    getchar();    return 0;}  

评论