花、耿、崔们丢不能估阶的余项只用主项“推得”的所谓“下限式”无效！
141不懂余项小于主项的1/2，胡说“丢不能估阶的余项只用主项“推得”的所谓“下限式”无效！”
例如100的主项是：
100/2/3*4/5*5/7=9.5238095238，
因为100的余项小于：
100/2/2/3*4/5*5/7=9.5238095238/2，
所以100的素对“下限式”大于：
100/2/2/3*4/5*5/7=9.5238095238/2，

估计主项：
N/2*1/3*...*9971/9973*...*(p-2)/p>40N^1/2,
(N>100000000,p是小于N^1/2的最大素数)
估计余项：
D(N)/(N/2*1/3*...*9971/9973*...*(p-2)/p)>1/2，
估计N以内的D(N)下限值：
N/4*1/3*...*9971/9973*...*(p-2)/p>20N^1/2,
根据剩余概率估计不小于100000000的每个偶数与200000个素数的差仍然是素数，
哥猜孪猜正确。

10^31以内的孪生素数个数：
10^31/2/3*...*(P-2)/p接近于2.74750556642637e+27，
2*10^31以内的孪生素数个数：
2*10^31/2/3*...*(P-2)/p接近于5.38931728330686e+27，
5.38931728330686*10^27-2.74750556642637*10^27=2.64181171688048*10^27,
(2.64181171688048*10^27)/(2.74750556642637*10^27)=0.9615309789223232,

2^10=16.4146635763/26.8603585795=0.61
2^11=26.8603585795/44.76726429909=0.6
2^12=44.76726429909/75.7599857369=0.5909090909
2^13=75.7599857369/129.8742612633=0.583333333
2^14=129.8742612633/225.1153861898=0.5769230769
2^15=225.1153861898/393.9519258321=0.5714285714,
2^16=393.9519258321/695.20928088=0.5666666667,
2^17=695.20928088/1235.927610453=0.5625,
2^18=1235.927610453/2211.659934495=0.5588235294
2^19=2211.659934495/3980.987882092=0.5555555556
2^20=3980.987882092,
2^21>3980.987882092*3/2=5971.481823138,(7203.692358072)

2^100的素对多于1.8^（100-5）,
2^200的素对多于1.8^（200-5）,
2^300的素对多于1.8^（300-5）,
2^400的素对多于1.8^（400-5）,
2^500的素对多于1.8^（500-5）,

计算证明哥德巴赫猜想要先根据剩余公式1/2*...计算出大偶数N的剩余数对近似值：
N/2/3*...*9971/9973*...*(p-2)/p>20N^1/2,（主项）
由于余项恒小于主项的1/2，所以N的素对下限值越来越大于：
N/2/2/3*...*9971/9973*...*(p-2)/p>10N^1/2,
例如：
2^20000/2/2/3*...*9971/9973*...*(p-2)/p>5*2^10001,
确定不小于2^20000的每个偶数存在的素数对越来越多于5*2^10001,

证明哥德巴赫猜想要根据剩余公式1/2*...计算出大偶数N的剩余数对近似值：
N/2/3*...*9971/9973*...*(p-2)/p>20N^1/2,（主项）
由于余项恒小于主项的1/2，所以N的素对下限值越来越大于：
N/2/2/3*...*9971/9973*...*(p-2)/p>10N^1/2,
例如：
2^20000/2/2/3*...*9971/9973*...*(p-2)/p>5*2^10001,
确定不小于2^20000的每个偶数存在的素数对越来越多于5*2^10001,
比例计算，精妙绝伦！

耿受明天的素对下限计算式 Ginf(N)=(M/log(M/2)- M/log(M) ) ——（单记）相对误差水平：
G(10^ 10 ) = 18200488 ;Ginf(M)≈ 17972409.5 δG(10^ 10 )≈-0.01253
G(10^ 11 ) = 149091160 ;Ginf(M)≈ 148114807.31 δG(10^ 11 )≈-0.006549
G(10^ 12 ) = 1243722370 ;Ginf(M)≈ 1241664436.16 δG(10^ 12 )≈-0.0016547
G(10^ 13 ) = 10533150855 ;Ginf(M)≈ 10558962596.08 δG(10^ 13 )≈ 0.002451
G(10^ 14 ) = 90350630388 ;Ginf(M)≈ 90889971838.51 δG(10^ 14 )≈ 0.005969
G(10^ 15 ) = 783538341852 ;Ginf(M)≈ 790594418665.97 δG(10^ 15 )≈ 0.0090054
从上面偶数的相对误差数据中可以看出，耿下限计算值增大的速度大于真值增大的速度，计算值在10^13 时已经大于真值，并且相对误差绝对值随10^n的指数n增大而逐渐增大。
正如我在9楼中预测的那样，耿计算式对10^ 15 的计算值的相对误差已经达到了0.009以上，因此耿计算式对10^ 16 的计算值的相对误差必然会达到0.01以上。（估计在0.012左右）
虽然说相对误差绝对值都不大，但是耿的下限计算式在大偶数时显然得出的不是下限计算值。

我使用计算式 Xi(M)=t1*c1*M/(logM)^2 ;( t1=1.358-log(M)^(2.045/3)*.03178 )；
c1——类同拉曼扭扬系数，但是只计算根号M内的素数。
计算偶数 M= 2^n 起的连续偶数的素对数量，在n=29处，计算值与真值处于相互纠缠，相对误差值很小，符号正负混杂。
而在n=30以后，基本处于下界状态，并且随n增大，下界计算值的相对误差绝对值逐渐增大，就是始终处于下降状态。
M= 2^29 = 536870912：
S( 536870912 ) = 975685 ;Xi(M)≈ 975604.11 δxi( 536870912 )≈-0.00008302
S( 536870914 ) = 1041368 ;Xi(M)≈ 1040644.39 δxi( 536870914 )≈-0.00069486
S( 536870916 ) = 2065478 ;Xi(M)≈ 2065985.14 δxi( 536870916 )≈ 0.00024546
S( 536870918 ) = 975861 ;Xi(M)≈ 975604.12 δxi( 536870918 )≈-0.00026336
S( 536870920 ) = 1334819 ;Xi(M)≈ 1335696.9 δxi( 536870920 )≈ 0.00065769
S( 536870922 ) = 2453728 ;Xi(M)≈ 2452947.45 δxi( 536870922 )≈-0.0003181
M=2^30 = 1073741824 ：
G(1073741824) = 1817111 ;Xi(M)≈ 1813876.74 δxi( 1073741824 )≈-0.001780
G(1073741826) = 3698190 ;Xi(M)≈ 3691398.42 δxi( 1073741826 )≈-0.001836
G(1073741828) = 1937221 ;Xi(M)≈ 1934801.87 δxi( 1073741828 )≈-0.001249
G(1073741830) = 2906799 ;Xi(M)≈ 2904091 δxi( 1073741830 )≈-0.000932
G(1073741832) = 3846703 ;Xi(M)≈ 3841150.69 δxi( 1073741832 )≈-0.001443
M=2^31 = 2147483648 ：
S( 2147483648 ) = 3390038 ;Xi(M)≈ 3380024.25 δxi( 2147483648 )≈-0.002954
S( 2147483650 ) = 5147510 ;Xi(M)≈ 5131813.66 δxi( 2147483650 )≈-0.003049
S( 2147483652 ) = 6897846 ;Xi(M)≈ 6878646.13 δxi( 2147483652 )≈-0.002783
S( 2147483654 ) = 3389472 ;Xi(M)≈ 3380024.26 δxi( 2147483654 )≈-0.002787
S( 2147483656 ) = 3615850 ;Xi(M)≈ 3605359.24 δxi( 2147483656 )≈-0.002901
M=2^32 = 4294967296 ：
S( 4294967296 ) = 6341424 ;Xi(M)≈ 6311717.92 δxi( 4294967296 )≈-0.004684
S( 4294967298 ) = 12679919 ;Xi(M)≈ 12623435.85 δxi( 4294967298 )≈-0.004455
S( 4294967300 ) = 9627145 ;Xi(M)≈ 9582937.21 δxi( 4294967300 )≈-0.004592
S( 4294967302 ) = 6424654 ;Xi(M)≈ 6394967.63 δxi( 4294967302 )≈-0.004621
S( 4294967304 ) = 12903730 ;Xi(M)≈ 12844900.15 δxi( 4294967304 )≈-0.004559
因此我可以说，计算式 Xi(M)=t1*c1*M/(logM)^2 在 M=2^31 起成为素对下界计算式是没有问题的，也是符合实际偶数的素对真实数量的。

我给出了计算偶数N的素数对近似和下限值公式：
p(N/ln（N/2）-N/ln（N）)=近似值,
p(N/ln（N/1.9）-N/ln（N）)=下限值,
有了这个正确计算素数对的计算公式，哥德巴赫猜想和孪生素数猜想迎刃而解。
例如2^10的孪生素数对多于：2*(2^10/ln(2^10/1.9)-2^10/ln(2^10))=30.15201930379，
偶数不小于2^10的素数对多于：(2^10/ln(2^10/1.9)-2^10/ln(2^10))>15,