|
[分享]孪素中项在孪素中项中的分拆-原有熊一兵介绍李明波加法猜想
2006年11月4日,天山草先生对李明波加法猜想给出如下电算结果:
'; 文件名: 孪中拆分展示
'; 计算大于等于 18 的“孪中数”拆分情况,给出每个拆分构成
Private Sub form_Click()
Open "F:\Microsoft Visual Studio\用户程序1\孪中拆分\ss.txt" For Input As 1 ';调入一千万以内的孪生素数表(最后一行末尾要补入一个任意数)
Open "F:\Microsoft Visual Studio\用户程序1\孪中拆分\dd简化2217.txt" For Output As 2 ';记录拆分结果
Dim s(100000) ';用于依次存放孪中偶数
Dim pp(10) As String
k = 0
For i = 1 To 5898 * 2 ';将孪生素数表调入 5898*2 行
For j = 1 To 5
k = k + 1
Input #1, s(k) ';读孪数的前一个,存入 s(i)
Input #1, w ';读孪数的后一个,舍弃
s(k) = s(k) + 1 ';化为孪中
Next j
Next i
For m = 4 To 58980 ';孪中数的序号
d = s(m)
n = d / 2
k = 0
For i = 2 To n
If s(i) > d / 2 Then GoTo 20
For j = i + 1 To n
If s(i) + s(j) > d Then GoTo 10
If s(i) + s(j) = d Then
k = k + 1
';a$ = "(" + Trim$(Str$(k)) + ") " + Trim$(Str$(d)) + "=" + Trim$(Str$(s(i))) + "+" + Trim$(Str$(s(j)))
' rint a$
' rint #2, a$ ’记录具体拆分情况,如不要,则只记录拆分数
GoTo 10
End If
Next j
10: Next i
20: Print " Z("; Trim$(Str$(m)); "): D("; Trim$(Str$(d)); ")="; k
Print #2, " Z("; Trim$(Str$(m)); "): D("; Trim$(Str$(d)); ")="; k
Next m
999: Close
End Sub
运行结果:
Z(4): D(18)= 1
Z(5): D(30)= 1
Z(6): D(42)= 1
Z(7): D(60)= 1
Z(8): D(72)= 2
Z(9): D(102)= 2
Z(10): D(108)= 1
……………………
如果看每个“孪中数”的具体拆分情况,运行结果为:
(1) 18=6+12
D(18)= 1
(1) 30=12+18
D(30)= 1
(1) 42=12+30
D(42)= 1
(1) 60=18+42
D(60)= 1
(1) 72=12+60
(2) 72=30+42
D(72)= 2
(1) 102=30+72
(2) 102=42+60
D(102)= 2
(1) 108=6+102
D(108)= 1
(1) 138=30+108
D(138)= 1
(1) 150=12+138
(2) 150=42+108
D(150)= 2
……………………
|
|