天文数的阶乘计算
在C语言里unsigned long int型的的整数的仅可表示0~4294967295之间的数,而12! = 479001600,13! = 6227020800。可见,用unsigned long int型的整数做阶乘运算时最多只能计算到12的阶乘。用函数double sqrt(double x)做开方算运算时,只能计算到16位有效数字。long double型的实数虽然可表示10-4931~10-4932的数,但其有效数字也只有18~19位。
我编写了一个“天文数字计算”程序突破了这一限制,可以把数字的长度扩充到无穷多位。除了能做加、减、乘、除、求模等基本运算外还可以做阶乘、乘方、开平方等运算。
例如:10000! = ?、(2002^2000)%9999 = 9394、2002的平方根的小数点后第10000位是4。
下面这个程序只是我用编写的“天文数字计算”里的一个计算。
因为在程序代码中使用了中文,所以这个程序如果不在中文DOS下运行,可能会出现乱码,但不影响程序的计算结果。
注:因为DOS能访问的内存有限,所以在DOS下可把计算结果扩充到30000多位。但如果用VC把它编绎成Windows程序,则可以计算到“真正的无穷多位”(与机子配置有关),一般计算到几十万位是没有问题的。
/* 此程序在TC2.0、TC3.0,BC,VC下都可编绎 */
#define M 20000 /* 结果位数,DOS能访问的内存有限,不要超过 30000 位*/
#define N (M+5)
main()
{
int Num;
reGISter int i,j,k,flag;
register unsigned int n,m,pc;
unsigned char str_n[5],result_0[N],result_1[N];
void fun_print_result(char *result,int flag);
int fun_mul(char *a,char *b,char *c,int flag);
printf("输入计算其阶乘的数:Num = ");
scanf("%d",&Num);
if(Num >= 100)
{
printf( "\n正在进行阶乘运算. . .\n");
printf( "已完成 00000 的阶乘,请等待. . .");
}
for (j=0;j
0条评论