|

楼主 |
发表于 2025-7-20 06:30
|
显示全部楼层
本帖最后由 ysr 于 2025-7-19 22:32 编辑
m=a(a-1),p=m-b,当a=2^100,b=12601~13801时,油:
p=1606938044258990275541962092339894951921974764381296132082631 p1=1606938044258990275541962092339894951921974764381296132082633
用时-76808.27秒,孪生素数对一个
代码如下:
Private Sub Command1_Click() '快速找孪生素数
Dim A, B
a1 = Text1: B1 = Text2
ts = Timer
A2 = ksm2(Val(a1))
s = 12601
Do While Val(s) <= B1
A = Trim(A2)
p = MPC(MbC(Trim(A), Trim(A)), Trim(A))
p = MPC(Trim(p), Trim(s))
c = fenjieyinzi0(Trim(p))
If InStr(c, "*") = 0 Then
p1 = MPC1(Trim(p), 2)
C1 = fenjieyinzi0(Trim(p1))
If InStr(C1, "*") = 0 Then
s2 = s2 + 1
s1 = s1 & "p=" & p & " p1=" & p1 & vbCrLf
End If
End If
s = s + 2
Loop
If Val(s2) > 0 Then
Text3 = s1 & "用时" & Timer - ts & "秒"
Else
Text3 = "wu jie " & "用时" & Timer - ts & "秒"
End If
End Sub
Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
End Sub |
|