我之前贴了一个大数运算之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
{
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;
}
来张图:
有问题可以随时交流!