大数的阶乘（C语言实现）

我之前贴了一个大数运算之100的阶乘（C语言实现），当时由于时间紧，没怎细看，也没优化就给传上来了，后来发现运行效率确实不行，主要因为使用了累计相加来实现乘法，怪不得被人踩了几下，今天又琢磨了一下，重新写了一个，不再使用累加代替相乘，感觉这个勉强要比前面的那个好，不过算法这个东西没有最好的，只有更好的，欢迎各位交流！

#include "stdafx.h"
#include "math.h"
#include "stdio.h"
#include "stdlib.h"

void dgCal(int a[],int pos)
{
	int i,cy=0;//进位值
	for(i=0;i<=pos;i++)//从0-pos逐个检查是否要进位
	{
		a[i]+=cy;
		if(a[i]<=9)//不需进位
		{
			cy=0;//进位制为0
		}
		else if(a[i]>9&&i<pos)//需进位，但未达到最高位
		{
			cy=a[i]/10;//进位值
			a[i]=a[i]%10;
		}
		else if(a[i]>9&&i>=pos)//需进位，已达到最高位
		{
			while(a[i]>9)//循环向前进位
			{
				cy=a[i]/10;
				a[i]=a[i]%10;
				i++;
				a[i]=cy;
			}
		}
	}
}

int _tmain(int argc, _TCHAR* argv[])
{
	int i,j;
	int *a;//结果数组
	int num=200;//要求的阶乘数
	int pos=0;//最高位下标

	int digit=0;//阶乘的位数
	double sum=0;
	for(i=1;i<=num;i++)
	{
		sum += log10((double)i);//计算阶乘的位数
	}
	digit=(int)sum+1;//阶乘的位数

	if(!(a=(int *)malloc( (digit+1) * sizeof(int) )))//分配内存
	{
		printf("内存分配失败！\n");
		return 0;
	}
	for(i=0;i<=digit;i++)
	{
		a[i]=0;//初始化
	}
	a[0]=1;//个位为1

	for(i=2;i<=num;i++)//2-num逐个与原来的结果相乘
	{
		for(j=digit;j>=0;j--)//查找最高位的下标
		{
			if(a[j]!=0)
			{
				pos=j;
				break;
			}
		}

		for(j=0;j<=pos;j++)
		{
			a[j]*=i;//逐个相乘
		}

		dgCal(a,pos);//乘完一轮后进位处理

	}

	for(j=digit;j>=0;j--)//查找最高位的下标
	{
		if(a[j]!=0)
		{
			pos=j;
			break;
		}
	}

	printf("200的阶乘：\n");
	printf("位数：%d\n",pos+1);
	for(i=pos;i>=0;i--)
	{
		printf("%d",a[i]);//打印结果
	}
	printf("\n");

	return 0;
}

来张图：

有问题可以随时交流！