超大数阶乘算法

本程序采用数组存放结果,故其计算上仍受到数组容量的限制。

如以对其进行改进采用链表存放中间结果及最后结果则可使计算任意大数的阶乘，只要你的内存足够大！

#include "stdafx.h"
#include "stdio.h"
#include "malloc.h"

const MAX=1000;
void pnext(int a[],int k)
{
 int *b, m=a[0],i,j,r,carry;
    b=(int *)malloc(sizeof(int)*(m+1));

 for(i=1;i<=m;i++) b[i]=a[i];
    for(j=1;j<k;j++)
 {
  for(carry=0,i=1;i<=m;i++)
  {
   r= (i<=a[0] ? a[i]+b[i] : a[i]) + carry;
   a[i]=r % 10; carry=r / 10;
  }
  if(carry)
   a[++m]=carry;
 }
 free(b);
 a[0]=m;

}
void showa(int *a,int k)
{
 int i;
 printf("% 4d !=",k);
    for(i=a[0];i>0;i--)
  printf("%d", a[i]);
  printf("/n/n");
 
}

int main(int argc, char* argv[])
{
    int a[MAX],n,k;
    printf("input n:");
    scanf("%d",&n);
 a[0]=1;a[1]=1; showa(a,1);
    for(k=2;k<=n;k++)
 {
  pnext(a,k);
  showa(a,k);
 }
 return 0;
}