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

愚工688 · 发表于 2017-7-7 13:04

偶数M的表法数值的变化是有规律性的，因此也是能够比较精确的进行计算的。

今天的日期是2017年7月7日，
继续以今天的日期作为随机数，计算偶数2017070700起连续50个偶数表为两个素数和的表法数值Sp(m)。

D( 2017070700 )= 8605750 Sp(m)= 8605704.95580 Δ(m)=-.00001 k(m)= 2.68613
D( 2017070702 )= 3204319 Sp(m)= 3204472.46332 Δ(m)= .00005 k(m)= 1.00022
D( 2017070704 )= 3769207 Sp(m)= 3769122.70231 Δ(m)=-.00002 k(m)= 1.17647
D( 2017070706 )= 8383404 Sp(m)= 8383869.03077 Δ(m)= .00006 k(m)= 2.61689
D( 2017070708 )= 3217272 Sp(m)= 3217159.13304 Δ(m)=-.00004 k(m)= 1.00418
D( 2017070710 )= 4272003 Sp(m)= 4271672.40866 Δ(m)=-8.00D-05 k(m)= 1.33333
D( 2017070712 )= 6417471 Sp(m)= 6418099.54268 Δ(m)= .0001    k(m)= 2.00331
D( 2017070714 )= 3265476 Sp(m)= 3266573.02486 Δ(m)= .00034 k(m)= 1.01961
D( 2017070718 )= 6990010 Sp(m)= 6990009.42371 Δ(m)= 0       k(m)= 2.18182
D( 2017070720 )= 5125036 Sp(m)= 5126006.91580 Δ(m)= .00019 k(m)= 1.6
D( 2017070722 )= 3203686 Sp(m)= 3203754.32555 Δ(m)= .00002 k(m)= 1
D( 2017070724 )= 6542816 Sp(m)= 6542107.34308 Δ(m)=-.00011 k(m)= 2.04201
D( 2017070726 )= 3557248 Sp(m)= 3559727.03545 Δ(m)= .0007 k(m)= 1.11111
D( 2017070728 )= 3228647 Sp(m)= 3227485.84868 Δ(m)=-.00036 k(m)= 1.00741
D( 2017070730 )= 8540346 Sp(m)= 8543344.90203 Δ(m)= .00035 k(m)= 2.66667
D( 2017070732 )= 3202151 Sp(m)= 3203754.34144 Δ(m)= .0005 k(m)= 1
D( 2017070734 )= 3846209 Sp(m)= 3844505.21354 Δ(m)=-.00044 k(m)= 1.2
D( 2017070736 )= 6409155 Sp(m)= 6407508.69558 Δ(m)=-.00026 k(m)= 2
D( 2017070738 )= 3373840 Sp(m)= 3374140.88382 Δ(m)= 9.00D-05 k(m)= 1.05318
D( 2017070740 )= 4553511 Sp(m)= 4556450.637    Δ(m)= .00065 k(m)= 1.42222
D( 2017070742 )= 6786656 Sp(m)= 6784420.99197 Δ(m)=-.00033 k(m)= 2.11765
D( 2017070744 )= 3660674 Sp(m)= 3661433.55485 Δ(m)= .00021 k(m)= 1.14286
D( 2017070746 )= 3290356 Sp(m)= 3289509.31891 Δ(m)=-.00026 k(m)= 1.02677
D( 2017070748 )= 8607527 Sp(m)= 8605199.67495 Δ(m)=-.00027 k(m)= 2.68597
D( 2017070750 )= 4271799 Sp(m)= 4271672.49337 Δ(m)=-.00003 k(m)= 1.33333
D( 2017070752 )= 3206274 Sp(m)= 3205074.85947 Δ(m)=-.00037 k(m)= 1.00041
D( 2017070754 )= 6408916 Sp(m)= 6407508.75276 Δ(m)=-.00022 k(m)= 2
D( 2017070756 )= 3342417 Sp(m)= 3342547.77917 Δ(m)= .00004 k(m)= 1.04332
D( 2017070758 )= 3204049 Sp(m)= 3203754.38273 Δ(m)=-9.00D-05 k(m)= 1
D( 2017070760 )= 8694268 Sp(m)= 8693228.27522 Δ(m)=-.00012 k(m)= 2.71345
D( 2017070762 )= 3952863 Sp(m)= 3954348.27453 Δ(m)= .00038 k(m)= 1.23429
D( 2017070764 )= 3202192 Sp(m)= 3203754.39226 Δ(m)= .00049 k(m)= 1
D( 2017070766 )= 6440065 Sp(m)= 6440034.21621 Δ(m)= 0       k(m)= 2.01015
D( 2017070768 )= 3319648 Sp(m)= 3320595.14609 Δ(m)= .00029 k(m)= 1.03647
D( 2017070770 )= 5176524 Sp(m)= 5177784.89179 Δ(m)= .00024 k(m)= 1.61616
D( 2017070772 )= 6407710 Sp(m)= 6408508.96334 Δ(m)= .00012 k(m)= 2.00031
D( 2017070774 )= 3482827 Sp(m)= 3484040.01345 Δ(m)= .00035 k(m)= 1.08749
D( 2017070776 )= 3842538 Sp(m)= 3844505.29359 Δ(m)= .00051 k(m)= 1.2
D( 2017070778 )= 6505883 Sp(m)= 6506085.88791 Δ(m)= .00003 k(m)= 2.03077
D( 2017070780 )= 4545345 Sp(m)= 4547041.9734 Δ(m)= .00037 k(m)= 1.41929
D( 2017070782 )= 3216106 Sp(m)= 3215994.40653 Δ(m)=-.00003 k(m)= 1.00382
D( 2017070784 )= 6423476 Sp(m)= 6420075.23823 Δ(m)=-.00053 k(m)= 2.00392
D( 2017070786 )= 3220983 Sp(m)= 3220183.93709 Δ(m)=-.00025 k(m)= 1.00513
D( 2017070788 )= 3215770 Sp(m)= 3215404.44649 Δ(m)=-.00011 k(m)= 1.00364
D( 2017070790 )= 10738234  Sp(m)= 10740205.33917  Δ(m)= .00018 k(m)= 3.35238
D( 2017070792 )= 3612612 Sp(m)= 3614043.93118 Δ(m)= .0004 k(m)= 1.12807
D( 2017070794 )= 3205126 Sp(m)= 3203754.43991 Δ(m)=-.00043 k(m)= 1
D( 2017070796 )= 6990110 Sp(m)= 6990009.69401 Δ(m)=-.00001 k(m)= 2.18182
D( 2017070798 )= 3234670 Sp(m)= 3236115.60229 Δ(m)= .00045 k(m)= 1.0101
--------------------------------------------------------------------------------
相对误差的统计计算：
2017070700 - 2017070798 :n= 50 ，μ= .00006 ，σx = .00029 ，δmin =-.00053 ，δmax = .0007

重生888 · 发表于 2017-7-7 16:01

我知道了，谢谢！

愚工688 · 发表于 2017-7-7 20:34

偶数M的表法数值的变化是有规律性的，因此也是能够比较精确的进行计算的。

我已经从3日到7日的日期×100为随机数，计算了5天随机数的连续50个偶数的相对误差情况，想必大家对我能够高精度的计算偶数的表为两个素数和的表法数的说法应该不会再怀疑了吧！

再补发一下1日、2日为随机数的连续50个偶数的相对误差数据，凑满连续的一个星期的数据就告一段落吧！

补以 2017-07-01 日的2017070100起连续的50个偶数的表法数计算与计算值相对误差的统计：

D( 2017070100 )= 8973871 Sp(m)= 8972599.55318  Δ(m)=-.00014 k(m)= 2.80065
D( 2017070102 )= 3201173 Sp(m)= 3203753.3408 Δ(m)= .00081 k(m)= 1
D( 2017070104 )= 3848752 Sp(m)= 3849725.01656  Δ(m)= .00025 k(m)= 1.20163
D( 2017070106 )= 6408681 Sp(m)= 6407506.6943 Δ(m)=-.00018 k(m)= 2
D( 2017070108 )= 3253276 Sp(m)= 3253041.86341  Δ(m)=-7.00D-05 k(m)= 1.01538
D( 2017070110 )= 4748419 Sp(m)= 4748528.05946  Δ(m)= .00002 k(m)= 1.48218
D( 2017070112 )= 6572315 Sp(m)= 6571801.75729  Δ(m)=-8.00D-05 k(m)= 2.05128
D( 2017070114 )= 3246670 Sp(m)= 3248159.85723  Δ(m)= .00046 k(m)= 1.01386
D( 2017070116 )= 3215494 Sp(m)= 3214121.49689  Δ(m)=-.00043 k(m)= 1.00324
D( 2017070118 )= 8248440 Sp(m)= 8249207.58167  Δ(m)= 9.00D-05 k(m)= 2.57486
D( 2017070120 )= 4657300 Sp(m)= 4660004.90092  Δ(m)= .00058 k(m)= 1.45455
D( 2017070122 )= 3203687 Sp(m)= 3203966.3851 Δ(m)= 9.00D-05 k(m)= 1.00007
D( 2017070124 )= 6408389 Sp(m)= 6407506.75148  Δ(m)=-.00014 k(m)= 2
D( 2017070126 )= 3203439 Sp(m)= 3203955.22585  Δ(m)= .00016 k(m)= 1.00006
D( 2017070128 )= 3416892 Sp(m)= 3417336.9409 Δ(m)= .00013 k(m)= 1.06667
D( 2017070130 )= 8540183 Sp(m)= 8543342.36072  Δ(m)= .00037 k(m)= 2.66667
D( 2017070132 )= 4270972 Sp(m)= 4271671.18459  Δ(m)= .00016 k(m)= 1.33333
D( 2017070134 )= 3419833 Sp(m)= 3419347.14927  Δ(m)=-.00014 k(m)= 1.06729
D( 2017070136 )= 6474675 Sp(m)= 6474277.5195 Δ(m)=-.00006 k(m)= 2.02084
D( 2017070138 )= 3217950 Sp(m)= 3219696.47384  Δ(m)= .00054 k(m)= 1.00498
D( 2017070140 )= 4269436 Sp(m)= 4271671.20154  Δ(m)= .00052 k(m)= 1.33333
D( 2017070142 )= 6601302 Sp(m)= 6600670.23527  Δ(m)=-.0001 k(m)= 2.06029
D( 2017070144 )= 3305581 Sp(m)= 3306476.34008  Δ(m)= .00027 k(m)= 1.03206
D( 2017070146 )= 4394433 Sp(m)= 4393718.96322  Δ(m)=-.00016 k(m)= 1.37143
D( 2017070148 )= 6409545 Sp(m)= 6409187.14761  Δ(m)=-.00006 k(m)= 2.00052
D( 2017070150 )= 4324155 Sp(m)= 4324407.90448  Δ(m)= .00006 k(m)= 1.34979
D( 2017070152 )= 3204654 Sp(m)= 3204352.36458  Δ(m)=-9.00D-05 k(m)= 1.00019
D( 2017070154 )= 7288672 Sp(m)= 7288530.51179  Δ(m)=-.00002 k(m)= 2.275
D( 2017070156 )= 3204307 Sp(m)= 3204042.78819  Δ(m)=-8.00D-05 k(m)= 1.00009
D( 2017070158 )= 3211786 Sp(m)= 3209631.8764 Δ(m)=-.00067 k(m)= 1.00183
D( 2017070160 )= 10248857  Sp(m)= 10252010.98534 Δ(m)= .00031 k(m)= 3.2
D( 2017070162 )= 3416825 Sp(m)= 3417336.9985 Δ(m)= .00015 k(m)= 1.06667
D( 2017070164 )= 3259858 Sp(m)= 3260809.55282  Δ(m)= .00029 k(m)= 1.01781
D( 2017070166 )= 6408715 Sp(m)= 6407506.8849 Δ(m)=-.00019 k(m)= 2
D( 2017070168 )= 3212080 Sp(m)= 3212530.85232  Δ(m)= .00014 k(m)= 1.00274
D( 2017070170 )= 4470253 Sp(m)= 4470801.80524  Δ(m)= .00012 k(m)= 1.39549
D( 2017070172 )= 7403427 Sp(m)= 7401804.46569  Δ(m)=-.00022 k(m)= 2.31035
D( 2017070174 )= 3845911 Sp(m)= 3844504.14619  Δ(m)=-.00037 k(m)= 1.2
D( 2017070176 )= 3690733 Sp(m)= 3691567.77092  Δ(m)= .00023 k(m)= 1.15226
D( 2017070178 )= 6407147 Sp(m)= 6407506.92302  Δ(m)= .00006 k(m)= 2
D( 2017070180 )= 4450150 Sp(m)= 4451890.88903  Δ(m)= .00039 k(m)= 1.38959
D( 2017070182 )= 3207112 Sp(m)= 3207414.9004 Δ(m)= 9.00D-05 k(m)= 1.00114
D( 2017070184 )= 6532635 Sp(m)= 6533144.3331 Δ(m)= 8.00D-05 k(m)= 2.03922
D( 2017070186 )= 3278701 Sp(m)= 3279522.35337  Δ(m)= .00025 k(m)= 1.02365
D( 2017070188 )= 3845939 Sp(m)= 3844504.17287  Δ(m)=-.00037 k(m)= 1.2
D( 2017070190 )= 8676239 Sp(m)= 8676726.28792  Δ(m)= .00006 k(m)= 2.7083
D( 2017070192 )= 3390133 Sp(m)= 3390215.32671  Δ(m)= .00002 k(m)= 1.0582
D( 2017070194 )= 3286305 Sp(m)= 3285901.01223  Δ(m)=-.00012 k(m)= 1.02564
D( 2017070196 )= 6836680 Sp(m)= 6834674.11221  Δ(m)=-.00029 k(m)= 2.13333
D( 2017070198 )= 3887660 Sp(m)= 3887013.27958  Δ(m)=-.00017 k(m)= 1.21327
-------------------------------------------------------------------------------
表法数计算值 Sp(m)的相对误差统计计算数据：
2017070100 - 2017070198 : n= 50 ，μ= .00005，σx = .00028 ，δmin =-.00067 ，δmax = .00081

重生888 · 发表于 2017-7-8 08:56

愚工先生计算很精确，比我的好，我的粗糙。但也比较好地反映了波动情况。如：尾数是10.  20的偶数，有两种加法，定为2. 而2. 4 .8. 14. 16. 22. 26. 28的尾数是有一种对称重复，所以是1.5；
104：有两种加法 1.       7 37 67 97
                                 97  67 37 7    有对称重复，取0.5

                     2.       13.  43.  73. 103
                                 91.  61.  31          没有对称重复，取1.    合起来1.5

还有您的D（100066）=590          100066=2*50033    那么K（m）该怎么算，谢谢！

愚工688 · 发表于 2017-7-8 10:20

若偶数M没有含有小于√（M-2)的奇素数，则K（m）=1;

M=?  100066
50033 + 50033 , 49943 + 50123 , 49937 + 50129 , …… ,  47 + 100019 , 23 + 100043 , 17 + 100049 ,
S( 100066 )= 590  ，Sp(m)= 615.6717  ,δ(m)= .044 ,K(m)= 1 ,r= 313
- Sp( 100066)=[( 100066/2- 2)/2]*( 1/ 3)*( 3/ 5)*( 5/ 7)*( 9/ 11)*( 11/ 13)*( 15/ 17)*( 17/ 19)*( 21/ 23)*( 27/ 29)*( 29/ 31)*( 35/ 37)*( 39/ 41)*( 41/ 43)*( 45/ 47)*( 51/ 53)*( 57/ 59)*( 59/ 61)*( 65/ 67)*( 69/ 71)*( 71/ 73)*( 77/ 79)*( 81/ 83)*( 87/ 89)*( 95/ 97)*( 99/ 101)*( 101/ 103)*( 105/ 107)*( 107/ 109)*( 111/ 113)*( 125/ 127)*( 129/ 131)*( 135/ 137)*( 137/ 139)*( 147/ 149)*( 149/ 151)*( 155/ 157)*( 161/ 163)*( 165/ 167)*( 171/ 173)*( 177/ 179)*( 179/ 181)*( 189/ 191)*( 191/ 193)*( 195/ 197)*( 197/ 199)*( 209/ 211)*( 221/ 223)*( 225/ 227)*( 227/ 229)*( 231/ 233)*( 237/ 239)*( 239/ 241)*( 249/ 251)*( 255/ 257)*( 261/ 263)*( 267/ 269)*( 269/ 271)*( 275/ 277)*( 279/ 281)*( 281/ 283)*( 291/ 293)*( 305/ 307)*( 309/ 311)*( 311/ 313)= 615.6717

愚工688 · 发表于 2017-7-8 10:25

偶数M的表为两个素数和的表法数值的变化是有规律性的，因此也是能够比较精确的进行计算的。

补发 2017-07-02 日数据，偶数20170702×100开始的连续50个偶数的表法数计算值与相对误差的统计

D( 2017070200 )= 4273103 Sp(m)= 4271671.32860 Δ(m)=-.00034 k(m)= 1.33333
D( 2017070202 )= 7688947 Sp(m)= 7689008.39911 Δ(m)= .00001 k(m)= 2.4
D( 2017070204 )= 3209287 Sp(m)= 3209119.92408 Δ(m)=-.00005 k(m)= 1.00168
D( 2017070206 )= 3204820 Sp(m)= 3203753.50598 Δ(m)=-.00033 k(m)= 1
D( 2017070208 )= 6428239 Sp(m)= 6427757.82286 Δ(m)=-7.00D-05 k(m)= 2.00632
D( 2017070210 )= 4523481 Sp(m)= 4522946.13506 Δ(m)=-.00012 k(m)= 1.41176
D( 2017070212 )= 3204498 Sp(m)= 3203753.51551 Δ(m)=-.00023 k(m)= 1
D( 2017070214 )= 6403541 Sp(m)= 6407507.03737 Δ(m)= .00062 k(m)= 2
D( 2017070216 )= 3868906 Sp(m)= 3870657.31621 Δ(m)= .00045 k(m)= 1.20816
D( 2017070218 )= 3234710 Sp(m)= 3235473.85697 Δ(m)= .00024 k(m)= 1.0099
D( 2017070220 )= 9521889 Sp(m)= 9520311.42871 Δ(m)=-.00017 k(m)= 2.97161
D( 2017070222 )= 3202004 Sp(m)= 3203753.53139 Δ(m)= .00055 k(m)= 1
D( 2017070224 )= 3497530 Sp(m)= 3500086.96555 Δ(m)= .00073 k(m)= 1.0925
D( 2017070226 )= 6453317 Sp(m)= 6454297.56515 Δ(m)= .00015 k(m)= 2.0146
D( 2017070228 )= 3205369 Sp(m)= 3203753.54092 Δ(m)=-.0005 k(m)= 1
D( 2017070230 )= 5603050 Sp(m)= 5601240.97761 Δ(m)=-.00032 k(m)= 1.74834
D( 2017070232 )= 6450401 Sp(m)= 6447305.27527 Δ(m)=-.00048 k(m)= 2.01242
D( 2017070234 )= 3355820 Sp(m)= 3355551.51821 Δ(m)=-8.00D-05 k(m)= 1.04738
D( 2017070236 )= 3280911 Sp(m)= 3281849.51914 Δ(m)= .00029 k(m)= 1.02438
D( 2017070238 )= 6718811 Sp(m)= 6716257.13240 Δ(m)=-.00038 k(m)= 2.09637
D( 2017070240 )= 4354265 Sp(m)= 4354600.57871 Δ(m)= 8.00D-05 k(m)= 1.35922
D( 2017070242 )= 3740379 Sp(m)= 3740061.04564 Δ(m)=-9.00D-05 k(m)= 1.1674
D( 2017070244 )= 7951243 Sp(m)= 7952197.42001 Δ(m)= .00012 k(m)= 2.48215
D( 2017070246 )= 3203915 Sp(m)= 3203753.56951 Δ(m)=-.00005 k(m)= 1
D( 2017070248 )= 3392617 Sp(m)= 3392209.66520 Δ(m)=-.00012 k(m)= 1.05882
D( 2017070250 )= 9317246 Sp(m)= 9320010.40252 Δ(m)= .0003 k(m)= 2.90909
D( 2017070252 )= 3204298 Sp(m)= 3203753.57904 Δ(m)=-.00017 k(m)= 1
D( 2017070254 )= 3204615 Sp(m)= 3204589.85306 Δ(m)=-.00001 k(m)= 1.00026
D( 2017070256 )= 6419933 Sp(m)= 6419853.04203 Δ(m)=-.00001 k(m)= 2.00385
D( 2017070258 )= 3853828 Sp(m)= 3853222.00313 Δ(m)=-.00016 k(m)= 1.20272
D( 2017070260 )= 4274149 Sp(m)= 4272820.87422 Δ(m)=-.00031 k(m)= 1.33369
D( 2017070262 )= 6410565 Sp(m)= 6407507.18985 Δ(m)=-.00048 k(m)= 2
D( 2017070264 )= 3797325 Sp(m)= 3797041.30146 Δ(m)=-7.00D-05 k(m)= 1.18519
D( 2017070266 )= 3223693 Sp(m)= 3224019.60921 Δ(m)= .0001 k(m)= 1.00633
D( 2017070268 )= 6468676 Sp(m)= 6468531.08709 Δ(m)=-.00002 k(m)= 2.01905
D( 2017070270 )= 4270198 Sp(m)= 4271671.47684 Δ(m)= .00035 k(m)= 1.33333
D( 2017070272 )= 3847021 Sp(m)= 3844504.33297 Δ(m)=-.00065 k(m)= 1.2
D( 2017070274 )= 6412043 Sp(m)= 6409678.52829 Δ(m)=-.00037 k(m)= 2.00068
D( 2017070276 )= 3585681 Sp(m)= 3584619.43179 Δ(m)=-.0003 k(m)= 1.11888
D( 2017070278 )= 3203757 Sp(m)= 3203753.62034 Δ(m)= 0 k(m)= 1
D( 2017070280 )= 8544005 Sp(m)= 8543342.99604 Δ(m)=-8.00D-05 k(m)= 2.66667
D( 2017070282 )= 3227259 Sp(m)= 3228202.19770 Δ(m)= .00029 k(m)= 1.00763
D( 2017070284 )= 3356956 Sp(m)= 3357753.39620 Δ(m)= .00024 k(m)= 1.04807
D( 2017070286 )= 9201145 Sp(m)= 9201656.17472 Δ(m)= .00006 k(m)= 2.87215
D( 2017070288 )= 3290074 Sp(m)= 3290594.62511 Δ(m)= .00016 k(m)= 1.02711
D( 2017070290 )= 4364284 Sp(m)= 4364604.19101 Δ(m)= 7.00D-05 k(m)= 1.36234
D( 2017070292 )= 6652669 Sp(m)= 6652167.60192 Δ(m)=-8.00D-05 k(m)= 2.07637
D( 2017070294 )= 3204912 Sp(m)= 3204189.58991 Δ(m)=-.00023 k(m)= 1.00014
D( 2017070296 )= 3204042 Sp(m)= 3203753.64893 Δ(m)=-9.00D-05 k(m)= 1
D( 2017070298 )= 6836083 Sp(m)= 6834674.45783 Δ(m)=-.00021 k(m)= 2.13333
-----------------------------------------------------------------------------
偶数表为两个素数和的表法数计算值 Sp(m)的相对误差统计计算数据：
2017070200 - 2017070298 : n= 50 ，μ=-.00004，σx = .00029 ，δmin =-.00065 ，δmax = .00073

愚工688 · 发表于 2017-7-8 16:34

偶数M表为两个素数和的表法数值的变化是有规律性的，因此也是能够比较精确的进行计算的。

下面将计算以日期×200或更大倍率的偶数的表为两个素数和的表法数计算值 Sp(m)以及它的相对误差统计.
今天的日期是2017年7月8日，
继续以今天的日期作为随机数，计算偶数20170708×200起连续50个偶数范围两个素数和的表法数值Sp(m)以及它的相对误差的统计计算数据。

D( 4034141600 )= 8447675 Sp(m)= 8451188.22656 Δ(m)= .00042 k(m)= 1.41037
D( 4034141602 )= 5998616 Sp(m)= 5999606.30233 Δ(m)= .00017 k(m)= 1.00124
D( 4034141604 )= 11984365  Sp(m)= 11984362.09995  Δ(m)= 0       k(m)= 2
D( 4034141606 )= 6734364 Sp(m)= 6735001.77357 Δ(m)= 9.00D-05 k(m)= 1.12396
D( 4034141608 )= 7187627 Sp(m)= 7190617.2671 Δ(m)= .00042 k(m)= 1.2
D( 4034141610 )= 15981536  Sp(m)= 15979149.49036  Δ(m)=-.00015 k(m)= 2.66667
D( 4034141612 )= 5993853 Sp(m)= 5993771.76562 Δ(m)=-.00001 k(m)= 1.00027
D( 4034141614 )= 5991929 Sp(m)= 5992181.06483 Δ(m)= .00004 k(m)= 1
D( 4034141616 )= 13075631  Sp(m)= 13075553.83638  Δ(m)=-.00001 k(m)= 2.1821
D( 4034141618 )= 6704089 Sp(m)= 6704538.261    Δ(m)= 7.00D-05 k(m)= 1.11888
D( 4034141620 )= 7992034 Sp(m)= 7991202.55852 Δ(m)=-.0001 k(m)= 1.3336
D( 4034141622 )= 14583011  Sp(m)= 14583787.18388  Δ(m)= .00005 k(m)= 2.4338
D( 4034141624 )= 6038915 Sp(m)= 6041110.9812 Δ(m)= .00036 k(m)= 1.00817
D( 4034141626 )= 6351098 Sp(m)= 6349678.4762 Δ(m)=-.00022 k(m)= 1.05966
D( 4034141628 )= 13979403  Sp(m)= 13980726.01748  Δ(m)= 9.00D-05 k(m)= 2.33316
D( 4034141630 )= 7987203 Sp(m)= 7989574.78479 Δ(m)= .0003 k(m)= 1.33333
D( 4034141632 )= 5990060 Sp(m)= 5992181.09156 Δ(m)= .00035 k(m)= 1
D( 4034141634 )= 11989458  Sp(m)= 11984362.18907  Δ(m)=-.00043 k(m)= 2
D( 4034141636 )= 7710476 Sp(m)= 7712472.66087 Δ(m)= .00026 k(m)= 1.28709
D( 4034141638 )= 6103131 Sp(m)= 6104355.07912 Δ(m)= .0002 k(m)= 1.01872
D( 4034141640 )= 16292048  Sp(m)= 16292466.2682 Δ(m)= .00003 k(m)= 2.71895
D( 4034141642 )= 5990810 Sp(m)= 5994778.49745 Δ(m)= .00066 k(m)= 1.00043
D( 4034141644 )= 6535856 Sp(m)= 6536924.84661 Δ(m)= .00016 k(m)= 1.09091
D( 4034141646 )= 11986155  Sp(m)= 11988272.29396  Δ(m)= .00018 k(m)= 2.00065
D( 4034141648 )= 5990810 Sp(m)= 5992181.11533 Δ(m)= .00023 k(m)= 1
D( 4034141650 )= 11363326  Sp(m)= 11363135.22011  Δ(m)=-.00002 k(m)= 1.89633
D( 4034141652 )= 11996087  Sp(m)= 11997431.33986  Δ(m)= .00011 k(m)= 2.00218
D( 4034141654 )= 5989482 Sp(m)= 5992181.12424 Δ(m)= .00045 k(m)= 1
D( 4034141656 )= 6015328 Sp(m)= 6013658.47892 Δ(m)=-.00028 k(m)= 1.00358
D( 4034141658 )= 11981415  Sp(m)= 11984362.26037  Δ(m)= .00025 k(m)= 2
D( 4034141660 )= 8003477 Sp(m)= 8003469.75698 Δ(m)= 0       k(m)= 1.33565
D( 4034141662 )= 6200262 Sp(m)= 6198808.07185 Δ(m)=-.00023 k(m)= 1.03448
D( 4034141664 )= 14735203  Sp(m)= 14731996.55661  Δ(m)=-.00022 k(m)= 2.45854
D( 4034141666 )= 6057248 Sp(m)= 6056400.52591 Δ(m)=-.00014 k(m)= 1.01072
D( 4034141668 )= 5991105 Sp(m)= 5992181.14504 Δ(m)= .00018 k(m)= 1
D( 4034141670 )= 17859161  Sp(m)= 17862148.71282  Δ(m)= .00017 k(m)= 2.98091
D( 4034141672 )= 7261841 Sp(m)= 7259251.5318 Δ(m)=-.00036 k(m)= 1.21145
D( 4034141674 )= 6382117 Sp(m)= 6381765.68718 Δ(m)=-.00006 k(m)= 1.06502
D( 4034141676 )= 11997994  Sp(m)= 11999906.23383  Δ(m)= .00016 k(m)= 2.00259
D( 4034141678 )= 7242399 Sp(m)= 7243274.86173 Δ(m)= .00012 k(m)= 1.20879
D( 4034141680 )= 7988597 Sp(m)= 7989574.88382 Δ(m)= .00012 k(m)= 1.33333
D( 4034141682 )= 11989390  Sp(m)= 11990774.51055  Δ(m)= .00012 k(m)= 2.00107
D( 4034141684 )= 6391068 Sp(m)= 6391659.91339 Δ(m)= 9.00D-05 k(m)= 1.06667
D( 4034141686 )= 6217220 Sp(m)= 6217011.63806 Δ(m)=-.00003 k(m)= 1.03752
D( 4034141688 )= 11985154  Sp(m)= 11986494.55026  Δ(m)= .00011 k(m)= 2.00036
D( 4034141690 )= 8163308 Sp(m)= 8165043.00361 Δ(m)= .00021 k(m)= 1.36262
D( 4034141692 )= 7191369 Sp(m)= 7190617.41682 Δ(m)=-.0001 k(m)= 1.2
D( 4034141694 )= 13612965  Sp(m)= 13612213.78534  Δ(m)=-.00006 k(m)= 2.27166
D( 4034141696 )= 6535171 Sp(m)= 6536924.93087 Δ(m)= .00027 k(m)= 1.09091
D( 4034141698 )= 5992169 Sp(m)= 5992409.83772 Δ(m)= .00004 k(m)= 1.00004
----------------------------------------------------------------------------------
表法数计算值的相对误差的统计计算：
4034141600 - 4034141698 :n= 50 ，μ= 0.00008，σx = .00021，δmin =-.00043 ，δmax = .00066

重生888 · 发表于 2017-7-8 17:31

谢谢网友，我的疑问基本没了。为您点赞，并顶起来！另外，这么大的数字实际素数对，是您程序做出来的，还是有数据库？望告知，谢谢！

愚工688 · 发表于 2017-7-8 22:32

本帖最后由愚工688 于 2017-7-9 02:23 编辑

重生888 发表于 2017-7-8 09:31
谢谢网友，我的疑问基本没了。为您点赞，并顶起来！另外，这么大的数字实际素数对，是您程序做出来的，还是 ...

没有数据库，即时计算的。

筛选大偶数的素数对真值使用的是网友黄博士给我的Fast Gn软件，
计算表法数与相对误差统计是我自己编的软件。
两者的结合就成为大偶数的表法数的高精度的计算与相对误差统计。（没有真值，就没有相对误差）
我自己编写的筛选素对程序，最多勉强筛选到百万级别的偶数，因为速度问题不可能处理大量的大偶数。

例如：今天日期加个0，的偶数的表法数计算，单独我的软件只能做到下面程度：（大约7秒左右）
Sp( 201707080 *)≈  574918 , Δ≈, k(m)= 1.41037
Sp( 201707082 *)≈  831258.7 , Δ≈, k(m)= 2.03922
Sp( 201707084 *)≈  407636.5 , Δ≈, k(m)= 1
Sp( 201707086 *)≈  490026.5 , Δ≈, k(m)= 1.20212
Sp( 201707088 *)≈  905858.9 , Δ≈, k(m)= 2.22222
Sp( 201707090 *)≈  592925.8 , Δ≈, k(m)= 1.45455
Sp( 201707092 *)≈  407985.2 , Δ≈, k(m)= 1.00086
start time :22:39:04, end time:22:39:11use time :0:0:07，
因为没有真值，就不能同时加上相对误差。故一般另外计算真值后，再手工计算相对误差。

使用软件FastGn得到：
201707080:12:2

G(201707080) = 574353
G(201707082) = 830506
G(201707084) = 407223
G(201707086) = 489912
G(201707088) = 905532
G(201707090) = 592549
G(201707092) = 408545
G(201707094) = 815169
G(201707096) = 406920
G(201707098) = 413341
G(201707100) = 1320382
G(201707102) = 427058

count = 12, algorithm = 2, working threads = 2, time use 0.036 sec

真值代入计算器计算：
Sp( 201707080 *)≈  574918 , Δ≈0.00098, k(m)= 1.41037
Sp( 201707082 *)≈  831258.7 , Δ≈0.000906, k(m)= 2.03922
Sp( 201707084 *)≈  407636.5 , Δ≈0.001015, k(m)= 1
Sp( 201707086 *)≈  490026.5 , Δ≈0.0002337, k(m)= 1.20212
Sp( 201707088 *)≈  905858.9 , Δ≈0.000361, k(m)= 2.22222
Sp( 201707090 *)≈  592925.8 , Δ≈0.0006359, k(m)= 1.45455
Sp( 201707092 *)≈  407985.2 , Δ≈-0.00137, k(m)= 1.00086

很明显，我的表法数计算式的计算大偶数的精度还是蛮高的。（小偶数没有这么高的精度，因为小偶数区域相对误差的分布范围比较大，无法预先对表法数计算值进行误差修正）

重生888 · 发表于 2017-7-9 10:14

谢谢愚工先生！您做的工作，不管别人相信不相信，反正我相信！与您交流受益匪浅！

		自动登录	找回密码
密码			注册