高精度计算大偶数表为两个素数和的表法数值的实例（以当天日期为随机数选择偶数）

重生888@

先生的数据非常多，也好。但是，我建议您，搞出一个式子，预计某某段的修正值，而且很准；将功不可没！

重生888@

我在您的数据里看到30整倍数的素数对是尾数10和20素数对的和，准确无误，所以我确定我的公式是对的。

愚工688

偶数M表为两个素数和的表法数值的变化是有规律性的，因此是能够比较精确的进行计算的。
今天的日期是2018年08月25日，继续以今天的日期作为随机数，计算更大一些的百亿级别的偶数20180825×2000起的连续偶数M表为两个素数和的表法数计算值Sp(m*)以及计算值的精度 jdz。
注：素对真值s(m)=G(M).系网友Ktprime 的高速筛选素对程序FastGn与我的素对计算程序对素对数量各自的表示形式。



G(40361650000) = 77936749；  Sp( 40361650000 *)≈  77916278.1 , jdz =sp(m)/s(m) ≈ 0.999737；
G(40361650002) = 103967418； Sp( 40361650002 *)≈ 103945890.9 , jdz =sp(m)/s(m) ≈ 0.999793；
G(40361650004) = 55685269；  Sp( 40361650004 *)≈  55674329.1 , jdz =sp(m)/s(m) ≈ 0.999804；
G(40361650006) = 48724048；  Sp( 40361650006 *)≈  48712464.6 , jdz =sp(m)/s(m) ≈ 0.999762；
G(40361650008) = 100663814； Sp( 40361650008 *)≈ 100631737.2 , jdz =sp(m)/s(m) ≈ 0.999681；
G(40361650010) = 74698256；  Sp( 40361650010 *)≈  74680638.8 , jdz =sp(m)/s(m) ≈ 0.999764；
G(40361650012) = 48704672；  Sp( 40361650012 *)≈  48697251.5 , jdz =sp(m)/s(m) ≈ 0.999848；
G(40361650014) = 116902104； Sp( 40361650014 *)≈ 116873403.6 , jdz =sp(m)/s(m) ≈ 0.999754；
G(40361650016) = 48893275；  Sp( 40361650016 *)≈  48884251.6 , jdz =sp(m)/s(m) ≈ 0.999815；
G(40361650018) = 48713683；  Sp( 40361650018 *)≈  48698959.3 , jdz =sp(m)/s(m) ≈ 0.999698；
G(40361650020) = 131497193； Sp( 40361650020 *)≈ 131462539.1 , jdz =sp(m)/s(m) ≈ 0.999736；
G(40361650022) = 50108932；  Sp( 40361650022 *)≈  50093159.3 , jdz =sp(m)/s(m) ≈ 0.999685；
G(40361650024) = 49415576；  Sp( 40361650024 *)≈  49403008.8 , jdz =sp(m)/s(m) ≈ 0.999746；
G(40361650026) = 97406348；  Sp( 40361650026 *)≈  97394503.1 , jdz =sp(m)/s(m) ≈ 0.999878；
G(40361650028) = 58474636；  Sp( 40361650028 *)≈  58461621.5 , jdz =sp(m)/s(m) ≈ 0.999777；
G(40361650030) = 75124650；  Sp( 40361650030 *)≈  75113197.3 , jdz =sp(m)/s(m) ≈ 0.999848；
G(40361650032) = 110881051； Sp( 40361650032 *)≈ 110855532.0 , jdz =sp(m)/s(m) ≈ 0.999770；
G(40361650034) = 49538025；  Sp( 40361650034 *)≈  49522628.7 , jdz =sp(m)/s(m) ≈ 0.999689；
G(40361650036) = 52377620；  Sp( 40361650036 *)≈  52359284.9 , jdz =sp(m)/s(m) ≈ 0.999650；
G(40361650038) = 98269750；  Sp( 40361650038 *)≈  98255153.8 , jdz =sp(m)/s(m) ≈ 0.999851；

素对数量计算式如下：
Sp( 40361650000 *) = 1/(1+ .15649 )*( 40361650000 /2 -2)*p(m) ≈ 77916278.1 , k(m)= 1.600014
Sp( 40361650002 *) = 1/(1+ .15649 )*( 40361650002 /2 -2)*p(m) ≈ 103945890.9 , k(m)= 2.134533
Sp( 40361650004 *) = 1/(1+ .15649 )*( 40361650004 /2 -2)*p(m) ≈ 55674329.1 , k(m)= 1.143275
Sp( 40361650006 *) = 1/(1+ .15649 )*( 40361650006 /2 -2)*p(m) ≈ 48712464.6 , k(m)= 1.000312
Sp( 40361650008 *) = 1/(1+ .15649 )*( 40361650008 /2 -2)*p(m) ≈ 100631737.2 , k(m)= 2.066477
Sp( 40361650010 *) = 1/(1+ .15649 )*( 40361650010 /2 -2)*p(m) ≈ 74680638.8 , k(m)= 1.53357
Sp( 40361650012 *) = 1/(1+ .15649 )*( 40361650012 /2 -2)*p(m) ≈ 48697251.5 , k(m)= 1
Sp( 40361650014 *) = 1/(1+ .15649 )*( 40361650014 /2 -2)*p(m) ≈ 116873403.6 , k(m)= 2.4
Sp( 40361650016 *) = 1/(1+ .15649 )*( 40361650016 /2 -2)*p(m) ≈ 48884251.6 , k(m)= 1.00384
Sp( 40361650018 *) = 1/(1+ .15649 )*( 40361650018 /2 -2)*p(m) ≈ 48698959.3 , k(m)= 1.000035
Sp( 40361650020 *) = 1/(1+ .15649 )*( 40361650020 /2 -2)*p(m) ≈ 131462539.1 , k(m)= 2.699588
Sp( 40361650022 *) = 1/(1+ .15649 )*( 40361650022 /2 -2)*p(m) ≈ 50093159.3 , k(m)= 1.028665
Sp( 40361650024 *) = 1/(1+ .15649 )*( 40361650024 /2 -2)*p(m) ≈ 49403008.8 , k(m)= 1.014493
Sp( 40361650026 *) = 1/(1+ .15649 )*( 40361650026 /2 -2)*p(m) ≈ 97394503.1 , k(m)= 2
Sp( 40361650028 *) = 1/(1+ .15649 )*( 40361650028 /2 -2)*p(m) ≈ 58461621.5 , k(m)= 1.200512
Sp( 40361650030 *) = 1/(1+ .15649 )*( 40361650030 /2 -2)*p(m) ≈ 75113197.3 , k(m)= 1.542452
Sp( 40361650032 *) = 1/(1+ .15649 )*( 40361650032 /2 -2)*p(m) ≈ 110855532 , k(m)= 2.276423
Sp( 40361650034 *) = 1/(1+ .15649 )*( 40361650034 /2 -2)*p(m) ≈ 49522628.7 , k(m)= 1.016949
Sp( 40361650036 *) = 1/(1+ .15649 )*( 40361650036 /2 -2)*p(m) ≈ 52359284.9 , k(m)= 1.0752
Sp( 40361650038 *) = 1/(1+ .15649 )*( 40361650038 /2 -2)*p(m) ≈ 98255153.8 , k(m)= 2.017673

start time =14:51:46,end time=15:03:33 ,

重生888@

我不懂最后一行意思，是公吗？

愚工688

偶数M表为两个素数和的表法数值的变化是有规律性的，因此是能够比较精确的进行计算的。
今天的日期是2018年08月30日，继续以今天的日期作为随机数，计算更大一些的百亿级别的偶数20180830×2000起的连续偶数M表为两个素数和的表法数计算值Sp(m*)以及计算值的精度 jdz。
注：素对真值s(m)=G(M).系网友Ktprime 的高速筛选素对程序FastGn与我的素对计算程序对素对数量各自的表示形式。

G(40361660000) = 65103082；Sp( 40361660000 *)≈  65092428.1 , jdz =sp(m)/s(m) ≈ 0.999836；
G(40361660002) = 48716513；Sp( 40361660002 *)≈  48711070.6 , jdz =sp(m)/s(m) ≈ 0.999888；
G(40361660004) = 97420084；Sp( 40361660004 *)≈  97394527.1 , jdz =sp(m)/s(m) ≈ 0.999738；
G(40361660006) = 48812520；Sp( 40361660006 *)≈  48807688.2 , jdz =sp(m)/s(m) ≈ 0.999901；
G(40361660008) = 49903673；Sp( 40361660008 *)≈  49894983.9 , jdz =sp(m)/s(m) ≈ 0.999826；
G(40361660010) =156087817；Sp( 40361660010 *)≈ 156044483.4 , jdz =sp(m)/s(m) ≈ 0.999722；
G(40361660012) = 52171755；Sp( 40361660012 *)≈  52163815.7 , jdz =sp(m)/s(m) ≈ 0.999848；
G(40361660014) = 53142407；Sp( 40361660014 *)≈  53124287.5 , jdz =sp(m)/s(m) ≈ 0.999659；
G(40361660016) = 98406076；Sp( 40361660016 *)≈  98378191.3 , jdz =sp(m)/s(m) ≈ 0.999717；
G(40361660018) = 48725589；Sp( 40361660018 *)≈  48726687.9 , jdz =sp(m)/s(m) ≈ 1.000023；
G(40361660020) = 74833977；Sp( 40361660020 *)≈  74816097.7 , jdz =sp(m)/s(m) ≈ 0.999761；
G(40361660022) = 98651551；Sp( 40361660022 *)≈  98627229.2 , jdz =sp(m)/s(m) ≈ 0.999753；
G(40361660024) = 62016986；Sp( 40361660024 *)≈  62002541.8 , jdz =sp(m)/s(m) ≈ 0.999767；
G(40361660026) = 48710007；Sp( 40361660026 *)≈  48697263.6 , jdz =sp(m)/s(m) ≈ 0.999738；
G(40361660028) = 99946969；Sp( 40361660028 *)≈  99922042.5 , jdz =sp(m)/s(m) ≈ 0.999751；
G(40361660030) = 64941632；Sp( 40361660030 *)≈  64933157.4 , jdz =sp(m)/s(m) ≈ 0.999869；
G(40361660032) = 54955123；Sp( 40361660032 *)≈  54942810.0 , jdz =sp(m)/s(m) ≈ 0.999776；
G(40361660034) = 97406700；Sp( 40361660034 *)≈  97394527.2 , jdz =sp(m)/s(m) ≈ 0.999875；
G(40361660036) = 49222703；Sp( 40361660036 *)≈  49209866.4 , jdz =sp(m)/s(m) ≈ 0.999739；
G(40361660038) = 59441395；Sp( 40361660038 *)≈  59427169.2 , jdz =sp(m)/s(m) ≈ 0.999761；
start time =14:15:30,end time=14:34:00 ,time use =

计算式如下：
Sp( 40361660000 *) = 1/(1+ .15649 )*( 40361660000 /2 -2)*p(m) ≈ 65092428.1 ,
Sp( 40361660002 *) = 1/(1+ .15649 )*( 40361660002 /2 -2)*p(m) ≈ 48711070.6 ,
Sp( 40361660004 *) = 1/(1+ .15649 )*( 40361660004 /2 -2)*p(m) ≈ 97394527.1 ,
Sp( 40361660006 *) = 1/(1+ .15649 )*( 40361660006 /2 -2)*p(m) ≈ 48807688.2 ,
Sp( 40361660008 *) = 1/(1+ .15649 )*( 40361660008 /2 -2)*p(m) ≈ 49894983.9 ,
Sp( 40361660010 *) = 1/(1+ .15649 )*( 40361660010 /2 -2)*p(m) ≈ 156044483.4 ,
Sp( 40361660012 *) = 1/(1+ .15649 )*( 40361660012 /2 -2)*p(m) ≈ 52163815.7 ,
Sp( 40361660014 *) = 1/(1+ .15649 )*( 40361660014 /2 -2)*p(m) ≈ 53124287.5 ,
Sp( 40361660016 *) = 1/(1+ .15649 )*( 40361660016 /2 -2)*p(m) ≈ 98378191.3 ,
Sp( 40361660018 *) = 1/(1+ .15649 )*( 40361660018 /2 -2)*p(m) ≈ 48726687.9 ,
Sp( 40361660020 *) = 1/(1+ .15649 )*( 40361660020 /2 -2)*p(m) ≈ 74816097.7 ,
Sp( 40361660022 *) = 1/(1+ .15649 )*( 40361660022 /2 -2)*p(m) ≈ 98627229.2 ,
Sp( 40361660024 *) = 1/(1+ .15649 )*( 40361660024 /2 -2)*p(m) ≈ 62002541.8 ,
Sp( 40361660026 *) = 1/(1+ .15649 )*( 40361660026 /2 -2)*p(m) ≈ 48697263.6 ,
Sp( 40361660028 *) = 1/(1+ .15649 )*( 40361660028 /2 -2)*p(m) ≈ 99922042.5 ,
Sp( 40361660030 *) = 1/(1+ .15649 )*( 40361660030 /2 -2)*p(m) ≈ 64933157.4 ,
Sp( 40361660032 *) = 1/(1+ .15649 )*( 40361660032 /2 -2)*p(m) ≈ 54942810 ,
Sp( 40361660034 *) = 1/(1+ .15649 )*( 40361660034 /2 -2)*p(m) ≈ 97394527.2 ,
Sp( 40361660036 *) = 1/(1+ .15649 )*( 40361660036 /2 -2)*p(m) ≈ 49209866.4 ,
Sp( 40361660038 *) = 1/(1+ .15649 )*( 40361660038 /2 -2)*p(m) ≈ 59427169.2 ,

从这些偶数的素对计算值的精度可以看到，已经有偶数的计算值精度大于1了，以后更大的4百亿级别偶数的计算仍然维持修正系数μ= .15649 ，那么出现计算值精度大于1的现象会逐渐增多。

愚工688

重生888@ 发表于 2018-8-26 16:09
我不懂最后一行意思，是公吗？

我对这些偶数的素对计算式的运行时间。
起始时间—— start time =14:51:46,
程序结束时间 —— end time=15:03:33 ,

重生888@

愚工688 发表于 2018-8-30 19:21
我对这些偶数的素对计算式的运行时间。
起始时间—— start time =14:51:46,
程序结束时间 —— end ti ...

谢谢好友的回复！如果是2020年8月30日，您的修正值该是多少？

愚工688

重生888@ 发表于 2018-8-31 02:46
谢谢好友的回复！如果是2020年8月30日，您的修正值该是多少？

看看实际计算式，没有变化，计算值精度也没有什么变化。
Sp( 40401660000 *) = 1/(1+ .15649 )*( 40401660000 /2 -2)*p(m) ≈ 141796609.3 , k(m)= 2.909147
Sp( 40401660002 *) = 1/(1+ .15649 )*( 40401660002 /2 -2)*p(m) ≈ 57389476.5 , k(m)= 1.177422
Sp( 40401660004 *) = 1/(1+ .15649 )*( 40401660004 /2 -2)*p(m) ≈ 49910557.5 , k(m)= 1.023982
Sp( 40401660006 *) = 1/(1+ .15649 )*( 40401660006 /2 -2)*p(m) ≈ 97483286.8 , k(m)= 2
Sp( 40401660008 *) = 1/(1+ .15649 )*( 40401660008 /2 -2)*p(m) ≈ 48762752.8 , k(m)= 1.000433
Sp( 40401660010 *) = 1/(1+ .15649 )*( 40401660010 /2 -2)*p(m) ≈ 67395852.6 , k(m)= 1.382716
Sp( 40401660012 *) = 1/(1+ .15649 )*( 40401660012 /2 -2)*p(m) ≈ 125000424.5 , k(m)= 2.564551
Sp( 40401660014 *) = 1/(1+ .15649 )*( 40401660014 /2 -2)*p(m) ≈ 48751774.7 , k(m)= 1.000208
Sp( 40401660016 *) = 1/(1+ .15649 )*( 40401660016 /2 -2)*p(m) ≈ 49784395 , k(m)= 1.021393
Sp( 40401660018 *) = 1/(1+ .15649 )*( 40401660018 /2 -2)*p(m) ≈ 98309644.7 , k(m)= 2.016954


G(40401660000) = 141832022； jdz =sp(m)/s(m) ≈ 0.999750；
G(40401660002) = 57409341；  jdz =sp(m)/s(m) ≈ 0.999654；
G(40401660004) = 49916594；  jdz =sp(m)/s(m) ≈ 0.999879；
G(40401660006) = 97512949；  jdz =sp(m)/s(m) ≈ 0.999696；
G(40401660008) = 48773320；  jdz =sp(m)/s(m) ≈ 0.999783；
G(40401660010) = 67409068；  jdz =sp(m)/s(m) ≈ 0.999804；
G(40401660012) = 125023806； jdz =sp(m)/s(m) ≈ 0.999813；
G(40401660014) = 48760534；  jdz =sp(m)/s(m) ≈ 0.999820；
G(40401660016) = 49802115；  jdz =sp(m)/s(m) ≈ 0.999644；
G(40401660018) = 98324976；  jdz =sp(m)/s(m) ≈ 0.999844；

重生888@

能不能再扩大，这个修正系数，对任何段的偶数都适应吗？谢谢！

愚工688

重生888@ 发表于 2018-8-31 12:19
能不能再扩大，这个修正系数，对任何段的偶数都适应吗？谢谢！

在一楼，我已经说清楚：μ系相对误差修正值，只适用一定范围的偶数区域。
能够修正一切偶数范围的计算值的相对误差的系数以达到高精度计算的，是不可能的，否则就直接是 个系数了。