我随意计算今天日期10倍的10个偶数的素数对，有计算值精度低于0.99的吗？欢迎验证

愚工688

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

  

  G( 202209270 ) = ?       ;Xi(M)≈ 1103774.81   δxi(M)≈?
  G( 202209272 ) = ?       ;Xi(M)≈ 412289.62    δxi(M)≈?
  G( 202209274 ) = ?       ;Xi(M)≈ 408251.13    δxi(M)≈?
  G( 202209276 ) = ?       ;Xi(M)≈ 816333.48    δxi(M)≈?
  G( 202209278 ) = ?       ;Xi(M)≈ 416640.66    δxi(M)≈?
  G( 202209280 ) = ?       ;Xi(M)≈ 737003.26    δxi(M)≈?
  G( 202209282 ) = ?       ;Xi(M)≈ 907037.24    δxi(M)≈?
  G( 202209284 ) = ?       ;Xi(M)≈ 435105.1     δxi(M)≈?
  G( 202209286 ) = ?       ;Xi(M)≈ 450480.42    δxi(M)≈?
  G( 202209288 ) = ?       ;Xi(M)≈ 870755.77    δxi(M)≈?
  time start =09:43:55, time end =09:44:02

  G( 2022092700 ) = ?      ;Xi(M)≈ 8682772.37   δxi(M)≈?
  G( 2022092702 ) = ?      ;Xi(M)≈ 4036455.44   δxi(M)≈?
  G( 2022092704 ) = ?      ;Xi(M)≈ 3242607.18   δxi(M)≈?
  G( 2022092706 ) = ?      ;Xi(M)≈ 6421633.81   δxi(M)≈?
  G( 2022092708 ) = ?      ;Xi(M)≈ 3422815.74   δxi(M)≈?
  G( 2022092710 ) = ?      ;Xi(M)≈ 5110164.57   δxi(M)≈?
  G( 2022092712 ) = ?      ;Xi(M)≈ 6433416.36   δxi(M)≈?
  G( 2022092714 ) = ?      ;Xi(M)≈ 3210816.92   δxi(M)≈?
  G( 2022092716 ) = ?      ;Xi(M)≈ 3946955.33   δxi(M)≈?
  G( 2022092718 ) = ?      ;Xi(M)≈ 6421633.85   δxi(M)≈?
  time start =09:44:19, time end =09:44:52

我懒得验证，至于计算值的计算精度牛不牛，大家可以验证试试：

时空伴随者

2022/9/27 10:28:21
G(202209270) = 1103799
G(202209272) = 411888
G(202209274) = 408156
G(202209276) = 815746
G(202209278) = 416733
G(202209280) = 737875
G(202209282) = 906020
G(202209284) = 435513
G(202209286) = 450594
G(202209288) = 869516
G(202209290) = 546695
用时: 318.06815 秒
2022/9/27 11:06:21
G(2022092700) = 8682983
G(2022092702) = 4037152
G(2022092704) = 3242484
G(2022092706) = 6421597
G(2022092708) = 3423316
G(2022092710) = 5110529
G(2022092712) = 6431460
G(2022092714) = 3212211
G(2022092716) = 3945564
G(2022092718) = 6422099
G(2022092720) = 4323417
用时: 4923.26115 秒

时空伴随者

2022-09-27 13:19:05
202209270 [2, 3, 5, 73, 92333]
202209272 [2, 2, 2, 101, 250259]
202209274 [2, 7607, 13291]
202209276 [2, 2, 3, 16850773]
202209278 [2, 71, 167, 8527]
202209280 [2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 5, 7, 7, 13, 31]
202209282 [2, 3, 3, 11, 1021259]
202209284 [2, 2, 23, 59, 37253]
202209286 [2, 19, 29, 281, 653]
202209288 [2, 2, 2, 3, 17, 495611]
202209290 [2, 5, 229, 88301]
用时 0.0009999275207519531 秒
2022-09-27 13:17:55
2022092700 [2, 2, 3, 5, 5, 73, 92333]
2022092702 [2, 7, 7, 7, 23, 128159]
2022092704 [2, 2, 2, 2, 2, 103, 613499]
2022092706 [2, 3, 337015451]
2022092708 [2, 2, 19, 149, 178567]
2022092710 [2, 5, 11, 11, 17, 197, 499]
2022092712 [2, 2, 2, 3, 3, 547, 51343]
2022092714 [2, 1011046357]
2022092716 [2, 2, 7, 43, 1679479]
2022092718 [2, 3, 337015453]
2022092720 [2, 2, 2, 2, 5, 101, 250259]
用时 0.0010001659393310547 秒

愚工688

时空伴随者 发表于 2022-9-27 05:10
2022/9/27 10:28:21
G(202209270) = 1103799
G(202209272) = 411888

你的真值都是正确的，就是筛选速度不快。基本上筛选百亿级别的偶数是可以的，千亿的偶数就勉强了，万亿的连续偶数就不行了。


  计算式： Xi(M)=t2*c1*M/(logM)^2   
        式中： 相对误差动态修正系数 t2=1.358-log(M)^(0.5)*.05484;
                 C1--类似拉曼扭杨系数，略作改进；（只计算√M内的素数）  

  G(202209270) = 1103799   ;Xi(M)≈ 1103774.81   infS(m)= 408166.73  jd(M) ≈? 0.999978；
  G(202209272) = 411888    ;Xi(M)≈ 412289.62    infS(m)= 408166.72  jd(M) ≈? 1.000975；
  G(202209274) = 408156    ;Xi(M)≈ 408251.13    infS(m)= 408166.74  jd(M) ≈? 1.000233；
  G(202209276) = 815746    ;Xi(M)≈ 816333.48    infS(m)= 408166.74  jd(M) ≈? 1.0007196；
  G(202209278) = 416733    ;Xi(M)≈ 416640.66    infS(m)= 408166.74  jd(M) ≈? 0.999779；
  G(202209280) = 737875    ;Xi(M)≈ 737003.26    infS(m)= 408166.74  jd(M) ≈? 0.998818；
  G(202209282) = 906020    ;Xi(M)≈ 907037.24    infS(m)= 408166.76  jd(M) ≈? 1.001122；
  G(202209284) = 435513    ;Xi(M)≈ 435105.1     infS(m)= 408166.77  jd(M) ≈? 0.999063；
  G(202209286) = 450594    ;Xi(M)≈ 450480.42    infS(m)= 408166.76  jd(M) ≈? 0.999747；
  G(202209288) = 869516    ;Xi(M)≈ 870755.77    infS(m)= 408166.77  jd(M) ≈? 1.001426；
  G(202209290) = 546695    ;Xi(M)≈ 546619.83    infS(m)= 408166.78  jd(M) ≈? 0.999863；
  G(202209292) = 420762    ;Xi(M)≈ 419828.67    infS(m)= 408166.76  jd(M) ≈? 0.997783；
  G(202209294) = 1004548   ;Xi(M)≈ 1004718.19   infS(m)= 408166.76  jd(M) ≈? 1.000169；
  G(202209296) = 410276    ;Xi(M)≈ 409889       infS(m)= 408166.78  jd(M) ≈? 0.999057；
  G(202209298) = 408504    ;Xi(M)≈ 408166.79    infS(m)= 408166.79  jd(M) ≈? 0.999175；
  time start =09:57:10, time end =09:57:21



  Xi(M)=t2*c1*M/(logM)^2   

  G(2022092700) = 8682983  ;Xi(M)≈ 8682772.37   infS(m)= 3210816.87         jd(M) ≈? 0.999976;
  G(2022092702) = 4037152  ;Xi(M)≈ 4036455.44   infS(m)= 3210816.83         jd(M) ≈? 0.999827;
  G(2022092704) = 3242484  ;Xi(M)≈ 3242607.18   infS(m)= 3210816.91         jd(M) ≈? 1.000038;
  G(2022092706) = 6421597  ;Xi(M)≈ 6421633.81   infS(m)= 3210816.9          jd(M) ≈? 1.000006;
  G(2022092708) = 3423316  ;Xi(M)≈ 3422815.74   infS(m)= 3210817.02         jd(M) ≈? 0.999854;
  G(2022092710) = 5110529  ;Xi(M)≈ 5110164.57   infS(m)= 3210816.76         jd(M) ≈? 0.999929;
  G(2022092712) = 6431460  ;Xi(M)≈ 6433416.36   infS(m)= 3210816.77         jd(M) ≈? 1.000304;
  G(2022092714) = 3212211  ;Xi(M)≈ 3210816.92   infS(m)= 3210816.92         jd(M) ≈? 0.999566;
  G(2022092716) = 3945564  ;Xi(M)≈ 3946955.33   infS(m)= 3210816.84         jd(M) ≈? 1.000353;
  G(2022092718) = 6422099  ;Xi(M)≈ 6421633.85   infS(m)= 3210816.92         jd(M) ≈? 0.999928;
  G(2022092720) = 4323417  ;Xi(M)≈ 4324332.52   infS(m)= 3210816.9          jd(M) ≈? 1.000212;
  G(2022092722) = 3564718  ;Xi(M)≈ 3564160.3    infS(m)= 3210816.82         jd(M) ≈? 0.999843;
  G(2022092724) = 6440134  ;Xi(M)≈ 6439929.2    infS(m)= 3210816.97         jd(M) ≈? 0.999968;
  G(2022092726) = 3211021  ;Xi(M)≈ 3211102.2    infS(m)= 3210817.02         jd(M) ≈? 1.000025;
  G(2022092728) = 3220806  ;Xi(M)≈ 3220838.93   infS(m)= 3210816.94         jd(M) ≈? 1.000010;
  time start =10:28:48, time end =10:29:40

对数计算式我只能控制相对误差的绝对值比较小，不能控制在下限的位置。连乘式就可以，因为连乘式的计算值是分区线性的。

愚工688

30亿的30个连续偶数的素对筛选就1秒多一点：

3000000000:30:2

G(3000000000) = 12224533
G(3000000002) = 4584282
G(3000000004) = 4783729
G(3000000006) = 9403062
G(3000000008) = 4662802
G(3000000010) = 7336233
G(3000000012) = 9403617
G(3000000014) = 5284363
G(3000000016) = 5000395
G(3000000018) = 9167945
G(3000000020) = 6112467
G(3000000022) = 5035637
G(3000000024) = 11170019
G(3000000026) = 4616381
G(3000000028) = 4583241
G(3000000030) = 13037886
G(3000000032) = 4803861
G(3000000034) = 4718646
G(3000000036) = 10552337
G(3000000038) = 5503450
G(3000000040) = 6204959
G(3000000042) = 10099510
G(3000000044) = 4601935
G(3000000046) = 4591669
G(3000000048) = 9168559
G(3000000050) = 6152020
G(3000000052) = 5502506
G(3000000054) = 9169787
G(3000000056) = 4588492
G(3000000058) = 5101569

count = 30, algorithm = 2, working threads = 2, time use 1.279 sec