《数论探秘》电子版

ysr

18361375335787470653与18361375335787510653之间有7对差大于200的素数对:
18361375335787483553/222
18361375335787487203/210
18361375335787488617/238
18361375335787494599/220
18361375335787497583/254
18361375335787499123/274
18361375335787505053/216
用时3249.008秒

18361375335787505053与18361375335787545053之间有6对差大于200的素数对:
18361375335787520819/208
18361375335787535629/218
18361375335787538501/210
18361375335787540417/246
18361375335787541051/210
18361375335787544369/250
用时3529.523秒

ysr

18361375335787544369与18361375335787614369之间有11对差大于200的素数对:
18361375335787550453/256
18361375335787559213/202
18361375335787562297/228
18361375335787566947/204
18361375335787568327/240
18361375335787570379/240
18361375335787573841/204
18361375335787574227/246
18361375335787591451/264
18361375335787600763/294
18361375335787601639/216
用时6470.105秒

ysr

36722750671574688447与36722750671574768447之间的素数开头差为1922760350154212639070最多有4个连续素数:
/36722750671574694761/1959483100825787333831/3882243450979999972901/5805003801134212611971
/36722750671574703677/1959483100825787342747/3882243450979999981817/5805003801134212620887
/36722750671574714737/1959483100825787353807/3882243450979999992877/5805003801134212631947
/36722750671574718241/1959483100825787357311/3882243450979999996381/5805003801134212635451
/36722750671574725519/1959483100825787364589/3882243450980000003659/5805003801134212642729
/36722750671574734489/1959483100825787373559/3882243450980000012629/5805003801134212651699
/36722750671574738539/1959483100825787377609/3882243450980000016679/5805003801134212655749
用时-59295.66秒

ysr

ysr 发表于 2022-10-19 22:59
36722750671574688447与36722750671574768447之间的素数开头差为1922760350154212639070最多有4个连续素数: ...

18361375335787601639与18361375335787681639之间有20对差大于200的素数对:
18361375335787621681/234
18361375335787629979/270
18361375335787630249/270
18361375335787631801/228
18361375335787632089/252
18361375335787636939/212
18361375335787640657/216
18361375335787644019/270
18361375335787645201/224
18361375335787649857/210
18361375335787660081/228
18361375335787666033/210
18361375335787666591/212
18361375335787668619/228
18361375335787669973/250
18361375335787670993/214
18361375335787672103/204
18361375335787673561/210
18361375335787677467/238
18361375335787678699/212
用时7220.544秒

ysr

'分解3535000647287053716310161223639238642910682574781777736223167的程序
Private Sub Command1_Click()
Dim A, B
A = Trim(Text1)
t = Timer
If Len(A) <= 10 Then
Text2 = fenjieyinzi(Trim(A))
Else
If Len(A) <= 13 Then
Text2 = fenjieyinzi0(Trim(A))
Else
B = zhengchuqy(MBBC(Trim(A)))
Print B
If InStr(MCC1(Trim(A), Trim(B)), "/") = 0 Then
Text2 = B & "*" & MCC1(Trim(A), Trim(B))
Else
s1 = 2
c = 2
Do While s1 <= 20 And InStr(MCC1(Trim(A), Trim(c)), "/") <> 0

c = zhengchuqy(MBBC(zhengchuqy(MCC(Trim(A), Trim(s1)))))
s1 = s1 + 2
Loop
B1 = zhengchuqy(MBBC(zhengchuqy(MCC(Trim(A), 4))))
Dim sa()
ReDim sa(1 To 10)
sa(1) = 13: sa(2) = 17: sa(3) = 19: sa(4) = 31: sa(5) = 61
sa(6) = 89: sa(7) = 107: sa(8) = 127: sa(9) = 521: sa(10) = 607
a1 = 1 & "/" & 2
sb = 1
s = 0
Do While Len(sb) <= Len(A) And InStr(a1, "/") <> 0 And s < 10
s = Val(s + 1)
sb = jsmssushu(Trim(sa(s)))
a1 = MCC1(Trim(A), Trim(sb))
Loop

If zhengchuqyushu(MCC(Trim(A), 4)) = 1 Then

C1 = zhengchuqy(MBBC(MPC(MbC(MPC1(Trim(B1), 1), MPC1(Trim(B1), 2)), zhengchuqy(MCC(Trim(A), 4)))))
Else

C1 = zhengchuqy(MBBC(MPC(MbC(MPC1(Trim(B1), 1), MPC1(Trim(B1), 1)), zhengchuqy(MCC(Trim(A), 4)))))
End If

If InStr(MCC1(Trim(A), Trim(c)), "/") = 0 Then
Text2 = c & "*" & MCC1(Trim(A), Trim(c))
Else
If InStr(MCC1(Trim(A), MPC1(MbC(MPC(MPC1(Trim(B1), 1), Trim(C1)), 2), 1)), "/") = 0 Then
Text2 = MPC1(MbC(MPC(MPC1(Trim(B1), 1), Trim(C1)), 2), 1) & "*" & MCC1(Trim(A), MPC1(MbC(MPC(MPC1(Trim(B1), 1), Trim(C1)), 2), 1))
Else
If InStr(a1, "/") = 0 Then
Text2 = sb & "*" & a1
Else
d = zhengchuqy(MCC(Trim(A), 4))
b3 = 1041
Do While b3 >= 1039 And b100 < 1
b4 = b3 * 853

n1 = zhengchuqy(MCC(Trim(d), Trim(b3)))
n2 = zhengchuqy(MBBC(Trim(n1)))
p1 = MPC1(MbC(Trim(n2), 2), 1)
r = zhengchuqyushu(MCC1(Trim(A), Trim(p1)))
q = zhengchuqy(MCC1(Trim(A), Trim(p1)))
q1 = MCC1(Trim(A), Trim(p1))
If MBJC(Trim(r), 0) = 0 Then
p2 = p1
Text2 = p2 & "*" & q
Else
r9 = p1
Print p1

p4 = r9


If zhengchuqyushu(MCC1(Trim(A), Trim(p4))) = 0 Then
Text2 = p4 & "*" & MCC1(Trim(A), Trim(p4))
Else
p5 = p1
Do While Len(p4) >= Len(p5) And zhengchuqyushu(MCC1(Trim(A), Trim(p4))) <> 0 And s12 <= 10
p4 = MPC(Trim(p4), 2)
s12 = s12 + 1
Loop
If zhengchuqyushu(MCC1(Trim(A), Trim(p4))) = 0 Then
Text2 = p4 & "*" & MCC1(Trim(A), Trim(p4))
Else
Dim x, y
x = Len(A) \ 5: y = Len(p4) \ 5
A = String(5 - Len(A) + 5 * x, "0") & A: p4 = String(5 - Len(p4) + 5 * y, "0") & p4
Print A
Dim aa() As String
ReDim aa(5 To 5 * x + 5)
Dim bb() As String
ReDim bb(5 To 5 * y + 5)
Dim i, j, js
For i = 5 To 5 * x + 5 Step 5
aa(i) = Mid(A, i - 4, 5)
Next
For j = 5 To 5 * y + 5 Step 5
bb(j) = Mid(p4, j - 4, 5)
Next
js = MCC1(qdqiandao0(aa(5) & aa(10) & aa(15)), qdqiandao0(bb(5) & bb(10)))
js17 = qdqiandao0(bb(5) & bb(10))
Print js17
s17 = Len(js17)
s19 = 42
Dim jsc()
ReDim jsc(1 To s19 + 1)
Do While s13 <= s19

s13 = s13 + 1
jsc(s13) = js17

If s13 = 33 Then
js6 = js17
Print js6
Else
js7 = js17

End If



js17 = Val(js17 - 1)

Loop
If y >= 2 Then
jss = jsc(33)
For i11 = 3 To y + 1
s19 = i11 * 5
jss = jss & bb(s19)
Next
Else
jss = js6
End If
jss4 = jss
Print r9, jss
If InStr(MCC1(qdqiandao0(Trim(A)), Trim(jss)), "/") = 0 Then
Text2 = jss & "*" & MCC1(qdqiandao0(Trim(A)), Trim(jss))
Else
p4 = r9
p = p4
q1 = MCC1(Trim(A), Trim(p))
q = zhengchuqy(Trim(q1))

B1 = zhengchuqy(MCC1(q & String(Len(p4), "0"), Trim(p)))
b6 = "100057991087214518562694508532"
b7 = "119124090279615294549011141875"
b8 = "99999999513536191490633492"
b12 = "828295": b13 = "830012": b14 = "831729"
b15 = "833446": b16 = "835163": b17 = "836880": b18 = "838597"
b19 = "840314": b20 = "842031": b21 = "843748"
b22 = "9823934683879298"
b4 = mbc2(Trim(B1), Trim(b6), 29)
b5 = mbc2(Trim(B1), Trim(b7), 29)
q5 = mbc2(Trim(p), Trim(b5), Len(p4))
q4 = mbc2(Trim(p), Trim(b4), Len(p4))
p6 = zhengchuqy(MCC1(Trim(A), Trim(q4)))
p7 = zhengchuqy(MCC1(Trim(A), Trim(q5)))
q1 = zhengchuqy(MCC1(Trim(A), Trim(p6)))
Do While s100 < 10 And b100 < 1
q2 = zhengchuqy(MCC1(Trim(A), Trim(p7)))
b9 = mbc2(zhengchuqy(MCC1(q1 & String(29, "0"), Trim(p6))), Trim(b6), 29)
b10 = mbc2(zhengchuqy(MCC1(q2 & String(6, "0"), Trim(p7))), Trim(b19), 6)

p8 = zhengchuqy(MCC1(Trim(A), mbc2(Trim(p6), Trim(b9), 29)))
p9 = zhengchuqy(MCC1(Trim(A), mbc2(Trim(p7), Trim(b10), 6)))
Print B1
Print p6
Print p7
Print p8
Print p9


A = qdqiandao0(Trim(A))
s21 = "9999999951353619149063" '22位的
s22 = "99999999513536191490633492" '26位的
s23 = "999999995135361914906334919" '27位的
s24 = "99956394676725327989068343" '26位的
s25 = "99993178279924191091877084380955816799674929190278759643399227469909582317774679409277779390361702881128085595674603314355069836029809524039751096000867814679064398130351092763713467281017813680758978657811271449078327588240759771011961998732872541429660632356218587832956076711272648205251375596651117006508796301486896750882609489302072424483319296843133847373293710768460052850309879" '386位的
jss = p9
q3 = MCC1(Trim(A), Trim(jss))
q4 = zhengchuqy(Trim(q3))
r13 = zhengchuqyushu(Trim(q3))
f7 = "295272111931977014044900 "
f4 = zhengchuqy(MCC1(q4 & String(Len(jss), "0"), Trim(jss)))
f2 = MbC(Trim(f4), Trim(s24))
f2 = Left(f2, Len(f2) - 26)
f6 = f2
b5 = MPC(Trim(q4), mbc2(Trim(f2), Trim(jss), Len(jss)))
b6 = MbC(Trim(r13), Trim(f2))
D1 = MPC(MbC(Trim(b5), Trim(b5)) & String(Len(jss), "0"), MbC(Trim(b6), 4))
D2 = zhengchuqy(MBBC(Trim(D1) & String(Len(jss), "0")))
f1 = zhengchuqy(MCC1(MPC1(Trim(b5) & String(Len(jss), "0"), Trim(D2)) & String(Len(jss) + 10, "0"), MbC(Trim(f2), 2)))
f5 = Len(f1) - Len(jss) - 10
f2 = Mid(f1, 1, f5)
f5 = MbC(Trim(f2), 98)
f5 = Left(f5, Len(f5) - 1)
x4 = mbc2(Trim(f2), Trim(f6), Len(jss))
Print f6
Print f2
Print x4
f3 = "283314068346029014014"
x3 = zhengchuqy(MCC1(MPC(MbC(Trim(q4), Trim(f2)), Trim(r13)), MPC1(Trim(jss), Trim(f2))))
jss = zhengchuqy(MCC1(Trim(A), MPC(Trim(q4), Trim(x3))))
'jss = MPC1(Trim(p9), Trim(f2))
jss1 = MPC1(Trim(jss4), Trim(f2))
jss2 = MPC(Trim(jss4), Trim(f2))
q5 = MPC(MbC(Trim(q4), Trim(q4)), Trim(A))
q6 = zhengchuqy(MBBC(Trim(q5)))
jss3 = MPC(Trim(q4), Trim(q6))
jss5 = zhengchuqy(MCC(Trim(jss4), 2))
jss3 = MPC1(Trim(jss3), Trim(jss5))
If Right(jss, 1) Mod 2 = 0 Then
jss = MPC(Trim(jss), 1)
Else
jss = jss
End If
Print jss
Print jss1
Print jss2
Print jss3
Print jss4
If Right(jss1, 1) Mod 2 = 0 Then
jss1 = MPC(Trim(jss1), 1)
Else
jss1 = jss1
End If
jss8 = jss
Do While InStr(MCC1(qdqiandao0(Trim(A)), Trim(jss)), "/") <> 0 And s15 <= 6000
s15 = s15 + 1
jss = MPC1(Trim(jss), 2)
Loop
If InStr(MCC1(Trim(A), Trim(jss)), "/") = 0 Then
Text2 = jss & "*" & MCC1(Trim(A), Trim(jss)) & "   中间数据X=" & x3
b100 = 1
Else
Do While InStr(MCC1(Trim(A), Trim(jss8)), "/") <> 0 And s16 <= 6500
s16 = s16 + 1
jss8 = MPC(Trim(jss8), 2)
Loop
If InStr(MCC1(Trim(A), Trim(jss8)), "/") = 0 Then
Text2 = jss8 & "*" & MCC1(Trim(A), Trim(jss8))
b100 = 1
Else
Text2 = "无法分解"
End If
End If
s100 = s100 + 1
p7 = MPC1(Trim(p7), 1 & String(Len(p7) - 1, "0"))
Loop



End If
End If
End If
End If
b3 = b3 - 2
Loop
End If
End If
End If
End If
End If
End If
Text3 = Timer - t & "秒"
End Sub

ysr

Private Function fenjieyinzi0(A As String) As String
x = 3
Do While InStr(MCC(Trim(A), 2), "/") = 0 And MBJC(Trim(A), 4) >= 0

  
If B = 0 Then
fenjieyinzi0 = fenjieyinzi0 & "2"
B = 1
Else
fenjieyinzi0 = fenjieyinzi0 & "*2"
End If
A = MCC(Trim(A), 2)
k = A
  
Loop

Do While MBJC(Trim(A), "0001") > 0
Do While MBJC(Trim(x), zhengchuqy(MBBC(Trim(A)))) <= 0
Do While InStr(MCC1(Trim(A), Trim(x)), "/") = 0 And MBJC(Trim(A), MbC(Trim(x), Trim(x))) >= 0
  
If B = 0 Then
fenjieyinzi0 = fenjieyinzi0 & x
B = 1
Else
fenjieyinzi0 = fenjieyinzi0 & "*" & x
End If
A = MCC1(Trim(A), Trim(x))
Loop
  
x = MPC1(Trim(x), 2)
Loop
  
k = A
A = 1
Loop
  
If B = 1 Then
fenjieyinzi0 = fenjieyinzi0 & "*" & k
Else
fenjieyinzi0 = "这是一个质数"
End If
End Function



Private Function fenjieyinzi(sa As String) As String
Dim x, A, B, k As String
A = Val(sa)

x = 3
If A <= 1 Or A > Int(A) Then
If A = 1 Then
fenjieyinzi = "它既不是质数，也不是合数"

Else
MsgBox "error"
End If
  
Else

Do While A / 2 = Int(A / 2) And A >= 4
  
If B = 0 Then
fenjieyinzi = fenjieyinzi & "2"
B = 1
Else
fenjieyinzi = fenjieyinzi & "*2"
End If
A = A / 2
k = A
  
Loop

Do While A > 1
Do While x <= Sqr(A)
Do While A / x = Int(A / x) And A >= x * x
  
If B = 0 Then
fenjieyinzi = fenjieyinzi & x
B = 1
Else
fenjieyinzi = fenjieyinzi & "*" & x
End If
A = A / x
Loop
  
x = x + 2
Loop
  
k = A
A = 1
Loop
  
If B = 1 Then
fenjieyinzi = fenjieyinzi & "*" & k
Else
fenjieyinzi = "这是一个质数"
End If
  
  
  
  

End If

End Function

ysr

18361375335787683699与18361375335787689699之间有1对差大于200的素数对:
18361375335787689061/210
用时3358.536秒

18361375335787689061与18361375335787699061之间有2对差大于200的素数对:
18361375335787692457/224
18361375335787698413/214
用时5665.412秒

ysr

18361375335787698413与18361375335787708413之间有4对差大于200的素数对:

18361375335787700599/210
18361375335787702249/210
18361375335787705757/216
18361375335787706149/380
用时5166.08700000001秒

ysr

Private Function qdqiandao0(sa As String) As String
  A = sa
  Do While Left(A, 1) = "0"
  A = Mid(A, 2)
  Loop
  If A = "" Then
  A = 0
  Else
  A = A
  End If
  qdqiandao0 = A
  End Function
  
   Private Function qdqd0(sa As String) As String
  A = sa
  Do While Left(A, 1) = "0"
  A = Mid(A, 2)
  Loop
  If A = "" Then
  A = 0
  Else
  A = A
  End If
  qdqd0 = A
  End Function

ysr

18361375335787726149与18361375335787786149之间有17对差大于200的素数对:
18361375335787726603/300
18361375335787732739/268
18361375335787736479/308
18361375335787742281/204
18361375335787747091/232
18361375335787748057/208
18361375335787754377/230
18361375335787756687/210
18361375335787759349/216
18361375335787765873/224
18361375335787766389/248
18361375335787768733/220
18361375335787770071/220
18361375335787770311/240
18361375335787775693/244
18361375335787782073/212
18361375335787785253/206
用时-79505.457秒