数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
楼主: vfbpgyfk

探索高精度计算素数对个数的弥合计算公式

[复制链接]
发表于 2023-3-13 19:17 | 显示全部楼层
yangchuanju
82楼愚公已经给出G(500000000000) = 655630055,除以4乘上3等于491722541.25,从哪里冒出一个645225314?  

那是愚工真值,我的计算值是645225314/4*3=483918985

点评

重生"/4*3"可能算翻了,杨的除以4乘上3实际上是消除波动系数1.3333;重生的645225314可能是按他的[N+F*N/ln(N)]/ln(N)^2算来的,645225314/4*3=483918985无意义。  发表于 2023-3-13 20:35
4*3是分母吗?如果是的话,应该用括号括起来。  发表于 2023-3-13 19:41
回复 支持 反对

使用道具 举报

发表于 2023-3-13 19:29 | 显示全部楼层

计算动态系数公式不变,没有人工干预,哪来的调低系数一说?

那就要看真值怎么说话?
回复 支持 反对

使用道具 举报

发表于 2023-3-14 06:52 | 显示全部楼层

4*3是分母吗?如果是的话,应该用括号括起来

不是哪个在前算那个吗? N/4*3不需要打括号;而真正要打括号的是N/(4*3)

点评

举个例子说吧,当N=12时,12/4*3=3*3=9,而12/(4*3)=12/12=1。所以,加不加括号很重要。如果不加括号,是否直接写3*N/4为好?或是0.75*N?  发表于 2023-3-14 07:17
回复 支持 反对

使用道具 举报

发表于 2023-3-14 07:09 | 显示全部楼层
本帖最后由 yangchuanju 于 2023-3-14 07:13 编辑

为验证重生系数是否适用于无穷大,下面先引用几组哥猜数据表:
网页A065577给出10的1-14次方的无序哥猜数(10的15次方哥猜数由愚公688计算):
A065577
1 2
2 6
3 28
4 127
5 810
6 5402
7 38807
8 291400
9 2274205
10 18200488
11 149091160
12 1243722370
13 10533150855
14 90350630388
10^15 = 783538341852

网页A006307给出2的0-40次方的无序哥猜数:
A006307-40个2^n哥猜数
0 0        21 7471
1 0        22 13705
2 1        23 24928
3 1        24 45746
4 2        25 83467
5 2        26 153850
6 5        27 283746
7 3        28 525236
8 8        29 975685
9 11        30 1817111
10 22        31 3390038
11 25        32 6341424
12 53        33 11891654
13 76        34 22336060
14 151        35 42034097
15 244        36 79287664
16 435        37 149711134
17 749        38 283277225
18 1314        39 536710100
19 2367        40 1018369893
20 4239      

以下2的41-50次方哥猜数由愚公688计算
41        1934814452
42        3680759328
43        7010898161
44        13369466800
45        25522944188
46        48776696083
47        93311971184
48        178680063951
49        342469661688
50        656978437719

网页A116979 和A062311分别给出12个p#和15个n!的哥猜数:
P#       A116979        N!        A062311
1        0 0        1        1 0
2        1 0        2        2 0
6        2 1        6        3 1
30        3 3        24        4 3
210        4 19        120        5 12
2310        5 114        720        6 39
30030        6 905        5040        7 184
510510        7 9493        40320        8 951
9699690        8 124180        362880        9 5531
223092870        9 2044847        3628800        10 38713
6469693230        10 43755729        39916800        11 346207
200560490130         11 1043468386        479001600        12 3130812
7420738134810         12 30309948241        6227020800        13 34444964
———————        ———————   87178291200        14 382437428
———————        ———————   1307674368000         15 4637235145

点评

是要我一个一个计算吗?  发表于 2023-3-14 16:37
看不懂杨先生发次帖子是啥意思?  发表于 2023-3-14 08:58
对于素数阶乘偶数来说,很快就与重生888@形成鲜明对比,一个是一成不变的刻舟求剑,一个是因时(偶数情况)而异,那个更符合实际,不言而喻。  发表于 2023-3-14 08:21
回复 支持 反对

使用道具 举报

发表于 2023-3-14 09:01 | 显示全部楼层
    6227020800  的素数对是34444964吗?如果是,是哪个计算的?我敢说他不对!

点评

我的筛选真值的软件只能够用于10^16以下,(10^16的素对只运行过一次,不敢再运行了,时间太长)所以一般偶数的范围控制在13位数,因为那软件可以自动筛选连续的偶数。14位数的偶数只能一个一个筛选。  发表于 2023-3-15 21:53
我的t2修正系数是有使用范围限制的,即t2≥1 。在t2<1的范围,哈-李计算式的值的相对误差依据很小,是不需要再进行修正了。(这个范围也是我目前无法计算到的范围。)我从来没有考虑过什么计算无穷大。  发表于 2023-3-15 21:44
愚公、那宝吉的计算式都不能用于无穷大偶数,你的计算式能不能用于无穷大,我打算试一试再说,现在不要你干任何事。  发表于 2023-3-14 17:18
我的计算值是: 34467665.7。  发表于 2023-3-14 11:05
回复 支持 反对

使用道具 举报

发表于 2023-3-14 09:12 | 显示全部楼层
本帖最后由 愚工688 于 2023-3-14 01:48 编辑
重生888@ 发表于 2023-3-13 11:06
我粗略计算:
D(11111111111110)=11405299316*52/51*78/77=11779957735
真值不小于11779957735/0.9 ...



G(11111111111110)=12000220328;1115.63s.
G(11111111111112)=17470455584;1994.64s.
G(11111111111114)=11172207831;
G(11111111111116)=8805605145;835.81s
G(11111111111118)=17496051995;1512.11sec.



百亿的没有计算过,时间太长。
回复 支持 反对

使用道具 举报

发表于 2023-3-14 09:54 | 显示全部楼层
本帖最后由 愚工688 于 2023-3-14 02:45 编辑

我只有使用连乘式计算的14位连续偶数的计算素对数据:

G(11111111111110)=12000220328;
Sp( 11111111111110 *)≈  12000127443.5 , Δ≈-0.00000774, k(m)= 1.37713
G(11111111111112)=17470455584;
Sp( 11111111111112 *)≈  17470524646.2 , Δ≈0.00000395, k(m)= 2.00491
G(11111111111114)=11172207831;
Sp( 11111111111114 *)≈  11172326401 , Δ≈0.0000106, k(m)= 1.28213
G(11111111111116)=8805605145;
Sp( 11111111111116 *)≈  8805574675.3 , Δ≈-0.00000346, k(m)= 1.01053
G(11111111111118)=17496051995;
Sp( 11111111111118 *)≈  17496043799.3 , Δ≈ -0.00000047, k(m)= 2.00784

Sp( 11111111111110 *) = 1/(1+ .17621 )*( 11111111111110 /2 -2)*p(m) ≈ 12000127443.5 , k(m)= 1.37713
Sp( 11111111111112 *) = 1/(1+ .17621 )*( 11111111111112 /2 -2)*p(m) ≈ 17470524646.2 , k(m)= 2.00491
Sp( 11111111111114 *) = 1/(1+ .17621 )*( 11111111111114 /2 -2)*p(m) ≈ 11172326401 , k(m)= 1.28213
Sp( 11111111111116 *) = 1/(1+ .17621 )*( 11111111111116 /2 -2)*p(m) ≈ 8805574675.3 , k(m)= 1.01053
Sp( 11111111111118 *) = 1/(1+ .17621 )*( 11111111111118 /2 -2)*p(m) ≈ 17496043799.3 , k(m)= 2.00784

点评

先生将比值,提高到0.99999... 为您点赞!  发表于 2023-3-14 10:35
回复 支持 反对

使用道具 举报

发表于 2023-3-14 10:44 | 显示全部楼层
本帖最后由 愚工688 于 2023-3-14 04:41 编辑
愚工688 发表于 2023-3-14 01:54
我只有使用连乘式计算的13位连续偶数的计算素对数据:

G(11111111111110)=12000220328;


这是因为相对误差修正系数取自的区域与计算的连续偶数区域重合,故计算精度高。下面我在计算12000000000000起始的连续偶数,看看计算精度有什么变化。


用了近2个小时的时间。

Sp( 12000000000000 *) = 1/(1+ .17621 )*( 12000000000000 /2 -2)*p(m) ≈ 24967794114.5 , k(m)= 2.66667
Sp( 12000000000002 *) = 1/(1+ .17621 )*( 12000000000002 /2 -2)*p(m) ≈ 11240581228    , k(m)= 1.20054
Sp( 12000000000004 *) = 1/(1+ .17621 )*( 12000000000004 /2 -2)*p(m) ≈ 9506967759      , k(m)= 1.01538
Sp( 12000000000006 *) = 1/(1+ .17621 )*( 12000000000006 /2 -2)*p(m) ≈ 19190753728.9 , k(m)= 2.04965
回复 支持 反对

使用道具 举报

发表于 2023-3-14 10:46 | 显示全部楼层
计算动态系数公式不变,没有人工干预,哪来的调低系数一说?  

精度怎么样?

点评

还没有筛选真值呢!要几个小时才能完成,并且只能一个个来。  发表于 2023-3-14 12:43
精度不错呀!  发表于 2023-3-14 11:19
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-3-14 11:21 | 显示全部楼层
我先把计算结果发上来,等真值贴上来后再论计算精度优劣,免得放马后炮,说不清楚是怎么搞出来的。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2025-6-27 00:07 , Processed in 0.087016 second(s), 15 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表