#include <stdio.h>
#include <malloc.h>
#define MAXN 1000
void pnext(int a[],int k)
 {
  int *b,m=a[0],j,i,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 write(int *a,int k)
 {
   int i;
   printf("%4d!=",k);
   for(i=a[0];i>0;i--)
    printf("%d",a[i]);
   printf("\n\n");

 }

  void main()
{
int a[MAXN],n,k;
printf("Enter the number n;");
scanf("%d",&n);
a[0]=1;
a[1]=1;
write(a,1);
for(k=2;k<=n;k++)
  {
    pnext(a,k);
    write(a,k);
    getchar();
  }
  getch();
}

问题：阶乘的计算！
   a[0]存储长整数N的M，已知4！=24，计算5！，可以对原来的24累加4次24后得到120。

运行后只能得到1！=1，2！=1，3！=1
 请问r=(i<a[0]?a[i]+b[i]:a[i])+carry;是不是有点问题，我改为r=(i<a[0]?a[i]:a[i]+b[i])+carry;后，只能得到正确的前3位的阶乘！后面的就不对了，

请高手指出其中的错误码，谢谢！

我加了两个标记发现
1.if (carry) {a[++m]=carry;printf("%-5d\n",carry);}
2.   for(i=a[0];i>0;i--)
    printf("%d  ",a[i]);
可见算法本身就有很多错误
不知为何楼主用如此复杂的算法，这样时间复杂度高而且会很容易出错的

谢谢!

其实只是一处有点错误!经过两个小时的研究终于搞定
r=(i<a[0]?a[i]+b[i]:a[i])+carry;

此处小了一个等于号:r=(i<=a[0]?a[i]+b[i]:a[i])+carry

别人问的一个题目!