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

愚工688 · 发表于 2017-7-3 19:55

本帖最后由愚工688 于 2018-9-10 07:17 编辑

高精度计算大偶数表为两个素数和的表法数值的实例

要想正确的计算出任意一个偶数表为两个素数和的表法数值我认为是不可能的。
但是对于哥德巴赫猜想的偶数表为两个素数和的表法数值，用一个计算式比较高精度的计算出其近似值还是能够做到的。
火车不是靠推的，牛皮不是靠吹的。
只有实际计算一下才会使人相信我说的:用一个计算式能够比较高精度的计算出偶数表为两个素数和的表法数的近似值。

今天的日期是2017-07-03，我就以日期这个随机数，从2017070300开始，计算50个连续的偶数，看看能够达到高精度的标准吗？

计算式： Sp(m*)=(A-2)P(m) /(1+μ)
=(A-2)×P(2·3·…·n·…·r)/(1+μ)
=(A-2)×P(2)×P(3)×…×P(n)×…×P(r)/(1+μ).
      =(A-2)×(1/2)×f(3)×…×f(n)×…×f(r)/(1+μ)；                {式3}
式中：3≤ n≤r；n是素数；μ系相对误差修正值，只适用一定范围的偶数区域。
f(n)=(n-1)/n, [jn=0时]；或f(n)=(n-2)/n,  [jn>0时] 。jn系A除以n时的余数。

{式3}经过数学变形，也可以用另一种形式表达:
Sp(m*)= (A-2)/[2(1+μ)]*π[(n-2)/n]*π[(k-1)/(k-2)]       {式4}
式中：3≤ n≤r；n是素数; k是偶数M含有的奇素数因子，k≤√(M-2)。

下面是这50个偶数的表法数计算值的相对误差统计计算数据

相对误差修正值μ=0.1406，适用范围15亿-25亿。)

其中：μ——区域内相对误差的平均值；
σx——对区域内表法数计算值的相对误差的统计计算的标准偏差。

D( 2017070300 )= 5197904 Sp(m)= 5199125.27165 Δ(m)= .00023
D( 2017070302 )= 3496132 Sp(m)= 3495003.99105 Δ(m)=-.00032
D( 2017070304 )= 6631608 Sp(m)= 6628455.85166 Δ(m)=-.00048
D( 2017070306 )= 3204009 Sp(m)= 3203919.67101 Δ(m)=-.00003
D( 2017070308 )= 3570308 Sp(m)= 3570546.13454 Δ(m)= 6.999999999999999D-05 ——我设定取5位小数，但是有时电脑会出来双精度值。
D( 2017070310 )= 8546996 Sp(m)= 8545893.45422 Δ(m)=-.00013
D( 2017070312 )= 3204253 Sp(m)= 3203753.67434 Δ(m)=-.00016
D( 2017070314 )= 3845448 Sp(m)= 3845568.32353 Δ(m)= .00003
D( 2017070316 )= 6647769 Sp(m)= 6646092.72651 Δ(m)=-.00025
D( 2017070318 )= 3312769 Sp(m)= 3314560.20226 Δ(m)= .00054
D( 2017070320 )= 4273722 Sp(m)= 4273960.70664 Δ(m)= .00006
D( 2017070322 )= 6407798 Sp(m)= 6407507.38045 Δ(m)=-.00005
D( 2017070324 )= 3392833 Sp(m)= 3392209.79301 Δ(m)=-.00018
D( 2017070326 )= 3206954 Sp(m)= 3205020.5031 Δ(m)=-5.999999999999999D-04
D( 2017070328 )= 8388271 Sp(m)= 8388009.68663 Δ(m)=-.00003
D( 2017070330 )= 5083030 Sp(m)= 5083111.74364 Δ(m)= .00002
D( 2017070332 )= 3417107 Sp(m)= 3417337.28652 Δ(m)= 6.999999999999999D-05
D( 2017070334 )= 6412179 Sp(m)= 6415036.80449 Δ(m)= .00045
D( 2017070336 )= 3201717 Sp(m)= 3203753.71246 Δ(m)= 6.400000000000001D-04
D( 2017070338 )= 3205713 Sp(m)= 3205860.06259 Δ(m)= .00005
D( 2017070340 )= 8621782 Sp(m)= 8620310.30648 Δ(m)=-.00017
D( 2017070342 )= 3844574 Sp(m)= 3844504.46639 Δ(m)=-.00002
D( 2017070344 )= 3259901 Sp(m)= 3257702.09423 Δ(m)=-.00067
D( 2017070346 )= 6417775 Sp(m)= 6416062.2063 Δ(m)=-.00027
D( 2017070348 )= 3203778 Sp(m)= 3203753.73152 Δ(m)=-.00001
D( 2017070350 )= 4429418 Sp(m)= 4429881.70724 Δ(m)= .0001
D( 2017070352 )= 7117121 Sp(m)= 7119452.75083 Δ(m)= .00033
D( 2017070354 )= 3496084 Sp(m)= 3495975.50504 Δ(m)=-.00003
D( 2017070356 )= 3844742 Sp(m)= 3846416.15054 Δ(m)= .00044
D( 2017070358 )= 6572300 Sp(m)= 6571802.55878 Δ(m)=-8.000000000000001D-05
D( 2017070360 )= 4270844 Sp(m)= 4271671.66744 Δ(m)= .00019
D( 2017070362 )= 3393692 Sp(m)= 3392209.85692 Δ(m)=-.00044
D( 2017070364 )= 6426294 Sp(m)= 6426634.40197 Δ(m)= .00005
D( 2017070366 )= 3546655 Sp(m)= 3546986.98429 Δ(m)= 9.000000000000001D-05
D( 2017070368 )= 3258242 Sp(m)= 3257943.82352 Δ(m)=-9.000000000000001D-05
D( 2017070370 )= 10275721 Sp(m)= 10274696.12092 Δ(m)=-.0001
D( 2017070372 )= 3204671 Sp(m)= 3204533.46148 Δ(m)=-.00004
D( 2017070374 )= 3580505 Sp(m)= 3582692.39111 Δ(m)= .00061
D( 2017070376 )= 6819858 Sp(m)= 6815898.14322 Δ(m)=-.00058
D( 2017070378 )= 3206183 Sp(m)= 3208720.83929 Δ(m)= .00079
D( 2017070380 )= 4710389 Sp(m)= 4709058.19109 Δ(m)=-.00028
D( 2017070382 )= 6543452 Sp(m)= 6544205.35055 Δ(m)= .00012
D( 2017070384 )= 3844154 Sp(m)= 3844504.54644 Δ(m)= 9.000000000000001D-05
D( 2017070386 )= 3205531 Sp(m)= 3204022.45048 Δ(m)=-.00047
D( 2017070388 )= 6458425 Sp(m)= 6458767.65083 Δ(m)= .00005
D( 2017070390 )= 4269097 Sp(m)= 4271671.73098 Δ(m)= 5.999999999999999D-04
D( 2017070392 )= 3313703 Sp(m)= 3312284.30389 Δ(m)=-.00043
D( 2017070394 )= 6472799 Sp(m)= 6472229.90825 Δ(m)=-9.000000000000001D-05
D( 2017070396 )= 3630533 Sp(m)= 3629525.01097 Δ(m)=-.00028
D( 2017070398 )= 3862515 Sp(m)= 3864845.86716 Δ(m)= 5.999999999999999D-04
------------------------------------------------------------------------------------------------------------------
2017070300 - 2017070398 : n= 50 ,μ= 0 ,σx= .00033 ,δmin=-.00067 ,δmax= .00079

偶数的表法数值的变化是有规律性的，因此也是能够比较精确的进行计算的。
看了我的偶数表为两个素数和的表法数计算值数据后，还会有人怀疑吗？

（注：把题目中的今天改为当天。2018-9-10）

愚工688 · 发表于 2026-9-20 20:10

今天日期的十倍、百倍的连续偶数的连乘式计算实例

Sp( 202609200 *)≈  1101030.8 , jd ≈ 0.99999；
Sp( 202609202 *)≈  408799.3 , jd ≈ 0.99992；
Sp( 202609204 *)≈  490559.1 , jd ≈ 0.99876；
Sp( 202609206 *)≈  850643 , jd ≈ 0.99946；
Sp( 202609208 *)≈  445211.3 , jd ≈ 0.99966；
start time =19:47:38,end time=19:47:44 ,time use =
202609200:6:2
真值：
G(202609200) = 1101037
G(202609202) = 408766
G(202609204) = 491170
G(202609206) = 851102
G(202609208) = 445362

1+1真值：
2026092000:6:2

G(2026092000) = 8662485
G(2026092002) = 3576466
G(2026092004) = 3216723
G(2026092006) = 6433739
G(2026092008) = 3216472
G(2026092010) = 4289770

Sp( 2026092000 *)≈  8659704.2 , jd ≈ 0.99968;
Sp( 2026092002 *)≈  3574976.2 , jd ≈ 0.99958;
Sp( 2026092004 *)≈  3215242.3 , jd ≈ 0.99954;
Sp( 2026092006 *)≈  6430484.6 , jd ≈ 0.99949;
Sp( 2026092008 *)≈  3215857.0 , jd ≈ 0.99981;
start time =19:48:00,end time=19:48:22 ,time use =

连乘式计算式
Sp( 202609200 *) = 1/(1+ .1264 )*( 202609200 /2 -2)*p(m) ≈ 1101030.8 , k(m)= 2.693329
Sp( 202609202 *) = 1/(1+ .1264 )*( 202609202 /2 -2)*p(m) ≈ 408799.3 , k(m)= 1
Sp( 202609204 *) = 1/(1+ .1264 )*( 202609204 /2 -2)*p(m) ≈ 490559.1 , k(m)= 1.2
Sp( 202609206 *) = 1/(1+ .1264 )*( 202609206 /2 -2)*p(m) ≈ 850643 , k(m)= 2.080833
Sp( 202609208 *) = 1/(1+ .1264 )*( 202609208 /2 -2)*p(m) ≈ 445211.3 , k(m)= 1.089071
Sp( 2026092000 *) = 1/(1+ .1411 )*( 2026092000 /2 -2)*p(m) ≈ 8659704.2 , k(m)= 2.693329
Sp( 2026092002 *) = 1/(1+ .1411 )*( 2026092002 /2 -2)*p(m) ≈ 3574976.2 , k(m)= 1.111884
Sp( 2026092004 *) = 1/(1+ .1411 )*( 2026092004 /2 -2)*p(m) ≈ 3215242.3 , k(m)= 1
Sp( 2026092006 *) = 1/(1+ .1411 )*( 2026092006 /2 -2)*p(m) ≈ 6430484.6 , k(m)= 2
Sp( 2026092008 *) = 1/(1+ .1411 )*( 2026092008 /2 -2)*p(m) ≈ 3215857 , k(m)= 1.000191

yangchuanju · 发表于 2026-9-11 05:52

顶起来，以便收集整理一下愚公的一些大偶数的哥猜数据；
也请吴代业温习温习您与愚公过去的交流贴！

愚工688 · 发表于 2023-6-16 21:20

重生888@ 发表于 2022-8-24 07:16
再次欣赏愚工先生高精度计算！

以今天日期的千倍的计算实例：

G(20230616000) = 44170095 ；inf( 20230616000 )≈  44139763.2 , jd ≈0.99931 , k(m)= 1.70364
G(20230616002) = 25928828 ；inf( 20230616002 )≈  25909144.1 , jd ≈0.99924 , k(m)= 1
G(20230616004) = 52021809 ；inf( 20230616004 )≈  51988184.3 , jd ≈0.99935 , k(m)= 2.00656
G(20230616006) = 26550802 ；inf( 20230616006 )≈  26541074.5 , jd ≈0.99963 , k(m)= 1.02439
G(20230616008) = 26436618 ；inf( 20230616008 )≈  26417166.6 , jd ≈0.99926 , k(m)= 1.01961
G(20230616010) = 77386937 ；inf( 20230616010 )≈  77336485.2 , jd ≈0.99935 , k(m)= 2.98491
G(20230616012) = 27672308 ；inf( 20230616012 )≈  27653511.6 , jd ≈0.99932 , k(m)= 1.06733

time start =20:57:28  ,time end =21:00:01 ,time use =

计算式：
inf( 20230616000 ) = 1/(1+ .1535 )*( 20230616000 /2 -2)*p(m) ≈ 44139763.2
inf( 20230616002 ) = 1/(1+ .1535 )*( 20230616002 /2 -2)*p(m) ≈ 25909144.1
inf( 20230616004 ) = 1/(1+ .1535 )*( 20230616004 /2 -2)*p(m) ≈ 51988184.3
inf( 20230616006 ) = 1/(1+ .1535 )*( 20230616006 /2 -2)*p(m) ≈ 26541074.5
inf( 20230616008 ) = 1/(1+ .1535 )*( 20230616008 /2 -2)*p(m) ≈ 26417166.6
inf( 20230616010 ) = 1/(1+ .1535 )*( 20230616010 /2 -2)*p(m) ≈ 77336485.2
inf( 20230616012 ) = 1/(1+ .1535 )*( 20230616012 /2 -2)*p(m) ≈ 27653511.6

重生888@ · 发表于 2022-8-24 15:16

再次欣赏愚工先生高精度计算！

愚工688 · 发表于 2022-8-22 12:35

yangchuanju 发表于 2022-8-21 02:45
愚公688老师一再强调：
高精度哥猜数的计算可采用下述公式：
Xi(M)=t2*c1*M/(logM)^2，（实则为Xi(M)=t ...

确实我使用类似哈代计算式时对拉曼扭杨系数只计算√M内的素数。
为什么呢，
1，计算M内的素数比只计算√M内的素数肯定需要更多的计算时间，而计算值相差很小，对计算式原来存在的相对误差几乎没有影响。
2,从素数筛选的埃拉托色尼筛法看，对筛选起作用的只是√M内的素数，大于√M内的素数不起筛选作用，因此哈代公式中计算到M内的全部素数没有必要。
3，使用与否哈李常数c=0.660161815…并不能提高计算值的精度，即使我电脑上运行哈李计算式时显示的极限是0.6601667也对计算值的精度没有什么影响。（运行的Basic程序是16位的，中精度显示是8位数）

因此我们采用计算式的目的是什么？
无非是得到比较高精度的计算值，那么为什么非得采用哈李常数c的值呢？

以今天日期的百倍为随机偶数计算的素对计算值：

   Xi(M)=t2*c1*M/(logM)^2
式中：t2=1.358-(log(M))^(.5)*0.05484;c1:只计算√M内素数的类似拉曼扭杨系数。



  G(2022082200) = 8563394  ;Xi(M)≈ 8562138.11 δxi(M)≈?-0.0001467；
  G(2022082202) = 3884529  ;Xi(M)≈ 3882829.96 δxi(M)≈?-0.0004374；
  G(2022082204) = 3434962  ;Xi(M)≈ 3433844.24 δxi(M)≈?-0.0003255；
  G(2022082206) = 6470850  ;Xi(M)≈ 6468476.5 δxi(M)≈?-0.0003667；
  G(2022082208) = 3224856  ;Xi(M)≈ 3226164.4 δxi(M)≈? 0.0004056；
  G(2022082210) = 4280169  ;Xi(M)≈ 4281069.08 δxi(M)≈? 0.0002103；
  G(2022082212) = 6425613  ;Xi(M)≈ 6425619.67 δxi(M)≈? 0.0000010；
  G(2022082214) = 3371151  ;Xi(M)≈ 3370748.72 δxi(M)≈?-0.0001192；
  G(2022082216) = 4326647  ;Xi(M)≈ 4325243.79 δxi(M)≈?-0.0003243；
  G(2022082218) = 7036101  ;Xi(M)≈ 7037391.62 δxi(M)≈? 0.0001835；
  G(2022082220) = 4610747  ;Xi(M)≈ 4612421.13 δxi(M)≈? 0.0003631；
  G(2022082222) = 3289104  ;Xi(M)≈ 3288532.62 δxi(M)≈?-0.0001737；
  G(2022082224) = 6578717  ;Xi(M)≈ 6578227.72 δxi(M)≈?-0.0000743；
  G(2022082226) = 3210588  ;Xi(M)≈ 3210801.76 δxi(M)≈? 0.0000667；
  G(2022082228) = 3219771  ;Xi(M)≈ 3220386.26 δxi(M)≈? 0.0001910；
  time start =12:06:44, time end =12:07:29

重生888@ · 发表于 2022-8-22 09:31

愚工688 发表于 2022-8-20 20:49
以今天日期的百倍为随机偶数计算偶数的素对数量的下界计算值与计算精度：

inf( 2022082000 )≈ 4537115 ...

G（2022082000）=4567849
D（2022082000）=4232402             D/G=0.926563       （由此知道有小因子）验证是17：
D=4232402*16/15=4514562             D/G=0.988333

4232402/2=2116201    （一份素对）
G（202208220）=10274717
D（202208220）=2116201*4=8464804       由误差过大，只有小因子，验证是7：
D（202208220）=8464804*6/5=10157764       D/G=0.988617

yangchuanju · 发表于 2022-8-21 10:45

愚工688 发表于 2021-11-27 15:54
不用说哈-李素对计算式的计算精度不高，那时你没有掌握哈-李素对计算式的计算的相对误差的变化的规律，不要 ...

愚公688老师一再强调：
高精度哥猜数的计算可采用下述公式：
Xi(M)=t2*c1*M/(logM)^2，（实则为Xi(M)=t2*c1*M/ln(M)^2）
相对误差动态修正系数 t2=1.358-log(M)^(0.5)*.05484;
C1--类似拉曼扭杨系数，略作改进；（只计算√M内的素数以提高计算速度）

公认的哈李对数计算式是：R2=2c*M/ln(M)^2*K，式中K——波动因子∏(p-2)/(p-1)，3≤p≥√M，p|M；
对于单计哥猜数计算式应改为：R1=c*M/ln(M)^2*K。
哈李对数计算式，适用于偶数无穷大；当偶数为有限数值时，计算值小于哥猜数真实值。
愚公计算式中的增加了一个修正系数t2，为的是使其计算值接近于哥猜数真实值；
愚公计算式中“类似拉曼扭杨系数”C1实则为哈李计算式中的c*K，至少笔者这样认为并一直如此处理涉及愚公的一些计算数据；
然而愚公又声称“只计算√M内的素数”，看来愚公的C1中的c并不是哈李常数c=0.660161815…；
查阅有关资料可知，哈李常数c是当p趋近于无穷大时连乘积∏[1-1/(p-1)^2]的极限值，
愚公计算式中的C1中的c只取到M平方根内的最大素数：
当p=3时，∏[1-1/(p-1)^2]=0.75；
当p=5时，∏[1-1/(p-1)^2]=0.703125；
当p=7时，∏[1-1/(p-1)^2]=0.68359375；
……
当p=97时，∏[1-1/(p-1)^2]=0.661377084547185；
当p=997时，∏[1-1/(p-1)^2]=0.6602457439708；
当p=9973时，∏[1-1/(p-1)^2]=0.660168296505504；
当p=99991时，∏[1-1/(p-1)^2]=0.660162345466730；
当p=999983时，∏[1-1/(p-1)^2]=0.660161860589839；
当p=9999991时，∏[1-1/(p-1)^2]=0.660161819715374；
当p=10570841（第100万个素数）时，∏[1-1/(p-1)^2]=0.660161819494655；
更多更大的∏[1-1/(p-1)^2]值笔者没有计算，相信连乘积的数值越来越小，并最终趋近于哈李常数c的值。

愚公声称如此处理是为了“提高计算速度”，看来使用变数c不如使用常数c计算速度更快；
要说是为了“提高计算精度”，倒是说得通，因为愚公的变数c要比哈李常数c大一些。
究竟是什么原因，请愚公老师给予解答！

愚工688 · 发表于 2022-8-20 20:49

以今天日期的百倍为随机偶数计算偶数的素对数量的下界计算值与计算精度：

inf( 2022082000 )≈  4537115.3 , jd ≈0.99327 ,infS(m) = 3190159.21 , k(m)= 1.42222
inf( 2022082002 )≈  6380318.4 , jd ≈0.99385 ,infS(m) = 3190159.22 , k(m)= 2
inf( 2022082004 )≈  3190159.2 , jd ≈0.99323 ,infS(m) = 3190159.22 , k(m)= 1
inf( 2022082006 )≈  3828191.1 , jd ≈0.99373 ,infS(m) = 3190159.22 , k(m)= 1.2
inf( 2022082008 )≈  6380318.5 , jd ≈0.99320 ,infS(m) = 3190159.22 , k(m)= 2
inf( 2022082010 )≈  4644414.1 , jd ≈0.99389 ,infS(m) = 3190159.23 , k(m)= 1.45586
inf( 2022082012 )≈  3300164.7 , jd ≈0.99308 ,infS(m) = 3190159.23 , k(m)= 1.03448
inf( 2022082014 )≈  6683461.2 , jd ≈0.99387 ,infS(m) = 3190159.23 , k(m)= 2.09502
inf( 2022082016 )≈  3190159.2 , jd ≈0.99364 ,infS(m) = 3190159.24 , k(m)= 1
inf( 2022082018 )≈  3549291.5 , jd ≈0.99312 ,infS(m) = 3190159.24 , k(m)= 1.11258
inf( 2022082020 )≈ 10208509.6 , jd ≈0.99356 ,infS(m) = 3190159.24 , k(m)= 3.2
inf( 2022082022 )≈  3229543.9 , jd ≈0.99320 ,infS(m) = 3190159.25 , k(m)= 1.01235
time start =20:25:04  ,time end =20:25:56 ,time use =

2022082000:12:2
素数对真值：
G(2022082000) = 4567849
G(2022082002) = 6419774
G(2022082004) = 3211911
G(2022082006) = 3852350
G(2022082008) = 6424014
G(2022082010) = 4672952
G(2022082012) = 3323164
G(2022082014) = 6724711
G(2022082016) = 3210573
G(2022082018) = 3573866
G(2022082020) = 10274717
G(2022082022) = 3251646

计算式：
inf( 2022082000 ) = 1/(1+ .148 )*( 2022082000 /2 -2)*p(m) ≈ 4537115.3
inf( 2022082002 ) = 1/(1+ .148 )*( 2022082002 /2 -2)*p(m) ≈ 6380318.4
inf( 2022082004 ) = 1/(1+ .148 )*( 2022082004 /2 -2)*p(m) ≈ 3190159.2
inf( 2022082006 ) = 1/(1+ .148 )*( 2022082006 /2 -2)*p(m) ≈ 3828191.1
inf( 2022082008 ) = 1/(1+ .148 )*( 2022082008 /2 -2)*p(m) ≈ 6380318.5
inf( 2022082010 ) = 1/(1+ .148 )*( 2022082010 /2 -2)*p(m) ≈ 4644414.1
inf( 2022082012 ) = 1/(1+ .148 )*( 2022082012 /2 -2)*p(m) ≈ 3300164.7
inf( 2022082014 ) = 1/(1+ .148 )*( 2022082014 /2 -2)*p(m) ≈ 6683461.2
inf( 2022082016 ) = 1/(1+ .148 )*( 2022082016 /2 -2)*p(m) ≈ 3190159.2
inf( 2022082018 ) = 1/(1+ .148 )*( 2022082018 /2 -2)*p(m) ≈ 3549291.5
inf( 2022082020 ) = 1/(1+ .148 )*( 2022082020 /2 -2)*p(m) ≈ 10208509.6
inf( 2022082022 ) = 1/(1+ .148 )*( 2022082022 /2 -2)*p(m) ≈ 3229543.9

愚工688 · 发表于 2022-5-15 16:57

使用连乘式计算3500亿的大偶数的素对数量，是比较化时间的：

G(350000000000) = 566240377；
Sp( 350000000000 *)≈  565994611.8 , Δ≈-0.000434, k(m)= 1.6
G(350000000002) = 353889363；
Sp( 350000000002 *)≈  353746632.4 , Δ≈-0.000403, k(m)= 1
G(350000000004) = 717784873；
Sp( 350000000004 *)≈  717457958.6 , Δ≈-0.000455, k(m)= 2.02817
G(350000000006) = 362950888；
Sp( 350000000006 *)≈  362817058.8 , Δ≈-0.000369, k(m)= 1.02564
G(350000000008) = 374985721；
Sp( 350000000008 *)≈  374813393.7 , Δ≈-0.000460, k(m)= 1.05955
G(350000000010) = 972138855
Sp( 350000000010 *)≈  971724198.1 , Δ≈-0.000427, k(m)= 2.74695
G(350000000012) = 395467221；
Sp( 350000000012 *)≈  395317744.4 , Δ≈-0.000378, k(m)= 1.11752
G(350000000014) = 424688917；
Sp( 350000000014 *)≈  424520117.7 , Δ≈-0.000397, k(m)= 1.20007
G(350000000016) = 707801727；
Sp( 350000000016 *)≈  707493264.7 , Δ≈-0.000436, k(m)= 2
G(350000000018) = 353914007；
Sp( 350000000018 *)≈  353748339 , Δ≈-0.000468, k(m)= 1
G(350000000020) = 539549822；
Sp( 350000000020 *)≈  539306176.5 , Δ≈-0.000452, k(m)= 1.52455
G(350000000022) = 734215313；
Sp( 350000000022 *)≈  733903026 , Δ≈-0.000425, k(m)= 2.07466
start time :09:52:32, end time:10:21:59use time :

计算式：（因子*p(m)的展开即是素数连乘式）
Sp( 350000000000 *) = 1/(1+ .16544 )*( 350000000000 /2 -2)*p(m) ≈ 565994611.8 , k(m)= 1.6
Sp( 350000000002 *) = 1/(1+ .16544 )*( 350000000002 /2 -2)*p(m) ≈ 353746632.4 , k(m)= 1
Sp( 350000000004 *) = 1/(1+ .16544 )*( 350000000004 /2 -2)*p(m) ≈ 717457958.6 , k(m)= 2.02817
Sp( 350000000006 *) = 1/(1+ .16544 )*( 350000000006 /2 -2)*p(m) ≈ 362817058.8 , k(m)= 1.02564
Sp( 350000000008 *) = 1/(1+ .16544 )*( 350000000008 /2 -2)*p(m) ≈ 374813393.7 , k(m)= 1.05955
Sp( 350000000010 *) = 1/(1+ .16544 )*( 350000000010 /2 -2)*p(m) ≈ 971724198.1 , k(m)= 2.74695
Sp( 350000000012 *) = 1/(1+ .16544 )*( 350000000012 /2 -2)*p(m) ≈ 395317744.4 , k(m)= 1.11752
Sp( 350000000014 *) = 1/(1+ .16544 )*( 350000000014 /2 -2)*p(m) ≈ 424520117.7 , k(m)= 1.20007
Sp( 350000000016 *) = 1/(1+ .16544 )*( 350000000016 /2 -2)*p(m) ≈ 707493264.7 , k(m)= 2
Sp( 350000000018 *) = 1/(1+ .16544 )*( 350000000018 /2 -2)*p(m) ≈ 353748339 , k(m)= 1
Sp( 350000000020 *) = 1/(1+ .16544 )*( 350000000020 /2 -2)*p(m) ≈ 539306176.5 , k(m)= 1.52455
Sp( 350000000022 *) = 1/(1+ .16544 )*( 350000000022 /2 -2)*p(m) ≈ 733903026 , k(m)= 2.07466

重生888@ · 发表于 2021-11-29 14:39

当然楼主这一次改进也很方便，方便与精度并存更好！不是道不同，而是略有差别！谢谢！

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

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

评分

点评

点评