《数论探秘》电子版

ysr

900内有3组蔡氏素数:
641,29999999999999429,29999999999999431,4n+2=30000000000000070
653,29999999999999429,29999999999999431,4n+2=30000000000000082
661,29999999999999429,29999999999999431,4n+2=30000000000000090
用时188.736秒

ysr

4000内有2组蔡氏素数:
3001,49999999999997081,49999999999997083,4n+2=50000000000000082
3917,49999999999996169,49999999999996171,4n+2=50000000000000086
用时567.269000000001秒

ysr

4000内有2组蔡氏素数:
3001,49999999999997081,49999999999997083,4n+2=50000000000000082
3917,49999999999996169,49999999999996171,4n+2=50000000000000086
用时136.774999999997秒，程序稍改速度变快。

ysr

8000内有5组蔡氏素数:
6701,5070602400912917605986812814869,5070602400912917605986812814871,4n+2=5070602400912917605986812821570
6709,5070602400912917605986812814869,5070602400912917605986812814871,4n+2=5070602400912917605986812821578
7321,5070602400912917605986812814257,5070602400912917605986812814259,4n+2=5070602400912917605986812821578
7901,5070602400912917605986812813681,5070602400912917605986812813683,4n+2=5070602400912917605986812821582
7333,5070602400912917605986812814257,5070602400912917605986812814259,4n+2=5070602400912917605986812821590
用时897.063999999997秒

程序算出来的31位的孪生素数，此法还是很有效的

ysr

Private Sub Command1_Click()
Dim a, B, q, m
Dim t As Double
t = Timer
a2 = Val(Text1)
m1 = Trim(Text4)
m2 = Trim(Text5)
m3 = MPC1(MbC(Trim(m2), 4), 66)
a1 = a2
If Right(a2, 1) Mod 2 = 0 Then
a2 = a2 + 1
Else
a2 = a2
End If
q = Val(Text2)
m = MPC1(MbC(Trim(m1), 4), 66)
Do While MBJC(Trim(m), Trim(m3)) <= 0
p1 = a2

Do While p1 <= q And p1 <= m

Do While InStr(fenjieyinzi0(Trim(p1)), "*") > 0
p1 = p1 + 2
Loop
a1 = 1
B1 = Sqr(Val(p1 - a1 ^ 2))
Do While InStr(B1, ".") <> 0 And p1 > Val(a1 ^ 2)
a1 = a1 + 1
B1 = Sqr(Abs(Val(p1 - a1 ^ 2)))
Loop
If InStr(B1, ".") = 0 Then
p2 = MPC(Trim(m), Trim(p1))



p3 = MPC1(Trim(p2), 2)


a = fenjieyinzi(Val(p1))
B = fenjieyinzi0(Trim(p2))
c = fenjieyinzi0(Trim(p3))


If InStr(a, "*") = 0 And InStr(B, "*") = 0 And InStr(c, "*") = 0 And InStr(B11, ".") = 0 And InStr(D12, ".") = 0 Then
s = s + 1
Print p1, p2, p3, p4
Text3 = Text3 & CStr(p1) & "," & CStr(p2) & "," & p3 & ",4n+2=" & m & vbCrLf

Else
s = s
End If
Else
p1 = p1
End If


p1 = Val(p1 + 2)


Loop
m = MPC1(Trim(m), 4)
Loop
Combo1 = q & "内有" & s & "组蔡氏素数:" & vbCrLf & Text3 & "用时" & Timer - t & "秒"

End Sub

ysr

146376311462000611635994205783与146376311462000611635994215783之间有8对差大于200的素数对:
146376311462000611635994207601/210
146376311462000611635994207861/210
146376311462000611635994209583/342
146376311462000611635994210379/216
146376311462000611635994211317/210
146376311462000611635994211569/222
146376311462000611635994212937/254
146376311462000611635994214681/202
用时2216.045秒

ysr

36,594,077,865,500,344,204,701,972,484-13830968^2-13830968-1=36,594,077,865,500,152,909,012,324,491,
36,594,077,865,500,152,909,012,324,491*4+3=146,376,311,462,000,611,636,049,297,967

ysr

1014000内有组蔡氏素数:
用时1674.714秒

ysr

146376311462000611635994214681与146376311462000611635994234681之间有12对差大于200的素数对:
146376311462000611635994217353/224
146376311462000611635994218611/342
146376311462000611635994220467/354
146376311462000611635994221811/248
146376311462000611635994225669/272
146376311462000611635994227221/270
146376311462000611635994227431/210
146376311462000611635994229357/204
146376311462000611635994229953/284
146376311462000611635994230313/222
146376311462000611635994233137/336
146376311462000611635994234471/222
用时4537.913秒

101004000内有组蔡氏素数:
用时1599.78秒