|

楼主 |
发表于 2019-1-8 13:08
|
显示全部楼层
公式0.660161815846869*∏((P-1)/(P-2))*(n前素数个数+n前2P的个数)^2/n, 当n是偶数时乘10/9,奇数时乘8/9。连乘积中的P是根号n前的素数,且能整除n。 这样计算出来的结果总累计值大于实际统计值,在10000以内平均多1.6个。10000以内统计累计值为2044424组解,而公式结算结果累计值为2060829组解,相差16405组解(绝对误差),相对误差0.8%.
如果用素数定理代替n前符合条件的元素个数,则计算结果偏差很大,10000以内公式计算值累计为1585948组解,公式为0.660161815846869*∏((P-1)/(P-2))*((n/2)/ln(n/2)+n/ln(n) )^2/n,当n是偶数时乘10/9,奇数时乘8/9,连乘积中的P是根号n前的素数,且能整除n。
把后边用素数定理代替符合条件元素的个数的公式,稍微处理变通一下可以用0.660161815846869*9/4*∏((P-1)/(P-2))*n/(ln(n)-ln2)^2,当n是偶数时乘10/9,奇数时乘8/9,连乘积中的P是根号n前的素数,且能整除n。
用此公式得到的解的组数累计值为1766563,偏差仍然很大,究其原因是素数定理不能很好的反应素数的个数。
|
|