阶乘和圆周率

无心人 · 发表于 2008-12-25 07:59:23

今天发现了个严重的问题
似乎编译的程序
没去掉#define TEST
导致连续两天的计算均没有存储中间结果
只有去掉#define TEST
才能每隔1000亿存储一次结果
哎
难道停止掉重新来?
还是赌10天内不停电？

gxqcn · 发表于 2008-12-25 08:01:18

那要看你们单位元旦放假期间是否要断电了。

无心人 · 发表于 2008-12-25 08:31:37

不断电
但电力可能停电啊

无心人 · 发表于 2008-12-25 09:11:42

重新启动，估计算到了37000亿
但在25000亿启动了，每天算1万亿左右

无心人 · 发表于 2009-1-6 07:54:39

目前进度：
15 7000 0000 0000
还没出现结果，已经是上个结果的10倍了
再等最多10天，到25万亿如果还没出现就结束掉
有点怀疑是否错过了

无心人 · 发表于 2009-1-12 16:41:22

目前计算到了22 8000 0000 0000
还没发现结果
期待第三方数据

暂时中断了

wayne · 发表于 2010-8-18 21:01:47

20# 无心人

计算N的阶乘的对数很容易，计算阶乘的位数也很容易。
1、设 f= LogGamma[N+1]/Log[10] ,即N的阶乘的底为10的对数。
2、设 L=Ceiling[n Log[10, n/E] + Log[10, 2 n Pi]/2] ，即N的阶乘的位数。

那么N的阶乘的前几位就是10^(f-L+1) .

上述过程都是瞬间即可完成。
无心人20楼的所有计算结果也是瞬间即可检验，可惜不能把上述过程逆着过来演算！
1154318938997的阶乘前100位：
31415926535841388545285400933820805520621988240910055245723955066982684225509763325500

wayne · 发表于 2010-8-18 21:03:21

可惜啊

wayne · 发表于 2010-8-19 11:38:59

俺四核4G内存的电脑跑了通宵，也是没找到结果，没能突破无心人的计算记录

wayne · 发表于 2010-8-19 11:47:06

问题本质就是解n，使得n的对数gamma函数的值的小数部分等于pi的对数值。

真不敢相信，除了搜索，就没其他招数了。

账号		自动登录	找回密码
密码			欢迎注册

[讨论] 阶乘和圆周率