《数论探秘》电子版

ysr

146376311462000611636049297967与146376311462000611636049307967之间有9对差大于200的素数对:
146376311462000611636049298973/242
146376311462000611636049299753/206
146376311462000611636049300647/216
146376311462000611636049301223/242
146376311462000611636049303069/208
146376311462000611636049303963/244
146376311462000611636049304677/238
146376311462000611636049306749/282
146376311462000611636049307223/390
用时6933.66秒

ysr

146376311462000611636049307223与146376311462000611636049327223之间有11对差大于200的素数对:
146376311462000611636049308363/374
146376311462000611636049308661/234
146376311462000611636049309713/216
146376311462000611636049311261/270
146376311462000611636049311609/260
146376311462000611636049313449/220
146376311462000611636049313787/240
146376311462000611636049316523/522
146376311462000611636049319469/290
146376311462000611636049320169/370
146376311462000611636049321459/252
用时10141.19秒

ysr

146376311462000611636049340701与146376311462000611636049360701之间有12对差大于200的素数对:
146376311462000611636049341813/224
146376311462000611636049342017/204
146376311462000611636049343503/270
146376311462000611636049345509/242
146376311462000611636049350063/270
146376311462000611636049353427/298
146376311462000611636049354639/262
146376311462000611636049355161/234
146376311462000611636049357371/330
146376311462000611636049359103/220
146376311462000611636049359873/216
146376311462000611636049360741/414
用时8929.452秒

9308429内有组蔡氏素数:
用时7171.52999999999秒

9309829内有组蔡氏素数:
用时2347.603秒

9364348内有组蔡氏素数:
用时2932.699秒

9368348内有组蔡氏素数:
用时11745.685秒
9449848与9469848之间的素数打头有0组差为2和2293498615990071511610820895302086940796564989168281123737588839386922876088484808070018553110125686554624和2的4生素数对: (用时6126.031秒）

ysr

传个图片看看，解方程的:

如下结果就是按此法算的：输入1:  a=1,  b=2,  c=3,  d=4，e=5;  输出结果1:    x1,2=-1.6369949192+ -2.4170266611,  x3,4=0.6369949192+ -0.6438375554.

不知道对不对，此方程求A的时候，根号里面的值是负的，所以变成了绝对值，结果按实数代入下一步了，不知道这样做对否？

ysr

Private Sub Command1_Click() '求由2个孪生素数对组成的4生素数组的程序
Dim A, B
A = Val(Text1)
a1 = A
q = Val(Text2)
ts = Timer
m = Sqr(q)
t = Trim(Text4)
If Right(A, 1) Mod 2 = 0 Then
A = A + 1
Else
A = A
End If
s = 0
a2 = A
b2 = a2 + 2
Do While a2 <= m
Do While InStr(fenjieyinzi0(Trim(a2)), "*") > 0 Or InStr(fenjieyinzi0(Trim(b2)), "*") > 0
a2 = a2 + 2
b2 = a2 + 2
Loop
a12 = 1
B12 = Sqr(Val(a2 - a12 ^ 2))
Do While InStr(B12, ".") <> 0 And a2 > Val(a12 ^ 2)
a12 = a12 + 1
B12 = Sqr(Abs(Val(a2 - a12 ^ 2)))
Loop
If InStr(B12, ".") = 0 Then
B1 = a2
b2 = a2 + 2
b3 = MPC1(Trim(b2), Trim(t))
b4 = MPC1(Trim(b3), 2)

C1 = fenjieyinzi0(Val(B1))
C2 = fenjieyinzi0(Val(b2))
c3 = fenjieyinzi0(Trim(b3))
D1 = fenjieyinzi0(Trim(b4))

If InStr(C1, "*") = 0 And InStr(D1, "*") = 0 And InStr(C2, "*") = 0 And InStr(c3, "*") = 0 Then
s = s + 1
Print B1, b2, b3, b4
Text3 = Text3 & "/" & B1 & "/" & b2 & "/" & b3 & "/" & b4 & vbCrLf
Else
s = s
End If
Else
a2 = a2
End If

a2 = a2 + 2

Loop
a2 = a2
s1 = s
Do While a2 <= q
b2 = a2 + 2
Do While InStr(fenjieyinzi0(Trim(a2)), "*") > 0 Or InStr(fenjieyinzi0(Trim(b2)), "*") > 0
a2 = a2 + 2
b2 = a2 + 2
Loop
a12 = 1
B12 = Sqr(Val(a2 - a12 ^ 2))
Do While InStr(B12, ".") <> 0 And a2 > Val(a12 ^ 2)
a12 = a12 + 1
B12 = Sqr(Abs(Val(a2 - a12 ^ 2)))
Loop
If InStr(B12, ".") = 0 Then
B1 = a2
b2 = a2 + 2
b3 = MPC1(Trim(b2), Trim(t))
b4 = MPC1(Trim(b3), 2)

C1 = fenjieyinzi0(Val(B1))
C2 = fenjieyinzi0(Val(b2))
c3 = fenjieyinzi0(Trim(b3))
D1 = fenjieyinzi0(Trim(b4))

If InStr(C1, "*") = 0 And InStr(D1, "*") = 0 And InStr(C2, "*") = 0 And InStr(c3, "*") = 0 Then
s1 = s1 + 1
Print B1, b2, b3, b4
Text3 = Text3 & "/" & B1 & "/" & b2 & "/" & b3 & "/" & b4 & vbCrLf
Else
s1 = s1
End If
Else
a2 = a2
End If
a2 = a2 + 2

Loop
Combo1 = a1 & "与" & q & "之间的素数打头有" & s1 & "组差为2和" & t & "和2的4生素数对: (用时" & Timer - ts & "秒）" & vbCrLf & Text3

End Sub

ysr

146376311462000611636049380187与146376311462000611636049390187之间有8对差大于200的素数对:
146376311462000611636049385061/258
146376311462000611636049385341/210
146376311462000611636049385649/246
146376311462000611636049386481/250
146376311462000611636049388157/204
146376311462000611636049388407/250
146376311462000611636049389289/216
146376311462000611636049389939/258
用时2260.223秒

ysr

36,594,077,865,500,344,204,701,972,484-13830967^2-13830967-1=36,594,077,865,500,152,909,039,986,427

36,594,077,865,500,152,909,039,986,427*4+3=146,376,311,462,000,611,636,159,945,711

ysr

146376311462000611636159945711与146376311462000611636159955711之间有10对差大于200的素数对:
146376311462000611636159946269/206
146376311462000611636159946869/348
146376311462000611636159948903/242
146376311462000611636159951069/276
146376311462000611636159951771/204
146376311462000611636159952051/280
146376311462000611636159953091/270
146376311462000611636159953821/244
146376311462000611636159954259/306
146376311462000611636159955591/258
用时2231.89秒

ysr

146376311462000611636159965587与146376311462000611636159975587之间有11对差大于200的素数对:
146376311462000611636159966193/256
146376311462000611636159966417/224
146376311462000611636159966657/240
146376311462000611636159968539/246
146376311462000611636159969241/220
146376311462000611636159970293/246
146376311462000611636159970671/278
146376311462000611636159972789/390
146376311462000611636159973201/264
146376311462000611636159973651/234
146376311462000611636159974311/222
用时16866.331秒

ysr

36,594,077,865,500,344,204,701,972,484-13830967^2-1=36,594,077,865,500,152,909,053,817,394

36,594,077,865,500,152,909,053,817,394*4+3=146,376,311,462,000,611,636,215,269,579