数学中国

用户名  找回密码
 注册
帖子
热搜: 活动 交友 discuz
楼主: ysr

朋友,请看一下重要启示!

[复制链接]
发表于 2023-2-10 10:00 | 显示全部楼层
10亿开始的10000个数字内的孪生素数:
  1. from sympy import primerange, isprime

  2. twins=[n for n in primerange(10**9, 10**9+10**4) if isprime(n + 2)]

  3. print(len(twins),twins)
复制代码



31 [1000000007, 1000000409, 1000000931, 1000001447, 1000001537, 1000001789, 1000001801, 1000002359, 1000002821, 1000003271, 1000003469, 1000003577, 1000003889, 1000004249, 1000004567, 1000004609, 1000004867, 1000004891, 1000005449, 1000005971, 1000006037, 1000006127, 1000006457, 1000006661, 1000007651, 1000007927, 1000008257, 1000008311, 1000009277, 1000009529, 1000009559]

点评

cz1
厉害  发表于 2023-2-10 11:54
ysr
我目前找到的最大的是106位的  发表于 2023-2-10 10:15
ysr
好的,我目前找到的最大的是106为的  发表于 2023-2-10 10:14
回复 支持 反对

使用道具 举报

发表于 2023-2-10 10:03 | 显示全部楼层
本帖最后由 Treenewbee 于 2023-2-10 10:27 编辑

Mathematica:

  1. Select[Range[10^9, 10^9 + 10^4], PrimeQ[#] && PrimeQ[# + 2] &]
复制代码

点评

ysr
这是python代码?我无法用,复制后无法工作,显示有错可能也不懂显示内容  发表于 2023-2-10 10:14
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-10 10:05 | 显示全部楼层
Treenewbee 发表于 2023-2-10 01:27
你需要算什么,我可以帮你优化下

Private Function zzxc(sa As String, sb As String) As String
'辗转相除法
Dim a, B, c, d, r
  a = Trim(sa)
  B = Trim(sb)
  If Len(a) < 10 And Len(B) < 10 Then
  
  If Val(a) > Val(B) Then
     c = a
     d = B
  Else
     c = B
     d = a
  End If
Do Until Val(c) Mod Val(d) = 0
     r = c Mod d
     c = d
     d = r
  Loop
  
  Else
  
  If MBJC(Trim(a), Trim(B)) >= 1 Then
  c = a
  d = B
  Else
  c = B
  d = a
  End If
  Do Until zhengchuqyushu(MCC1(Trim(c), Trim(d))) = 0
  r = zhengchuqyushu(MCC1(Trim(c), Trim(d)))
  c = d
  d = r
  Loop
  End If

  
  zzxc = d
  
End Function

Private Function qniyuan(sa As String, sb As String) As String
'求乘法的逆元
Dim n, p, a, B, c, d, r
  n = Trim(sa)
  p = Trim(sb)
  a = 1
  B = 0
  c = 0
  d = 1
  If Len(n) < 10 And Len(p) < 10 Then
  
  If Val(n) > Val(p) Then
     m = n
     q = p
     s1 = 1
  Else
     m = p
     q = n
     s1 = 0
  End If
Do Until Val(m) Mod Val(q) = 0
    s = m \ q
     r = m Mod q
     s1 = s1 + 1
     If s1 Mod 2 = 1 Then
     a = a
     B = a * s + B
     c = c
     d = c * s + d
     Else
     B = B
     a = a + B * s
     d = d
     c = c + d * s
     End If
     m = q
     q = r
  Loop
  If Val(a + B * m) = p Then
  B = B
  a = a + B * (m - 1)
  d = d
  c = c + d * (m - 1)
  Else
  If Val(B + a * m) = p Then
  a = a
  B = B + a * m
  c = c
  d = d + c * m
  Else
  B = B
  a = a + B * (m - 1)
  d = d
  c = c + d * (m - 1)
  End If
  End If
x = (a + B) Mod p
  Y = (c + d) Mod n
  
  
  Else
  
  If MBJC(Trim(n), Trim(p)) >= 1 Then
  m = n
  q = p
  s1 = 1
  Else
  m = p
  q = n
  s1 = 0
  End If
  Do Until zhengchuqyushu(MCC1(Trim(m), Trim(q))) = 0
  s = zhengchuqy(MCC1(Trim(m), Trim(q)))
  r = zhengchuqyushu(MCC1(Trim(m), Trim(q)))
  s1 = s1 + 1
  If s1 Mod 2 = 1 Then
  a = a
  B = MPC1(MbC(Trim(a), Trim(s)), Trim(B))
  c = c
  d = MPC1(MbC(Trim(c), Trim(s)), Trim(d))
  Else
  B = B
  a = MPC1(Trim(a), MbC(Trim(B), Trim(s)))
  d = d
  c = MPC1(Trim(c), MbC(Trim(d), Trim(s)))
  End If
  
  m = q
  q = r
  Loop
  
  If MPC1(Trim(a), MbC(Trim(B), Trim(m))) = p Then
  B = B
  a = MPC1(Trim(a), MbC(Trim(B), MPC(Trim(m), 1)))
  d = d
  c = MPC1(Trim(c), MbC(Trim(d), MPC(Trim(m), 1)))
  Else
  If MPC1(Trim(B), MbC(Trim(a), Trim(m))) = p Then
  a = a
  B = MPC1(Trim(B), MbC(Trim(a), Trim(m)))
  c = c
  d = MPC1(Trim(d), MbC(Trim(c), Trim(m)))
  Else
  B = B
  a = MPC1(Trim(a), MbC(Trim(B), MPC(Trim(m), 1)))
  d = d
  c = MPC1(Trim(c), MbC(Trim(d), MPC(Trim(m), 1)))
  End If
  End If
Do While Left(a, 1) = "0"
    a = Mid(a, 2)
Loop
  
  End If
  
  qniyuan = a
End Function

Private Function qksmimo(sa As String, sb As String, sc As String) As String
’快速幂模
Dim c, e, n, d
c = Trim(sa)
e = Trim(sb)
n = Trim(sc)
d = 1
If Len(c) < 5 And Len(e) < 5 And Len(n) < 5 Then
c = Val(c): n = Val(n)
Do While e > 0
If Right(e, 1) Mod 2 = 0 Then
c = c * c Mod n
e = e / 2

Else
d = d * c Mod n
e = e - 1
End If
Loop
Else
c = c
Do While MBJC(Trim(e), 1) >= 0
If Right(e, 1) Mod 2 = 0 Then
c = zhengchuqyushu(MCC1(MbC(Trim(c), Trim(c)), Trim(n)))
e = zhengchuqy(MCC1(Trim(e), 2))
Else
d = zhengchuqyushu(MCC1(MbC(Trim(c), Trim(d)), Trim(n)))
e = MPC(Trim(e), 1)
End If
Loop
End If

qksmimo = d
End Function

Private Function fenjieyinzi(sa As String) As String
Dim a, B
Dim x As String
x = sa
B1 = Sqr(Val(x)) / 2
If InStr(B1, ".") = 0 Then
B = B1
Else
B = Left(B1, InStr(B1, ".") - 1)
End If
If x = 3 Or x = 2 Then
a = True
Else
If Right(x, 1) Mod 2 = 0 Then
a = False
Else
For I = 3 To 2 * B + 1 Step 2
b2 = x / I
If InStr(b2, ".") = 0 Then
a = False
Exit For

Else: a = True

End If
Next
End If
End If
If a = True Then
fenjieyinzi = "这是个素数"
Else
fenjieyinzi = "2*2"
End If


End Function
Private Function fenjieyinzi0(sa As String) As String
'判断素数
Dim a, n
n = Trim(sa)
If Len(n) < 11 Then
fenjieyinzi0 = fenjieyinzi(Trim(n))
Else
n1 = MPC(Trim(n), 1)
a = 123
'a为明文
a1 = zzxc(Trim(n), Trim(a))
If Val(a1) > 1 Then
fenjieyinzi0 = a1 & "*"
Else
c = 999
'c为公约
Do While zzxc(Trim(n1), Trim(c)) > 1
c = Val(c - 1)
Loop
d = qniyuan(Trim(c), Trim(n1))
'd为逆元为私钥
a2 = qksmimo(Trim(a), Trim(c), Trim(n))
'a2为密文
a3 = qksmimo(Trim(a2), Trim(d), Trim(n))
If MBJC(Trim(a3), Trim(a)) = 0 Then
fenjieyinzi0 = "这是素数有"
Else
fenjieyinzi0 = "2*2"
End If
End If
End If


End Function

回复 支持 反对

使用道具 举报

发表于 2023-2-10 10:26 | 显示全部楼层
  1. from sympy import isprime,nextprime
  2. p=10**123
  3. while p<10**123+10**5:
  4.     p=nextprime(p)
  5.     if isprime(p+2):
  6.         print(p)
复制代码


1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000074011
1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000088897
1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000091357

点评

cz1
赞  发表于 2023-2-10 11:55
ysr
您的程序速度快我比不了,复制了代码也用不了  发表于 2023-2-10 11:00
ysr
能找到差大于1550的相邻素数吗?  发表于 2023-2-10 10:59
ysr
能算几十万位的素数吗?  发表于 2023-2-10 10:58
ysr
好的,谢谢!  发表于 2023-2-10 10:47
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-10 10:58 | 显示全部楼层
Treenewbee 发表于 2023-2-10 02:26
100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 ...

Private Sub Command1_Click()
'求差大于200的相邻素数对的程序
Dim a, B, q
Dim t As Double
t = Timer
a = Trim(Text1)
q = Trim(Text2)
a1 = a
If Right(a, 1) Mod 2 = 0 Then
a = MPC1(Trim(a), 1)
Else
a = a
End If
a2 = a
Do While MBJC(Trim(a2), Trim(q)) <= 0

Do While InStr(fenjieyinzi0(Trim(a2)), "*") > 0
a2 = MPC1(Trim(a2), 2)
Loop
s2 = s2 + 1
If s2 = 1 Then
a3 = 0
Else
a3 = MPC(Trim(a2), Trim(a4))
End If
If Val(a3) > 200 Then
Text3 = Text3 & a2 & "/" & a3 & vbCrLf
s1 = s1 + 1
Else
s1 = s1
End If
a4 = a2
a2 = MPC1(Trim(a2), 2)
Loop


Combo1 = a1 & "与" & q & "之间有" & Val(s1) & "对差大于200的素数对:" & vbCrLf & Text3 & "用时" & Timer - t & "秒"






End Sub
回复 支持 反对

使用道具 举报

 楼主| 发表于 2023-2-10 11:09 | 显示全部楼层
ysr 发表于 2023-2-10 02:58
Private Sub Command1_Click()
'求差大于200的相邻素数对的程序
Dim a, B, q


刚才有领导来检查工作,我的运行的程序都没有了,工作人员给关闭了可能是运行中的程序关闭就消失了,那个找大孪生素数的程序运行了两天了也没有结束呢,没有出来结果。
回复 支持 反对

使用道具 举报

发表于 2023-2-10 11:21 | 显示全部楼层
能找到差大于1550的相邻素数吗?

  1. from sympy import nextprime
  2. n=876129783640337274493518
  3. p1=nextprime(n)
  4. p2=nextprime(p1)
  5. print(p1,p2,p2-p1)
复制代码


876129783640337274493519 876129783640337274495071 1552

点评

cz1
厉害了  发表于 2023-2-10 11:55
ysr
我的程序速度太慢  发表于 2023-2-10 11:34
ysr
876129783640337274493519~876129783640337274495071之间的素数有2个:(用时92.94531秒) 876129783640337274493519 876129783640337274495071  发表于 2023-2-10 11:33
ysr
好的,这个本论坛好像有人发过了,还能找到更大的吗?  发表于 2023-2-10 11:23
回复 支持 反对

使用道具 举报

发表于 2023-2-10 11:51 | 显示全部楼层
[99000033712000137027043172639, 99000033712000137027043174639, 2000]

[8731133256746730441012679218808000172688990926807, 8731133256746730441012679218808000172688990929807, 3000]

[15848645757237170332266306025344629951337383026446270255661, 15848645757237170332266306025344629951337383026446270259661, 4000]

[2787894262996051954686348019177002923466704960301632922923673701702143609638833244607, 2787894262996051954686348019177002923466704960301632922923673701702143609638833250607, 6000]

点评

ysr
99000033712000137027043172639~99000033712000137027043174639之间的素数有2个:(用时183.3906秒) 99000033712000137027043172639 99000033712000137027043174639  发表于 2023-2-10 12:18
ysr
哇,厉害!大大点赞!  发表于 2023-2-10 12:11
cz1
第一高手!  发表于 2023-2-10 11:56
回复 支持 反对

使用道具 举报

发表于 2023-2-10 11:54 | 显示全部楼层
[420517575233115220178844026853036232825701025431059360625902934237070595672284391371291677111958273103698978716908840985606999571661690196653,
420517575233115220178844026853036232825701025431059360625902934237070595672284391371291677111958273103698978716908840985606999571661690205541,
8888]

点评

ysr
哇,厉害!给你大大点赞!  发表于 2023-2-10 12:11
回复 支持 反对

使用道具 举报

发表于 2023-2-10 12:03 | 显示全部楼层
[951895138203100127916737476292130673069912407655034360374674907451915559257242306072366345648136065623290646203791991319822831899374526583698117148686038334669411410010108480912448126277269022405660338562895665755407319195436581538533626564258408530085796520848596096952779594754601690394294192314860034484791195037659285312486949839743585798783298269791570495228796729532766418251125169591057715567160285742111507032666142917306422252659062561446917562741167573982353072096374046895462399965748510717288597970677160434518098590344549560229299010020966982552562121978448985746322982995564229991357326870161501485419686017788369607680703902961412762266210375967055394214610600921731442394506874692106657990429288773602418246717943918578490383639866735784356224775156089891189697758271184377925985683392432675230395145724788731541637333256283403458088191164140479127711331365043041950124583421728181486456945698297609431545452542355815731880906966651869369797291552253303908490143908589618078919246634615769319350139460156097994614630647951612485653646067072704347155478969833165042837484328597311825348869709250620568070429431799844813,
951895138203100127916737476292130673069912407655034360374674907451915559257242306072366345648136065623290646203791991319822831899374526583698117148686038334669411410010108480912448126277269022405660338562895665755407319195436581538533626564258408530085796520848596096952779594754601690394294192314860034484791195037659285312486949839743585798783298269791570495228796729532766418251125169591057715567160285742111507032666142917306422252659062561446917562741167573982353072096374046895462399965748510717288597970677160434518098590344549560229299010020966982552562121978448985746322982995564229991357326870161501485419686017788369607680703902961412762266210375967055394214610600921731442394506874692106657990429288773602418246717943918578490383639866735784356224775156089891189697758271184377925985683392432675230395145724788731541637333256283403458088191164140479127711331365043041950124583421728181486456945698297609431545452542355815731880906966651869369797291552253303908490143908589618078919246634615769319350139460156097994614630647951612485653646067072704347155478969833165042837484328597311825348869709250620568070429431799911479,
66666]
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

LaTEX棰勮杈撳叆銆€鏁欑▼銆€绗﹀彿搴�銆€鍔犺鍐呮爣绛�銆€鍔犺闂存爣绛�銆€
瀵瑰簲鐨� LaTEX 鏁堟灉锛�

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2025-7-25 18:01 , Processed in 0.570677 second(s), 22 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表
\frac{\square}{\square}\sqrt{\square}\square_{\baguet}^{\baguet}\overarc{\square}\ \dot{\baguet}\left(\square\right)\binom{\square}{\square}\begin{cases}\square\\\square\end{cases}\ \begin{bmatrix}\square&\square\\\square&\square\end{bmatrix}\to\Rightarrow\mapsto\alpha\ \theta\ \pi\times\div\pm\because\angle\ \infty
\frac{\square}{\square}\sqrt{\square}\sqrt[\baguet]{\square}\square_{\baguet}\square^{\baguet}\square_{\baguet}^{\baguet}\sum_{\baguet}^{\baguet}\prod_{\baguet}^{\baguet}\coprod_{\baguet}^{\baguet}\int_{\baguet}^{\baguet}\lim_{\baguet}\lim_{\baguet}^{\baguet}\bigcup_{\baguet}^{\baguet}\bigcap_{\baguet}^{\baguet}\bigwedge_{\baguet}^{\baguet}\bigvee_{\baguet}^{\baguet}
\underline{\square}\overline{\square}\overrightarrow{\square}\overleftarrow{\square}\overleftrightarrow{\square}\underrightarrow{\square}\underleftarrow{\square}\underleftrightarrow{\square}\dot{\baguet}\hat{\baguet}\vec{\baguet}\tilde{\baguet}
\left(\square\right)\left[\square\right]\left\{\square\right\}\left|\square\right|\left\langle\square\right\rangle\left\lVert\square\right\rVert\left\lfloor\square\right\rfloor\left\lceil\square\right\rceil\binom{\square}{\square}\boxed{\square}
\begin{cases}\square\\\square\end{cases}\begin{matrix}\square&\square\\\square&\square\end{matrix}\begin{pmatrix}\square&\square\\\square&\square\end{pmatrix}\begin{bmatrix}\square&\square\\\square&\square\end{bmatrix}\begin{Bmatrix}\square&\square\\\square&\square\end{Bmatrix}\begin{vmatrix}\square&\square\\\square&\square\end{vmatrix}\begin{Vmatrix}\square&\square\\\square&\square\end{Vmatrix}\begin{array}{l|l}\square&\square\\\hline\square&\square\end{array}
\to\gets\leftrightarrow\nearrow\searrow\downarrow\uparrow\updownarrow\swarrow\nwarrow\Leftarrow\Rightarrow\Leftrightarrow\rightharpoonup\rightharpoondown\impliedby\implies\Longleftrightarrow\leftharpoonup\leftharpoondown\longleftarrow\longrightarrow\longleftrightarrow\Uparrow\Downarrow\Updownarrow\hookleftarrow\hookrightarrow\mapsto
\alpha\beta\gamma\Gamma\delta\Delta\epsilon\varepsilon\zeta\eta\theta\Theta\iota\kappa\varkappa\lambda\Lambda\mu\nu\xi\Xi\pi\Pi\varpi\rho\varrho\sigma\Sigma\tau\upsilon\Upsilon\phi\Phi\varphi\chi\psi\Psi\omega\Omega\digamma\vartheta\varsigma\mathbb{C}\mathbb{H}\mathbb{N}\mathbb{P}\mathbb{Q}\mathbb{R}\mathbb{Z}\Re\Im\aleph\partial\nabla
\times\cdot\ast\div\pm\mp\circ\backslash\oplus\ominus\otimes\odot\bullet\varnothing\neq\equiv\not\equiv\sim\approx\simeq\cong\geq\leq\ll\gg\succ\prec\in\ni\cup\cap\subset\supset\not\subset\not\supset\notin\not\ni\subseteq\supseteq\nsubseteq\nsupseteq\sqsubset\sqsupset\sqsubseteq\sqsupseteq\sqcap\sqcup\wedge\vee\neg\forall\exists\nexists\uplus\bigsqcup\bigodot\bigotimes\bigoplus\biguplus\bigcap\bigcup\bigvee\bigwedge
\because\therefore\angle\parallel\perp\top\nparallel\measuredangle\sphericalangle\diamond\diamondsuit\doteq\propto\infty\bowtie\square\smile\frown\bigtriangledown\triangle\triangleleft\triangleright\bigcirc \wr\amalg\models\preceq\mid\nmid\vdash\dashv\nless\ngtr\ldots\cdots\vdots\ddots\surd\ell\flat\sharp\natural\wp\clubsuit\heartsuit\spadesuit\oint\lfloor\rfloor\lceil\rceil\lbrace\rbrace\lbrack\rbrack\vert\hbar\aleph\dagger\ddagger

MathQuill输入:

Latex代码输入: