数学建模社区-数学中国

标题: 论给定区间素数的分布规律公式 [打印本页]

作者: tysh670407    时间: 2015-7-24 10:34
标题: 论给定区间素数的分布规律公式
本帖最后由 tysh670407 于 2016-6-6 07:35 编辑 + L5 u) ?% Y; o' U% [$ o

7 V# S6 c2 O8 V$ M
论给定区间素数的分布规律公式5 ?6 e6 u! f! m6 m8 ?$ j2 `
: D1 M6 n7 }  _- A' m
田永胜
(内蒙古自治区  吉兰泰  750333
        摘要:通过对自然数按照一定方向旋转排列,找到了自然数的等势区间并集,并对每个区间的素数分布情况进行研究,给出了在给定区间内素数的分布定理、公式及推论。
        关键词 自然数;螺旋排列;给定区间;素数分布;规律;
        引言
        自然数沿数轴方向排列时,素数的分布没有规律可循;当把自然数按一定的方向旋转排列时,素数的分布就变得有规律。下面揭示它的分布规律。
        1 自然数的排列规律
        首先,按逆时针方向把自然数进行排列,如下图:
                                                                自数数排列图20150423.png
$ _& O% o" V* p; X5 {7 y) c' l) H2 K: U9 Y
自然数螺旋排列图
        从上图可以看出,自然数集合N+也可以由一连串连续区间的并集组成,[1]∪(19]∪(9,25]∪(25,49]∪(49,81]∪(81,121]∪……∪((2x-3)^2,(2x-1)^2]…。并且,每个区间的最大数都是奇数(2x-1)的平方。
3 u/ u# e  n* }/ P* W, j$ {
        2 素数分布定理和公式0 G5 ?0 Z6 _0 G5 P
        首先,来研究每一区间数字的素数分布情况:
第一区间只有自然数1,素数个数为0
第二区间为(19],有8个数字,其中素数有4个,所占比例为 4/8=0.5
第三区间为(1025],有16个数字,其中素数有5个,所占比例为 5/16=0.3125
第四区间为(2549],有24个数字,其中素数有6个,所占比例为 6/24=0. 25;以此类推。
其次,再来看每一个区间的素数分布与区间内的数有什么内在规律。1在中心,不是素数;在区间(19]有8个自然数,最大数是9,求9的自然对数的倒数,1/ln90.455,与该区间实际素数所占比例接近;乘以总数8,值约等于3.64,取整数后为4,与该区间实际素数个数相同。在区间(1025]有16个自然数,最大数是25,求25的自然对数的倒数, 1/ln250.311,与区间内实际素数所占比例0.3125很接近,乘以总数16,值约等于4.97,取整数后为5,与该区间实际素数个数相同。在区间(2549]有24个自然数,最大数是49,求49的自然对数的倒数, 1/ln490.2569,与区间内实际素数所占比例0. 25很接近,乘以总数24,值约等于6.16,取整数后为6,与该区间实际素数个数相同。以此类推,如素数分布规律表所示。
5 x; K: C- y8 ?                                                          素数分布规律表20150423.png : N  o1 U! z* K- E+ g' X
素数分布规律表
        由上表可以看出,在第2到第8区间,实际素数个数与理论素数个数相等,其他的区间实际素数个数在理论素数个数左右波动,每个区间实际素数的所占比例和理论素数分布密度非常接近。
下面,给出素数分布定理的一般形式。
定理
        设x为自然数,在给定区间((2x-3)^2,(2x-1)^2]内,素数的分布密度公式为
1/ln2x-1^2
        给定区间内自然数的个数为  
2x-1^2-2x-3^2=8x-8
        用π(x)表示给定区间内的素数个数,则给定区间素数个数与自然数的个数之间存在如下线性关系
π(x)=( 8x-8)/ ln2x-1^2
        若用Sn表示n圈内素数的总和,则
, X  @- n. c  Y7 r                                                         123.PNG + g6 i! ?1 Q- x
        推论1 在区间((2x-3)^2,(2x-1)^2]内,只有有限个素数,当x趋向无穷大时,素数也趋向无穷大,即 + r9 i7 [+ P8 j7 `% o# n

- h: X6 p4 t; o7 M% ~; I                                                          456.PNG
        接着,再来看每一个区间的孪生素数的分布情况:在区间(19]内有2357两对孪生素数,在区间(925]内有11131719两对孪生素数,在区间(2549]内有29314143两对孪生素数,在区间(4981]内有59617173两对孪生素数,在区间(81121]内有101103107109两对孪生素数,在区间(121169]内有137139149151两对孪生素数,在区间(169225]内有179181191193两对孪生素数,每一区间内被小于或等于(2x-1)的素数约去后,都有两对孪生素数。因此,得出推论在每一个区间至少有两对孪生素数。
8 ?! P7 }  s/ P/ b9 H' g
+ B0 i7 ?& n9 k' N; d! o
        推论2 在区间((2x-3)^2,(2x-1)^2]内至少有两对孪生素数。当x趋向无穷时,孪生素数也趋向无穷。

( U7 i! W, H3 c+ v2 U' K
4 v9 D- }& E, j: E3 c
5 k7 T. X4 _6 i1 z+ W, V
( E% ?0 v. V+ Z/ {' e' W/ l
        推论3
在区间((2x-3^2,(2x-1^2]内,实际素数个数总是在理论素数个数左右波动,即它们的比值在1左右波动,当x取有限数值时,所有区间实际素数与理论素数之比(π(x)/(8x-8/ln(2x-1)^2))的平均值趋向1。当x取无穷大时,无穷区间实际素数与理论素数之比(π(x)/(8x-8/ln(2x-1)^2))的平均值等于1。即x∞时,
{π(1)/ [8×1-8/ln(2×1-1)^2]+ π(2)/ [8×2-8/ln(2×2-1)^2]+
π(3)/[8×3-8/ln(2×3-1)^2]+…+π(x)/ [8x-8/ln(2x-1)^2]}/1+2+3+…+x=1

( B: p% W0 R0 I1 ~: R; B( G9 l& |6 X( \  B- |1 {# Q
$ \4 a8 L5 M* _5 T! ^# l1 m

( u- _  C# U. h1 v0 p, a# y, p/ j1 Y) y& B! w

% A. p5 ^- p( Q. Y8 Y5 b+ ?  e8 p
3 C4 }5 ^! ?& h' e7 c  I

% T' P; o; a; _7 e; Z% X6 n, @3 T+ M

8 o& V9 |7 a& ^) R& a0 k$ K$ R/ {  L- H: L  O: E. J/ b- }" B) A

$ q& s0 Z* E% L5 q8 U7 e9 \0 d
作者: tysh670407    时间: 2015-8-7 15:51
本帖最后由 tysh670407 于 2016-6-2 09:22 编辑
, w' [, u( z( p3 W0 V) S4 Q, C- ?! ~5 P

* M2 y6 G) D" W: F: V0 b) z10031至10040区间素数分布情况,实际素数在理论素数左右波动,实际素数总数40541,理论素数总数40515,差值26,占比26/40515=0.000641737。单个区间最大偏差(4111-4052)/4052=0.0145607。见附图。
# T% I4 u) V! z6 h# h- {5 D. p) B. o# E

: [2 W6 l% v7 U2 O# z
8 l1 L  e9 ~" @+ A4 w
# O5 o7 [5 a# L
+ ^- ^. }$ G7 g0 O. P& e
% p; G1 |3 e: X1 w6 H% v2 c
; \; [, Z  i0 p! s, {

10031-10040区间素数分布.PNG (16.25 KB, 下载次数: 338)

10031-10040区间素数分布.PNG


作者: tysh670407    时间: 2015-8-10 21:05
本帖最后由 tysh670407 于 2016-6-2 09:30 编辑
2 F9 u; m. p# B
8 R/ g' f3 }* C( |, |6 [1021至1030区间素数分布情况,实际素数在理论素数左右波动,实际素数总数5374,理论素数总数5374.单个区间最大偏差(562-537)/537=0.04655493,见附图。
6 ~. h* J' O/ ^2 V4 Q- f- b  H% z2 H1 x% `
0 V& ?/ z9 f( f) U

# u0 h( Q/ C! \3 {/ |" c
& y6 J- b! X9 P0 T/ Q& `

1021-1030素数分布表.PNG (15.42 KB, 下载次数: 326)

1021-1030素数分布表.PNG


作者: tysh670407    时间: 2015-8-20 11:12
本帖最后由 tysh670407 于 2015-9-1 21:50 编辑
6 ]: F4 y# P" g) R, J2 J
+ `. v: |; i6 ^5 K给定区间素数的分布问题,其实就是素数的筛选问题,我们不可能一下子解决所有的素数的筛选问题,只能分区间来进行,这样问题就好分析了。自然数集可以表示为一系列连续区间的并集,分别对这些连续区间的素数分布情况进行研究,就会发现一些规律和公式。5 \( {6 Z. ], i8 f- F; y2 V
对于给定区间((2x-3)^2,(2x-1)^2],素数的多少可以用小于或等于2x+1的素数去约,剩余的就是素数,用素数筛选后素数所占的比例是多少呢?
1 V2 E; `; A9 L, z, s+ u设给定区间数的比例为1,
) q( z, y: H+ Z& F被2除后余1-1/2=1/2,
+ d2 W4 A/ ?2 D5 a1 V被3除后余1/2-1/2*1/3=1/3,/ a- ?0 Q2 ]8 L0 S' a- N9 I0 P
被5除后余1/3-1/3×1/5=4/15,
" q3 o5 x/ p) r( L" {3 P" W/ H/ R# Q被7除后余4/15-4/15×1/7=24/105;
6 q8 P5 L% @" D9 i( Z被11除后余24/105-24/105×1/11=240/1155;
( J0 D7 P  X8 ]3 v( k9 R6 c  [被13除后余240/1155-240/1155×1/13=2880/15015;
# J, _9 w( U) s& L8 R9 T1 O被17除后余2880/15015-2880/15015×1/17=46080/255255;! @; k# J" G- R1 Y
被19除后余46080/255255-46080/255255×1/19=829440/4849845;
0 Q( B/ t: ]0 v' P; b: h& g8 I被23除后余829440/4849845-829440/4849845×1/23=18247680/111546435;! j9 ~; L/ Q" i; Y: E: X, e
以此类推。8 d5 N# q4 G8 f/ K

4 D( I5 N: |. @; \- U* {, ~& T第2区间(1,9],自然数有8个,被2约后剩余1/2,8×1/2=4;5 V. l# X  F! }' D! G( k
第3区间(9,25],自然数有16个,被3约后剩余1/3,16×1/3=5.33,取整为5;+ |& A  W5 C$ T# Y7 t: e
第4区间(25,49],自然数有24个,被5约后剩余4/15,24×4/15=6.4,取整为6;
% V7 O% `# F8 L1 G# u第5区间(49,81],自然数有32个,被7约后剩余24/105,32×24/105=7.31,取整为7;
* D6 I4 R: G) S2 V第6区间(81,121],自然数有40个,被11约后剩余240/1155,40×240/1155=8.31,取整为8;& _5 w% u9 \% [; Z
第7区间(121,169],自然数有48个,被13约后剩余2880/15015,48×2880/15015=9.21,取整为9;第8区间(169,225],自然数有56个,被13约后剩余2880/15015,56×2880/15015=10.74,取整为10;
6 e1 L: _( x6 h- Y第9区间(225,289],自然数有64个,被17约后剩余46080/255255,64×46080/255255=11.55,取整为11;1 D" O6 I' v) V% i# Z6 J& {
第10区间(289,361],自然数有72个,被19约后剩余829440/4849845,72×829440/4849845=12.31,取整为12;
2 J3 s, p: P  z0 U0 w  以此类推。9 h7 z7 j/ E. B. [- o; U; _
2 v. {$ i8 G4 Z8 y- j

% Q& x* q) t& X8 L4 i8 l
作者: tysh670407    时间: 2015-9-1 21:34
本帖最后由 tysh670407 于 2018-6-29 10:10 编辑 9 b4 |6 a) j) u3 z8 D6 b

8 l$ r3 {, |; d    突然发现素数是有限的,当对素数的分布密度1/ln(2x-1)^2求极限时,即x→∞时,lim1/ln(2x-1)^2=0,也就是说,无穷区间的素数概率为零,素数也为零,即( 8x-8)/ ln(2x-1)^2=0,因此,推论1的极限应该是0,而不是∞,所以得出1到∞区间的素数的总和是有限的结论。' W1 e/ L, X0 A4 ~' K
6 ?1 X# H3 ]" i* ^

作者: tysh670407    时间: 2015-9-1 21:48
又一个数学悖论,素数真的是无穷的吗?2 P; e* @2 U9 e" E+ N7 V$ d

作者: tysh670407    时间: 2015-10-8 21:50
本帖最后由 tysh670407 于 2018-6-29 10:11 编辑 # {" C" H: Q2 G3 x3 h
5 r( ]8 w& k1 y4 ^
想了一个月,突然想到,如果引入无穷小的概念,就可以解决这个问题。
/ q% O: Q- o0 w# [/ L7 s
5 l" j5 O# k* B! X, S. |$ a    因为1位于无穷个区间的中心,数字按照等角螺旋进行排列,无穷大的倒数自然就是无穷小了。如果我们用符号⊙表示无穷小,那么" Z7 y3 M7 w" L0 y( ~8 b+ u% T
lim1/ln(2x-1)^2的极限值就等于⊙。
% E' x- G! s% h( P; D5 R
0 f) p, w2 [* w; s6 W% J    这个结果验证了我们常说的一句话:宇宙其大无外,其小无内% q# ^; q: v' |9 {! F3 V+ W

2 k5 v+ R' |2 E# f
作者: tysh670407    时间: 2016-6-2 09:55
发张JPG格式给定区间素数分布表。8 R, G$ f2 b, z) b2 V  t

给定区间素数分布规律表.jpg (237.85 KB, 下载次数: 338)

给定区间素数分布规律表.jpg


作者: tysh670407    时间: 2016-6-2 11:24
    1793年,德国数学家高斯由直觉看出了素数的分布的渐近分布定律:π(x)=∫dx/lnx,从2-x的积分。
  H2 S) F( t) V( {2 _高斯和勒让德还猜想极限(x→∞)limπ(x)/(x/lnx)=1,这个猜想就是著名的素数定理。1849年,俄国数学家切比雪夫发表博士论文,在假定π(x)/(x/lnx)极限存在的前提下证明了(x→∞)limπ(x)/(x/lnx)=1。' Q* b% p% c% Y

) T5 h9 L) ~7 r) m5 I    其实,前辈的猜想是对的,高斯是以1000为单位的区间,看出了素数的分布的渐近分布定律,而本人通过研究相邻两个奇数平方之间的区间素数分布,给出了区间素数分布的公式,实际与理论素数之比在1左右波动,如果把所有区间的比值平均,平均值趋向数值1。
9 T9 f2 r0 `0 r) H4 e. t
作者: tysh670407    时间: 2016-6-2 11:35
本帖最后由 tysh670407 于 2016-6-2 11:39 编辑
/ `3 A2 z6 r. q: n/ }) n
' ^+ i* W, h& t# l6 y# ?给定区间实际素数与理论素数之比的平均值趋向1。图中三个平均值平均后为1.000894746。( p+ b, j+ t0 H! F0 o+ ^% Q1 F
1 H, ]- d& S& ?! X$ C! {

给定区间实际素数与理论素数之比平均值1.JPG (530.38 KB, 下载次数: 358)

给定区间实际素数与理论素数之比平均值1.JPG


作者: tysh670407    时间: 2016-6-2 11:35
本帖最后由 tysh670407 于 2016-6-6 07:57 编辑 . R- ?8 L0 g- N* {  I# H% X4 X

" |; O2 _2 \6 k+ R5 _! w/ E        从1到无穷区间实际素数与理论素数之比(π(x)/((8x-8)/ln(2x-1)^2))的平均值等于1。: ]7 `% g* ~: ?3 t6 b9 {& b& `
  M4 B6 t" k' {
推论3 在区间((2x-3^2,(2x-1^2]内,实际素数个数总是在理论素数个数左右波动,即它们的比值在1左右波动,当x取有限数值时,所有区间实际素数与理论素数之比(π(x)/(8x-8/ln(2x-1)^2))的平均值趋向1。当x取无穷大时,无穷区间实际素数与理论素数之比(π(x)/(8x-8/ln(2x-1)^2))的平均值等于1。即x∞时,
{π(1)/ [8×1-8/ln(2×1-1)^2]+ π(2)/ [8×2-8/ln(2×2-1)^2]+
π(3)/[8×3-8/ln(2×3-1)^2]+…+π(x)/ [8x-8/ln(2x-1)^2]}/1+2+3+…+x=1
( c/ z5 S- C& G2 V: Y
9 m) ^9 N8 J9 q$ \

7 h6 q' G2 o% I4 Y) Z3 |
作者: 店铺买饺子    时间: 2016-11-16 17:09
真厉害 好好学习
0 R& u0 e8 q, j6 \. T
作者: tysh670407    时间: 2018-5-21 14:53
好久没有来论坛了,谢谢各位同仁前辈的关爱,如果大家觉得对,请应用到自己的论文中。# @; L$ g3 `. u

作者: 551987369    时间: 2018-6-6 17:50
无用的,你在一百万以内做的很精确,可是跑的十亿,一百亿,你的这些规律将被无情的砸烂,误差大到不能接受.我在10^10做了几个比较精确的近似公式,到了10^20,白做了,没有半点用处.
! q1 @+ ?' n. z* v6 E4 M我们不用 瞎忙了,现成的黎曼公式无法超越.虽然大数仍需修补,但目前没有精度比它更高的公式了.
, C6 o9 D) q8 X9 N1 C, U% A看你的文章,好像是民科,没有官方语言符号.虽然以前我也是民科,但我现在转而研究消化前人的成果,不说证明了什么,也不说发现了什么.静静的读书才是好事.
3 Z; L4 \. e0 \3 p1 h
作者: tysh670407    时间: 2018-6-12 14:53
551987369 发表于 2018-6-6 17:50
, |, V: q2 f0 a* f+ C: l无用的,你在一百万以内做的很精确,可是跑的十亿,一百亿,你的这些规律将被无情的砸烂,误差大到不能接受.我在 ...

+ v  z: y" ]. W( G4 [& ~' D* l感谢前辈的关注,黎曼也只是对高斯的猜想,作了一些修改,如果没有继承和修订,科学就不会发展,如果你能对学生的这些研究成果,做一些研究,你会有所收获的,高斯以1000为间隔区间,而本人只是对区间扩大为相邻两个奇数的平方之间,对高斯素数定理作了一些修订而已,科学是不需要别人承认的,只需要把你的成果告诉大家就行了。
2 A5 K/ o5 |* i" F# v
作者: tysh670407    时间: 2018-6-19 15:34
551987369 发表于 2018-6-6 17:50
3 X( f# U  i+ r- V9 S9 u无用的,你在一百万以内做的很精确,可是跑的十亿,一百亿,你的这些规律将被无情的砸烂,误差大到不能接受.我在 ...

0 p5 |$ B1 W1 Y- L# }5 z: J0 ~20021、30031、32741区间实际和理论素数对比表,最大数42.9亿。10个区间累计误差小于0.3%。3 s( N0 {+ G6 B" [+ k' k
20021、30031、32741区间实际与理论素数对比图.JPG
  t+ d* {0 ?; e2 Y0 H: C! w
( j5 {6 L! N* B! I+ u$ q9 A, O
作者: tysh670407    时间: 2018-6-22 10:15
本帖最后由 tysh670407 于 2018-6-22 10:18 编辑 9 S* G5 W3 ~2 e) Y

3 c1 l7 W& L3 V1 Z7 S1 e0 v1 d: [自己用EXCEL和其他软件做的一个统计表,能力有限。作了1000间隔区间和10000区间间隔的统计表。手头软件只能计算最大数为42.9亿的素数。2 ]0 {# S- H0 n6 D6 U+ u" o
2-100区间实际素数与理论素数总数对比图.jpg : i$ z' \! m4 Z1 Z
论相邻两个奇数平方之间的素数个数001.jpg
8 [: f  i; p# W* q 论相邻两个奇数平方之间的素数个数002.jpg - N2 |$ {7 |4 N! _. y5 w
论相邻两个奇数平方之间的素数个数003.jpg - [. e9 T6 C+ f- ?/ F
论相邻两个奇数平方之间的素数个数004.jpg , R% b0 C* d+ Q8 U+ ^8 Q" ]
论相邻两个奇数平方之间的素数个数005.jpg
' M6 `- d! b9 x) V/ {& r5 x: Z从表上可以看出,数越大,误差越小,比值越接近1。
6 _  P' Z; X" Q9 f
' E( E3 X9 j  |6 {7 W
作者: 数学1+1    时间: 2018-6-25 10:29
本帖最后由 数学1+1 于 2018-6-25 10:48 编辑
1 G' W4 W* |) `4 O0 m8 h; A: x. F* J1 r3 B* T, s
定理:命 A≥0 ,M≥3 .记在 A 与 A+M 之间的素数个数为 π(A;M) .则+ u/ m& s) H! I- w0 h. P! F* t
             π(A;M)≤(2M/log M)[1+o(log logM/log M)] .
; @% x, q& `- \; k  L; U% Y+ X这里与o有关的常数与 A 及 M 无关.7 b5 C, C( f6 z
问题 :求
+ s9 a7 D& h. S' V6 ]6 N0 X9 o             π(A;M)≥?: n1 k4 d% C0 N4 b

7 p% c1 d# ?! x- u
, o1 B+ [9 ~/ U
作者: tysh670407    时间: 2018-6-29 10:07
本帖最后由 tysh670407 于 2018-6-29 10:09 编辑 / |# X( v/ Y1 M3 G. o5 y
. v7 ^% t4 @' x1 c! i1 [/ F
做了一张相邻两个奇数平方之间的素数个数与素数定理误差统计计算表,供大家研究。
% J. e3 S& N9 m7 t7 ^( D9 f  }3 m9 D( S; W4 I
相邻两个奇数平方之间的素数个数与素数定理误差比较表.JPG
# t, G/ S: s' a" }' _8 Q  
1 D, i9 k- _0 A: X0 {  C+ k' Y 相邻两个奇数平方之间的素数个数与素数定理误差比较表001.JPG % j1 ^! H/ `! I5 R

$ a" N' S. q2 E. s$ K9 p& v, |7 L8 ]$ v; F5 v: U( z; z% z, z6 o6 i3 u
相邻两个奇数平方之间的素数个数与素数定理误差比较表002.JPG
4 B6 S0 A% c' m# M$ e7 _4 A5 G
  |) H7 ?4 v# f0 H* k/ D 相邻两个奇数平方之间的素数个数与素数定理误差比较表003.JPG 8 Y6 J& `8 \: |& \$ }# Z
  8 Z7 W9 w$ k! R/ L$ n" r
相邻两个奇数平方之间的素数个数与素数定理误差比较表004.JPG
% N2 B9 i* c9 Z$ ?' Z, \! ~# x  \7 ]: ^- B! X6 J
从表中可以看出,素数定理的误差远远大于相邻两个奇数平方之间的素数个数的误差,误差率总体呈下降趋势,但依然比相邻两个奇数平方之间的素数个数的误差率大几千倍。$ G, S: B+ c" z4 `5 \+ w: \3 S6 M: N
4.98/0.0015=3320) `0 ~/ W) I* {0 I; B2 Q# a2 V

- m" _# B8 T9 y, s! d$ F. X# p6 K3 t' H  c9 j" _  B





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5