|

楼主 |
发表于 2011-1-16 16:46
|
显示全部楼层
[原创]素数个数的公式和第m个素数出现的范围的探索
[这个贴子最后由ysr在 2011/01/24 11:38am 第 5 次编辑]
准确素数个数公式推导过称:由已知素数P及P以内的全体素数,可用如下公式求出P^2之内的素数个数,其中i为P以内的素数个数,素数个数=((P^2+1)/2)*(1-1/3)*(1-1/5)*……*(1-1/P)+i-1,减1是为了去掉数字1,如下为结果与实际比较:
9之内,公式=3+1-1=3,实际=3,连乘积=2/3,
25之内,公式=6.9+2-1=7.9,实际=8,连乘积=8/15,
49之内,公式=11+3-1=13,实际=14,连乘积=48/105,
121之内,公式=25+4-1=28,实际=29,连乘积=480/1155,
……
利用李联忠老师的部分等式,上述公式可等效为:
x/(a*lnx)+2x^(1/2)/lnx,其中a=0.908402461为 Π(1-1/(P-1))/lnx的倒数,其中x>=1000,是平均值不是极限值,
由李老师的等式,当x→∞时,limπ(x)=limx/2*∏(1-1/p)得,上限公式π(x)=2x/lnx,此公式远大于实际个数的波峰值,乘以欧拉常数可得接近实际上限的公式π(x)=2γx/lnx,此公式较实用,
除以欧拉常数得略大于实际波峰值的公式π(x)=x/lnx/γ,
没有严格证明,验证如下:当x=2000时,2x/lnx=526,2γx/lnx=303,x/lnx/γ=459,实际=303
当x>=2000时,上述公式的关系:
x/lnx<x/(a*lnx)+2x^(1/2)/lnx<=实际<=2γx/lnx<x/lnx/γ,
第m个素数出现范围的公式推导:
由于m=x/lnx,素数概率为1/lnx,则素数f=m*lnx,由于lnx,为lnm的延续,近似等于lnm,故f=m*lnm,素数概率的上限为1/lnx/γ,素数出现的上限为f=m*lnm/γ,
上下限的平均值为:
f=m(lnm+1),f误差范围为+-2f^(1/2)+1
由于素数必为非平方数,由非平方数公式(n+x)^2-(2x+1),其中x为平方根的整数部分,
则得f误差范围为f+-2f^(1/2)+1,此范围较接近实际,
更大的范围为f1与f2之间,其中f1=m*lnx,f2=m*lnm/γ,范围太大不实用. |
|