偶数素对数量计算式——连乘式的数学理论依据

愚工688

连乘式的理论依据是概率的乘法定理。

运用到的数学定理
【相互独立事件同时发生的概率】两个相互独立的事件同时发生的事件记作A·B事件，则A·B事件的概率等于事A与事件B发生的概率的积，即P（A·B）=P(A)·P(B)。
        一般地，如果事件A1、A2、…An相互独立，那么这n个事件同时发生的概率等于每个事件发生的概率的积，
  即  P(A1·A2·…·An)=P(A1)·P(A2)·P(An).
-摘自《高中数理化概念公式定理手册》189页  上海远东出版社  ISBN 7-80613-324-0. 98年12月第一版

而在自然数中除以任意二个素数的余数变化是相互独立的，因此除以素数的余数变化是独立事件。
比如：
在单数列中除以任意奇素数的余数是以该素数的值为周期循环变化的；
同样在偶数列中除以任意奇素数的余数是以该素数的值为周期循环变化的；
它们并没有因为除以2的余数不同而对于除以其它素数时的余数的周期性变化而改变。
这就是自然数中除以任意素数的余数变化具有独立性的体现。

偶数2A分成的两个数表示为A±x，则x与A在除以√（2A）内的素数时的余数对应关系是A±x构成为素对的必要条件。
由已知A可以确定 A除以这些素数时的余数：j2、j3、j5、j7、……；
那么当x值除以这些素数时的余数同时满足不等于j2，不等于j3与3-j3，不等于j5与5-j5，不等于j7与7-j7，……时，这样的x值使得A±x不能被2，3，5，7 ，……这些素数整除而成为素数对，即{1+1}的解。
这样的x值的数量依据概率的乘法定理有：
Sp(m)= (A-2) P(2·3·…·n·…·r)
     =(A-2)*P(2)P(3)…P(n)…P(r)
     =(A-2)*(1/2)*f(3)*…*f(n)*…*f(r). --------------{式3}----这是人们通常称为“连乘式”的素对数量计算式。
式中：3≤ n≤r；n是素数。f(n)=(n-1)/n, [jn=0时]；或f(n)=(n-2)/n, [jn>0时] ；jn系A除以n时的余数。

实例：
例：偶数908，其√（908-2）内的最大素数是29，其半值A= 454，其分成两个素数对A±x的变量x的取值区间[0，A-3]中含有的整数为( 908/2- 2)个，
因此，其构成素对的x值的计算式是：
Sp( 908)=[( 908/2- 2)/2]*( 1/ 3)*( 3/ 5)*( 5/ 7)*( 9/ 11)*( 11/ 13)*( 15/ 17)*( 17/ 19)*( 21/ 23)*( 27/ 29)= 15

具体到每一步的含义：
1/2——[0，A-3]中满足除以2的余数不等于j2的数的发生概率；
( 1/ 3)—— [0，A-3]中满足除以3的余数不等于j3与（3-j3）的数的发生概率；
( 3/ 5)—— [0，A-3]中满足除以5的余数不等于j5与（5-j5）的数的发生概率；
( 5/ 7)—— [0，A-3]中满足除以7的余数不等于j7与（7-j7）的数的发生概率；
……
这里的j2,j3,…,jn,…，jr系偶数半值A除以素数2，3，…，n，…，r时的余数。

因此依据概率的独立事件的乘法定理：
在自然数[0，A-3]区域中除以素数2，3，…，n，…，r时余数同时满足不等于j2、j3及(3－j3)、j5及(5－j5)、…、jr及(r-jr)的x值的分布概率P(m)，
有P(m)=P(2·3·5·…·n·…·r)
      =P(2)P(3)…P(n)…P(r).
即有
Sp( 908)=( 908/2- 2)*P(m)=[( 908/2- 2)/2]*( 1/ 3)*( 3/ 5)*( 5/ 7)*( 9/ 11)*( 11/ 13)*( 15/ 17)*( 17/ 19)*( 21/ 23)*( 27/ 29)= 15
A= 454 ，
x= : 33 , 45 , 87 , 117 , 123 , 147 , 177 , 255 , 273 , 297 , 303 , 315 , 357 , 375 , 423 ,
把x代人到A±x中，得到全部素对：
[ 908 = ]  421 + 487  409 + 499  367 + 541  337 + 571  331 + 577  307 + 601  277 + 631  199 + 709  181 + 727  157 + 751  151 + 757  139 + 769  97 + 811  79 + 829  31 + 877
M= 908        S(m)= 15    S1(m)= 15   Sp(m)= 15.00   δ(m)≈ 0     δ1(m)≈  0    r= 29

连乘式的计算值Sp(m)与实际上不能被√M内素数整除的素对数量 S1(m)的相似程度是蛮高的。看下例的数据折线图：



图中的k(m),体现了偶数素对数量的波动程度的幅度。k(m)=π[(p1-1)/(p1-2)]，
在图中可以看到波动系数k(m)的波动与S1(m),Sp(m)二条折线的波动规律基本是同步的。
   
波动系数来源于上述的（式3）的变形：
{式3}经过数学变形，也可以用另一种形式表达:
        Sp(m)= (A-2)/2*π[(n-2)/n]*π[(p1-1)/(p1-2)] -------------------{式4}
        式中：3≤ n≤r；n是素数;
        p1是偶数半值A所含的奇素数因子. K(m)= π[(p1-1)/(p1-2)], K(m)为偶数M的素因子系数，也可以称为表法数波动系数。

而偶数全部素对S(m)的图形因为S(m)≥S1(m)，并且它们的差距很小，故S(m)图形的波动基本与S1(m)图形类似。

愚工688

连乘式的素对计算值的相对误差平均水平，在一亿以下的范围内是小于哈代计算式的相对误差平均水平的。
而在50000以上，连乘式的计算值的相对误差中值将逐渐离开0位，向0.20区域偏移；
而对应的哈代计算式，在偶数趋大时的相对误差将逐渐的趋向于0位。这就是数学家把哈代公式称作渐进式的因素。这个趋于0的趋势，是连乘式的计算值所不具备的优点。

连乘式的优点，是贴切的符合埃拉脱色尼筛法，对于素数的判断定理，也是比较直观的。如1#的实例。

在偶数趋大的各个阶段，连乘式计算值的相对误差的统计计算的样本数据如下：
相对误差δ(m)的统计计算:
M=[ 6 , 100 ]         r= 7    n= 48    μ=-.2418  σχ= .2292  δ(min)=-.625  δ(max)= .3429
M=[ 6 , 10000 ]       r= 97   n= 4998  μ=-.075   σχ= .0736  δ(min)=-.625  δ(max)= .3429
M=[ 10002 , 20000 ]   r= 139  n= 5000  μ=-.0315  σχ= .0361  δ(min)=-.1603 δ(max)= .1017
M=[ 20002 , 30000 ]   r= 173  n= 5000  μ=-.0100  σχ= .0288  δ(min)=-.1145 δ(max)= .1245  
M=[ 30002 , 40000 ]   r= 199  n= 5000  μ=-.0037  σχ= .0263  δ(min)=-.1034 δ(max)= .1101
M=[ 40002 , 50000 ]   r= 223  n= 5000  μ= .005   σχ= .0253  δ(min)=-.1021 δ(max)= .1131
M=[ 50002 , 60000 ]   r= 241  n= 5000  μ= .0082  σχ= .0219  δ(min)=-.0688 δ(max)= .1064
M=[ 60002 , 70000 ]   r= 263  n= 5000  μ= .0139  σχ= .0213  δ(min)=-.0681 δ(max)= .0993
M=[ 70002 , 80000 ]   r= 281  n= 5000  μ= .0145  σχ= .0202  δ(min)=-.051  δ(max)= .1006  
M=[ 80002 , 90000 ]   r= 293  n= 5000  μ= .0129  σχ= .0196  δ(min)=-.0597 δ(max)= .0976
M=[ 90002 , 100000 ]  r= 313  n= 5000  μ= .0218  σχ= .0174  δ(min)=-.038  δ(max)= .112
M=[ 100002 , 110000 ] r= 331  n= 5000  μ= .0233  σχ= .017   δ(min)=-.0381 δ(max)= .0906


15万——1亿的偶数样本的相对误差E(m)的统计计算数据如下：

[ 150002 , 150100 ]   :   n= 50    μ= .0316   σχ= .0135   δ(min)= .0004  δ(max)= .0589
[ 500002 , 500100 ]   :   n= 50    μ= .0536   σχ= .0084   δ(min)= .0359  δ(max)= .0698
[ 1000002 , 1000100 ] :   n= 50    μ= .0691   σχ= .0069   δ(min)= .0508  δ(max)= .0879
[ 2000002 , 2000100 ] :   n= 50    μ= .0804   σχ= .0063   δ(min)= .0641  δ(max)= .0951
[ 3000002 , 3000100 ] :   n= 50    μ= .0825   σχ= .005    δ(min)= .0706  δ(max)= .0923
[ 4000002 , 4000100 ] :   n= 50    μ= .0874   σχ= .004    δ(min)= .0792  δ(max)= .0967
[ 5000002 , 5000100 ] :   n= 50    μ= .0923   σχ= .004    δ(min)= .0841  δ(max)= .1012
[ 6000002 , 6000100 ] :   n= 50    μ= .0897   σχ= .0041   δ(min)= .0799  δ(max)= .0999
[ 7000002 , 7000100 ] :   n= 50    μ= .096    σχ= .0033   δ(min)= .0884  δ(max)= .1027
[ 8000002 , 8000100 ] :   n= 50    μ= .0935   σχ= .0035   δ(min)= .0842  δ(max)= .1003
[9000002 -  9000100 ] :   n= 50    μ= .0959   σχ= .0026   δ(min)= .0893  δ(max)= .1007
[10000000 - 10000100] :   n= 51    μ= .10032  σχ= .00256  δ(min)= .09543 δ(max)= .10503
[20000002 - 20000100] :   n= 50    μ= .1046   σχ= .0022   δ(min)= .0969  δ(max)= .1094


1亿-500亿的取样样本的相对误差的统计计算数据：
（标准偏差的通用符号为σx ，μ-样本平均值）

100000000 -   100000098 : n= 50 μ= .1192  σx= .0013  δ(min)= .1156  δ(max)= .1224
1000000000 - 1000000098 : n= 50 μ= .1368  σx= .0004  δ(min)= .1356  δ(max)= .138
2000000000 - 2000000098 : n= 50 μ= .1406  σx= .0003  δ(min)= .1399  δ(max)= .141
3000000000 - 3000000098 : n= 50 μ= .1431  σx= .0002  δ(min)= .1425  δ(max)= .1435
4000000000 - 4000000098 : n= 50 μ= .1449  σx= .0003  δ(min)= .1441  δ(max)= .1456
5000000000 - 5000000098 : n= 50 μ= .1462  σx= .0003  δ(min)= .1456  δ(max)= .1468
5999999990 - 6000000088 : n= 50 μ= .1471  σx= .0002  δ(min)= .1466  δ(max)= .1474  
8000000000 - 8000000050 : n= 26 μ= .1486  σx= .0002  δ(min)= .1481  δ(max)= .1490
10000000000-10000000098 : n= 50 μ= .1494  σx= .0002  δ(min)= .1491  δ(max)= .1497
15000000000-15000000098 : n= 50 μ= .15159 σx= .00014 δ(min)= .1511  δ(max)= .15185
20000000002-20000000100 : n= 50 μ= .15281 σx= .00011 δ(min)= .1525  δ(max)= .15307
30000000002-30000000100 : n= 50 μ= .15494 σx= .0001  δ(min)= .15474 δ(max)= .15519
40000000002-40000000100 : n= 50 μ= .15614 σx= .00008 δ(min)= .1559  δ(max)= .15637  
50000000002-50000000100 : n= 50 μ= .1571  σx= .0001  δ(min)= .1569  δ(max)= .1573
70000000000 - 70000000048 : n= 25 μ= .158689 σx = .000061  δmin = .158571 δmax = .158863
80000000000 - 80000000048 : n= 25 μ= .159080 σx = .000052  δmin = .158896 δmax = .159196
100000000000-100000000048 : n= 25 μ= .160175 σx = .000049  δmin = .16005  δmax = .16026
200000000000-200000000048 : n= 25 μ= .162808 σx = .000041  δmin = .16272  δmax = .16289
400000000000-400000000038 : n= 20 μ= .16544  σx = .000024  δmin = .165403 δmax = .165486

愚工688

相对误差修正系数的运用原理

由相对误差的定义知：
相对误差  Δ=（计算值-真值）/真值
整理可得  真值=计算值/（1+Δ） {式A}

式A显示了计算值及相对误差与真值的关系。
而我们使用素对计算式计算实际的真值，是不可能预先知道偶数的相对误差Δ的。但是通过样本区域的偶数的相对误差的统计计算，我们可以知道样本区域相对误差的平均值μ是与实际的相对误差值很接近的，
因此，使用平均值μ代替{式A}的真实的相对误差值Δ，将能够大幅度的提高连乘式的计算值的计算精度。
同时由于临近的样本区域的相对误差的平均值μ变化比较小，因此可以把一个样本区域的相对误差的均值μ扩大使用到一个比较大的范围内的偶数，修正后的计算值精度不会下降很多。


由于连乘式计算值的相对误差变化具有规律性，即样本区域的最大值与最小值趋近，标准偏差比较小。
因此当我们对偶数趋大时的素对计算值的相对误差趋大的现象可以进行一定的修正。
采用的具有修正系数的计算式如下：
计算式： 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时的余数。
更多的计算实例数据，可以参见我的帖子：高精度计算大偶数表为两个素数和的表法数值的实例（……）
  http://www.mathchina.com/bbs/for ... id=59160&extra=

愚工688

用inf( m )=Sp( m )/(1+μ) 来计算1亿-100亿偶数的素对数量下界：
（这里的μ=0.1502  ，系125亿样本小区域的素对计算值相对误差值进行统计计算得到的相对误差均值μ)

G(100000000) = 291400,inf( 100000000 )≈  283684.9 , Δ≈-0.026, k(m)= 1.33333
G(100000002) = 464621,inf( 100000002 )≈  451550.5 , Δ≈-0.028, k(m)= 2.12231
G(100000004) = 247582,inf( 100000004 )≈  240749.6 , Δ≈-0.028, k(m)= 1.13154
G(100000006) = 218966,inf( 100000006 )≈  213198.8 , Δ≈-0.026, k(m)= 1.00204
G(100000008) = 437717,inf( 100000008 )≈  425527.4 , Δ≈-0.028, k(m)= 2
G(100000010) = 323687,inf( 100000010 )≈  315205.5 , Δ≈-0.026, k(m)= 1.48148
G(100000012) = 263241,inf( 100000012 )≈  255316.5 , Δ≈-0.03, k(m)= 1.2

G(1000000000) = 2274205,inf( 1000000000 )≈  2245348.2 , Δ≈-0.013, k(m)= 1.33333
G(1000000002) = 3496205,inf( 1000000002 )≈  3454562.6 , Δ≈-0.012, k(m)= 2.05139
G(1000000004) = 1747858,inf( 1000000004 )≈  1727191 , Δ≈-0.012, k(m)= 1.02564
G(1000000006) = 1704301,inf( 1000000006 )≈  1684011.2 , Δ≈-0.012, k(m)= 1
G(1000000008) = 4151660,inf( 1000000008 )≈  4104122.7 , Δ≈-0.011, k(m)= 2.43711
G(1000000010) = 2422662,inf( 1000000010 )≈  2395038.2 , Δ≈-0.011, k(m)= 1.42222
G(1000000012) = 1960129,inf( 1000000012 )≈  1937861 , Δ≈-0.011, k(m)= 1.15074

G(2000000000) = 4238417,inf( 2000000000 )≈  4203544.7 , Δ≈-0.0082, k(m)= 1.33333
G(2000000002) = 4897539,inf( 2000000002 )≈  4855431.3 , Δ≈-0.0086, k(m)= 1.54011
G(2000000004) = 6519934,inf( 2000000004 )≈  6467330.1 , Δ≈-0.0081, k(m)= 2.05139
G(2000000006) = 3342074,inf( 2000000006 )≈  3313613.9 , Δ≈-0.0085, k(m)= 1.05105
G(2000000008) = 3261215,inf( 2000000008 )≈  3233495.9 , Δ≈-0.0085, k(m)= 1.02564
G(2000000010) = 8478380,inf( 2000000010 )≈  8407089.5 , Δ≈-0.0084, k(m)= 2.66667
G(2000000012) = 3180443,inf( 2000000012 )≈  3152658.5 , Δ≈-0.0087, k(m)= 1

G(4000000000) = 7930427, inf( 4000000000 )≈  7891735.0 , Δ≈-0.0049, k(m)= 1.33333
G(4000000002) = 11887591,inf( 4000000002 )≈ 11837602.6 , Δ≈-0.0042, k(m)= 2
G(4000000004) = 9156520, inf( 4000000004 )≈  9115760.4 , Δ≈-0.0045, k(m)= 1.54014
G(4000000006) = 6404412, inf( 4000000006 )≈  6373721.7 , Δ≈-0.0048, k(m)= 1.07686
G(4000000008) = 12198479,inf( 4000000008 )≈ 12141765.9 , Δ≈-0.0046, k(m)= 2.05139
G(4000000010) = 7926931, inf( 4000000010 )≈  7892524.4 , Δ≈-0.0043, k(m)= 1.33347
G(4000000012) = 6249883, inf( 4000000012 )≈  6220979.0 , Δ≈-0.0046, k(m)= 1.05105

G(6000000000) = 22899781,inf( 6000000000 )≈ 22831687.7 , Δ≈-0.0030, k(m)= 2.66667
G(6000000002) = 8585981 ,inf( 6000000002 )≈  8563011.4 , Δ≈-0.0027, k(m)= 1.00013
G(6000000004) = 8588030 ,inf( 6000000004 )≈  8561882.9 , Δ≈-0.0030, k(m)= 1
G(6000000006) = 26447626,inf( 6000000006 )≈ 26372932.4 , Δ≈-0.0028, k(m)= 3.08027
G(6000000008) = 8957244 ,inf( 6000000008 )≈  8934138.7 , Δ≈-0.0026, k(m)= 1.04348
G(6000000010) = 11446102,inf( 6000000010 )≈ 11415843.9 , Δ≈-0.0026, k(m)= 1.33333
G(6000000012) = 17617549,inf( 6000000012 )≈ 17563755.4 , Δ≈-0.0031, k(m)= 2.05139

G(8000000000) = 14862150,inf( 8000000000 )≈  14841172.2 , Δ≈-0.0014, k(m)= 1.33333
G(8000000002) = 11485548,inf( 8000000002 )≈  11469257.9 , Δ≈-0.0014, k(m)= 1.0304
G(8000000004) = 22296318,inf( 8000000004 )≈  22261758.3 , Δ≈-0.0016, k(m)= 2
G(8000000006) = 11146652,inf( 8000000006 )≈  11131349.1 , Δ≈-0.0014, k(m)= 1.00004
G(8000000008) = 17167422,inf( 8000000008 )≈  17143070.4 , Δ≈-0.0014, k(m)= 1.54014
G(8000000010) = 29840750,inf( 8000000010 )≈  29801550.7 , Δ≈-0.0013, k(m)= 2.67738
G(8000000012) = 11998604,inf( 8000000012 )≈  11986401.1 , Δ≈-0.0010, k(m)= 1.07686

G(10000000000) = 18200488,inf( 10000000000 )≈  18189357.2 , Δ≈-0.00061, k(m)= 1.333
G(10000000002) = 27302893,inf( 10000000002 )≈  27284035.8 , Δ≈-0.00069, k(m)= 2
G(10000000004) = 13655366,inf( 10000000004 )≈  13642017.9 , Δ≈-0.00098, k(m)= 1
G(10000000006) = 13742400,inf( 10000000006 )≈  13734820.7 , Δ≈-0.00055, k(m)= 1.0068
G(10000000008) = 27563979,inf( 10000000008 )≈  27543883.7 , Δ≈-0.00073, k(m)= 2.019
G(10000000010) = 28031513,inf( 10000000010 )≈  28014088.3 , Δ≈-0.00062, k(m)= 2.0535
G(10000000012) = 13654956,inf( 10000000012 )≈  13644784.5 , Δ≈-0.00074, k(m)= 1.0002

显然，只要修正系数μ的取值不要离开统计计算的区域过远，那么素对下界计算值的计算精度就能够得到保障。

yangchuanju

愚工688 发表于 2022-5-28 19:28
用inf( m )=Sp( m )/(1+μ) 来计算1亿-100亿偶数的素对数量下界：
（这里的μ=0.1502  ，系125亿样本小 ...

愚工688老师：
您好！
老师曾计算并提供给10^10-10^16的单计哥猜数，学生一度转载过：
G( 10^10) = 18200488; ( 2.539 sec)
G( 10^11) = 149091160; (38.146 sec)
G( 10^12 ) = 1243722370 (625.098 sec)
G( 10^13 ) = 10533150855, ;(1090.54 sec，约18分10秒)
G( 10^14 ) = 90350630388 ;(12740.44 sec，约3个半小时多一点 ) ;
G( 10^15 ) = 783538341852 ;(169664.44 sec，约 47.13h ）
G(10^16)=68526741128787 (单记法)—— 挂机计算，忘记记录时间了。
其中的10^10-10^14数据与A065577完全相同，10^16哥猜数是第一次见到的！

白新岭先生也曾引用过愚工老师的这些数据，并与它的计算数据进行过比对：
10^n    偶数为10的n次方(理论值)    愚工688给的数据单记    比值(理论/实际）
10      36548552          18200488          1.004054177
11      299158483         149091160         1.003273712
12      2494079818        1243722370        1.002667427
13      21112797736       10533150855       1.002207128
14      1.8104035318900000E+11     9.0350630388000000E+10    1.001876536
15      1.5696113224770000E+12     7.8353834185200000E+11    1.001617431
16      1.3738923405990000E+13     6.8526741128787000E+12    1.002449787
其中白新岭的理论值为双计，愚工的10^16哥猜数为单计6.8526741128787000E+12，有效数字14位，实为6.8526741128787*10^13；

与10^1-10^15的哥猜数进行综合比对，14位数字68526741128787之中可能多了一个数字，实际哥猜数应为13位数字。
学生相信愚工老师花费百多个小时（估计数）的计算结果是不会出错的，可能是发帖录入时多写了一个数字：
或许是411...中多了一个1，或许是28787之中多了一个7或8，或许是其它，学生无法分清。
请愚工老师查对（复核）一下10^16的哥猜数究竟是多少！

敬候佳音！
                  您的学生  yangchuanju

愚工688

yangchuanju 发表于 2022-6-15 00:38
愚工688老师：
您好！
老师曾计算并提供给10^10-10^16的单计哥猜数，学生一度转载过：

谢谢，确实如此，在手工记录时10^16的素对数量时发生了错误，多出了一位数。把它删去吧！没有兴趣再重新计算一遍。

  G( 10^10 ) = 18200488         ;Xi(M)≈ 18179890.52          δxi(M)≈-0.001132  
  G( 10^11 ) = 149091160        ;Xi(M)≈ 148486029.78         δxi(M)≈-0.004059  
  G( 10^12 ) = 1243722370       ;Xi(M)≈ 1233556241.87        δxi(M)≈-0.008174  
  G( 10^13 ) = 10533150855,     ;Xi(M)≈ 10395227871.57       δxi(M)≈-0.013094
  G( 10^14 ) = 90350630388      ;Xi(M)≈ 88673642506.88       δxi(M)≈-0.018561
  G( 10^15 ) = 783538341852     ;Xi(M)≈ 764388083252.93      δxi(M)≈-0.024441
   G(10^16)   =68526741128787 (单记法)——明显多了一位数。
可以从下面两个偶数的素对计算时间看出，偶数大10倍，用时要10倍多，故在原始的记录中能否找到少一位的记录，不再重新计算了。——没有找到，只能猜测多写了一个“1”。
G( 10^14 ) = 90350630388  ;(12740.44 sec，约3个半小时多一点 ) ;     
  G( 10^15 ) = 783538341852 ;(169664.44 sec，约 47.13h ）

因为使用的这个高速计算程序：Gpartiton 没有文本输出（也许是我不熟悉使用其功能），只是在屏幕上面显示运算结果，手工纪录下来，眼睛一花出错的可能性比较大。  

重生888@

令10^16=N
D(N)=5/6*(N+F*N/lnN)/(lnN)^2                F=3.3556
        =6639661079555
G(N)=68526741128787或6852674128787
1. D/G=0.096891
2. D/G=0.968915
明显是第二个是对的！

大傻8888888

yangchuanju 发表于 2022-6-15 08:38
愚工688老师：
您好！
老师曾计算并提供给10^10-10^16的单计哥猜数，学生一度转载过：

yangchuanju先生：您好！
您在这个帖子提到白新岭先生也曾引用过愚工老师的这些数据，并与它的计算数据进行过比对。其中 偶数为10的n次方(理论值)是用什么方法和公式得出的，看起来比哈李公式要准确，不知白新岭先生是怎样做到的，愿闻其详，谢谢啦 ！

大傻8888888

yangchuanju 发表于 2022-6-15 08:38
愚工688老师：
您好！
老师曾计算并提供给10^10-10^16的单计哥猜数，学生一度转载过：

       白新岭在其《k生素数的数量公式》帖子中一般只提供计算结果，不提供计算公式是因为他怕别人抢夺他辛勤努力得来的成果，我也就不强人所难了。不过光有计算公式，如果不能证明其成立，这样的公式不过是个猜想而已，有些看起来比较符合实际值，不过是凑数罢了。不过也有可能是我吃不到葡萄就说葡萄是酸的。还是八仙过海各显神通好！

重生888@

大傻8888888 发表于 2022-6-22 23:50
白新岭在其《k生素数的数量公式》帖子中一般只提供计算结果，不提供计算公式是因为他怕别人抢夺他 ...

怕人抢夺，言不由衷；把人用，人家还不用呢！见过随便加个系数，能计算连续偶数吗？只把别人拿来而已！