|
本帖最后由 愚工688 于 2022-9-16 04:05 编辑
时空先生71楼的素对真值数据是正确的:
我使用连乘式的计算:
G(1000000364) = 1913317;Sp( 1000000364 *)≈ 1914805.0 , jdz =sp(m)/G(m) ≈ 1.00078;
G(1000000366) = 1704613;Sp( 1000000366 *)≈ 1703862.0 , jdz =sp(m)/G(m) ≈ 0.999559;
G(1000000368) = 3424584;Sp( 1000000368 *)≈ 3425199.6 , jdz =sp(m)/G(m) ≈ 1.000180;
G(1000000370) = 2270076;Sp( 1000000370 *)≈ 2271816.0 , jdz =sp(m)/G(m) ≈ 1.000766;
G(1000000372) = 2044767;Sp( 1000000372 *)≈ 2045640.8 , jdz =sp(m)/G(m) ≈ 1.000427;
G(1000000374) = 3405467;Sp( 1000000374 *)≈ 3407724.1 , jdz =sp(m)/G(m) ≈ 1.000663;
G(1000000376) = 1710698;Sp( 1000000376 *)≈ 1711434.8 , jdz =sp(m)/G(m) ≈ 1.000431;
G(1000000378) = 1867330;Sp( 1000000378 *)≈ 1866752.3 , jdz =sp(m)/G(m) ≈ 0.999690;
G(1000000380) = 4545539;Sp( 1000000380 *)≈ 4547480.1 , jdz =sp(m)/G(m) ≈ 0.999987;
G(1000000382) = 1703897;Sp( 1000000382 *)≈ 1704021.2 , jdz =sp(m)/G(m) ≈ 1.000073;
G(1000000384) = 1828147;Sp( 1000000384 *)≈ 1828326.6 , jdz =sp(m)/G(m) ≈ 1.000098;
G(1000000386) = 4592599;Sp( 1000000386 *)≈ 4594412.1 , jdz =sp(m)/G(m) ≈ 1.000395;
G(1000000388) = 1821491;Sp( 1000000388 *)≈ 1820001.0 , jdz =sp(m)/G(m) ≈ 0.999182;
G(1000000390) = 2270835;Sp( 1000000390 *)≈ 2271816.1 , jdz =sp(m)/G(m) ≈ 1.000432;
start time =09:33:02,end time=09:33:43 ,time use =
计算式:
Sp( 1000000364 *) = 1/(1+ .1368 )*( 1000000364 /2 -2)*p(m) ≈ 1914805.0 , k(m)= 1.123803
Sp( 1000000366 *) = 1/(1+ .1368 )*( 1000000366 /2 -2)*p(m) ≈ 1703862.0 , k(m)= 1
Sp( 1000000368 *) = 1/(1+ .1368 )*( 1000000368 /2 -2)*p(m) ≈ 3425199.6 , k(m)= 2.010256
Sp( 1000000370 *) = 1/(1+ .1368 )*( 1000000370 /2 -2)*p(m) ≈ 2271816 , k(m)= 1.333333
Sp( 1000000372 *) = 1/(1+ .1368 )*( 1000000372 /2 -2)*p(m) ≈ 2045640.8 , k(m)= 1.200591
Sp( 1000000374 *) = 1/(1+ .1368 )*( 1000000374 /2 -2)*p(m) ≈ 3407724.1 , k(m)= 2
Sp( 1000000376 *) = 1/(1+ .1368 )*( 1000000376 /2 -2)*p(m) ≈ 1711434.8 , k(m)= 1.004444
Sp( 1000000378 *) = 1/(1+ .1368 )*( 1000000378 /2 -2)*p(m) ≈ 1866752.3 , k(m)= 1.095601
Sp( 1000000380 *) = 1/(1+ .1368 )*( 1000000380 /2 -2)*p(m) ≈ 4547480.1 , k(m)= 2.668925
Sp( 1000000382 *) = 1/(1+ .1368 )*( 1000000382 /2 -2)*p(m) ≈ 1704021.2 , k(m)= 1.000093
Sp( 1000000384 *) = 1/(1+ .1368 )*( 1000000384 /2 -2)*p(m) ≈ 1828326.6 , k(m)= 1.073048
Sp( 1000000386 *) = 1/(1+ .1368 )*( 1000000386 /2 -2)*p(m) ≈ 4594412.1 , k(m)= 2.696469
Sp( 1000000388 *) = 1/(1+ .1368 )*( 1000000388 /2 -2)*p(m) ≈ 1820001.0 , k(m)= 1.068162
Sp( 1000000390 *) = 1/(1+ .1368 )*( 1000000390 /2 -2)*p(m) ≈ 2271816.1 , k(m)= 1.333333
Sp( 1000000392 *) = 1/(1+ .1368 )*( 1000000392 /2 -2)*p(m) ≈ 3407724.1 , k(m)= 2
使用修正系数,能够使得连乘式计算值的相对误差绝对值比较小,由于偶数的素对数量中存在不可计算的部分,(即小素数小于√(M-2)的素对),因此呈现精度在1上下变化时就是达到使用修正系数极限了,如果要控制相对误差呈现负值的情况,就需要把修正系数的μ值适当调大一些,以得到偶数素对的下界计算值。
计算的下界计算值:
G(1000000390) = 2270835;inf( 1000000390 )≈ 2264247.4 , Δ≈,infS(m) = 1698185.53 , jd(m)= 0.997099;
G(1000000392) = 3407509;inf( 1000000392 )≈ 3396371.1 , Δ≈,infS(m) = 1698185.53 , jd(m)= 0.996731;
G(1000000394) = 1804040;inf( 1000000394 )≈ 1797021.7 , Δ≈,infS(m) = 1698185.54 , jd(m)= 0.996110;
G(1000000396) = 1736785;inf( 1000000396 )≈ 1731483.3 , Δ≈,infS(m) = 1698185.54 , jd(m)= 0.996947;
G(1000000398) = 3405959;inf( 1000000398 )≈ 3396901.7 , Δ≈,infS(m) = 1698185.54 , jd(m)= 0.997341;
G(1000000400) = 2886129;inf( 1000000400 )≈ 2877079.2 , Δ≈,infS(m) = 1698185.55 , jd(m)= 0.996864;
G(1000000402) = 1704261;inf( 1000000402 )≈ 1698185.6 , Δ≈,infS(m) = 1698185.55 , jd(m)= 0.996435;
G(1000000404) = 3718140;inf( 1000000404 )≈ 3705132.1 , Δ≈,infS(m) = 1698185.55 , jd(m)= 0.996501;
time start =11:37:17 ,time end =11:37:40 ,time use =
计算式:
inf( 1000000390 ) = 1/(1+ .1406 )*( 1000000390 /2 -2)*p(m) ≈ 2264247.4
inf( 1000000392 ) = 1/(1+ .1406 )*( 1000000392 /2 -2)*p(m) ≈ 3396371.1
inf( 1000000394 ) = 1/(1+ .1406 )*( 1000000394 /2 -2)*p(m) ≈ 1797021.7
inf( 1000000396 ) = 1/(1+ .1406 )*( 1000000396 /2 -2)*p(m) ≈ 1731483.3
inf( 1000000398 ) = 1/(1+ .1406 )*( 1000000398 /2 -2)*p(m) ≈ 3396901.7
inf( 1000000400 ) = 1/(1+ .1406 )*( 1000000400 /2 -2)*p(m) ≈ 2877079.2
inf( 1000000402 ) = 1/(1+ .1406 )*( 1000000402 /2 -2)*p(m) ≈ 1698185.6
inf( 1000000404 ) = 1/(1+ .1406 )*( 1000000404 /2 -2)*p(m) ≈ 3705132.1
|
|