求大整数的阶乘（hdu1042）

呵呵！刚刚在chinaunix注册了，这是第一篇日志，大部分借鉴的是hhccbb的算法。写的不好的，希望大牛们谅解。现在转入正题。

hdu1042要求求n！，而n的值最大到10000。属于规模型解题范围，10000！的位数为35670。

（一）首先，了解n！的阶乘范围：

   由于int有符号型可表示的最大数为2^31-1，为2147483647。位数为10位。13！＝6227020800，就已经超过了int型表示范围。

   所以，需要借鉴更大存储空间的类型，可以选用__int64类型。或者利用数组存储，联系高精度算法的设计。

（二）n！（n在两位数时）的基本求法：

   从网上借鉴的算法，要求n！，可以用数组a[]来存放n！的每位数字。例如：4！＝24，可以用a[1]=4,a[2]=2,a[0]=2(a[0]存放位数)；5！＝120，可以用a[1]=0,a[2]=2,a[3]=1,a[0]=3;

   故算法初始化可以使a[0]=1,a[1]=1来表示初值1。然后利用循环，数组a[]依次乘以2,3,4,....,n。n！的低位存放在a[]的前部，高位往后写。

   自己写的一个函数如下：

void fac(int n,short *a)//n!的求法  200!---375位
{//a[0]存放n!的位数，a[1],...,a[a[0]]逆序存放每位数字。
//若于高精度连用，注意将a[1],...,a[a[0]]转置。
/*for(i=1,k=a[0];i<k;i++,k--)  {temp=a[i];a[i]=a[k];a[k]=temp;}*/
 int i,j,k=1;
  a[1]=1;
  for(i=2;i<maxsize;i++) a[i]=0;
  for(i=2;i<=n;i++)
  {
   for(j=1;j<=k;j++) a[j]*=i;
   for(j=1;j<=k;j++)
   {
  if(a[j]>=10)
  {if(j==k) k++;//k计位数。
   a[j+1]+=a[j]/10;
   a[j]%=10;
  }
   }//for j
  }//for i
  a[0]=k;
}

说明：（1）a[0]存放n！的位数；（2）a[1],a[2],...,a[a[0]]逆序存放n!各位数字。

（三）n！（n很大时，处于三四位时）的求解：

    假设要求10000！（35670位数），若采用标题（二）的算法，利用每个数组元素来存放每位数字的话，显然造成了空间的浪费，而循环也使得算法（二）对10000!不适合。

    所以可以考虑用每个数组元素来存放n！的每个10位数字，只需要把算法（二）的a[]类型设置为__int64型，将进制扩大为10^10。

    hdu1042的代码设计如下：

#include<stdio.h>
//#include<math.h>
#define I 10000000000
#define maxsize 35670

//10000! 位数：35660

__int64 a[maxsize];

void fac2(int n,__int64 *a)
{int i,j;__int64 k;
 if(n==0||n==1)
 {a[0]=1;a[1]=1;return;
 }
 a[0]=1;a[1]=1;k=1;
 for(i=2;i<maxsize;i++) a[i]=0;
 for(i=2;i<=n;i++)
 {
  for(j=1;j<=k;j++) if(a[j]) a[j]*=i;
  for(j=1;j<=k;j++)
  {
   if(a[j]>=I)
   {
   if(j==k) k++;
   a[j+1]+=a[j]/I;
   a[j]%=I;
   }
  }//for j
 }//for i
 a[0]=k;
}

int main()
{int n;__int64 i,j;
 while(scanf("%d",&n)!=EOF)
//for(n=3;n<=50;n++)
 {fac2(n,a);
 //printf("%d!=/n",n);
  printf("%I64d",a[a[0]]);

  for(i=a[0]-1;i>=1;i--)
  {
    //若a[i]不足9位数，前面补零。
   if(a[i]<10) printf("000000000");
   else if(a[i]<100) printf("00000000");
   else if(a[i]<1000) printf("0000000");
   else if(a[i]<10000) printf("000000");
   else if(a[i]<100000) printf("00000");
   else if(a[i]<1000000) printf("0000");
   else if(a[i]<10000000) printf("000");
   else if(a[i]<100000000) printf("00");
   else if(a[i]<1000000000) printf("0");
   printf("%I64d",a[i]);
  }
  printf("/n");
 }
return 1;
}注意n!的输出，由于除a[a[0]]外前面的每个数组元素位数不一定为10位。所以这些元素输出时，需要前面补零凑够10位。

   if(a[i]<10) printf("000000000");
   else if(a[i]<100) printf("00000000");
   else if(a[i]<1000) printf("0000000");
   else if(a[i]<10000) printf("000000");
   else if(a[i]<100000) printf("00000");
   else if(a[i]<1000000) printf("0000");
   else if(a[i]<10000000) printf("000");
   else if(a[i]<100000000) printf("00");
   else if(a[i]<1000000000) printf("0");

此部分的设计，也可以先求出a[i]的位数，然后根据位数补足10位。

 for(j=1;j<=10-digit(a[i]);j++) printf("0");
   printf("%I64d",a[i]);//若a[i]不足10位数，前面补零。

digit()定义如下：

int digit(__int64 a)
{int num=0;
 while(a) {a/=10;num++;}
 return num;
}