一个程序拿出来show show，求n的阶乘保留小数点后16位，在P41.7G下计算100万的阶乘大概是1秒。大家给点意见

[ywhbn]

// 在VC++6.0下写的
// 求n的阶乘保留小数点后16位，1.0e+300 进制
// 代码简单就不注释了，在P41.7G，256M内存下计算100万的阶乘大概是1秒

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

// 读浮点数的指数
#define get_exp(p) ( atoi(p + strlen(p) - 3) )

main()
{
time_t
stime,
etime;
unsigned char
pr[50];
unsigned long int
n,
re;
unsigned long int register
i;
double
r;

while(1)
{
printf("input n [<=1.0e+8 , 0 -- exit] :");
scanf("%ld", &n); // n不是整数会死循环
if(!n)
break;
r = 1.0;
re = 0;
// 开始计算
time(&stime);
for(i = 2; i <= n; i++)
{
r = r * i;
if(r >= 1.0e+300)
{
r /= 1.0e+300;
re += 300;
}
}
time(&etime);
sprintf(pr, " %d! = %.20e", n, r);
sprintf(pr+strlen(pr)-3, "%ld", get_exp(pr)+re);
printf(" %s \n ", pr);
printf(" use time: %.4f s\n", difftime(etime, stime));
}
}

[masse]

……汗

[mathe]

可是你的代码太慢了:)

[dot99]

去算法版研究一下先:D

[imRainman]

不错，有创意。

可以用在不要求精确结果的地方~

[mathe]

#define M_PI 3.141592653589793
double v;
while(1)
{
printf("input n [<=1.0e+8 , 0 -- exit] :");
scanf("%ld", &n); // n不是整数会死循环
if(!n)
break;
if(n>1000){
int m=n+1;
int intpart;
v=log(2*PI*m)/2.0+m*log(m)-m+1/(12*m)+1/(288.0*m*m)-139/(51840.0*m*m*m)-571/(2488320.0*m*m*m));
v/=log(10.0);
intpart=(int)v;v-=intpart;
v*=log(10.0);v=exp(v);
printf("%fe+%d\n",v,intpart);
}else{
v=1.0;
for(i=1;i<=n;i++)v*=i;
printf("%f\n",v);
}

[mathe]

弄错了一步，m=n就可以了，不需要m=n+1.

[ywhbn]

计算的就是这几句：
for(i = 2; i <= n; i++)
{
r = r * i;
if(r >= 1.0e+300)
{
r /= 1.0e+300;
re += 300;
}
}
结果保存在 r和re 里

[clin003]

不客气了 ，先坐了先

[idler]

怎么又发了。。。
100万的精度是10位有效数字

[aladar]

看到1e300进制就吐了，你拿浮点玩来着？精度太低了吧。。。还100万呢-_-!!