|
楼主 |
发表于 2018-2-20 16:59
|
显示全部楼层
本帖最后由 愚工688 于 2018-2-25 04:55 编辑
偶数M表为两个素数和的数量的表法数值Sp(m)的变化是有规律的,因此是能够比较高精度的进行计算的。
今天的日期是2018年02月20日,狗年初五,继续以今天的日期作为随机数,计算百亿级别的偶数20180220×1000起的连续偶数M表为两个素数和的表法数值Sp(m)以及计算值的精度 jdz。
注:素对真值s(m)=G(M).两个不同程序的各自表示形式。
G(20180220000) = 69256004;Sp( 20180220000 *)≈ 69231573.8 , jdz =sp(m)/s(m) ≈ 0.999647,
G(20180220002) = 25860023;Sp( 20180220002 *)≈ 25855267.7 , jdz =sp(m)/s(m) ≈ 0.999816,
G(20180220004) = 31704579;Sp( 20180220004 *)≈ 31691799.7 , jdz =sp(m)/s(m) ≈ 0.999597,
G(20180220006) = 53025651;Sp( 20180220006 *)≈ 53003408.8 , jdz =sp(m)/s(m) ≈ 0.999581,
G(20180220008) = 29027902;Sp( 20180220008 *)≈ 29012205.3 , jdz =sp(m)/s(m) ≈ 0.999459,
G(20180220010) = 36640473;Sp( 20180220010 *)≈ 36623960.7 , jdz =sp(m)/s(m) ≈ 0.999542,
G(20180220012) = 52110271;Sp( 20180220012 *)≈ 52093576.4 , jdz =sp(m)/s(m) ≈ 0.999680,
G(20180220014) = 25866591;Sp( 20180220014 *)≈ 25855267.7 , jdz =sp(m)/s(m) ≈ 0.999562,
G(20180220016) = 28736696;Sp( 20180220016 *)≈ 28728075.2 , jdz =sp(m)/s(m) ≈ 0.999700,
G(20180220018) = 62082973;Sp( 20180220018 *)≈ 62052642.5 , jdz =sp(m)/s(m) ≈ 0.999511,
G(20180220020) = 34630330;Sp( 20180220020 *)≈ 34613876.0 , jdz =sp(m)/s(m) ≈ 0.999525,
G(20180220022) = 27591157;Sp( 20180220022 *)≈ 27578952.2 , jdz =sp(m)/s(m) ≈ 0.999558,
G(20180220024) = 55316810;Sp( 20180220024 *)≈ 55291936.5 , jdz =sp(m)/s(m) ≈ 0.999550,
G(20180220026) = 25869622;Sp( 20180220026 *)≈ 25855267.7 , jdz =sp(m)/s(m) ≈ 0.999445,
G(20180220028) = 25864426;Sp( 20180220028 *)≈ 25855267.7 , jdz =sp(m)/s(m) ≈ 0.999646,
G(20180220030) = 68977499;Sp( 20180220030 *)≈ 68947380.5 , jdz =sp(m)/s(m) ≈ 0.999563,
G(20180220032) = 31035882;Sp( 20180220032 *)≈ 31026321.2 , jdz =sp(m)/s(m) ≈ 0.999692,
G(20180220034) = 28589281;Sp( 20180220034 *)≈ 28576447.3 , jdz =sp(m)/s(m) ≈ 0.999551,
G(20180220036) = 54783585;Sp( 20180220036 *)≈ 54752331.6 , jdz =sp(m)/s(m) ≈ 0.999430,
G(20180220038) = 29184922;Sp( 20180220038 *)≈ 29170303.4 , jdz =sp(m)/s(m) ≈ 0.999499,
G(20180220040) = 35677013;Sp( 20180220040 *)≈ 35662438.2 , jdz =sp(m)/s(m) ≈ 0.999591,
G(20180220042) = 51760134;Sp( 20180220042 *)≈ 51740408.7 , jdz =sp(m)/s(m) ≈ 0.999619,
G(20180220044) = 25864971;Sp( 20180220044 *)≈ 25855267.7 , jdz =sp(m)/s(m) ≈ 0.999625,
G(20180220046) = 31040389;Sp( 20180220046 *)≈ 31026321.3 , jdz =sp(m)/s(m) ≈ 0.999547,
G(20180220048) = 52877916;Sp( 20180220048 *)≈ 52859658.5 , jdz =sp(m)/s(m) ≈ 0.999655,
表法数计算式摘录:
Sp( 20180220000 *) = 1/(1+ .1533 )*( 20180220000 /2 -2)*p(m) ≈ 69231573.8 ,
Sp( 20180220002 *) = 1/(1+ .1533 )*( 20180220002 /2 -2)*p(m) ≈ 25855267.7 ,
Sp( 20180220004 *) = 1/(1+ .1533 )*( 20180220004 /2 -2)*p(m) ≈ 31691799.7 ,
Sp( 20180220006 *) = 1/(1+ .1533 )*( 20180220006 /2 -2)*p(m) ≈ 53003408.8 ,
……
Sp( 20180220044 *) = 1/(1+ .1533 )*( 20180220044 /2 -2)*p(m) ≈ 25855267.7 ,
Sp( 20180220046 *) = 1/(1+ .1533 )*( 20180220046 /2 -2)*p(m) ≈ 31026321.3 ,
Sp( 20180220048 *) = 1/(1+ .1533 )*( 20180220048 /2 -2)*p(m) ≈ 52859658.5 ,
式中:
p(m)= π[(n-2)/n]*π[(p-1)/(p-2)] ; 3≤ n≤r;n是素数; p是偶数M含有的奇素数因子,p≤√(M-2)。
与294# 的偶数相同,计算值也同样具有比较高的精度。
|
|