|
[watermark] 与平方根数关联的数论公式
数内素数个数约为0.5(√N)π(√N)或(√N)∏[h/(h-1)]
一个公式可直观探察素数多少,另一公式可明确递增特性,有特色。
数N内素数个数的符号:π(N),(各参数)连乘的符号:∏,平方根符号:√,
几种求数内素数个数公式的转换:
由素数定理推出的公式:π(√N)≈(√N)/Ln(√N)≈(√N)∏[(p-1)/p],
π(N)≈N/LnN=(√N)(√N)/[2Ln(√N)]≈0.5(√N)π(√N),
π(N)≈N∏[(P-1)/P]≈(√N)(√N)/∏[P/(P-1)]≈(√N)∏[h/(h-1)],
递减方式公式与递增方式公式的关系。
(√N)(√N)``2``3``4``5``6```````√N-2```√N-1``1
----------·-·-·-·-·-·....·-----·-----≡--
(√N)·2....3..4..5..6..7.......√N-1....√N...1
``N````2```/`/`/``/````√N
-----·-·....·----≡-----
..2...././././..√N.....1
N```2``4``5``10`````P-1``√N`4``6``8``9````````h
--·-·-·-·--·..·---≡--·-·-·-·-·..·---
2...3..5..7..11......P....1..3..5..7..8.......h-1
例如:
100``2``4``6``800````````10``4``6``8``9``10
---·-·-·-==---=22.86==--·-·-·-·-·---
.2...3..5..7..35..........1..3..5..7..8...9
可直观探察素数的公式: 0.5(√N)π(√N)
以平方根为单位表示数中素数的多少,约等于平方根内素数
个数的一半。举例说明:
100内素数的个数约等于10内素数个数4的一半。20+5个。
400内素数的个数约等于20内素数个数8的一半。80-2个。
1000内素数的个数约等于31.6内素数个数11的一半。173-5个。
10000内素数的个数约等于100内素数个数25的一半。1250-20个。
用100来深入一下平方根内素数个数的一半的细节。
把100分在10·10的行列中,|100/2|个偶数占{2,4,6,8,10}的5列,
|50/3|个数占{3,9}的1.6列,|34/5|与|28/7|个数占{5}的列,
没被占用的列有{1,7}两列,作为筛除掉合数用到的4个素数还需补解上。
素数个数的一半是指:平方根的后半区域的素数个数再加上1个。
最后,作为筛除掉合数用到的一些素数还需补在解上。
用平方根的后半区域的素数个数求解数内素数个数的下限解最合适。
明确递增特性的素数个数公式:(√N)∏[h/(h-1)]
h表示小于(√N)的所有合数。“4,6,8,9,10,12,14,15,..,"
素数个数约等于的平方根数与(诸合数与合数减一的比)的连乘积。
``````````````4``6``8``9``10``12``14``15`````````h
π(N)≈(√N)·-·-·-·-·---·--·-·--......·---
..............3..5..7..8...9..11..13..14.......(h-1)
递增特性:以平方根数为基底数,每多一项,都要增加分母数份的一份。
举例说明:
100内素数的个数约等于把10个,加10/3个素数,得13.4。加13/5个素数,
得16个。加16/7,得18.3。 加18.3/8,得20.6。
加22.6/9,得24.3。
````````````10`13.3`16``18.3``20.6
π(100)≈10+--+----+---·---·----
.............3...5...7....8....9
π(100)≈10+3.3333+2.6666+2.2857+2.2857+2.2857=22.857
实际分布:个数大于两个平方根数,素数个数与合数个数同步增大。
.....10+2+3+2+2+3+2+1=25
2,11,23|31|41|53|61|71|83|97|
3,13,29|37|43|59|67|73|89|
5,17,..|..|47|.....|79|
7,19,..|
青岛 王新宇
2009.9.26
分析数论单筛公式,双筛公式
老文章:“20041221分析数论单筛公式,双筛公式 ”
现在,用一个具体数来分析单筛,双筛公式,看看如何精确它。
本文用到我的“首创筛法确定实际解下限的方法”一文的内容。
求素数的个数的方法,
取偶数N为“962”,(说明,本文中素数的个数是指奇素数的个数)
1,先求出给定数以内的不包括“1”的奇数的总个数,符号“Q”。
(962/2)-1==480,表示962中,不包括“1”的奇数的总个数有480个。
2,按照奇素数分类求出,没有重合的,实际真实奇合数的数值。
480含有奇素因子“3,5”的类型,
含素因子3的奇合数的个数有(480-3)/3==159,其他有(480-159-1)=321个,
含素因子5的奇合数的个数有(320-5)/5==63,其他有(321-63-1)=357个。
480没含有奇素因子的奇素数类型,素数符号“P”
利用各级素数的互素表求,等于[P,(961/P)]区间的互素数个数。
含素因子7的奇合数的个数,每30有8个数,区间为[7,137.2],有36个互素数。
其他小素数的奇合数,利用在[P,(961/P)]区间的互素数表等于素数表求:
11至87.4之间有19个素数,含素因子11的奇合数的个数==19,
13至73.93之间有16个素数,含素因子13的奇合数的个数==16,
17至56.6之间有10个素数,含素因子17的奇合数的个数==10,
23至41.8之间有5个,23·{23,29,31,37,41},含23的奇合数的个数=5。
29至33.2之间有2个,29·{29,31},含素因子29的奇合数的个数==2。
31仅自己1个。31·,含素因子31的奇合数的个数==1。
下面是按照素因子分类的实际真实奇合数的数值。
奇素因子=3````5````7```11``13``17``19``23``29``31..类型
奇合数==159...63...36..19..16..10...8..5...2...1...个数
3,分级求解:求出与3,与5,与7,...,与31,都互素的数的个数,
素数个数==不含1的奇数的个数连减各类型素因子的奇合数的个数。
S==480..-159.-63.-36.-19.-16.-10..-8..-5..-2..-1
各级解..,321,258,222,203,187,177,169,164,162,161
连减公式,通过分级求解,很容易改写成,分式连乘公式。
````````960`321`258`222`203`187`177`169`164`162`161
S(962)==-—·—·--·--·--·--·--·--·--·--·-—
.........2..480.321.258.222.203.187.177.169.164.162
把各项的分母变换成素因子,分子也恒等变换,
即改成((素因子减“1数”)/素因子)形式,
```````960``2.`4.`6`10.05`11.98`16`18.1`22.3`28.6`30.8
S(962)=-—·-·-·-·---·----·--·---·---·----·---—
........2...3..5..7..11...13....17..19...23...29...31
各项“1数”为
|0.994.|0.984|0.972|0.95|1.02|0.9|0.86|0.7|0.36|0.19|
....3....5.....7....11....13...17..19...23..29....31
数据说明:单筛公式,“1数”取“1”,会有误差。
把素因子素数取N开方数以内的素数,
解值误差是很大,但适当去掉后面的项,解接近正确解。
对素数的连减公式的{连减系列数},进行恒等变换,求解。
S==480--{159+63+36+19+16+10+8+5+2+1}
===480--{159+62+34+16+12.+7+4..+2
..........+1.+2.+3++4.+3.+4+5..+1}
===480--160--64-37-20-15-11-9-3
各级解为320,256,160,64,37,20,15,11,9,3
将连减公式“S==480-160-64-37-20-15-11-9-3”化为连乘公式。
````````960`320`256`219`199`184`173`164`161
s(961)==-—·—·--·—·—·--·--·--·—==161
.........2..480.320.256.219.199.184.173.164
把各项的分母变换成素因子,分子也恒等变换,
即改成((素因子减“1数”)/素因子)形式,
```````960``2.`4.`6`10``12``16``18``22.58
S(962)=-—·-·-·-·-·--·--·--·---—
........2...3..5..7.11..13..17..19..23.
其中最后的素因子分数项停止在“22/23”这一级:即
误差项等效于后面几项,适当去掉后面的项,解接近正确解。
实际解大于公式解,可作为上限公式用,
单筛公式就是素数个数的上限公式(3)的主项,最后再加上附加项。
````````961````2``4```6```10```12```16```18```22```28``30
π(961)=——·—·—·—·—·—·—·—·—-·—·—
.........2.....3..5...7...11...13...17...19...23...29..31
各级解为480)320)256)219)199.7)184.3)173)164.3)157)152)147+11
各素因子减数为160,64,37,20,15,11,9,7,5,5
实际解........159..63..36..19..16..10..8..5..2..1.
...............-1..-1..-1..-1..+1..-1.-1.-2.-3.-4
(3)的附加项,是用“+√N”改善了误差,理论上应加√N内素数个数。
由奇数的个数起逐级减小,逐级接近到素数个数的真值。
进一步的研究
逐级减少是素数个数的上限公式,可变换成素数个数的下限公式,
根据,“N”除最大素因子~√N,少了末项素因子,各分子可左移一项;
``````````````高P-1``√ N``4``6`10`12`16`18`22`28`
π(N)>√ N∏(———)=——·-·-·-·-·-·-·-·-·..(10)
................P......2...3..5..7.11.13.17.19.23
````````√961````4````6```10```12```16```18```22```28```30
π(961)=———·—-·—-·—-·—-·—-·—-·—-·—-·—
...........1.....3....5....7...11...13...17...19...23...29
各级解为.31)41.33.)49.6)70.85)77.2)95.1)100.7)116.6)142)147+11
素数个数的下限公式,由数的开方数起,逐级增加到素数个数的真值,
两个公式的解是同一个数,因为是同一个公式的变换,都有附加项,
素数的个数肯定在“上限,下限”两公式解的中间,
逐级求解,因为双向解是内夹解,两解相交点,必是正确解。
双向求解是动态解决多变误差参数的好方法,误差会相补抵消。
待续
青岛 王新宇
2004.12.21
20041222qdxinyu分析数论双筛公式
下面,用一个具体数来分析双筛公式,看看如何精确它。
本文用到我的“首创筛法确定实际解下限的方法”一文的内容。
求对称素数的个数的方法,
求对称素数的个数的公式与求孪生素数的成对的对数的公式
是一样的。故先求孪生素数的成对的对数。
取偶数N为“962”,(说明,本文中素数的个数是指奇素数的个数)
1,先得知道给定数以内的不包括“2”的素数的总个数,符号“S”。
S(962)=161,表示962中,不包括“2”的奇素数的总个数有161个。
2,按照奇素数分类求出,实际真实的伴生奇合数的数值。
伴生奇合数是我新定义的一种奇合数,简称“伴数”
它是与的奇合数。
求解31的平方数以内的各级伴数的方法另介绍,先用结果,
按照奇素数分类求出伴数,奇素数按渐大的顺序,求到顺序数。
``在“(6n+1),(6n-1)”数系中,
(1)“3”素数级的伴数个数是零,没有伴数。
(2)其他素数级的伴数个数:
“5”素数级的伴树的集合(略),伴数个数有44个,
“7”素数级的伴树的集合(略),伴数个数有18个,
11·{19;23,29,37;41,59,61;67;83,}“11”素数级伴数有9个,
13·{13,17;29;31,37,43,47;53;59;73,}“13”素数级伴数有10个,
17·{23,29,37;47,}“17”伴数有4个,伴素数389,491,631,797
19·{19,23;31;37,}“19”伴数有4个,伴素数361,439,589,701
23·{37;41,}“23”级伴数有2个,伴素数853,941
29·{29,}“29”级伴数有1个,伴素数29·29-2==839
31·{空级}“31”级伴数是零,没有伴数,
有了伴数,很容易得到连减公式,连乘公式.
孪生素数个数==奇素数个数连减各级伴数。孪生素数符号"L"
S==480-159--63--36--19--16--10--8--5--2--1==161
L==161--0---44--18---9--10--4---4--2--1--0==69
各级解:161,117,.99,.90,.80,.76,.72,70,69,..=69
连减公式,通过分级求解,很容易改写成,分式连乘公式。
各级解={161,117,99,90,80,76,72,70.69,69,}
分式连乘公式
```````````117``99`90`80`76`72`70`69`69
L(962)=161·--·--·-·-·-·-·-·-·-
...........161.117.99.90.80.76.72.70.69
对应的级.3...5...7..11.13.17.19.23.29.31
公式表示:“961以内的孪生素数有69个,”等于实际数。
把各项的分母变换成各(级素数-1),分子也恒等变换,
即改成((级素数减“2数”)/(级素数-1))形式,
````````````3.`5.`9``10.7`15.2`17.05`21.4`27.2``30
L(962)=161·-·-·-·----·---·----·---·---·--
............4..6.10...12...16....18...22...28...30
各项“2数与1的差数”为
`````````|1.093|0.923|1..|1.3|0.8|0.95|0.6|0.8|0.|
数据说明:双筛公式,“2数与1的差数”取“1”,会有误差。
把级素数取N开方数以内的素数,
解值误差是很大,但适当去掉后面的项,解接近正确解。
对孪生素数的连减公式的{连减系列数},进行恒等变换,求解。
L==161--0---44--18---9--10---4---4----2----1---0==69
===161--0---40--21--10-7.5---5.1-4.3--3--1.1---0==69
各级解===={161,121,100,90,..82.5,77.4,73,70.,69}
各级缩比为..4,..6,..10,12,..16.1,18..,24,64,..1
分式连乘公式:代入各级缩比,23/24=21/22,63/64=27.56/28
```````161``3``5``9``11``15`17``21.1`27.56``1
L(962)=-—·-·-·-·-·--·--·----·—--·—
........1...4..6.10..12..16.18..22....28....1
其中最后的级素数项停止在“27/28”;少“29/30”项。表示
误差项等效于后面的项,适当去掉后面的项,解接近正确解。
实际解大于公式解,可作为上限公式用。孪生素数个数转换为
孪生素数对数,要乘“1/2”,等效于级素数含“3”
孪生素数对数公式如下:
S个素数包含的孪生素数的个数:其中,P为起始部分的素数;
`````````(P-2)`````1``3``5``9`11`15`17`21
L(S)~S∏(---)==S·-·-·-·-·-·-·-·-...孪生素数对数
.........(P-1).....2..4..6.10.12.16.18.22
由奇素数的个数起逐级减小,逐级接近到孪生素数对数的真值。
进一步的研究,将素数个数公式代入:
`````````(P-2)``N```(P-1)``(P-2)``N```(P-2)
L(S)~S∏(---)==--∏(---)∏(---)==--∏(---)
.........(P-1)..2...(.P.)..(P-1)..2...(.P.)
上面是奇数逐级减少到孪生素数对数的上限公式,可变换成下限公式,
根据,“N”除最大素因子~√N,少了末项素因子,各分子又可左移一项;
``````√ N```高P-2``√ N``3``5``9`11`15`17```√ N``(奇合数)
L(S)~----∏(——-)=——·-·-·-·-·-·-..=----∏——-----)
........2......P......2...3..5..7.11.13.17....2...(奇合数-2)
这就是半开方数逐级增加到孪生素数对数的下限公式。
孪生素数对数的下限公式,由数的半开方数起,逐级增加到解的真值,
两个公式的解是同一个数,因为是同一个公式的变换,
孪生素数成对的对数肯定在“上限,下限”两公式解的中间,
逐级求解,因为双向解是内夹解,两解相交点,必是正确解。
双向求解是动态解决多变误差参数的好方法,误差会相补抵消。
下面,用一个具体数来介绍求对称素数的个数的方法,
相对于偶数的两个对称位置都是素数的情形,有两种类型。
第一种类型是:与开方数以内的素数相对称的位置是素数,少,
第二种类型是:与开方数以外的素数相对称的位置是素数,多。
双筛公式:不含第一种类型的解,只计算占主体解的第二种类型解。
`````````1``3-r3``5-r5``7-r7``11-r11```````P-rP``````p-rp
G(x)~x·-·----·----·----·------·....·----·..·-----
.........2...3.....5.....7......11...........P.........p
表示x大约有G(x)个主体对称的素数。与首尾对称的素数没计入。
其中:P表示不大于x开方数的诸有效素数,p为P中的最大的有效素数。
(注意rP的P是下角标 , 不是数) r3,r5,...rp为对应于P的删除比例,
是x素因子的,rP选1; 不是x素因子的,rP选2 ;
大素数时,按实际的删除比例修正。
偶数=962,由962=2*13*37知道,962含开方数以内的素因子13,
``````````````1```1```3```5``9`12`15`17`21`27`29
G(962)===962·--·--·--·-·----·-·-·-·-·-·--
..............2...3...5...7.11.13.17.19.23.29.31
逐级减少的解为480,160,96,68,56,51,45,41,37,34,32,
表示偶数962约有32个对称的素数, 事实正好是32个对称的素数,如下:
43,919,79,883,103,859,109,853,139,823,151,811,
193,769,211,751,229,733,261,701,271,691,331,631,
349,613,421,541,439,523,463,499,
公式中:分母是2和13的分子减一,其他项的分子减二。
其中,分子与分母的差不大于2。分母是小于该偶数开方数的诸素数。
对称素数的个数的公式与孪生素数对数的公式区别是:分数项的
分子增大1个数的项随着素因子的增多而增多。即有:素因子增大系数。
对称素数的个数等于素数个数和孪生素数对数之间,
用素因子增大系数确定的比例上。
各种比例数的数量按2的幂递增;
2种,4种,8种,16,32,64,128,256,....
若只考虑解的有无,则只考虑对称素数的个数的最小解就可以了。
青岛 王新宇
2004.12.22
20041223qdxinyu分析数论双筛公式(续)
检查了一下,更正一下伴数的数据,上面公式的
L==161--0---44--18---9--10--4---4--2--1--0==69 应该是
L==161--0---41--20---9--10--5---4--2--1--0==69
虽然中间参数都该相应纠正,但原计算公式最后的解值不变。
其孪生素数如下所示;为69个。
...,..3,..5,..7,.11,.13,.17,.19,
.29,.31,.41,.43,.59,.61,.71,.73,
101,103,107,109,137,139,149,151,
179,181,191,193,197,199,227,229,
239,241,269,271,281,283,311,313,
347,349,419,421,431,433,461,463,
521,523,569,571,599,601,617,619,
641,643,659,661,809,811,821,823,
827,829,857,859,881,883,
对称素数的公式与孪生素数对数的公式解的值的区别是:
对称素数有:素因子增大系数。
相对于偶数的两个对称位置都是素数的情形,有两种类型的解。
第一种类型是:与开方数以内的素数相对称的位置的解,少,
第二种类型是:与开方数以外的素数相对称的位置的解,多。
双筛公式:不含第一种类型的解,只计算占主体解的第二种类型解。
偶数962的双筛公式的解,不含{3,5,7,11,13,17,19,29,31},少九个,
孪生素数个数减九个,为30对,
偶数962含素因子13, 素因子增大系数为“(13-1)/(13-2)” ,
数论的解常默认要取整运算,
偶数962对称素数的对数等于(30+2)对是合理的,正确。
对称素数就是符合哥德巴赫猜想的两素数和等于偶数的表达式的个数。
对称素数大于1就是哥德巴赫猜想的证明。
待续
青岛 王新宇
2004.12.23
以上均为草稿,希望谅解.
[/watermark] |
|