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

vfbpgyfk

高精度计算素数对个数，是众人向往和追求及奋斗的目标，下面是经近一个月来的探索结果。敬请评议、质疑和赐教。谢谢！

白新岭

那宝吉先生对哥德巴赫猜想中的，偶数的素数对数量，精度的不断追求，探讨，精益求精。锲而不舍的精神值的我学习。
       那老师在vfp编程方面给予我的帮助，念念不忘。使我在数据验证方面有了一把强有力的工具。
        

重生888@

没有真值作参考，屁精度倒做不成！（当然，随便加个系数，也能算精度！）

重生888@

看不懂表中内容，休要乱说，尊重点客观事实好吗

我并不是乱说：你没有真值，哪来动态平均素数对？
要不来个13位偶数：
1111111111112
1111111111114
......
......
1111111111130
试试计算精度？

vfbpgyfk

没有金刚钻，不揽磁器活，你来看看，计算13位51个连续偶数的类偶数平均素数对个数用时3秒。从2023-3-6 08：36发贴，到现在谁也编造不出来这么多的数据。兴好，还有空余表格，那就把因式分解也贴上来吧。
发完贴，看了一下时间，09：02减去08：36=26分钟。
*******************************************************************
回复愚工688：
不应该呀？其实编因式分解程序很简单：
1、先对能够被2整除的进行分解，先设GD=N，进入条件循环，满足条件就退出循环。循环内判断MOD(GD,2)=0，如果=0，GD=GD/2,记下或累计分解字符串。否则，跳出条件循环。
2、对奇数进行分解。继续使用GD，并对GD开平方，余下方法同分解2。全部结束奇数分解后，要判断商的GD是不大于1或√GD<3，如果符合条件，就需要将GD加入因子行列。
3、如果不使用素数表作因式分解，因需要判断素数，则需要时间要长一些。特别对那种大偶数，且是GD/2就是素数，或是最后一个因子比较大，费了半天劲，竟然是个大素数。

重生888@

经与先生数据对比，已和我的数据相差无几,先生的动态新方法不错！不愧为工程师，衷心祝贺！
D（1111111111112）=1059609250
D（1111111111114）=1011445194
D（1111111111116）=2141883940
D（1111111111140）=2697187184

你的数据略大，也许精度高些！究竟谁高，只能求助高手愚工先生！

重生888@

实践出真知呀！  
都是相互学习的结果！在前怎么不是金刚钻？学我分类，反推愚工计算成果，难道忘记了？

重生888@

能够借鉴或是取长补短，善于潜心分析研究，归纳总结出一套独具一格的方式方法，就是在促进科学技术的进步，就是在解决实质性的问题。  发表于 2023-3-7 08:34

别人陈赞你，能否有点谦虚之心？

愚工688

我通常计算到11位的偶数，12位及以上的偶数需要花费的时间太长了，一般不计算。

例：

G(20180125000) = 41387273；Sp( 20180125000 *)≈  41370609.8 , jdz =sp(m)/s(m) ≈0.999597；
G(20180125002) = 51734336；Sp( 20180125002 *)≈  51714552.8 , jdz =sp(m)/s(m) ≈0.999618；
G(20180125004) = 28213467；Sp( 20180125004 *)≈  28206010.9 , jdz =sp(m)/s(m) ≈0.999736；
G(20180125006) = 25874986；Sp( 20180125006 *)≈  25855510.0 , jdz =sp(m)/s(m) ≈0.999247；
G(20180125008) = 53212163；Sp( 20180125008 *)≈  53188477.7 , jdz =sp(m)/s(m) ≈0.999555；
G(20180125010) = 34488610；Sp( 20180125010 *)≈  34474013.3 , jdz =sp(m)/s(m) ≈0.999577；
G(20180125012) = 26043997；Sp( 20180125012 *)≈  26029036.9 , jdz =sp(m)/s(m) ≈0.999426；
G(20180125014) = 63464123；Sp( 20180125014 *)≈  63432184.5 , jdz =sp(m)/s(m) ≈0.999497；
G(20180125016) = 25869177；Sp( 20180125016 *)≈  25855510.0 , jdz =sp(m)/s(m) ≈0.999472；
G(20180125018) = 26137004；Sp( 20180125018 *)≈  26122143.6 , jdz =sp(m)/s(m) ≈0.999431；
G(20180125020) = 83264234；Sp( 20180125020 *)≈  83231258.0 , jdz =sp(m)/s(m) ≈0.999604；
G(20180125022) = 26141818；Sp( 20180125022 *)≈  26136363.0 , jdz =sp(m)/s(m) ≈0.999791；
G(20180125024) = 26182685; Sp( 20180125024 *)≈  26174906.5 , jdz =sp(m)/s(m) ≈0.999703；
G(20180125026) = 55388635; Sp( 20180125026 *)≈  55369615.2 , jdz =sp(m)/s(m) ≈0.999657；
G(20180125028) = 31648432; Sp( 20180125028 *)≈  31634976.9 , jdz =sp(m)/s(m) ≈0.999575；
G(20180125030) = 37628253; Sp( 20180125030 *)≈  37608014.6 , jdz =sp(m)/s(m) ≈0.999462；
G(20180125032) = 53042918; Sp( 20180125032 *)≈  53026876.4 , jdz =sp(m)/s(m) ≈0.999698；
G(20180125034) = 25873522; Sp( 20180125034 *)≈  25859639.6 , jdz =sp(m)/s(m) ≈0.999463；
G(20180125036) = 27599978; Sp( 20180125036 *)≈  27588642.7 , jdz =sp(m)/s(m) ≈0.999589；
G(20180125038) = 51785236; Sp( 20180125038 *)≈  51760886.0 , jdz =sp(m)/s(m) ≈0.999530；
G(20180125040) = 34490848; Sp( 20180125040 *)≈  34478550.9 , jdz =sp(m)/s(m) ≈0.999643；
G(20180125042) = 35475312; Sp( 20180125042 *)≈  35462757.1 , jdz =sp(m)/s(m) ≈0.999646；
G(20180125044) = 51729098; Sp( 20180125044 *)≈  51711020.0 , jdz =sp(m)/s(m) ≈0.999651；
G(20180125046) = 25867808; Sp( 20180125046 *)≈  25855510.0 , jdz =sp(m)/s(m) ≈0.999525；
G(20180125048) = 25868050; Sp( 20180125048 *)≈  25855510.0 , jdz =sp(m)/s(m) ≈0.999515；
以上计算了25个偶数的表法数，全部计算值的精度在0.999以上，余下5个偶数不再计算，精度不会低于0.999；
G(20180125050) = 68979731; Sp( 20180125050 *)≈  68948026.7 , jdz =sp(m)/s(m) ≈
G(20180125052) = 26525287; Sp( 20180125052 *)≈  26518471.8 , jdz =sp(m)/s(m) ≈
G(20180125054) = 25867347; Sp( 20180125054 *)≈  25855510.0 , jdz =sp(m)/s(m) ≈
G(20180125056) = 70100894; Sp( 20180125056 *)≈  70077653.8 , jdz =sp(m)/s(m) ≈
G(20180125058) = 26099275; Sp( 20180125058 *)≈  26097150.3 , jdz =sp(m)/s(m) ≈

start time =18:46:48,end time=19:00:44 ,time use =

G(20180210000) = 34482209；Sp( 20180210000 *)≈  34473673.1 , jdz =sp(m)/s(m) ≈0.99975；
G(20180210002) = 25913942；Sp( 20180210002 *)≈  25906251.4 , jdz =sp(m)/s(m) ≈0.99970；
G(20180210004) = 56188319；Sp( 20180210004 *)≈  56159409.6 , jdz =sp(m)/s(m) ≈0.99949；
G(20180210006) = 28738147；Sp( 20180210006 *)≈  28728061.0 , jdz =sp(m)/s(m) ≈0.99965；
G(20180210008) = 31041710；Sp( 20180210008 *)≈  31026305.8 , jdz =sp(m)/s(m) ≈0.99950;
G(20180210010) = 69977385；Sp( 20180210010 *)≈  69946583.2 , jdz =sp(m)/s(m) ≈0.99956;
G(20180210012) = 25869852；Sp( 20180210012 *)≈  25857625.4 , jdz =sp(m)/s(m) ≈0.99953;
G(20180210014) = 25972971；Sp( 20180210014 *)≈  25957926.9 , jdz =sp(m)/s(m) ≈0.99942;
G(20180210016) = 51730048；Sp( 20180210016 *)≈  51710509.7 , jdz =sp(m)/s(m) ≈0.99962;
G(20180210018) = 25861279；Sp( 20180210018 *)≈  25855254.9 , jdz =sp(m)/s(m) ≈0.99977;

表法数计算式：
Sp( 20180210000 *) = 1/(1+ .1533 )*( 20180210000 /2 -2)*p(m) ≈ 34473673.1 , k(m)= 1.333333
Sp( 20180210002 *) = 1/(1+ .1533 )*( 20180210002 /2 -2)*p(m) ≈ 25906251.4 , k(m)= 1.001972
Sp( 20180210004 *) = 1/(1+ .1533 )*( 20180210004 /2 -2)*p(m) ≈ 56159409.6 , k(m)= 2.172069
Sp( 20180210006 *) = 1/(1+ .1533 )*( 20180210006 /2 -2)*p(m) ≈ 28728061.0 , k(m)= 1.111111
Sp( 20180210008 *) = 1/(1+ .1533 )*( 20180210008 /2 -2)*p(m) ≈ 31026305.8 , k(m)= 1.2
Sp( 20180210010 *) = 1/(1+ .1533 )*( 20180210010 /2 -2)*p(m) ≈ 69946583.2 , k(m)= 2.705314
Sp( 20180210012 *) = 1/(1+ .1533 )*( 20180210012 /2 -2)*p(m) ≈ 25857625.4 , k(m)= 1.000092
Sp( 20180210014 *) = 1/(1+ .1533 )*( 20180210014 /2 -2)*p(m) ≈ 25957926.9 , k(m)= 1.003971
Sp( 20180210016 *) = 1/(1+ .1533 )*( 20180210016 /2 -2)*p(m) ≈ 51710509.7 , k(m)= 2
Sp( 20180210018 *) = 1/(1+ .1533 )*( 20180210018 /2 -2)*p(m) ≈ 25855254.9 , k(m)= 1

愚工688

10个13位的偶数的素数对数量的计算，用了近84分钟。
计算值的精度（ jd ）：

G(1111111111110)= 4144343488 ；Sp( 1111111111110 *)≈  4143428568.6 , jd  ≈ 0.99978；
G(1111111111112)= 1077281340 ；Sp( 1111111111112 *)≈  1077020774.2 , jd  ≈ 0.99976；
G(1111111111114)= 1028319123 ；Sp( 1111111111114 *)≈  1028090715.9 , jd  ≈ 0.99978；
G(1111111111116)= 2177524305 ；Sp( 1111111111116 *)≈  2177079425.9 , jd  ≈ 0.99980；
G(1111111111118)= 1028290769 ；Sp( 1111111111118 *)≈  1028065284.5 , jd  ≈ 0.99978；
G(1111111111120)= 1371883867 ；Sp( 1111111111120 *)≈  1371576988.1 , jd  ≈ 0.99978；
G(1111111111122)= 2056565991 ；Sp( 1111111111122 *)≈  2056130569.0 , jd  ≈ 0.99979；
G(1111111111124)= 1242965181 ；Sp( 1111111111124 *)≈  1242683292.8 , jd  ≈ 0.99977；
G(1111111111126)= 1072884644 ；Sp( 1111111111126 *)≈  1071616188.7 , jd  ≈ 0.99882；
G(1111111111128)= 2091397902 ；Sp( 1111111111128 *)≈  2090980239.6 , jd  ≈ 0.99980；
start time =11:18:11,end time=12:41:37 ,time use =

计算式：
Sp( 1111111111110 *) = 1/(1+ .16935 )*( 1111111111110 /2 -2)*p(m) ≈ 4143428568.6 , k(m)= 4.030317
Sp( 1111111111112 *) = 1/(1+ .16935 )*( 1111111111112 /2 -2)*p(m) ≈ 1077020774.2 , k(m)= 1.047619
Sp( 1111111111114 *) = 1/(1+ .16935 )*( 1111111111114 /2 -2)*p(m) ≈ 1028090715.9 , k(m)= 1.000025
Sp( 1111111111116 *) = 1/(1+ .16935 )*( 1111111111116 /2 -2)*p(m) ≈ 2177079425.9 , k(m)= 2.117647
Sp( 1111111111118 *) = 1/(1+ .16935 )*( 1111111111118 /2 -2)*p(m) ≈ 1028065284.5 , k(m)= 1
Sp( 1111111111120 *) = 1/(1+ .16935 )*( 1111111111120 /2 -2)*p(m) ≈ 1371576988.1 , k(m)= 1.334134
Sp( 1111111111122 *) = 1/(1+ .16935 )*( 1111111111122 /2 -2)*p(m) ≈ 2056130569 , k(m)= 2
Sp( 1111111111124 *) = 1/(1+ .16935 )*( 1111111111124 /2 -2)*p(m) ≈ 1242683292.8 , k(m)= 1.208759
Sp( 1111111111126 *) = 1/(1+ .16935 )*( 1111111111126 /2 -2)*p(m) ≈ 1071616188.7 , k(m)= 1.042362
Sp( 1111111111128 *) = 1/(1+ .16935 )*( 1111111111128 /2 -2)*p(m) ≈ 2090980239.6 , k(m)= 2.033898