数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
楼主: ysr

求数列的通项公式

[复制链接]
 楼主| 发表于 2026-4-13 07:37 | 显示全部楼层
4034918759与4047918759之间的素数开头最多有公差为6469693230的10个连续素数:
/4035549169/10505242399/16974935629/23444628859/29914322089/36384015319/42853708549/49323401779/55793095009/62262788239
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 15:04 | 显示全部楼层
ysr 发表于 2026-4-9 07:11
现在开始从4047918759到7985202243看看,能否找到那个62项的等差素数数列

单位停电了,这个算了好几天的程序白费劲了,来电后还要重新开始计算呢
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 15:17 | 显示全部楼层
有来电了,重新开始吧!不定啥时候就又会停电了,没法,多开几个程序,分段计算把!
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 17:44 | 显示全部楼层
Private Function fenjieyinzi(sa As String) As String
Dim x, a, b
x = sa
b = Int(Sqr(Val(x)) / 2)
If x = 3 Or x = 2 Then
a = True
Else
If x Mod 2 = 0 Then
a = False
Else

For I = 3 To 2 * b + 1 Step 2
If x Mod I = 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
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 18:05 | 显示全部楼层
本帖最后由 ysr 于 2026-4-13 13:42 编辑

Private Sub Command1_Click()
Dim a, b
a = Val(Text1)
a1 = a
Q = Val(Text2)
r = Val(Text4)
t = Timer
m = Sqr(Q)
If Right(a, 1) Mod 2 = 0 Then
a = a + 1
Else
a = a
End If

Do While InStr(fenjieyinzi(Val(a)), "*") <> 0
a = a + 2
Loop

Do While a <= Q
b = a
d = fenjieyinzi(Val(b))
s = 0
s3 = ""
Do While s >= 0 And InStr(d, "*") = 0
s4 = s + 1
If s4 Mod 20 = 0 Then
s3 = s3 & "/" & b & vbCrLf
Else
s3 = s3 & "/" & b
End If
s = s + 1
s2 = s * (s + 1)
b = b + r
d = fenjieyinzi(Val(b))
Loop
If s >= 10 Then
Text3 = Text3 & s3 & vbCrLf
End If
If s > 1 Then
s5 = s5 & "/" & s
Else
s5 = s5
End If
Print s5
a = a + 2
s2 = 0
Loop
s6 = Mid(s5, 2)
     s7 = max(Trim(s6))
Combo1 = a1 & "与" & Q & "之间的素数开头最多有公差为" & r & "的" & s7 & "个连续素数:用时" & Timer - t  & "秒 " & vbCrLf & Text3

     

End Sub

Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
Text4 = ""
Combo1 = ""
Form1.Cls
End Sub
Private Function max(ByVal lists As String) As String
Dim temp As Long
Dim a() As String
a = Split(lists, "/")
Dim b As Long

For b = 0 To UBound(a)
If temp < a(b) Then temp = a(b)
Next
max = temp
End Function

Private Function fenjieyinzi(sa As String) As String
Dim x, a, b
x =  Val(sa)
b = Int(Sqr(Val(x)) / 2)
If x = 3 Or x = 2 Then
a = True
Else
If Val(x / 2) - Int(Val(x / 2)) = 0 Then
a = False
Else

For i = 3 To 2 * b + 1 Step 2
If Val(x / i) - Int(Val(x / i)) = 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
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 19:07 | 显示全部楼层
分段计算:7985202243-4047918759=3,937,283,484
3,937,283,484/4=984,320,871
984,320,871+1=984,320,872
4047918759+984,320,872=5,032,239,631
5,032,239,631+984,320,872=6,016,560,503
6,016,560,503+984,320,872=7,000,881,375
7,000,881,375+984,320,872=7,985,202,247
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 19:20 | 显示全部楼层
4049918759与4050018759之间的素数开头最多有公差为6469693230的8个连续素数:用时1829.266秒
合0.5小时
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 20:30 | 显示全部楼层
本帖最后由 ysr 于 2026-4-13 12:33 编辑
ysr 发表于 2026-4-13 11:20
4049918759与4050018759之间的素数开头最多有公差为6469693230的8个连续素数:用时1829.266秒
合0.5小时


就这个速度,从4047918759算到5,032,239,631至少205天,那还怎么能算完?明天重新做个程序,分成更小的段,分成最多两天时间的计算量好了
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 21:00 | 显示全部楼层
4050018759与4050218759之间的素数开头最多有公差为6469693230的7个连续素数:用时3069.391秒
范围小了
回复 支持 反对

使用道具 举报

 楼主| 发表于 2026-4-13 21:15 | 显示全部楼层
本帖最后由 ysr 于 2026-4-13 13:36 编辑

4047918759与4049918759之间的素数开头最多有公差为6469693230的8个连续素数:用时13571.02秒
合3.77小时,这个程序稍快一点,(按这个速度计算,从4047918759算到5,032,239,631至少77天)代码如下:

Private Sub Command1_Click()
Dim a, b
a = Val(Text1)
a1 = a
Q = Val(Text2)
r = Val(Text4)
t = Timer
m = Sqr(Q)
If Right(a, 1) Mod 2 = 0 Then
a = a + 1
Else
a = a
End If

Do While InStr(fenjieyinzi(Val(a)), "*") <> 0
a = a + 2
Loop

Do While a <= Q
b = a
d = fenjieyinzi(Val(b))
s = 0
s3 = ""
Do While s >= 0 And InStr(d, "*") = 0
s4 = s + 1
If s4 Mod 20 = 0 Then
s3 = s3 & "/" & b & vbCrLf
Else
s3 = s3 & "/" & b
End If
s = s + 1
s2 = s * (s + 1)
b = b + r
d = fenjieyinzi(Val(b))
Loop
If s >= 10 Then
Text3 = Text3 & s3 & vbCrLf
End If
If s > 1 Then
s5 = s5 & "/" & s
Else
s5 = s5
End If
Print s5
a = a + 2
s2 = 0
Loop
s6 = Mid(s5, 2)
     s7 = max(Trim(s6))
Combo1 = a1 & "与" & Q & "之间的素数开头最多有公差为" & r & "的" & s7 & "个连续素数:用时" & Timer - t & "秒" & vbCrLf & Text3

     

End Sub

Private Sub Command2_Click()
Text1 = ""
Text2 = ""
Text3 = ""
Text4 = ""
Combo1 = ""
Form1.Cls
End Sub
Private Function max(ByVal lists As String) As String
Dim temp As Long
Dim a() As String
a = Split(lists, "/")
Dim b As Long

For b = 0 To UBound(a)
If temp < a(b) Then temp = a(b)
Next
max = temp
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
回复 支持 反对

使用道具 举报

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

本版积分规则

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

GMT+8, 2026-4-30 08:03 , Processed in 0.111575 second(s), 14 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表